yjwx.m

信号实验常用的simulink模型和mfile,可直接在matlan下运行。

function [ret,x0,str,ts,xts]=yjwx(t,x,u,flag);
%YJWX	is the M-file description of the SIMULINK system named YJWX.
%	The block-diagram can be displayed by typing: YJWX.
%
%	SYS=YJWX(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 YJWX 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 YJWX with a FLAG of zero:
%	[SIZES]=YJWX([],[],[],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',[21,60,469,432])
     open_system(sys)
end;
set_param(sys,'algorithm',     'RK-45')
set_param(sys,'Start time',    '0.0')
set_param(sys,'Stop time',     '100')
set_param(sys,'Min step size', '0.0001')
set_param(sys,'Max step size', '0.1')
set_param(sys,'Relative error','1e-3')
set_param(sys,'Return vars',   '')

add_block('built-in/Step Fcn',[sys,'/','阶跃输入'])
set_param([sys,'/','阶跃输入'],...
		'Font Name','宋体',...
		'position',[35,250,70,290])


%     Subsystem  '方波信号'.

new_system([sys,'/','方波信号'])
set_param([sys,'/','方波信号'],'Location',[158,441,759,682])

add_block('built-in/Outport',[sys,'/','方波信号/out_1'])
set_param([sys,'/','方波信号/out_1'],...
		'hide name',0,...
		'position',[560,105,580,125])

add_block('built-in/Gain',[sys,'/','方波信号/Gain'])
set_param([sys,'/','方波信号/Gain'],...
		'hide name',0,...
		'Gain','ht',...
		'position',[510,102,535,128])

add_block('built-in/Product',[sys,'/','方波信号/Product'])
set_param([sys,'/','方波信号/Product'],...
		'hide name',0,...
		'position',[455,105,480,125])

add_block('built-in/Logical Operator',[sys,'/',['方波信号/Logical',13,'Operator1']])
set_param([sys,'/',['方波信号/Logical',13,'Operator1']],...
		'Operator','NOT',...
		'Number of Input Ports','1',...
		'position',[355,38,385,62])

add_block('built-in/Logical Operator',[sys,'/',['方波信号/Logical',13,'Operator']])
set_param([sys,'/',['方波信号/Logical',13,'Operator']],...
		'Operator','XOR',...
		'position',[285,38,315,62])

add_block('built-in/Constant',[sys,'/','方波信号/Constant1'])
set_param([sys,'/','方波信号/Constant1'],...
		'Value','stt',...
		'position',[295,149,320,171])

add_block('built-in/Clock',[sys,'/','方波信号/Clock1'])
set_param([sys,'/','方波信号/Clock1'],...
		'hide name',0,...
		'position',[295,100,315,120])

add_block('built-in/Relational Operator',[sys,'/',['方波信号/Relational',13,'operator']])
set_param([sys,'/',['方波信号/Relational',13,'operator']],...
		'hide name',0,...
		'position',[370,102,400,133])

add_block('built-in/Sum',[sys,'/','方波信号/Sum1'])
set_param([sys,'/','方波信号/Sum1'],...
		'hide name',0,...
		'inputs','+-',...
		'position',[130,120,150,140])

add_block('built-in/Sum',[sys,'/','方波信号/Sum'])
set_param([sys,'/','方波信号/Sum'],...
		'hide name',0,...
		'inputs','+-',...
		'position',[135,35,155,55])

add_block('built-in/Constant',[sys,'/','方波信号/Constant'])
set_param([sys,'/','方波信号/Constant'],...
		'position',[25,30,45,50])

add_block('built-in/Unit Delay',[sys,'/','方波信号/Unit Delay'])
set_param([sys,'/','方波信号/Unit Delay'],...
		'orientation',2,...
		'Sample time','[Ts,st1]',...
		'position',[145,75,195,95])

add_block('built-in/Unit Delay',[sys,'/','方波信号/Unit Delay1'])
set_param([sys,'/','方波信号/Unit Delay1'],...
		'orientation',2,...
		'Sample time','[Ts,st2]',...
		'x0','ini',...
		'position',[140,165,190,185])
add_line([sys,'/','方波信号'],[390,50;430,50;430,110;450,110])
add_line([sys,'/','方波信号'],[320,50;350,50])
add_line([sys,'/','方波信号'],[325,160;345,160;345,125;365,125])
add_line([sys,'/','方波信号'],[540,115;555,115])
add_line([sys,'/','方波信号'],[485,115;505,115])
add_line([sys,'/','方波信号'],[405,120;450,120])
add_line([sys,'/','方波信号'],[320,110;365,110])
add_line([sys,'/','方波信号'],[140,85;105,85;105,50;130,50])
add_line([sys,'/','方波信号'],[160,45;230,45;230,85;200,85])
add_line([sys,'/','方波信号'],[50,40;130,40])
add_line([sys,'/','方波信号'],[155,130;210,130;210,175;195,175])
add_line([sys,'/','方波信号'],[135,175;100,175;100,135;125,135])
add_line([sys,'/','方波信号'],[50,40;70,40;70,125;125,125])
add_line([sys,'/','方波信号'],[160,45;280,45])
add_line([sys,'/','方波信号'],[155,130;245,130;245,55;280,55])
set_param([sys,'/','方波信号'],...
		'Mask Display','plot(0,0,100,100,[90,75,75,60,60,35,35,20,20,10],[20,20,80,80,20,20,80,80,20,20])',...
		'Mask Type','Pulse Generator')
set_param([sys,'/','方波信号'],...
		'Mask Dialogue','Pulse Generator.|Pulse period (secs):|Pulse width:|Pulse height:|Pulse start time:')
set_param([sys,'/','方波信号'],...
		'Mask Translate','Ts=@1; du=@2; ht=@3; stt=@4; ini=ones(length(ht),1); st1=rem(stt,Ts); st2=rem(stt+du,Ts);if(Ts<=1.2*du),ini=zeros(length(ht),1);end;')
set_param([sys,'/','方波信号'],...
		'Mask Help','Pulse generator which ensures pulse\ntransitions are hit. Provides a vector of pulses when the height is entered as a vector.\nUnmask to see how it works.',...
		'Mask Entries','20\/10\/1\/0\/')


%     Finished composite block '方波信号'.

set_param([sys,'/','方波信号'],...
		'Font Name','宋体',...
		'position',[35,151,80,189])

add_block('built-in/Scope',[sys,'/','高通系统的响应'])
set_param([sys,'/','高通系统的响应'],...
		'Font Name','宋体',...
		'Vgain','1.500000',...
		'Hgain','20.000000',...
		'Vmax','3.000000',...
		'Hmax','40.000000',...
		'Window',[324,143,639,289])
open_system([sys,'/','高通系统的响应'])
set_param([sys,'/','高通系统的响应'],...
		'position',[385,148,415,192])

add_block('built-in/Scope',[sys,'/','低通系统的响应'])
set_param([sys,'/','低通系统的响应'],...
		'Font Name','宋体',...
		'Vgain','1.500000',...
		'Hgain','20.000000',...
		'Vmax','3.000000',...
		'Hmax','40.000000',...
		'Window',[323,290,639,451])
open_system([sys,'/','低通系统的响应'])
set_param([sys,'/','低通系统的响应'],...
		'position',[390,294,420,336])

add_block('built-in/Transfer Fcn',[sys,'/','一阶低通'])
set_param([sys,'/','一阶低通'],...
		'Font Name','宋体',...
		'Denominator','[R*C 1]')
set_param([sys,'/','一阶低通'],...
		'Mask Display','plot([0 2 2 2],[2 2 1.5 2.5],[2 6],[2.5 2.5],[2 6],[1.5 1.5],[6 6],[2.5 1.5],[6 11],[2 2],[8 8 6 10],[2 0 0 0],[6 10],[-1 -1],[8 8 7 9],[-1 -3 -3 -3]);',...
		'Mask Type','积分')
set_param([sys,'/','一阶低通'],...
		'Mask Dialogue','积分|电容量（法）：|电阻值（欧姆）：',...
		'Mask Translate','C=@1;R=@2;',...
		'Mask Entries','.1\/10\/',...
		'position',[185,290,255,340])


%     Subsystem  '频谱分析仪1'.

new_system([sys,'/','频谱分析仪1'])
set_param([sys,'/','频谱分析仪1'],'Location',[0,0,371,185])

add_block('built-in/Mux',[sys,'/','频谱分析仪1/Mux'])
set_param([sys,'/','频谱分析仪1/Mux'],...
		'inputs','2',...
		'position',[145,71,175,104])

add_block('built-in/Inport',[sys,'/','频谱分析仪1/system output'])
set_param([sys,'/','频谱分析仪1/system output'],...
		'Port','2',...
		'position',[55,110,75,130])

add_block('built-in/S-Function',[sys,'/',['频谱分析仪1/S-function block',13,'which calls M-file']])
set_param([sys,'/',['频谱分析仪1/S-function block',13,'which calls M-file']],...
		'function name','sfuntf',...
		'parameters','fftpts,npts,HowOften,offset,ts,0',...
		'position',[250,72,295,108])

add_block('built-in/Inport',[sys,'/','频谱分析仪1/system input'])
set_param([sys,'/','频谱分析仪1/system input'],...
		'position',[55,30,75,50])
add_line([sys,'/','频谱分析仪1'],[80,120;110,120;110,95;140,95])
add_line([sys,'/','频谱分析仪1'],[80,40;110,40;110,80;140,80])
add_line([sys,'/','频谱分析仪1'],[180,90;245,90])
set_param([sys,'/','频谱分析仪1'],...
		'Mask Display','plot(0,0,100,100,[11,93,93,11,11],[90,90,40,40,90],[90,80,72,64,55,45,41,36,20],[50,50,52,59,67,75,78,80,80])',...
		'Mask Type','Spectrum analyzer.')
set_param([sys,'/','频谱分析仪1'],...
		'Mask Dialogue','Spectrum analyser using graph window.\nHook input 1 to system input.\nHook input 2 to system output.|Length of buffer:|Number of points for fft:|Plot after how many points:|Sample time:')
set_param([sys,'/','频谱分析仪1'],...
		'Mask Translate','npts=@1; fftpts=@2; HowOften=@3; ts=@4(1); if length(@4) > 1, offset = @4(2), else, offset = 0; end')
set_param([sys,'/','频谱分析仪1'],...
		'Mask Help','Hook the system''s input to the first port of this block and the system''s output to second port. Results are displayed as the frequency response of this system. Use signals with high harmonic content for best results. ')
set_param([sys,'/','频谱分析仪1'],...
		'Mask Entries','512\/512\/64\/0.05\/')


%     Finished composite block '频谱分析仪1'.

set_param([sys,'/','频谱分析仪1'],...
		'Font Name','宋体',...
		'Drop Shadow',4,...
		'position',[390,219,420,266])

add_block('built-in/Transfer Fcn',[sys,'/','一阶高通'])
set_param([sys,'/','一阶高通'],...
		'Font Name','宋体',...
		'Numerator','[R*C 0]',...
		'Denominator','[R*C 1]')
set_param([sys,'/','一阶高通'],...
		'Mask Display','plot([0 2 2 2],[2 2 0 4],[3 3 3 8],[0 4 2 2],[6 6 5.5 6.5],[2 0 0 0],[5.5 5.5],[0 -4],[6.5 6.5],[0 -4],[5.5 6.5 6 6],[-4 -4 -4 -5.5],[5.5 6.5],[-5.5 -5.5]);',...
		'Mask Type','微分')
set_param([sys,'/','一阶高通'],...
		'Mask Dialogue','微分|电容量（法）：|电阻值（欧姆）：',...
		'Mask Translate','C=@1;R=@2;',...
		'Mask Entries','0.1\/10\/',...
		'position',[185,142,255,198])


%     Subsystem  '频谱分析仪'.

new_system([sys,'/','频谱分析仪'])
set_param([sys,'/','频谱分析仪'],'Location',[0,0,371,185])

add_block('built-in/Mux',[sys,'/','频谱分析仪/Mux'])
set_param([sys,'/','频谱分析仪/Mux'],...
		'inputs','2',...
		'position',[145,71,175,104])

add_block('built-in/Inport',[sys,'/','频谱分析仪/system output'])
set_param([sys,'/','频谱分析仪/system output'],...
		'Port','2',...
		'position',[55,110,75,130])

add_block('built-in/S-Function',[sys,'/',['频谱分析仪/S-function block',13,'which calls M-file']])
set_param([sys,'/',['频谱分析仪/S-function block',13,'which calls M-file']],...
		'function name','sfuntf',...
		'parameters','fftpts,npts,HowOften,offset,ts,0',...
		'position',[250,72,295,108])

add_block('built-in/Inport',[sys,'/','频谱分析仪/system input'])
set_param([sys,'/','频谱分析仪/system input'],...
		'position',[55,30,75,50])
add_line([sys,'/','频谱分析仪'],[80,120;110,120;110,95;140,95])
add_line([sys,'/','频谱分析仪'],[80,40;110,40;110,80;140,80])
add_line([sys,'/','频谱分析仪'],[180,90;245,90])
set_param([sys,'/','频谱分析仪'],...
		'Mask Display','plot(0,0,100,100,[11,93,93,11,11],[90,90,40,40,90],[90,80,72,64,55,45,41,36,20],[50,50,52,59,67,75,78,80,80])',...
		'Mask Type','Spectrum analyzer.')
set_param([sys,'/','频谱分析仪'],...
		'Mask Dialogue','Spectrum analyser using graph window.\nHook input 1 to system input.\nHook input 2 to system output.|Length of buffer:|Number of points for fft:|Plot after how many points:|Sample time:')
set_param([sys,'/','频谱分析仪'],...
		'Mask Translate','npts=@1; fftpts=@2; HowOften=@3; ts=@4(1); if length(@4) > 1, offset = @4(2), else, offset = 0; end')
set_param([sys,'/','频谱分析仪'],...
		'Mask Help','Hook the system''s input to the first port of this block and the system''s output to second port. Results are displayed as the frequency response of this system. Use signals with high harmonic content for best results. ')
set_param([sys,'/','频谱分析仪'],...
		'Mask Entries','512\/512\/64\/0.05\/')


%     Finished composite block '频谱分析仪'.

set_param([sys,'/','频谱分析仪'],...
		'Font Name','宋体',...
		'Drop Shadow',4,...
		'position',[385,79,415,126])

add_block('built-in/Scope',[sys,'/','输入信号'])
set_param([sys,'/','输入信号'],...
		'Font Name','宋体',...
		'Vgain','1.500000',...
		'Hgain','20.000000',...
		'Vmax','3.000000',...
		'Hmax','40.000000',...
		'Window',[324,3,640,143])
open_system([sys,'/','输入信号'])
set_param([sys,'/','输入信号'],...
		'position',[385,14,415,56])
add_line(sys,[85,170;180,170])
add_line(sys,[135,170;135,315;180,315])
add_line(sys,[135,170;135,35;380,35])
add_line(sys,[260,170;380,170])
add_line(sys,[135,90;380,90])
add_line(sys,[315,170;315,115;380,115])
add_line(sys,[260,315;385,315])
add_line(sys,[135,230;385,230])
add_line(sys,[320,315;320,255;385,255])

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