com_cdmo.m

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

function [ret,x0,str,ts,xts]=com_cdmo(t,x,u,flag);
%COM_CDMO	is the M-file description of the SIMULINK system named COM_CDMO.
%	The block-diagram can be displayed by typing: COM_CDMO.
%
%	SYS=COM_CDMO(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 COM_CDMO 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 COM_CDMO with a FLAG of zero:
%	[SIZES]=COM_CDMO([],[],[],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',[82,154,523,684])
     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', '10')
set_param(sys,'Relative error','1e-3')
set_param(sys,'Return vars',   '')


%     Subsystem  'Demonstration blocks'.

new_system([sys,'/','Demonstration blocks'])
set_param([sys,'/','Demonstration blocks'],'Location',[175,552,335,697])
set_param([sys,'/','Demonstration blocks'],...
		'Mask Display','.')


%     Finished composite block 'Demonstration blocks'.

set_param([sys,'/','Demonstration blocks'],...
		'orientation',1,...
		'BackGround',7,...
		'position',[50,498,105,521])


%     Subsystem  ['SSB-AM',13,'passband'].

new_system([sys,'/',['SSB-AM',13,'passband']])
set_param([sys,'/',['SSB-AM',13,'passband']],'Location',[76,558,518,713])


%     Subsystem  ['SSB-AM',13,'passband/QAM discr-filt'].

new_system([sys,'/',['SSB-AM',13,'passband/QAM discr-filt']])
set_param([sys,'/',['SSB-AM',13,'passband/QAM discr-filt']],'Location',[302,478,834,743])

add_block('built-in/Clock',[sys,'/',['SSB-AM',13,'passband/QAM discr-filt/Clock']])
set_param([sys,'/',['SSB-AM',13,'passband/QAM discr-filt/Clock']],...
		'position',[80,85,100,105])

add_block('built-in/Gain',[sys,'/',['SSB-AM',13,'passband/QAM discr-filt/Gain']])
set_param([sys,'/',['SSB-AM',13,'passband/QAM discr-filt/Gain']],...
		'Gain','Fc',...
		'position',[130,76,180,114])

add_block('built-in/Inport',[sys,'/',['SSB-AM',13,'passband/QAM discr-filt/in_2']])
set_param([sys,'/',['SSB-AM',13,'passband/QAM discr-filt/in_2']],...
		'Port','2',...
		'position',[160,145,180,165])

add_block('built-in/Inport',[sys,'/',['SSB-AM',13,'passband/QAM discr-filt/in_1']])
set_param([sys,'/',['SSB-AM',13,'passband/QAM discr-filt/in_1']],...
		'position',[160,30,180,50])

add_block('built-in/Outport',[sys,'/',['SSB-AM',13,'passband/QAM discr-filt/out_1']])
set_param([sys,'/',['SSB-AM',13,'passband/QAM discr-filt/out_1']],...
		'position',[490,90,510,110])

add_block('built-in/Sum',[sys,'/',['SSB-AM',13,'passband/QAM discr-filt/Sum2']])
set_param([sys,'/',['SSB-AM',13,'passband/QAM discr-filt/Sum2']],...
		'position',[430,90,450,110])

add_block('built-in/Product',[sys,'/',['SSB-AM',13,'passband/QAM discr-filt/Product1']])
set_param([sys,'/',['SSB-AM',13,'passband/QAM discr-filt/Product1']],...
		'position',[355,112,380,143])

add_block('built-in/Product',[sys,'/',['SSB-AM',13,'passband/QAM discr-filt/Product']])
set_param([sys,'/',['SSB-AM',13,'passband/QAM discr-filt/Product']],...
		'position',[355,47,380,78])

add_block('built-in/Note',[sys,'/',['SSB-AM',13,'passband/QAM discr-filt/Ref: Hambley: An Intro to Comm. Systems. pp72']])
set_param([sys,'/',['SSB-AM',13,'passband/QAM discr-filt/Ref: Hambley: An Intro to Comm. Systems. pp72']],...
		'position',[250,305,255,310])

add_block('built-in/Fcn',[sys,'/',['SSB-AM',13,'passband/QAM discr-filt/Fcn1']])
set_param([sys,'/',['SSB-AM',13,'passband/QAM discr-filt/Fcn1']],...
		'Expr','sin(u[1]+Ph)',...
		'position',[230,105,325,135])

add_block('built-in/Fcn',[sys,'/',['SSB-AM',13,'passband/QAM discr-filt/Fcn']])
set_param([sys,'/',['SSB-AM',13,'passband/QAM discr-filt/Fcn']],...
		'Expr','cos(u[1]+Ph)',...
		'position',[230,55,325,85])
add_line([sys,'/',['SSB-AM',13,'passband/QAM discr-filt']],[185,95;200,95;200,70;225,70])
add_line([sys,'/',['SSB-AM',13,'passband/QAM discr-filt']],[200,95;200,120;225,120])
add_line([sys,'/',['SSB-AM',13,'passband/QAM discr-filt']],[105,95;125,95])
add_line([sys,'/',['SSB-AM',13,'passband/QAM discr-filt']],[385,130;400,130;400,105;425,105])
add_line([sys,'/',['SSB-AM',13,'passband/QAM discr-filt']],[385,65;400,65;400,95;425,95])
add_line([sys,'/',['SSB-AM',13,'passband/QAM discr-filt']],[185,155;330,155;330,135;350,135])
add_line([sys,'/',['SSB-AM',13,'passband/QAM discr-filt']],[330,120;350,120])
add_line([sys,'/',['SSB-AM',13,'passband/QAM discr-filt']],[185,40;335,40;335,55;350,55])
add_line([sys,'/',['SSB-AM',13,'passband/QAM discr-filt']],[330,70;350,70])
add_line([sys,'/',['SSB-AM',13,'passband/QAM discr-filt']],[455,100;485,100])
set_param([sys,'/',['SSB-AM',13,'passband/QAM discr-filt']],...
		'Mask Display','QM DSB\nAM',...
		'Mask Type','Quadrature Multiplexed DSB-SC AM')
set_param([sys,'/',['SSB-AM',13,'passband/QAM discr-filt']],...
		'Mask Dialogue','Quadrature multiplexed double-sideband suppressed-carrier amplitude modulation:|Amplitude scale factor:|Carrier frequency (Hz):|Initial phase (rad):|Sample time:')
set_param([sys,'/',['SSB-AM',13,'passband/QAM discr-filt']],...
		'Mask Translate','Am=@1;Fc=@2*2*pi;Ph=@3;ts=@4;')
set_param([sys,'/',['SSB-AM',13,'passband/QAM discr-filt']],...
		'Mask Help','The Quadratic Multiplexed Double-Sideband-Supressed-Carrier AM block causes the amplitude of the sinusoidal carrier to depend on the input message signal, resulting in amplitude modulation. There are two input signal being modulated at the same time by using two orthogonal sinusoidal carriers.')
set_param([sys,'/',['SSB-AM',13,'passband/QAM discr-filt']],...
		'Mask Entries','1\/Fc\/Ph\/ts\/')


%     Finished composite block ['SSB-AM',13,'passband/QAM discr-filt'].

set_param([sys,'/',['SSB-AM',13,'passband/QAM discr-filt']],...
		'ForeGround',3,...
		'position',[280,45,360,85])

add_block('built-in/Outport',[sys,'/',['SSB-AM',13,'passband/out_1']])
set_param([sys,'/',['SSB-AM',13,'passband/out_1']],...
		'position',[400,55,420,75])

add_block('built-in/Note',[sys,'/',['SSB-AM',13,'passband/Ref: Hambley: An Intro to Comm. Systems. pp82']])
set_param([sys,'/',['SSB-AM',13,'passband/Ref: Hambley: An Intro to Comm. Systems. pp82']],...
		'position',[355,290,360,295])

add_block('built-in/Inport',[sys,'/',['SSB-AM',13,'passband/in_1']])
set_param([sys,'/',['SSB-AM',13,'passband/in_1']],...
		'position',[30,20,50,40])

add_block('built-in/Filter',[sys,'/',['SSB-AM',13,'passband/Filter']])
set_param([sys,'/',['SSB-AM',13,'passband/Filter']],...
		'Numerator','num',...
		'Denominator','den',...
		'Sample time','ts',...
		'position',[95,66,150,104])

add_block('built-in/Gain',[sys,'/',['SSB-AM',13,'passband/Gain']])
set_param([sys,'/',['SSB-AM',13,'passband/Gain']],...
		'Gain','k',...
		'position',[185,72,215,98])
add_line([sys,'/',['SSB-AM',13,'passband']],[365,65;395,65])
add_line([sys,'/',['SSB-AM',13,'passband']],[55,30;255,30;255,55;275,55])
add_line([sys,'/',['SSB-AM',13,'passband']],[55,30;65,30;65,85;90,85])
add_line([sys,'/',['SSB-AM',13,'passband']],[155,85;180,85])
add_line([sys,'/',['SSB-AM',13,'passband']],[220,85;240,85;240,75;275,75])
set_param([sys,'/',['SSB-AM',13,'passband']],...
		'Mask Display','plot(0,0,100,100,x,y);SSB AM',...
		'Mask Type','Passband SSB-AM')
set_param([sys,'/',['SSB-AM',13,'passband']],...
		'Mask Dialogue','Modulate the input signal using \nsingle-sideband amplitude modulation\nmethod with Hilbert transform filter.|Carrier frequency (Hz):|Initial Phase (rad):|Bandwidth of the input signal (Hz):|Time delay for Hilbert transform filter (sec):|Sample time (sec):|"upper" sideband or "lower" sideband:')
set_param([sys,'/',['SSB-AM',13,'passband']],...
		'Mask Translate','Fc=@1;Ph=@2;ban=@3;dly=@4;ts=@5;sd=@6;if ban>1/ts, ban=1/ts-ts;end;[num,den]=hilbiir(ts,dly,ban);pi2=2*pi;pih=pi/2;k=-1;if sd(1)==''l'',k=1;end;[x,y]=moduicon(1);')
set_param([sys,'/',['SSB-AM',13,'passband']],...
		'Mask Help','This block modulates the input signal using the SSB-AM method. The block uses Hilbert transfer filter. Hilbert transform filter is a non-causal filter. A positive time delay is requested in the implementation. You can specify the modulation to be an USB or LSB modulation.')
set_param([sys,'/',['SSB-AM',13,'passband']],...
		'Mask Entries','100/2/pi\/0\/1/.01/10\/.01*10\/.01\/''upper''\/')


%     Finished composite block ['SSB-AM',13,'passband'].

set_param([sys,'/',['SSB-AM',13,'passband']],...
		'position',[50,348,130,392])


%     Subsystem  ['FDM',13,'passband'].

new_system([sys,'/',['FDM',13,'passband']])
set_param([sys,'/',['FDM',13,'passband']],'Location',[30,490,489,666])

add_block('built-in/Gain',[sys,'/',['FDM',13,'passband/Sensitivity']])
set_param([sys,'/',['FDM',13,'passband/Sensitivity']],...
		'Gain','2',...
		'position',[75,25,115,55])

add_block('built-in/Product',[sys,'/',['FDM',13,'passband/Multiplier',13,'phase detector']])
set_param([sys,'/',['FDM',13,'passband/Multiplier',13,'phase detector']],...
		'position',[190,35,215,55])

add_block('built-in/Filter',[sys,'/',['FDM',13,'passband/Filter']])
set_param([sys,'/',['FDM',13,'passband/Filter']],...
		'Numerator','num',...
		'Denominator','den',...
		'Sample time','ts',...
		'position',[275,26,330,64])


%     Subsystem  ['FDM',13,'passband/Discrete-time',13,'VCO'].

new_system([sys,'/',['FDM',13,'passband/Discrete-time',13,'VCO']])
set_param([sys,'/',['FDM',13,'passband/Discrete-time',13,'VCO']],'Location',[93,129,873,276])

add_block('built-in/Outport',[sys,'/',['FDM',13,'passband/Discrete-time',13,'VCO/out_1']])
set_param([sys,'/',['FDM',13,'passband/Discrete-time',13,'VCO/out_1']],...
		'position',[690,35,710,55])

add_block('built-in/Gain',[sys,'/',['FDM',13,'passband/Discrete-time',13,'VCO/Sensitivity3']])
set_param([sys,'/',['FDM',13,'passband/Discrete-time',13,'VCO/Sensitivity3']],...
		'Gain','Ac',...
		'position',[605,30,655,60])

add_block('built-in/Fcn',[sys,'/',['FDM',13,'passband/Discrete-time',13,'VCO/sin']])
set_param([sys,'/',['FDM',13,'passband/Discrete-time',13,'VCO/sin']],...
		'Expr','cos(u[1])',...
		'position',[510,32,585,58])

add_block('built-in/Constant',[sys,'/',['FDM',13,'passband/Discrete-time',13,'VCO/Phase']])
set_param([sys,'/',['FDM',13,'passband/Discrete-time',13,'VCO/Phase']],...
		'Value','Ph',...
		'position',[400,100,420,120])

add_block('built-in/Sum',[sys,'/',['FDM',13,'passband/Discrete-time',13,'VCO/Sum']])
set_param([sys,'/',['FDM',13,'passband/Discrete-time',13,'VCO/Sum']],...
		'inputs','+++',...
		'position',[465,27,485,63])

add_block('built-in/Gain',[sys,'/',['FDM',13,'passband/Discrete-time',13,'VCO/Sensitivity2']])
set_param([sys,'/',['FDM',13,'passband/Discrete-time',13,'VCO/Sensitivity2']],...
		'Gain','Fc',...
		'position',[255,75,300,105])