com_ddmo.m

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

function [ret,x0,str,ts,xts]=com_ddmo(t,x,u,flag);
%COM_DDMO	is the M-file description of the SIMULINK system named COM_DDMO.
%	The block-diagram can be displayed by typing: COM_DDMO.
%
%	SYS=COM_DDMO(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_DDMO 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_DDMO with a FLAG of zero:
%	[SIZES]=COM_DDMO([],[],[],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',[178,178,705,718])
     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',   '')

add_block('built-in/Note',[sys,'/','(passband simulation)'])
set_param([sys,'/','(passband simulation)'],...
		'ForeGround',2,...
		'position',[240,20,245,25])


%     Subsystem  'MASK mod'.

new_system([sys,'/','MASK mod'])
set_param([sys,'/','MASK mod'],'Location',[652,470,1010,578])

add_block('built-in/Outport',[sys,'/','MASK mod/out_1'])
set_param([sys,'/','MASK mod/out_1'],...
		'position',[300,25,320,45])

add_block('built-in/Inport',[sys,'/','MASK mod/in_1'])
set_param([sys,'/','MASK mod/in_1'],...
		'position',[30,25,50,45])


%     Subsystem  'MASK mod/MASK map1'.

new_system([sys,'/','MASK mod/MASK map1'])
set_param([sys,'/','MASK mod/MASK map1'],'Location',[520,162,902,335])

add_block('built-in/Outport',[sys,'/','MASK mod/MASK map1/out_1'])
set_param([sys,'/','MASK mod/MASK map1/out_1'],...
		'position',[315,65,335,85])

add_block('built-in/Zero-Order Hold',[sys,'/',['MASK mod/MASK map1/Zero-Order',13,'Hold']])
set_param([sys,'/',['MASK mod/MASK map1/Zero-Order',13,'Hold']],...
		'Sample time','td',...
		'position',[100,24,135,56])

add_block('built-in/Inport',[sys,'/','MASK mod/MASK map1/in_1'])
set_param([sys,'/','MASK mod/MASK map1/in_1'],...
		'position',[40,30,60,50])

add_block('built-in/Gain',[sys,'/','MASK mod/MASK map1/Gain'])
set_param([sys,'/','MASK mod/MASK map1/Gain'],...
		'Gain','Kc',...
		'position',[180,25,220,55])

add_block('built-in/Sum',[sys,'/','MASK mod/MASK map1/Sum1'])
set_param([sys,'/','MASK mod/MASK map1/Sum1'],...
		'inputs','+-',...
		'position',[250,63,270,87])

add_block('built-in/Constant',[sys,'/','MASK mod/MASK map1/Gain diff'])
set_param([sys,'/','MASK mod/MASK map1/Gain diff'],...
		'Value','Ac',...
		'position',[180,96,215,124])
add_line([sys,'/','MASK mod/MASK map1'],[225,40;230,40;230,70;245,70])
add_line([sys,'/','MASK mod/MASK map1'],[220,110;230,110;230,80;245,80])
add_line([sys,'/','MASK mod/MASK map1'],[275,75;310,75])
add_line([sys,'/','MASK mod/MASK map1'],[65,40;95,40])
add_line([sys,'/','MASK mod/MASK map1'],[140,40;175,40])
set_param([sys,'/','MASK mod/MASK map1'],...
		'Mask Display','MASK\nmap',...
		'Mask Type','MASK map')
set_param([sys,'/','MASK mod/MASK map1'],...
		'Mask Dialogue','Map a digital symbol onto an analog\nsignal in the range [-1, 1] using M-ary\namplitude shift keying modulation mapping.|M-ary number (input integers in range [0, M-1]):|Input symbol interval and offset (sec):')
set_param([sys,'/','MASK mod/MASK map1'],...
		'Mask Translate','N=@1;Ac=1;Kc=Ac*2/(N-1);td=@2;')
set_param([sys,'/','MASK mod/MASK map1'],...
		'Mask Help','This block maps an input digit signal onto a signal in the range [-1, 1]. The coded signal can be modulated using analog\nmodulation blocks. Input digits are in the range [0, M-1], where M is the M-ary number. The input symbol interval can be a two-element vector with the second element being offset. The default value offset is 0. Use the MODMAP function to see the constellation.')
set_param([sys,'/','MASK mod/MASK map1'],...
		'Mask Entries','M\/td\/')


%     Finished composite block 'MASK mod/MASK map1'.

set_param([sys,'/','MASK mod/MASK map1'],...
		'position',[80,13,160,57])


%     Subsystem  ['MASK mod/DSB-SC AM',13,'passband'].

new_system([sys,'/',['MASK mod/DSB-SC AM',13,'passband']])
set_param([sys,'/',['MASK mod/DSB-SC AM',13,'passband']],'Location',[103,459,555,641])

add_block('built-in/Zero-Order Hold',[sys,'/',['MASK mod/DSB-SC AM',13,'passband/Zero-Order',13,'Hold',13,'timing control']])
set_param([sys,'/',['MASK mod/DSB-SC AM',13,'passband/Zero-Order',13,'Hold',13,'timing control']],...
		'Sample time','ts',...
		'position',[50,219,85,251])

add_block('built-in/Gain',[sys,'/',['MASK mod/DSB-SC AM',13,'passband/Gain']])
set_param([sys,'/',['MASK mod/DSB-SC AM',13,'passband/Gain']],...
		'Gain','Fc',...
		'position',[70,56,120,94])

add_block('built-in/Clock',[sys,'/',['MASK mod/DSB-SC AM',13,'passband/Clock']])
set_param([sys,'/',['MASK mod/DSB-SC AM',13,'passband/Clock']],...
		'position',[20,65,40,85])

add_block('built-in/Product',[sys,'/',['MASK mod/DSB-SC AM',13,'passband/Product']])
set_param([sys,'/',['MASK mod/DSB-SC AM',13,'passband/Product']],...
		'position',[340,42,365,73])

add_block('built-in/Outport',[sys,'/',['MASK mod/DSB-SC AM',13,'passband/out_1']])
set_param([sys,'/',['MASK mod/DSB-SC AM',13,'passband/out_1']],...
		'position',[380,50,400,70])

add_block('built-in/Inport',[sys,'/',['MASK mod/DSB-SC AM',13,'passband/in_1']])
set_param([sys,'/',['MASK mod/DSB-SC AM',13,'passband/in_1']],...
		'position',[205,20,225,40])

add_block('built-in/Constant',[sys,'/',['MASK mod/DSB-SC AM',13,'passband/Constant1']])
set_param([sys,'/',['MASK mod/DSB-SC AM',13,'passband/Constant1']],...
		'Value','Ph',...
		'position',[60,115,100,145])

add_block('built-in/Fcn',[sys,'/',['MASK mod/DSB-SC AM',13,'passband/Fcn']])
set_param([sys,'/',['MASK mod/DSB-SC AM',13,'passband/Fcn']],...
		'Expr','cos(u[1])',...
		'position',[205,85,265,115])

add_block('built-in/Sum',[sys,'/',['MASK mod/DSB-SC AM',13,'passband/Sum2']])
set_param([sys,'/',['MASK mod/DSB-SC AM',13,'passband/Sum2']],...
		'position',[170,90,190,110])
add_line([sys,'/',['MASK mod/DSB-SC AM',13,'passband']],[270,100;270,65;335,65])
add_line([sys,'/',['MASK mod/DSB-SC AM',13,'passband']],[370,60;375,60])
add_line([sys,'/',['MASK mod/DSB-SC AM',13,'passband']],[230,30;270,30;270,50;335,50])
add_line([sys,'/',['MASK mod/DSB-SC AM',13,'passband']],[45,75;65,75])
add_line([sys,'/',['MASK mod/DSB-SC AM',13,'passband']],[125,75;150,75;150,95;165,95])
add_line([sys,'/',['MASK mod/DSB-SC AM',13,'passband']],[105,130;130,130;130,105;165,105])
add_line([sys,'/',['MASK mod/DSB-SC AM',13,'passband']],[195,100;200,100])
set_param([sys,'/',['MASK mod/DSB-SC AM',13,'passband']],...
		'Mask Display','plot(0,0,100,100,x,y);DSB AM',...
		'Mask Type','Passband DSB-SC AM')
set_param([sys,'/',['MASK mod/DSB-SC AM',13,'passband']],...
		'Mask Dialogue','Modulate the input signal using\nDSB-SC AM method.|Carrier frequency (Hz):|Initial Phase (rad):|Sample time (sec):')
set_param([sys,'/',['MASK mod/DSB-SC AM',13,'passband']],...
		'Mask Translate','Fc=@1;Fc=Fc*2*pi;Ph=@2;ts=@3;[x,y]=moduicon(1);')
set_param([sys,'/',['MASK mod/DSB-SC AM',13,'passband']],...
		'Mask Help','This block modulates the input signal using the double-sideband suppressed-carrier amplitude modulation. In general, the carrier frequency should be much higher than the message signal frequency.')
set_param([sys,'/',['MASK mod/DSB-SC AM',13,'passband']],...
		'Mask Entries','Fc\/Ph\/ts\/')


%     Finished composite block ['MASK mod/DSB-SC AM',13,'passband'].

set_param([sys,'/',['MASK mod/DSB-SC AM',13,'passband']],...
		'ForeGround',4,...
		'position',[195,13,275,57])
add_line([sys,'/','MASK mod'],[280,35;295,35])
add_line([sys,'/','MASK mod'],[55,35;75,35])
add_line([sys,'/','MASK mod'],[165,35;190,35])
set_param([sys,'/','MASK mod'],...
		'Mask Display','plot(0,0,100,100,x,y);MASK',...
		'Mask Type','Passband MASK modulation')
set_param([sys,'/','MASK mod'],...
		'Mask Dialogue','Modulate the input digital signal using M-ary amplitude shift keying method.|M-ary number (input integer in range [0, M-1]):|Input symbol interval (sec):|Carrier frequency (Hz):|Carrier initial phase (rad):|Sample time (sec):')
set_param([sys,'/','MASK mod'],...
		'Mask Translate','M=@1;td=@2;Fc=@3;Ph=@4;ts=@5;[x,y]=moduicon(5);')
set_param([sys,'/','MASK mod'],...
		'Mask Help','This block modulates a digital signal using MASK method. The inputs are integers in the range [0, M-1], where M is the M-ary number. When the symbol interval is a two-element vector, the second element is the offset (default 0). The block keeps the input value at the offset time point for the entire interval.')
set_param([sys,'/','MASK mod'],...
		'Mask Entries','8\/.1\/100/2/pi\/0\/.01\/')


%     Finished composite block 'MASK mod'.

set_param([sys,'/','MASK mod'],...
		'position',[40,48,120,92])


%     Subsystem  'MASK demod'.

new_system([sys,'/','MASK demod'])
set_param([sys,'/','MASK demod'],'Location',[482,504,899,622])

add_block('built-in/Inport',[sys,'/','MASK demod/in_1'])
set_param([sys,'/','MASK demod/in_1'],...
		'position',[30,30,50,50])

add_block('built-in/Outport',[sys,'/','MASK demod/out_1'])
set_param([sys,'/','MASK demod/out_1'],...
		'position',[360,30,380,50])


%     Subsystem  ['MASK demod/DSB-SC-ADM',13,'passband'].

new_system([sys,'/',['MASK demod/DSB-SC-ADM',13,'passband']])
set_param([sys,'/',['MASK demod/DSB-SC-ADM',13,'passband']],'Location',[588,690,716,690])


%     Subsystem  ['MASK demod/DSB-SC-ADM',13,'passband/Dis module',13,'integrator'].

new_system([sys,'/',['MASK demod/DSB-SC-ADM',13,'passband/Dis module',13,'integrator']])
set_param([sys,'/',['MASK demod/DSB-SC-ADM',13,'passband/Dis module',13,'integrator']],'Location',[40,403,515,507])

add_block('built-in/Outport',[sys,'/',['MASK demod/DSB-SC-ADM',13,'passband/Dis module',13,'integrator/out_1']])
set_param([sys,'/',['MASK demod/DSB-SC-ADM',13,'passband/Dis module',13,'integrator/out_1']],...
		'position',[445,30,465,50])

add_block('built-in/Unit Delay',[sys,'/',['MASK demod/DSB-SC-ADM',13,'passband/Dis module',13,'integrator/Unit Delay']])
set_param([sys,'/',['MASK demod/DSB-SC-ADM',13,'passband/Dis module',13,'integrator/Unit Delay']],...
		'Sample time','ts',...
		'x0','init',...
		'position',[345,30,395,50])