f14n.m

本书是电子通信类的本科、研究生辅助教材

function [ret,x0,str]=f14n(t,x,u,flag);
%F14N	is the M-file description of the SIMULINK system named F14N.
%	The block-diagram can be displayed by typing: F14N.
%
%	SYS=F14N(T,X,U,FLAG) returns depending on FLAG certain
%	system values given time point, T, current state vector, X,
%	and input vector, U.
%	FLAG is used to indicate the type of output to be returned in SYS.
%
%	Setting FLAG=1 causes F14N to return state derivitives, FLAG=2
%	discrete states, FLAG=3 system outputs and FLAG=4 next sample
%	time. For more information and other options see SFUNC.
%
%	Calling F14N with a FLAG of zero:
%	[SIZES]=F14N([],[],[],0),  returns a vector, SIZES, which
%	contains the sizes of the state vector and other parameters.
%		SIZES(1) number of states
%		SIZES(2) number of discrete states
%		SIZES(3) number of outputs
%		SIZES(4) number of inputs.
%	For the definition of other parameters in SIZES, see SFUNC.
%	See also, TRIM, LINMOD, LINSIM, EULER, RK23, RK45, ADAMS, GEAR.

%	Copyright (c) 1990-94 by The MathWorks, Inc.

% Note: This M-file is only used for saving graphical information;
%       after the model is loaded into memory an internal model
%       representation is used.

% the system will take on the name of this mfile: 
sys = mfilename;
new_system(sys)
simver(1.3)
if(0 == (nargin + nargout))
     set_param(sys,'Location',[2,13,650,456])
     open_system(sys)
end;
set_param(sys,'algorithm',		'Linear')
set_param(sys,'Start time',	'0.0')
set_param(sys,'Stop time',		'2')
set_param(sys,'Min step size',	'1e-8')
set_param(sys,'Max step size',	'.01')
set_param(sys,'Relative error','1e-4')
set_param(sys,'Return vars',	'')

add_block('built-in/Inport',[sys,'/',['Noise',13,'Input']])
set_param([sys,'/',['Noise',13,'Input']],...
		'Font Number',720899,...
		'Drop Shadow',720896,...
		'Port','3',...
		'position',[180,280,200,300])


%     Subsystem  ['Dryden Wind',13,'Gust Models'].

new_system([sys,'/',['Dryden Wind',13,'Gust Models']])
set_param([sys,'/',['Dryden Wind',13,'Gust Models']],'Location',[115,54,581,332])

add_block('built-in/Transfer Fcn',[sys,'/',['Dryden Wind',13,'Gust Models/Q-gust model']])
set_param([sys,'/',['Dryden Wind',13,'Gust Models/Q-gust model']],...
		'Font Number',720899,...
		'Drop Shadow',720896,...
		'Numerator','pi/(4*b)*[1 0]',...
		'Denominator','[1 pi*Vto/(4*b)]',...
		'position',[150,160,290,210])

add_block('built-in/Outport',[sys,'/',['Dryden Wind',13,'Gust Models/wGust',13,'(ft//sec)']])
set_param([sys,'/',['Dryden Wind',13,'Gust Models/wGust',13,'(ft//sec)']],...
		'Font Number',720899,...
		'Drop Shadow',720896,...
		'Port','1',...
		'position',[375,37,400,63])

add_block('built-in/Outport',[sys,'/',['Dryden Wind',13,'Gust Models/qGust',13,'(rad//sec)']])
set_param([sys,'/',['Dryden Wind',13,'Gust Models/qGust',13,'(rad//sec)']],...
		'Font Number',720899,...
		'Drop Shadow',720896,...
		'Port','2',...
		'position',[375,172,400,198])

add_block('built-in/Transfer Fcn',[sys,'/',['Dryden Wind',13,'Gust Models/W-gust model']])
set_param([sys,'/',['Dryden Wind',13,'Gust Models/W-gust model']],...
		'Font Number',720899,...
		'Drop Shadow',720896,...
		'Numerator','Swg/sqrt(a^3)*[sqrt(3)*a,1]',...
		'Denominator','[1 (2/a) 1/a^2]',...
		'position',[110,22,330,78])

add_block('built-in/Inport',[sys,'/',['Dryden Wind',13,'Gust Models/Noise',13,'Input']])
set_param([sys,'/',['Dryden Wind',13,'Gust Models/Noise',13,'Input']],...
		'Font Number',720899,...
		'Drop Shadow',720896,...
		'Port','1',...
		'position',[35,37,60,63])
add_line([sys,'/',['Dryden Wind',13,'Gust Models']],[335,50;370,50])
add_line([sys,'/',['Dryden Wind',13,'Gust Models']],[350,50;350,115;115,115;115,185;145,185])
add_line([sys,'/',['Dryden Wind',13,'Gust Models']],[295,185;370,185])
add_line([sys,'/',['Dryden Wind',13,'Gust Models']],[65,50;105,50])


%     Finished composite block ['Dryden Wind',13,'Gust Models'].

set_param([sys,'/',['Dryden Wind',13,'Gust Models']],...
		'Font Number',720899,...
		'Drop Shadow',720896,...
		'position',[240,262,270,313])

add_block('built-in/Gain',[sys,'/','Gain1'])
set_param([sys,'/','Gain1'],...
		'hide name',0,...
		'Font Number',720899,...
		'Drop Shadow',720896,...
		'Gain','Mq',...
		'position',[360,305,400,345])

add_block('built-in/Gain',[sys,'/','Gain2'])
set_param([sys,'/','Gain2'],...
		'hide name',0,...
		'Font Number',720899,...
		'Drop Shadow',720896,...
		'Gain','Mw',...
		'position',[360,253,400,297])

add_block('built-in/Gain',[sys,'/','Gain'])
set_param([sys,'/','Gain'],...
		'hide name',0,...
		'Font Number',720899,...
		'Drop Shadow',720896,...
		'Gain','Zw',...
		'position',[360,206,400,244])


%     Subsystem  ['Command',13,'Generator'].

new_system([sys,'/',['Command',13,'Generator']])
set_param([sys,'/',['Command',13,'Generator']],'Location',[62,227,509,447])

add_block('built-in/Outport',[sys,'/',['Command',13,'Generator/alpha',13,'command']])
set_param([sys,'/',['Command',13,'Generator/alpha',13,'command']],...
		'Font Number',720899,...
		'Drop Shadow',720896,...
		'Port','1',...
		'position',[375,50,395,70])

add_block('built-in/Transfer Fcn',[sys,'/',['Command',13,'Generator/Slow Pole']])
set_param([sys,'/',['Command',13,'Generator/Slow Pole']],...
		'Font Number',720899,...
		'Drop Shadow',720896,...
		'Numerator','[0.01]',...
		'Denominator','[ 1 0.01]',...
		'position',[170,22,250,78])

add_block('built-in/Constant',[sys,'/',['Command',13,'Generator/Constant']])
set_param([sys,'/',['Command',13,'Generator/Constant']],...
		'Font Number',720899,...
		'Drop Shadow',720896,...
		'Value','1',...
		'position',[20,37,45,63])

add_block('built-in/Transfer Fcn',[sys,'/',['Command',13,'Generator/Fast Pole']])
set_param([sys,'/',['Command',13,'Generator/Fast Pole']],...
		'Font Number',720899,...
		'Drop Shadow',720896,...
		'Numerator','[426.435]',...
		'Denominator','[ 1 426.435]',...
		'position',[170,96,250,154])

add_block('built-in/Sum',[sys,'/',['Command',13,'Generator/Sum1']])
set_param([sys,'/',['Command',13,'Generator/Sum1']],...
		'hide name',0,...
		'Font Number',720899,...
		'Drop Shadow',720896,...
		'inputs','-+',...
		'position',[295,40,320,75])

add_block('built-in/Gain',[sys,'/',['Command',13,'Generator/Gain1']])
set_param([sys,'/',['Command',13,'Generator/Gain1']],...
		'hide name',0,...
		'Font Number',720899,...
		'Drop Shadow',720896,...
		'Gain','0.0349',...
		'position',[65,21,120,79])

add_block('built-in/Note',[sys,'/',['Command',13,'Generator/Second Order Command Generator Model']])
set_param([sys,'/',['Command',13,'Generator/Second Order Command Generator Model']],...
		'Font Number',720897,...
		'Drop Shadow',720896,...
		'position',[210,180,211,181])
add_line([sys,'/',['Command',13,'Generator']],[325,60;370,60])
add_line([sys,'/',['Command',13,'Generator']],[255,50;290,50])
add_line([sys,'/',['Command',13,'Generator']],[255,125;265,125;265,65;290,65])
add_line([sys,'/',['Command',13,'Generator']],[50,50;60,50])
add_line([sys,'/',['Command',13,'Generator']],[125,50;165,50])
add_line([sys,'/',['Command',13,'Generator']],[140,50;140,125;165,125])


%     Finished composite block ['Command',13,'Generator'].

set_param([sys,'/',['Command',13,'Generator']],...
		'Font Number',720899,...
		'Drop Shadow',720896,...
		'position',[35,25,65,75])

add_block('built-in/Outport',[sys,'/','Nz Pilot (g)'])
set_param([sys,'/','Nz Pilot (g)'],...
		'Font Number',720899,...
		'Drop Shadow',720896,...
		'Port','4',...
		'position',[740,120,760,140])

add_block('built-in/Sum',[sys,'/','Sum1'])
set_param([sys,'/','Sum1'],...
		'orientation',1,...
		'hide name',0,...
		'Font Number',720899,...
		'Drop Shadow',720896,...
		'inputs','++',...
		'position',[85,95,120,115])

add_block('built-in/Outport',[sys,'/','alpha (rad)'])
set_param([sys,'/','alpha (rad)'],...
		'Font Number',720899,...
		'Drop Shadow',720896,...
		'Port','1',...
		'position',[755,200,775,220])

add_block('built-in/Inport',[sys,'/',['alpha',13,'command']])
set_param([sys,'/',['alpha',13,'command']],...
		'orientation',2,...
		'Font Number',720899,...
		'Drop Shadow',720896,...
		'Port','1',...
		'position',[140,40,160,60])

add_block('built-in/Note',[sys,'/','w'])
set_param([sys,'/','w'],...
		'Font Number',720899,...
		'Drop Shadow',720896,...
		'position',[645,195,646,196])

add_block('built-in/Gain',[sys,'/','Gain5'])
set_param([sys,'/','Gain5'],...
		'Font Number',720899,...
		'Drop Shadow',720896,...
		'Gain','1/Uo',...
		'position',[655,186,705,234])


%     Subsystem  ['Nz pilot',13,'calculation'].

new_system([sys,'/',['Nz pilot',13,'calculation']])
set_param([sys,'/',['Nz pilot',13,'calculation']],'Location',[50,50,587,397])

add_block('built-in/Outport',[sys,'/',['Nz pilot',13,'calculation/Pilot',13,'g force (g)']])
set_param([sys,'/',['Nz pilot',13,'calculation/Pilot',13,'g force (g)']],...
		'Font Number',720899,...
		'Drop Shadow',720896,...
		'Port','1',...
		'position',[470,115,490,135])

add_block('built-in/Sum',[sys,'/',['Nz pilot',13,'calculation/Sum1']])
set_param([sys,'/',['Nz pilot',13,'calculation/Sum1']],...
		'Font Number',720899,...
		'Drop Shadow',720896,...
		'inputs','-++',...
		'position',[340,95,360,155])

add_block('built-in/Inport',[sys,'/',['Nz pilot',13,'calculation/wDot']])
set_param([sys,'/',['Nz pilot',13,'calculation/wDot']],...
		'Font Number',720899,...
		'Drop Shadow',720896,...
		'Port','1',...
		'position',[25,35,45,55])

add_block('built-in/Product',[sys,'/',['Nz pilot',13,'calculation/Product']])
set_param([sys,'/',['Nz pilot',13,'calculation/Product']],...
		'Font Number',720899,...
		'Drop Shadow',720896,...
		'inputs','2',...
		'position',[160,242,185,273])

add_block('built-in/Note',[sys,'/',['Nz pilot',13,'calculation/Nz pilot calculation']])
set_param([sys,'/',['Nz pilot',13,'calculation/Nz pilot calculation']],...
		'Font Number',720899,...
		'Drop Shadow',720896,...
		'position',[369,270,370,271])

add_block('built-in/Inport',[sys,'/',['Nz pilot',13,'calculation/qDot']])
set_param([sys,'/',['Nz pilot',13,'calculation/qDot']],...
		'Font Number',720899,...
		'Drop Shadow',720896,...
		'Port','2',...
		'position',[25,115,45,135])

add_block('built-in/Gain',[sys,'/',['Nz pilot',13,'calculation/Gain2']])
set_param([sys,'/',['Nz pilot',13,'calculation/Gain2']],...
		'Font Number',720899,...
		'Drop Shadow',720896,...
		'Gain','1/32.2',...
		'position',[390,102,440,148])

add_block('built-in/Gain',[sys,'/',['Nz pilot',13,'calculation/Gain1']])
set_param([sys,'/',['Nz pilot',13,'calculation/Gain1']],...
		'Font Number',720899,...
		'Drop Shadow',720896,...
		'Gain','22.8',...
		'position',[215,102,265,148])

add_block('built-in/Constant',[sys,'/',['Nz pilot',13,'calculation/Constant']])
set_param([sys,'/',['Nz pilot',13,'calculation/Constant']],...
		'Font Number',720899,...
		'Drop Shadow',720896,...
		'Value','Uo',...
		'position',[25,252,50,278])

add_block('built-in/Note',[sys,'/',['Nz pilot',13,'calculation/Nz = (dq//dt*22.8 + q*Uo - dw//dt)//g']])
set_param([sys,'/',['Nz pilot',13,'calculation/Nz = (dq//dt*22.8 + q*Uo - dw//dt)//g']],...
		'Font Number',720899,...
		'Drop Shadow',720896,...
		'position',[368,295,369,296])

add_block('built-in/Inport',[sys,'/',['Nz pilot',13,'calculation/q']])
set_param([sys,'/',['Nz pilot',13,'calculation/q']],...
		'Font Number',720899,...
		'Drop Shadow',720896,...
		'Port','3',...
		'position',[25,190,45,210])
add_line([sys,'/',['Nz pilot',13,'calculation']],[190,260;290,260;290,145;335,145])
add_line([sys,'/',['Nz pilot',13,'calculation']],[270,125;335,125])
add_line([sys,'/',['Nz pilot',13,'calculation']],[50,45;290,45;290,105;335,105])
add_line([sys,'/',['Nz pilot',13,'calculation']],[365,125;385,125])
add_line([sys,'/',['Nz pilot',13,'calculation']],[445,125;465,125])
add_line([sys,'/',['Nz pilot',13,'calculation']],[55,265;155,265])
add_line([sys,'/',['Nz pilot',13,'calculation']],[50,125;210,125])
add_line([sys,'/',['Nz pilot',13,'calculation']],[50,200;105,200;105,250;155,250])


%     Finished composite block ['Nz pilot',13,'calculation'].

set_param([sys,'/',['Nz pilot',13,'calculation']],...
		'Font Number',720899,...
		'Drop Shadow',720896,...
		'position',[655,32,690,88])

add_block('built-in/Clock',[sys,'/','Clock'])
set_param([sys,'/','Clock'],...
		'Font Number',720899,...
		'Drop Shadow',720896,...
		'position',[310,42,335,68])

add_block('built-in/Sum',[sys,'/','Sum'])
set_param([sys,'/','Sum'],...
		'hide name',0,...
		'Font Number',720899,...
		'Drop Shadow',720896,...
		'inputs','++',...
		'position',[450,280,470,320])

add_block('built-in/Note',[sys,'/','q'])
set_param([sys,'/','q'],...
		'Font Number',720899,...
		'Drop Shadow',720896,...
		'position',[640,250,641,251])


%     Subsystem  ['Aircraft Dynamics',13,'Model'].

new_system([sys,'/',['Aircraft Dynamics',13,'Model']])
set_param([sys,'/',['Aircraft Dynamics',13,'Model']],'Location',[41,8,640,480])

add_block('built-in/Integrator',[sys,'/',['Aircraft Dynamics',13,'Model/Integrator1']])
set_param([sys,'/',['Aircraft Dynamics',13,'Model/Integrator1']],...
		'Font Number',720897,...
		'Drop Shadow',720896,...
		'Initial','0',...
		'position',[425,315,445,335])

add_block('built-in/Sum',[sys,'/',['Aircraft Dynamics',13,'Model/Sum3']])
set_param([sys,'/',['Aircraft Dynamics',13,'Model/Sum3']],...
		'Font Number',720897,...
		'Drop Shadow',720896,...
		'inputs','++',...
		'position',[360,315,380,335])

add_block('built-in/Gain',[sys,'/',['Aircraft Dynamics',13,'Model/Gain8']])
set_param([sys,'/',['Aircraft Dynamics',13,'Model/Gain8']],...
		'orientation',2,...
		'Font Number',720899,...
		'Drop Shadow',720896,...
		'Gain','Zw',...
		'position',[380,360,420,400])

add_block('built-in/Gain',[sys,'/',['Aircraft Dynamics',13,'Model/Gain7']])
set_param([sys,'/',['Aircraft Dynamics',13,'Model/Gain7']],...
		'orientation',2,...
		'Font Number',720899,...
		'Drop Shadow',720896,...
		'Gain','Zw',...
		'position',[355,20,395,60])

add_block('built-in/Sum',[sys,'/',['Aircraft Dynamics',13,'Model/Sum']])
set_param([sys,'/',['Aircraft Dynamics',13,'Model/Sum']],...
		'Font Number',720897,...
		'Drop Shadow',720896,...
		'inputs','++',...
		'position',[330,95,350,115])

add_block('built-in/Integrator',[sys,'/',['Aircraft Dynamics',13,'Model/Integrator']])
set_param([sys,'/',['Aircraft Dynamics',13,'Model/Integrator']],...
		'Font Number',720897,...
		'Drop Shadow',720896,...
		'Initial','0',...
		'position',[405,95,425,115])

add_block('built-in/Inport',[sys,'/',['Aircraft Dynamics',13,'Model/Rotary Gust',13,'qGust (rad//sec)']])
set_param([sys,'/',['Aircraft Dynamics',13,'Model/Rotary Gust',13,'qGust (rad//sec)']],...