chou.m

一个带界面的信号与系统实验系统

function [ret,x0,str,ts,xts]=chou(t,x,u,flag);
%CHOU	is the M-file description of the SIMULINK system named CHOU.
%	The block-diagram can be displayed by typing: CHOU.
%
%	SYS=CHOU(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 CHOU 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 CHOU with a FLAG of zero:
%	[SIZES]=CHOU([],[],[],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',[4,42,623,456])
     open_system(sys)
end;
set_param(sys,'algorithm',     'RK-45')
set_param(sys,'Start time',    '0.0')
set_param(sys,'Stop time',     '9999')
set_param(sys,'Min step size', '0.000000001')
set_param(sys,'Max step size', '0.001')
set_param(sys,'Relative error','1e-5')
set_param(sys,'Return vars',   '')

add_block('built-in/Note',[sys,'/','连续时间抽样实验'])
set_param([sys,'/','连续时间抽样实验'],...
		'Font Name','黑体',...
		'Font Angle','italic',...
		'Font Size',15,...
		'position',[125,25,130,30])


%     Subsystem  ['抽样脉冲',13,'（抽样函数）'].

new_system([sys,'/',['抽样脉冲',13,'（抽样函数）']])
set_param([sys,'/',['抽样脉冲',13,'（抽样函数）']],'Location',[158,441,759,682])

add_block('built-in/Outport',[sys,'/',['抽样脉冲',13,'（抽样函数）/out_1']])
set_param([sys,'/',['抽样脉冲',13,'（抽样函数）/out_1']],...
		'hide name',0,...
		'position',[560,105,580,125])

add_block('built-in/Gain',[sys,'/',['抽样脉冲',13,'（抽样函数）/Gain']])
set_param([sys,'/',['抽样脉冲',13,'（抽样函数）/Gain']],...
		'hide name',0,...
		'Gain','ht',...
		'position',[510,102,535,128])

add_block('built-in/Product',[sys,'/',['抽样脉冲',13,'（抽样函数）/Product']])
set_param([sys,'/',['抽样脉冲',13,'（抽样函数）/Product']],...
		'hide name',0,...
		'position',[455,105,480,125])

add_block('built-in/Logical Operator',[sys,'/',['抽样脉冲',13,'（抽样函数）/Logical',13,'Operator1']])
set_param([sys,'/',['抽样脉冲',13,'（抽样函数）/Logical',13,'Operator1']],...
		'Operator','NOT',...
		'Number of Input Ports','1',...
		'position',[355,38,385,62])

add_block('built-in/Logical Operator',[sys,'/',['抽样脉冲',13,'（抽样函数）/Logical',13,'Operator']])
set_param([sys,'/',['抽样脉冲',13,'（抽样函数）/Logical',13,'Operator']],...
		'Operator','XOR',...
		'position',[285,38,315,62])

add_block('built-in/Constant',[sys,'/',['抽样脉冲',13,'（抽样函数）/Constant1']])
set_param([sys,'/',['抽样脉冲',13,'（抽样函数）/Constant1']],...
		'Value','stt',...
		'position',[295,149,320,171])

add_block('built-in/Clock',[sys,'/',['抽样脉冲',13,'（抽样函数）/Clock1']])
set_param([sys,'/',['抽样脉冲',13,'（抽样函数）/Clock1']],...
		'hide name',0,...
		'position',[295,100,315,120])

add_block('built-in/Relational Operator',[sys,'/',['抽样脉冲',13,'（抽样函数）/Relational',13,'operator']])
set_param([sys,'/',['抽样脉冲',13,'（抽样函数）/Relational',13,'operator']],...
		'hide name',0,...
		'position',[370,102,400,133])

add_block('built-in/Sum',[sys,'/',['抽样脉冲',13,'（抽样函数）/Sum1']])
set_param([sys,'/',['抽样脉冲',13,'（抽样函数）/Sum1']],...
		'hide name',0,...
		'inputs','+-',...
		'position',[130,120,150,140])

add_block('built-in/Sum',[sys,'/',['抽样脉冲',13,'（抽样函数）/Sum']])
set_param([sys,'/',['抽样脉冲',13,'（抽样函数）/Sum']],...
		'hide name',0,...
		'inputs','+-',...
		'position',[135,35,155,55])

add_block('built-in/Constant',[sys,'/',['抽样脉冲',13,'（抽样函数）/Constant']])
set_param([sys,'/',['抽样脉冲',13,'（抽样函数）/Constant']],...
		'position',[25,30,45,50])

add_block('built-in/Unit Delay',[sys,'/',['抽样脉冲',13,'（抽样函数）/Unit Delay']])
set_param([sys,'/',['抽样脉冲',13,'（抽样函数）/Unit Delay']],...
		'orientation',2,...
		'Sample time','[Ts,st1]',...
		'position',[145,75,195,95])

add_block('built-in/Unit Delay',[sys,'/',['抽样脉冲',13,'（抽样函数）/Unit Delay1']])
set_param([sys,'/',['抽样脉冲',13,'（抽样函数）/Unit Delay1']],...
		'orientation',2,...
		'Sample time','[Ts,st2]',...
		'x0','ini',...
		'position',[140,165,190,185])
add_line([sys,'/',['抽样脉冲',13,'（抽样函数）']],[390,50;430,50;430,110;450,110])
add_line([sys,'/',['抽样脉冲',13,'（抽样函数）']],[320,50;350,50])
add_line([sys,'/',['抽样脉冲',13,'（抽样函数）']],[325,160;345,160;345,125;365,125])
add_line([sys,'/',['抽样脉冲',13,'（抽样函数）']],[540,115;555,115])
add_line([sys,'/',['抽样脉冲',13,'（抽样函数）']],[485,115;505,115])
add_line([sys,'/',['抽样脉冲',13,'（抽样函数）']],[405,120;450,120])
add_line([sys,'/',['抽样脉冲',13,'（抽样函数）']],[320,110;365,110])
add_line([sys,'/',['抽样脉冲',13,'（抽样函数）']],[140,85;105,85;105,50;130,50])
add_line([sys,'/',['抽样脉冲',13,'（抽样函数）']],[160,45;230,45;230,85;200,85])
add_line([sys,'/',['抽样脉冲',13,'（抽样函数）']],[50,40;130,40])
add_line([sys,'/',['抽样脉冲',13,'（抽样函数）']],[155,130;210,130;210,175;195,175])
add_line([sys,'/',['抽样脉冲',13,'（抽样函数）']],[135,175;100,175;100,135;125,135])
add_line([sys,'/',['抽样脉冲',13,'（抽样函数）']],[50,40;70,40;70,125;125,125])
add_line([sys,'/',['抽样脉冲',13,'（抽样函数）']],[160,45;280,45])
add_line([sys,'/',['抽样脉冲',13,'（抽样函数）']],[155,130;245,130;245,55;280,55])
set_param([sys,'/',['抽样脉冲',13,'（抽样函数）']],...
		'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,'/',['抽样脉冲',13,'（抽样函数）']],...
		'Mask Dialogue','Pulse Generator.|Pulse period (secs):|Pulse width:|Pulse height:|Pulse start time:')
set_param([sys,'/',['抽样脉冲',13,'（抽样函数）']],...
		'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,'/',['抽样脉冲',13,'（抽样函数）']],...
		'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.')
set_param([sys,'/',['抽样脉冲',13,'（抽样函数）']],...
		'Mask Entries','pi/100\/0.001\/5\/0.0005\/')


%     Finished composite block ['抽样脉冲',13,'（抽样函数）'].

set_param([sys,'/',['抽样脉冲',13,'（抽样函数）']],...
		'Font Name','Arial',...
		'Font Size',11,...
		'position',[55,257,100,293])


%     Subsystem  '乘法器'.

new_system([sys,'/','乘法器'])
set_param([sys,'/','乘法器'],'Location',[33,479,253,634])

add_block('built-in/Product',[sys,'/','乘法器/Product'])
set_param([sys,'/','乘法器/Product'],...
		'position',[65,55,90,75])

add_block('built-in/Sum',[sys,'/','乘法器/Sum'])
set_param([sys,'/','乘法器/Sum'],...
		'inputs','+',...
		'position',[125,55,145,75])

add_block('built-in/Outport',[sys,'/','乘法器/out_1'])
set_param([sys,'/','乘法器/out_1'],...
		'position',[185,55,205,75])

add_block('built-in/Inport',[sys,'/','乘法器/in_2'])
set_param([sys,'/','乘法器/in_2'],...
		'Port','2',...
		'position',[15,90,35,110])

add_block('built-in/Inport',[sys,'/','乘法器/in_1'])
set_param([sys,'/','乘法器/in_1'],...
		'position',[15,25,35,45])
add_line([sys,'/','乘法器'],[40,100;60,70])
add_line([sys,'/','乘法器'],[40,35;60,60])
add_line([sys,'/','乘法器'],[95,65;120,65])
add_line([sys,'/','乘法器'],[150,65;180,65])
set_param([sys,'/','乘法器'],...
		'Mask Display','X',...
		'Mask Type','Dot Product',...