详细的螺旋桨matlab模型

很详细的船舶推进系统模型，包括桨距角控制、主轴转速控制等。

function u_out = comb_n_ref（u）
% 生成过程转速和桨距角指令
% This function calculates the command setpoint signals n_com and
% theta_com for the engines and propeller under chosen mode:
% （Economy or Manouvering）

% ============================================================
%    Inputs :: 
%	   u（1）     : Mode （1: Economy 2: Manouvering）
%	   u（2）	: Port Handle at Center            PH_C  !! active in Mode 1 and 2.
% ------------------------------------------------------------
%    Outputs:: 
%          n_com    : reference revolution per minut 
%	   theta_com: reference pitch angle          
% ------------------------------------------------------------
%   Date : 10-03-1998
%   By   : Roozbeh Izadi-Zamanabadi
% ============================================================

Mode = u（1）;

PH_C = max（min（u（2）10） -10）;

Xp_PH_C = [-10   -7  -6   -4.6  -3];
Yp_PH_C = [-0.41 -0.41 -0.41 -0.22 -0.055];
Xp_PH_C = [Xp_PH_C -2.99   -0.01     0    0.9   1.6 1.8   2   3.6   5.3   6    6.2  7   7.1  7.5  8    8.5   9    9.5  10 ];
Yp_PH_C = [Yp_PH_C -0.0549 -0.0549 -0.055 0.075 0.2 0.18 0.24 0.395 0.62 0.72 0.71 0.87 0.91 0.94 0.94 0.94 0.94 0.94 0.94];

Xn_PH_C = [-10   -9.8   -8   -6.2 -5.65 -0.6  -0.5  4.3   5.3  6   6.2   7    8.9  10];
Yn_PH_C = [30.62 30.62 28.57 25.2 23.55 23.55 23.75 23.75 26.4 27.5 27.7 28.99 30 30.5]/2.45;


%% ----------------------------------------------------------
%% 		Calculation of setpoint commands for 
%%		Economy operational mode （Mode = 1）
%% ----------------------------------------------------------

if Mode == 1  
%% -- Calculation of setpoint command for shaft speed ------
	if abs（PH_C）<= 4
		n_com = 9.7;
	else
		n_com = 0.4667*abs（PH_C）+7.833;
	end
%% -- Calculation of setpoint command for pitch angle ------ 
	if PH_C<= （-5.43）
		theta_com = -0.39;
	elseif （PH_C> -5.43） & （PH_C <= -2）
		theta_com = 0.1137*PH_C + 0.2274;
	elseif （PH_C> -2） & （PH_C <= 0）
		theta_com = 0;
	elseif （PH_C> 0） & （PH_C <= 8）
		theta_com = 0.1138*PH_C;
	else
		theta_com = 0.91;
	end

elseif Mode == 2
%% ----------------------------------------------------------
%% 		Calculation of setpoint commands for 
%%		Manouvering operational mode （Mode = 2）
%% ----------------------------------------------------------
%% -- Calculation of setpoint command for shaft speed ------ 

	%% ---------------- Shaft speed  ------------------%%
		n_com = interp1（Xn_PH_CYn_PH_CPH_C‘linear‘）;

%% -- Calculation of setpoint command for pitch angle ------ 

	%% ---------------- Propeller angle  ------------------%%
		theta_com = interp1（Xp_PH_CYp_PH_CPH_C‘linear‘）;


%% ************************************************************
%% 		The third mode is removed for the ferry
%% ************************************************************ 
%elseif Mode == 3
%% ----------------------------------------------------------
%% 		Calculation of setpoint commands for 
%%		Separate operationa

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2018-12-13 15:04  one_engine00\
     文件        3416  2018-11-13 09:09  one_engine00\comb_n.m
     文件        8150  2018-11-13 09:07  one_engine00\compdata.asv
     文件        8160  2018-11-13 09:07  one_engine00\compdata.m
     文件       64019  2001-08-09 20:25  one_engine00\complex.mdl
     文件      105847  2018-11-13 09:16  one_engine00\complex1.mdl
     文件       60321  2001-08-09 20:34  one_engine00\compsim.mdl
     文件        1659  2018-11-13 09:10  one_engine00\cpp88new.m
     文件        3276  1998-03-24 14:15  one_engine00\dmi_tow.m
     文件        4597  1998-04-08 11:35  one_engine00\gen_extf.m
     文件        7804  2001-08-09 17:28  one_engine00\gen_extf.mdl
     文件        2619  2018-11-13 09:12  one_engine00\gen_nois.m
     文件        2208  2018-11-13 09:13  one_engine00\gen_ref.m
     文件        1182  1998-03-11 16:26  one_engine00\limit_y.m
     文件        2890  2001-08-09 20:31  one_engine00\miss_det.m
     文件       27920  2018-11-13 08:55  one_engine00\n_p_ref.mat
     文件      122130  2018-12-13 15:04  one_engine00\noise_E.mat
     文件      258414  2018-12-13 15:04  one_engine00\np_ref.mat
     文件        1408  2001-08-09 18:28  one_engine00\overld.m
     文件         421  1998-04-30 10:38  one_engine00\plotdata.m
     文件         677  2001-08-09 20:54  one_engine00\ru_func.m
     文件        2062  2001-08-10 11:53  one_engine00\savedata.m
     文件       63892  2018-12-13 15:04  one_engine00\T_extern.mat
     文件        2283  2001-08-09 15:02  one_engine00\TQ_const.m
     文件       16111  2018-12-13 14:58  one_engine00\模型说明.docx