sysm4.txt

与上一个上传源码属同一问题的不同方面

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

add_block('built-in/Note',[sys,'/','同步解调'])
set_param([sys,'/','同步解调'],...
		'Font Name','宋体',...
		'Font Size',12,...
		'position',[276,311,281,316])

add_block('built-in/Note',[sys,'/','导频信号'])
set_param([sys,'/','导频信号'],...
		'Font Name','宋体',...
		'Font Size',12,...
		'position',[216,276,221,281])

add_block('built-in/Note',[sys,'/','上边带 SSB'])
set_param([sys,'/','上边带 SSB'],...
		'Font Name','宋体',...
		'Font Size',12,...
		'position',[201,246,206,251])

add_block('built-in/Note',[sys,'/','发送端'])
set_param([sys,'/','发送端'],...
		'Font Name','宋体',...
		'Font Size',12,...
		'position',[151,201,156,206])

add_block('built-in/Note',[sys,'/','接收端'])
set_param([sys,'/','接收端'],...
		'Font Name','宋体',...
		'Font Size',12,...
		'position',[146,281,151,286])

add_block('built-in/Scope',[sys,'/','解调信号'])
set_param([sys,'/','解调信号'],...
		'Font Name','宋体',...
		'Font Size',12,...
		'Vgain','1.000000',...
		'Hgain','1.000000',...
		'Vmax','2.000000',...
		'Hmax','2.000000',...
		'Window',[157,149,488,463],...
		'position',[465,260,495,290])

add_block('built-in/Product',[sys,'/','乘法器1'])
set_param([sys,'/','乘法器1'],...
		'Font Name','宋体',...
		'Font Size',12,...
		'Mask Display','plot([0 1],[1 0],[0 1],[0 1]);',...
		'Mask Type','乘法器',...
		'Mask Dialogue','乘法器',...
		'position',[260,264,290,286])

add_block('built-in/Zero-Pole',[sys,'/',['Butterworth',13,'LP 滤波器']])
set_param([sys,'/',['Butterworth',13,'LP 滤波器']],...
		'Font Name','宋体',...
		'Font Size',12,...
		'Zeros','z',...
		'Poles','p',...
		'Gain','k',...
		'Mask Display','plot(20*log10(w),20*log10(abs(h)))',...
		'Mask Type','buttap')
set_param([sys,'/',['Butterworth',13,'LP 滤波器']],...
		'Mask Dialogue','Analog Butterworth low pass filter.\n\n(Requires Signal Processing Toolbox)|Cutoff frequency (rad/sec)|Order')
set_param([sys,'/',['Butterworth',13,'LP 滤波器']],...
		'Mask Translate','[z,p,k]=buttap(@2); z=@1*z; p=@1*p; k=real(prod(-p)); [num,den]=zp2tf(z,p,k); w=logspace(log10(@1/10),log10(10*@1),20); h=freqs(num,den,w); ')
set_param([sys,'/',['Butterworth',13,'LP 滤波器']],...
		'Mask Help','Butterworth low pass filter using MATLAB''s butter command.',...
		'Mask Entries','15\/4\/',...
		'position',[345,253,405,297])

add_block('built-in/Sum',[sys,'/','加法器'])
set_param([sys,'/','加法器'],...
		'Font Name','宋体',...
		'Font Size',12,...
		'position',[365,165,385,185])

add_block('built-in/Sum',[sys,'/','减法器'])
set_param([sys,'/','减法器'],...
		'Font Name','宋体',...
		'Font Size',12,...
		'inputs','+-',...
		'position',[365,75,385,95])

add_block('built-in/Product',[sys,'/','乘法器'])
set_param([sys,'/','乘法器'],...
		'ForeGround',3,...
		'Font Name','宋体',...
		'Font Size',12,...
		'Mask Display','plot([0 1],[1 0],[0 1],[0 1]);',...
		'Mask Type','乘法器',...
		'Mask Dialogue','乘法器',...
		'position',[285,69,315,91])

add_block('built-in/Product',[sys,'/','乘法器2'])
set_param([sys,'/','乘法器2'],...
		'ForeGround',3,...
		'Font Name','宋体',...
		'Font Size',12,...
		'Mask Display','plot([0 1],[1 0],[0 1],[0 1]);',...
		'Mask Type','乘法器',...
		'Mask Dialogue','乘法器',...
		'position',[285,169,315,191])

add_block('built-in/Transport Delay',[sys,'/',['移',13,'相']])
set_param([sys,'/',['移',13,'相']],...
		'orientation',1,...
		'Font Name','宋体',...
		'Font Size',12,...
		'Delay Time','0.25/(W/(2*pi))',...
		'Mask Display','90度',...
		'Mask Type','90度移相')
set_param([sys,'/',['移',13,'相']],...
		'Mask Dialogue','90度移相|输入信号的角频率（Rad/s）：',...
		'Mask Translate','W=@1;',...
		'Mask Entries','10\/',...
		'position',[217,120,263,150])

add_block('built-in/Transport Delay',[sys,'/',['移',13,'相',13,'1']])
set_param([sys,'/',['移',13,'相',13,'1']],...
		'orientation',1,...
		'move name',0,...
		'Font Name','宋体',...
		'Font Size',12,...
		'Delay Time','0.25/(W/(2*pi))',...
		'Mask Display','90度',...
		'Mask Type','90度移相')
set_param([sys,'/',['移',13,'相',13,'1']],...
		'Mask Dialogue','90度移相|输入信号的角频率（Rad/s）：',...
		'Mask Translate','W=@1;',...
		'Mask Entries','100\/',...
		'position',[162,120,208,150])

add_block('built-in/Scope',[sys,'/','下边带信号'])
set_param([sys,'/','下边带信号'],...
		'ForeGround',3,...
		'Font Name','宋体',...
		'Font Size',12,...
		'Vgain','2.000000',...
		'Hgain','1.000000',...
		'Vmax','4.000000',...
		'Hmax','2.000000',...
		'Window',[69,158,400,472])
set_param([sys,'/','下边带信号'],...
		'position',[465,160,495,190])

add_block('built-in/Scope',[sys,'/','上边带信号'])
set_param([sys,'/','上边带信号'],...
		'ForeGround',3,...
		'Font Name','宋体',...
		'Font Size',12,...
		'Vgain','2.000000',...
		'Hgain','1.000000',...
		'Vmax','4.000000',...
		'Hmax','2.000000',...
		'Window',[72,40,403,354])
set_param([sys,'/','上边带信号'],...
		'position',[465,70,495,100])

add_block('built-in/Note',[sys,'/',['cos ( Wc * t + p )',13,'Wc 载波角频率',13,'p 载波初始相位']])
set_param([sys,'/',['cos ( Wc * t + p )',13,'Wc 载波角频率',13,'p 载波初始相位']],...
		'ForeGround',3,...
		'Font Name','宋体',...
		'Font Size',12,...
		'position',[76,156,81,161])

add_block('built-in/Note',[sys,'/','已调 SSB 信号'])
set_param([sys,'/','已调 SSB 信号'],...
		'ForeGround',3,...
		'Font Name','宋体',...
		'Font Size',12,...
		'position',[416,41,421,46])

add_block('built-in/Note',[sys,'/','X(t)'])
set_param([sys,'/','X(t)'],...
		'ForeGround',3,...
		'Font Name','宋体',...
		'Font Size',12,...
		'position',[141,51,146,56])

add_block('built-in/Note',[sys,'/','调制信号'])
set_param([sys,'/','调制信号'],...
		'ForeGround',3,...
		'Font Name','宋体',...
		'Font Size',12,...
		'position',[81,31,86,36])

add_block('built-in/Sine Wave',[sys,'/','载波'])
set_param([sys,'/','载波'],...
		'ForeGround',3,...
		'Font Name','宋体',...
		'Font Size',12,...
		'frequency','100',...
		'position',[70,125,90,145])

add_block('built-in/Signal Generator',[sys,'/',['信号发生器',13,'频率（弧度//秒）']])
set_param([sys,'/',['信号发生器',13,'频率（弧度//秒）']],...
		'ForeGround',3,...
		'Font Name','宋体',...
		'Font Size',12,...
		'Peak','1.000000',...
		'Peak Range','5.000000',...
		'Freq','10.000000',...
		'Freq Range','50.000000')
set_param([sys,'/',['信号发生器',13,'频率（弧度//秒）']],...
		'Wave','Sin',...
		'Units','Rads',...
		'position',[60,58,105,92])

add_block('built-in/Note',[sys,'/',['对 SSB 信号，无法',13,'直接从其中提取参考',13,'信号，故发送端在发',13,'射 SSB 信号的同时，',13,'还须发射受到一定程',13,'度抑制的载波信号。',13,'（称为导频信号）']])
set_param([sys,'/',['对 SSB 信号，无法',13,'直接从其中提取参考',13,'信号，故发送端在发',13,'射 SSB 信号的同时，',13,'还须发射受到一定程',13,'度抑制的载波信号。',13,'（称为导频信号）']],...
		'ForeGround',3,...
		'Font Name','宋体')
set_param([sys,'/',['对 SSB 信号，无法',13,'直接从其中提取参考',13,'信号，故发送端在发',13,'射 SSB 信号的同时，',13,'还须发射受到一定程',13,'度抑制的载波信号。',13,'（称为导频信号）']],...
		'Font Size',12,...
		'position',[66,216,71,221])

add_block('built-in/Note',[sys,'/','同步调制//解调 AM 系统 （ 单边带 SSB ）—— 移相法'])
set_param([sys,'/','同步调制//解调 AM 系统 （ 单边带 SSB ）—— 移相法'],...
		'Font Name','宋体',...
		'Font Size',12,...
		'position',[266,11,271,16])
add_line(sys,[295,275;340,275])
add_line(sys,[410,275;460,275])
add_line(sys,[110,75;240,75;240,115])
add_line(sys,[95,135;130,135;130,85;280,85])
add_line(sys,[240,75;280,75])
add_line(sys,[185,155;185,185;280,185])
add_line(sys,[240,155;240,175;280,175])
add_line(sys,[185,85;185,115])
add_line(sys,[320,80;360,80])
add_line(sys,[320,180;360,180])
add_line(sys,[390,85;460,85])
add_line(sys,[390,175;460,175])
add_line(sys,[335,80;335,170;360,170])
add_line(sys,[345,180;345,90;360,90])
add_line(sys,[425,85;425,225;155,225;155,270;255,270])
add_line(sys,[130,135;130,280;255,280])

drawnow

% Return any arguments.
if (nargin | nargout)
	% Must use feval here to access system in memory
	if (nargin > 3)
		if (flag == 0)
			eval(['[ret,x0,str,ts,xts]=',sys,'(t,x,u,flag);'])
		else
			eval(['ret =', sys,'(t,x,u,flag);'])
		end
	else
		[ret,x0,str,ts,xts] = feval(sys);
	end
else
	drawnow % Flash up the model and execute load callback
end