one_engine53.zip

one_engine53.zip

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 operational mode （Mode =

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件        2062  2001-08-10 11:53  savedata.m
     文件       68047  2001-08-14 16:31  complex1.mdl
     文件        2604  2001-08-14 16:36  gen_nois.m
     文件         677  2001-08-09 20:54  ru_func.m
     文件       60321  2001-08-09 20:34  compsim.mdl
     文件        2890  2001-08-09 20:31  miss_det.m
     文件       64019  2001-08-09 20:25  complex.mdl
     文件        2245  2001-08-09 20:21  gen_ref.m
     文件        1408  2001-08-09 18:28  overld.m
     文件        7916  2001-08-09 17:42  compdata.m
     文件        7804  2001-08-09 17:28  gen_extf.mdl
     文件        2283  2001-08-09 15:02  TQ_const.m
     文件         421  1998-04-30 10:38  plotdata.m
     文件        4597  1998-04-08 11:35  gen_extf.m
     文件        3276  1998-03-24 14:15  dmi_tow.m
     文件        3390  1998-03-12 12:28  comb_n.m
     文件        1182  1998-03-11 16:26  limit_y.m
     文件        1650  1998-03-11 15:49  cpp88new.m
     文件       19968  2001-08-14 16:47  readme.1st.doc