f14c.m

数字通信第四版原书的例程

		'Port','3',...
		'position',[160,293,185,317])
add_line([sys,'/',['Aircraft Dynamics ',13,'Model']],[430,130;505,130])
add_line([sys,'/',['Aircraft Dynamics ',13,'Model']],[465,130;465,190;430,255;350,255])
add_line([sys,'/',['Aircraft Dynamics ',13,'Model']],[430,305;510,305])
add_line([sys,'/',['Aircraft Dynamics ',13,'Model']],[465,305;465,255;430,195;350,195])
add_line([sys,'/',['Aircraft Dynamics ',13,'Model']],[175,375;235,375;235,320;265,320])
add_line([sys,'/',['Aircraft Dynamics ',13,'Model']],[170,35;235,35;235,115;265,115])
add_line([sys,'/',['Aircraft Dynamics ',13,'Model']],[300,195;235,195;235,145;265,145])
add_line([sys,'/',['Aircraft Dynamics ',13,'Model']],[190,305;265,305])
add_line([sys,'/',['Aircraft Dynamics ',13,'Model']],[295,305;360,305])
add_line([sys,'/',['Aircraft Dynamics ',13,'Model']],[295,130;360,130])
add_line([sys,'/',['Aircraft Dynamics ',13,'Model']],[300,255;235,255;235,290;265,290])
add_line([sys,'/',['Aircraft Dynamics ',13,'Model']],[185,130;265,130])
add_line([sys,'/',['Aircraft Dynamics ',13,'Model']],[65,35;125,35])
add_line([sys,'/',['Aircraft Dynamics ',13,'Model']],[90,35;90,375;130,375])


%     Finished composite block ['Aircraft Dynamics ',13,'Model'].

set_param([sys,'/',['Aircraft Dynamics ',13,'Model']],...
		'Font Number',720899,...
		'Drop Shadow',720896,...
		'position',[530,200,560,250])

add_block('built-in/Note',[sys,'/','q'])
set_param([sys,'/','q'],...
		'Font Number',720899,...
		'Drop Shadow',720896,...
		'position',[590,220,591,221])

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/Clock',[sys,'/','Clock'])
set_param([sys,'/','Clock'],...
		'Font Number',720899,...
		'Drop Shadow',720896,...
		'position',[310,42,335,68])


%     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/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/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/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/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/Inport',[sys,'/',['Nz pilot ',13,'calculation/q']])
set_param([sys,'/',['Nz pilot ',13,'calculation/q']],...
		'Font Number',720899,...
		'Drop Shadow',720896,...
		'Port','2',...
		'position',[25,115,45,135])

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/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/Inport',[sys,'/',['Nz pilot ',13,'calculation/w']])
set_param([sys,'/',['Nz pilot ',13,'calculation/w']],...
		'Font Number',720899,...
		'Drop Shadow',720896,...
		'Port','1',...
		'position',[25,35,45,55])

add_block('built-in/Derivative',[sys,'/',['Nz pilot ',13,'calculation/Derivative']])
set_param([sys,'/',['Nz pilot ',13,'calculation/Derivative']],...
		'Font Number',720899,...
		'Drop Shadow',720896,...
		'position',[125,32,165,58])

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/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/Derivative',[sys,'/',['Nz pilot ',13,'calculation/Derivative1']])
set_param([sys,'/',['Nz pilot ',13,'calculation/Derivative1']],...
		'Font Number',720899,...
		'Drop Shadow',720896,...
		'position',[125,112,165,138])
add_line([sys,'/',['Nz pilot ',13,'calculation']],[50,125;85,125;85,250;155,250])
add_line([sys,'/',['Nz pilot ',13,'calculation']],[85,125;120,125])
add_line([sys,'/',['Nz pilot ',13,'calculation']],[55,265;155,265])
add_line([sys,'/',['Nz pilot ',13,'calculation']],[445,125;465,125])
add_line([sys,'/',['Nz pilot ',13,'calculation']],[365,125;385,125])
add_line([sys,'/',['Nz pilot ',13,'calculation']],[170,45;290,45;290,105;335,105])
add_line([sys,'/',['Nz pilot ',13,'calculation']],[50,45;120,45])
add_line([sys,'/',['Nz pilot ',13,'calculation']],[170,125;210,125])
add_line([sys,'/',['Nz pilot ',13,'calculation']],[270,125;335,125])
add_line([sys,'/',['Nz pilot ',13,'calculation']],[190,260;290,260;290,145;335,145])


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

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

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

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

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/Outport',[sys,'/','alpha (rad)'])
set_param([sys,'/','alpha (rad)'],...
		'Font Number',720899,...
		'Drop Shadow',720896,...
		'Port','1',...
		'position',[755,205,775,225])

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,'/','Nz Pilot (g)'])
set_param([sys,'/','Nz Pilot (g)'],...
		'Font Number',720899,...
		'Drop Shadow',720896,...
		'Port','4',...
		'position',[740,120,760,140])


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

new_system([sys,'/',['Command ',13,'Generator']])
set_param([sys,'/',['Command ',13,'Generator']],'Location',[86,67,533,300])

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_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/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/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/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/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/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_line([sys,'/',['Command ',13,'Generator']],[125,50;165,50])
add_line([sys,'/',['Command ',13,'Generator']],[140,50;140,125;165,125])
add_line([sys,'/',['Command ',13,'Generator']],[50,50;60,50])
add_line([sys,'/',['Command ',13,'Generator']],[255,125;265,125;265,65;290,65])
add_line([sys,'/',['Command ',13,'Generator']],[255,50;290,50])
add_line([sys,'/',['Command ',13,'Generator']],[325,60;370,60])


%     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/Gain',[sys,'/','Gain'])
set_param([sys,'/','Gain'],...
		'hide name',0,...
		'Font Number',720899,...
		'Drop Shadow',720896,...
		'Gain','Zw',...
		'position',[360,206,400,244])

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,'/','Gain1'])
set_param([sys,'/','Gain1'],...
		'hide name',0,...
		'Font Number',720899,...
		'Drop Shadow',720896,...
		'Gain','Mq',...
		'position',[360,305,400,345])


%     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/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_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/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/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/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_line([sys,'/',['Dryden Wind ',13,'Gust Models']],[65,50;105,50])
add_line([sys,'/',['Dryden Wind ',13,'Gust Models']],[295,185;370,185])
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])


%     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/Inport',[sys,'/',['Noise',13,'Input']])
set_param([sys,'/',['Noise',13,'Input']],...
		'Font Number',720899,...
		'Drop Shadow',720896,...
		'Port','3',...
		'position',[180,280,200,300])
add_line(sys,[340,55;365,55])
add_line(sys,[70,50;95,50;95,90])
add_line(sys,[185,165;215,165;215,195;225,195])
add_line(sys,[565,240;630,240;630,70;650,70])
add_line(sys,[630,240;630,360;105,360;105,155;145,155])
add_line(sys,[630,300;755,300])
add_line(sys,[740,300;740,355;770,355])
add_line(sys,[400,415;355,415])
add_line(sys,[205,290;235,290])
add_line(sys,[145,220;125,220;125,165;145,165])
add_line(sys,[710,215;750,215])
add_line(sys,[720,215;720,375;225,375;225,395;180,395])
add_line(sys,[150,405;85,405;85,145;145,145])
add_line(sys,[245,415;180,415])
add_line(sys,[105,120;105,135;145,135])
add_line(sys,[565,215;650,215])
add_line(sys,[610,215;610,45;650,45])
add_line(sys,[275,275;355,275])
add_line(sys,[320,275;320,225;355,225])
add_line(sys,[475,300;490,300;490,240;525,240])
add_line(sys,[185,140;280,140])
add_line(sys,[405,325;425,325;425,310;445,310])
add_line(sys,[405,275;425,275;425,290;445,290])
add_line(sys,[275,300;320,300;320,325;355,325])
add_line(sys,[405,225;525,225])
add_line(sys,[350,140;510,140;510,210;525,210])
add_line(sys,[135,50;110,50;110,90])
add_line(sys,[695,60;725,60])
add_line(sys,[705,60;705,130;735,130])
% Return any arguments.
if (nargin | nargout)
	% Must use feval here to access system in memory
	if (nargin > 3)
		if (flag == 0)
			eval(['[ret,x0,xstr]=',sys,'(t,x,u,flag);'])
		else
			eval(['ret =', sys,'(t,x,u,flag);'])
		end
	else
		[ret,x0,str] = feval(sys);
	end
end