tstpll.m

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

function [ret,x0,str,ts,xts]=tstpll(t,x,u,flag);
%TSTPLL	is the M-file description of the SIMULINK system named TSTPLL.
%	The block-diagram can be displayed by typing: TSTPLL.
%
%	SYS=TSTPLL(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 TSTPLL to return state derivatives, FLAG=2
%	discrete states, FLAG=3 system outputs and FLAG=4 next sample
%	time. For more information and other options see SFUNC.
%
%	Calling TSTPLL with a FLAG of zero:
%	[SIZES]=TSTPLL([],[],[],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
%		SIZES(5) number of roots (currently unsupported)
%		SIZES(6) direct feedthrough flag
%		SIZES(7) number of sample times
%
%	For the definition of other parameters in SIZES, see SFUNC.
%	See also, TRIM, LINMOD, LINSIM, EULER, RK23, RK45, ADAMS, GEAR.

% 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',[84,271,719,701])
     open_system(sys)
end;
set_param(sys,'algorithm',     'RK-45')
set_param(sys,'Start time',    '0.0')
set_param(sys,'Stop time',     '999999')
set_param(sys,'Min step size', '0.0001')
set_param(sys,'Max step size', '.005')
set_param(sys,'Relative error','1e-3')
set_param(sys,'Return vars',   '')

add_block('built-in/Note',[sys,'/','The curve color in the graphical scope matches the color of the blocks in the block diagram'])
set_param([sys,'/','The curve color in the graphical scope matches the color of the blocks in the block diagram'],...
		'position',[295,391,300,396])


%     Subsystem  ['Charge pump',13,'PLL'].

new_system([sys,'/',['Charge pump',13,'PLL']])
set_param([sys,'/',['Charge pump',13,'PLL']],'Location',[156,174,638,416])

add_block('built-in/Outport',[sys,'/',['Charge pump',13,'PLL/out_2']])
set_param([sys,'/',['Charge pump',13,'PLL/out_2']],...
		'Port','2',...
		'position',[320,25,340,45])

add_block('built-in/Outport',[sys,'/',['Charge pump',13,'PLL/out_3']])
set_param([sys,'/',['Charge pump',13,'PLL/out_3']],...
		'Port','3',...
		'position',[370,205,390,225])

add_block('built-in/Transfer Fcn',[sys,'/',['Charge pump',13,'PLL/Transfer Fcn']])
set_param([sys,'/',['Charge pump',13,'PLL/Transfer Fcn']],...
		'Numerator','Fn',...
		'Denominator','Fd',...
		'position',[240,59,290,101])


%     Subsystem  ['Charge pump',13,'PLL/Voltage-controlled',13,'oscilator'].

new_system([sys,'/',['Charge pump',13,'PLL/Voltage-controlled',13,'oscilator']])
set_param([sys,'/',['Charge pump',13,'PLL/Voltage-controlled',13,'oscilator']],'Location',[422,82,872,213])

add_block('built-in/Outport',[sys,'/',['Charge pump',13,'PLL/Voltage-controlled',13,'oscilator/out_1']])
set_param([sys,'/',['Charge pump',13,'PLL/Voltage-controlled',13,'oscilator/out_1']],...
		'position',[410,30,430,50])

add_block('built-in/Gain',[sys,'/',['Charge pump',13,'PLL/Voltage-controlled',13,'oscilator/Sensitivity']])
set_param([sys,'/',['Charge pump',13,'PLL/Voltage-controlled',13,'oscilator/Sensitivity']],...
		'Gain','Kc',...
		'position',[65,21,105,49])

add_block('built-in/Inport',[sys,'/',['Charge pump',13,'PLL/Voltage-controlled',13,'oscilator/in_1']])
set_param([sys,'/',['Charge pump',13,'PLL/Voltage-controlled',13,'oscilator/in_1']],...
		'position',[20,25,40,45])

add_block('built-in/Sum',[sys,'/',['Charge pump',13,'PLL/Voltage-controlled',13,'oscilator/Sum']])
set_param([sys,'/',['Charge pump',13,'PLL/Voltage-controlled',13,'oscilator/Sum']],...
		'position',[150,30,170,50])


%     Subsystem  ['Charge pump',13,'PLL/Voltage-controlled',13,'oscilator/Module',13,'Integrator'].

new_system([sys,'/',['Charge pump',13,'PLL/Voltage-controlled',13,'oscilator/Module',13,'Integrator']])
set_param([sys,'/',['Charge pump',13,'PLL/Voltage-controlled',13,'oscilator/Module',13,'Integrator']],'Location',[529,321,954,458])

add_block('built-in/Outport',[sys,'/',['Charge pump',13,'PLL/Voltage-controlled',13,'oscilator/Module',13,'Integrator/out_1']])
set_param([sys,'/',['Charge pump',13,'PLL/Voltage-controlled',13,'oscilator/Module',13,'Integrator/out_1']],...
		'position',[370,30,390,50])

add_block('built-in/Sum',[sys,'/',['Charge pump',13,'PLL/Voltage-controlled',13,'oscilator/Module',13,'Integrator/Sum']])
set_param([sys,'/',['Charge pump',13,'PLL/Voltage-controlled',13,'oscilator/Module',13,'Integrator/Sum']],...
		'position',[195,30,215,50])

add_block('built-in/Constant',[sys,'/',['Charge pump',13,'PLL/Voltage-controlled',13,'oscilator/Module',13,'Integrator/one']])
set_param([sys,'/',['Charge pump',13,'PLL/Voltage-controlled',13,'oscilator/Module',13,'Integrator/one']],...
		'orientation',3,...
		'move name',0,...
		'position',[70,55,90,75])

add_block('built-in/Reset Integrator',[sys,'/',['Charge pump',13,'PLL/Voltage-controlled',13,'oscilator/Module',13,'Integrator/Reset',13,'integrator']])
set_param([sys,'/',['Charge pump',13,'PLL/Voltage-controlled',13,'oscilator/Module',13,'Integrator/Reset',13,'integrator']],...
		'position',[110,19,150,51])

add_block('built-in/Fcn',[sys,'/',['Charge pump',13,'PLL/Voltage-controlled',13,'oscilator/Module',13,'Integrator/module']])
set_param([sys,'/',['Charge pump',13,'PLL/Voltage-controlled',13,'oscilator/Module',13,'Integrator/module']],...
		'Expr','rem(u[1],modu)',...
		'position',[235,29,325,51])

add_block('built-in/Memory',[sys,'/',['Charge pump',13,'PLL/Voltage-controlled',13,'oscilator/Module',13,'Integrator/Memory']])
set_param([sys,'/',['Charge pump',13,'PLL/Voltage-controlled',13,'oscilator/Module',13,'Integrator/Memory']],...
		'orientation',2,...
		'x0','init',...
		'position',[210,75,260,105])

add_block('built-in/Inport',[sys,'/',['Charge pump',13,'PLL/Voltage-controlled',13,'oscilator/Module',13,'Integrator/in_1']])
set_param([sys,'/',['Charge pump',13,'PLL/Voltage-controlled',13,'oscilator/Module',13,'Integrator/in_1']],...
		'position',[35,15,55,35])
add_line([sys,'/',['Charge pump',13,'PLL/Voltage-controlled',13,'oscilator/Module',13,'Integrator']],[330,40;365,40])
add_line([sys,'/',['Charge pump',13,'PLL/Voltage-controlled',13,'oscilator/Module',13,'Integrator']],[340,40;340,90;265,90])
add_line([sys,'/',['Charge pump',13,'PLL/Voltage-controlled',13,'oscilator/Module',13,'Integrator']],[80,50;80,35;105,35])
add_line([sys,'/',['Charge pump',13,'PLL/Voltage-controlled',13,'oscilator/Module',13,'Integrator']],[60,25;105,25])
add_line([sys,'/',['Charge pump',13,'PLL/Voltage-controlled',13,'oscilator/Module',13,'Integrator']],[220,40;230,40])
add_line([sys,'/',['Charge pump',13,'PLL/Voltage-controlled',13,'oscilator/Module',13,'Integrator']],[155,35;190,35])
add_line([sys,'/',['Charge pump',13,'PLL/Voltage-controlled',13,'oscilator/Module',13,'Integrator']],[205,90;180,90;190,45])
set_param([sys,'/',['Charge pump',13,'PLL/Voltage-controlled',13,'oscilator/Module',13,'Integrator']],...
		'Mask Display','Mod\nint',...
		'Mask Type','Module Integrator')
set_param([sys,'/',['Charge pump',13,'PLL/Voltage-controlled',13,'oscilator/Module',13,'Integrator']],...
		'Mask Dialogue','Module the given number when state exceed limit|Module bound:|Initial value:',...
		'Mask Translate','init=@2; modu=@1;')
set_param([sys,'/',['Charge pump',13,'PLL/Voltage-controlled',13,'oscilator/Module',13,'Integrator']],...
		'Mask Help','The absolute value of this integrator will not exceed the specified module bound. The state value is set to be zero when the module bound is reached. This block can only be used for scale value calculation.')
set_param([sys,'/',['Charge pump',13,'PLL/Voltage-controlled',13,'oscilator/Module',13,'Integrator']],...
		'Mask Entries','1\/0\/')


%     Finished composite block ['Charge pump',13,'PLL/Voltage-controlled',13,'oscilator/Module',13,'Integrator'].

set_param([sys,'/',['Charge pump',13,'PLL/Voltage-controlled',13,'oscilator/Module',13,'Integrator']],...
		'position',[195,24,250,56])

add_block('built-in/Fcn',[sys,'/',['Charge pump',13,'PLL/Voltage-controlled',13,'oscilator/sin']])
set_param([sys,'/',['Charge pump',13,'PLL/Voltage-controlled',13,'oscilator/sin']],...
		'Expr','Ac*sin(u[1]*pi2+Ph)',...
		'position',[270,25,395,55])

add_block('built-in/Constant',[sys,'/',['Charge pump',13,'PLL/Voltage-controlled',13,'oscilator/Carrier',13,'frequency']])
set_param([sys,'/',['Charge pump',13,'PLL/Voltage-controlled',13,'oscilator/Carrier',13,'frequency']],...
		'Value','Fc',...
		'position',[90,75,110,95])
add_line([sys,'/',['Charge pump',13,'PLL/Voltage-controlled',13,'oscilator']],[255,40;265,40])
add_line([sys,'/',['Charge pump',13,'PLL/Voltage-controlled',13,'oscilator']],[45,35;60,35])
add_line([sys,'/',['Charge pump',13,'PLL/Voltage-controlled',13,'oscilator']],[110,35;145,35])
add_line([sys,'/',['Charge pump',13,'PLL/Voltage-controlled',13,'oscilator']],[115,85;125,85;125,45;145,45])
add_line([sys,'/',['Charge pump',13,'PLL/Voltage-controlled',13,'oscilator']],[175,40;190,40])
add_line([sys,'/',['Charge pump',13,'PLL/Voltage-controlled',13,'oscilator']],[400,40;405,40])
set_param([sys,'/',['Charge pump',13,'PLL/Voltage-controlled',13,'oscilator']],...
		'Mask Display','VCO',...
		'Mask Type','VCO')
set_param([sys,'/',['Charge pump',13,'PLL/Voltage-controlled',13,'oscilator']],...
		'Mask Dialogue','Voltage-Controlled Oscilator:|Amplitude constant:|Carier frequency (Hz):|Sensitivity constant:|Initial phase (rad):')
set_param([sys,'/',['Charge pump',13,'PLL/Voltage-controlled',13,'oscilator']],...
		'Mask Translate','Ac=@1;Fc=@2;Kc=@3;Ph=@4;pi2=2*pi;')
set_param([sys,'/',['Charge pump',13,'PLL/Voltage-controlled',13,'oscilator']],...
		'Mask Help','The VCO block cause the frequency of the sinusoidal carrier to depend on the input signal voltage, resulting in frequency of the output change depending on the amplitude of the input signal.')
set_param([sys,'/',['Charge pump',13,'PLL/Voltage-controlled',13,'oscilator']],...
		'Mask Entries','Ac\/Fc\/1\/Ph\/')


%     Finished composite block ['Charge pump',13,'PLL/Voltage-controlled',13,'oscilator'].

set_param([sys,'/',['Charge pump',13,'PLL/Voltage-controlled',13,'oscilator']],...
		'orientation',2,...
		'position',[165,137,220,173])


%     Subsystem  ['Charge pump',13,'PLL/Phase//freq',13,'detector'].

new_system([sys,'/',['Charge pump',13,'PLL/Phase//freq',13,'detector']])
set_param([sys,'/',['Charge pump',13,'PLL/Phase//freq',13,'detector']],'Location',[14,335,990,748])

add_block('built-in/Constant',[sys,'/',['Charge pump',13,'PLL/Phase//freq',13,'detector/Constant1']])
set_param([sys,'/',['Charge pump',13,'PLL/Phase//freq',13,'detector/Constant1']],...
		'Value','[Thre; Thre]',...
		'position',[130,272,205,298])

add_block('built-in/Combinatorial Logic',[sys,'/',['Charge pump',13,'PLL/Phase//freq',13,'detector/phase',13,' flip-flop']])
set_param([sys,'/',['Charge pump',13,'PLL/Phase//freq',13,'detector/phase',13,' flip-flop']],...
		'Truth Table','lo',...
		'position',[440,125,495,165])

add_block('built-in/Mux',[sys,'/',['Charge pump',13,'PLL/Phase//freq',13,'detector/Mux2']])
set_param([sys,'/',['Charge pump',13,'PLL/Phase//freq',13,'detector/Mux2']],...
		'inputs','[2 2]',...
		'position',[390,125,420,165])

add_block('built-in/Memory',[sys,'/',['Charge pump',13,'PLL/Phase//freq',13,'detector/Memory']])
set_param([sys,'/',['Charge pump',13,'PLL/Phase//freq',13,'detector/Memory']],...
		'orientation',2,...
		'position',[425,205,465,235])

add_block('built-in/Outport',[sys,'/',['Charge pump',13,'PLL/Phase//freq',13,'detector/out_1']])
set_param([sys,'/',['Charge pump',13,'PLL/Phase//freq',13,'detector/out_1']],...
		'position',[940,145,960,165])

add_block('built-in/Gain',[sys,'/',['Charge pump',13,'PLL/Phase//freq',13,'detector/Gain1']])
set_param([sys,'/',['Charge pump',13,'PLL/Phase//freq',13,'detector/Gain1']],...
		'Gain','Cu',...
		'position',[900,145,920,165])