tstrsbnv.m

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

function [ret,x0,str,ts,xts]=tstrsbnv(t,x,u,flag);
%TSTRSBNV	is the M-file description of the SIMULINK system named TSTRSBNV.
%	The block-diagram can be displayed by typing: TSTRSBNV.
%
%	SYS=TSTRSBNV(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 TSTRSBNV 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 TSTRSBNV with a FLAG of zero:
%	[SIZES]=TSTRSBNV([],[],[],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',[39,71,792,341])
     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', '1')
set_param(sys,'Relative error','1e-3')
set_param(sys,'Return vars',   '')
set_param(sys,'Load callback','N_rs_demo=15;K_rs_demo=9;M_rs_demo=4;GWN=.05;');


%     Subsystem  ['Pulses divide',13,'same sample time'].

new_system([sys,'/',['Pulses divide',13,'same sample time']])
set_param([sys,'/',['Pulses divide',13,'same sample time']],'Location',[55,192,283,286])

add_block('built-in/Outport',[sys,'/',['Pulses divide',13,'same sample time/Outport']])
set_param([sys,'/',['Pulses divide',13,'same sample time/Outport']],...
		'position',[160,35,180,55])

add_block('built-in/S-Function',[sys,'/',['Pulses divide',13,'same sample time/S-function1']])
set_param([sys,'/',['Pulses divide',13,'same sample time/S-function1']],...
		'function name','homopuls',...
		'parameters','samp_time, out_divid, offset',...
		'position',[40,31,105,59])
add_line([sys,'/',['Pulses divide',13,'same sample time']],[110,45;155,45])
set_param([sys,'/',['Pulses divide',13,'same sample time']],...
		'Mask Display','plot(0,0,100,100,x,y);Pulse',...
		'Mask Type','Pulse generator')
set_param([sys,'/',['Pulses divide',13,'same sample time']],...
		'Mask Dialogue','Vector pulses with sample rate of each pulse as Sample_time / Divider.|Sample time (scalar, sec):|Divider (integer vector):|Offset (same length as Divider):')
set_param([sys,'/',['Pulses divide',13,'same sample time']],...
		'Mask Translate','samp_time=@1;out_divid=@2;offset=@3;[x,y]=srcsicon(2);')
set_param([sys,'/',['Pulses divide',13,'same sample time']],...
		'Mask Help','This block generates a vector pulse. The output length is the same as the divider vector length. The ith element of the output has the sample time Sample_time/Divider(i). All elements in divider must be integers.')
set_param([sys,'/',['Pulses divide',13,'same sample time']],...
		'Mask Entries','1\/N_rs_demo*M_rs_demo\/.5/N_rs_demo\/')


%     Finished composite block ['Pulses divide',13,'same sample time'].

set_param([sys,'/',['Pulses divide',13,'same sample time']],...
		'hide name',0,...
		'position',[250,138,330,182])


%     Subsystem  ['Binary sequence',13,'RS decode'].

new_system([sys,'/',['Binary sequence',13,'RS decode']])
set_param([sys,'/',['Binary sequence',13,'RS decode']],'Location',[108,572,804,731])

add_block('built-in/Outport',[sys,'/',['Binary sequence',13,'RS decode/out_2']])
set_param([sys,'/',['Binary sequence',13,'RS decode/out_2']],...
		'Port','2',...
		'position',[480,55,500,75])


%     Subsystem  ['Binary sequence',13,'RS decode/Rising edge',13,'detector'].

new_system([sys,'/',['Binary sequence',13,'RS decode/Rising edge',13,'detector']])
set_param([sys,'/',['Binary sequence',13,'RS decode/Rising edge',13,'detector']],'Location',[282,116,671,290])

add_block('built-in/Relational Operator',[sys,'/',['Binary sequence',13,'RS decode/Rising edge',13,'detector/Relational',13,'Operator']])
set_param([sys,'/',['Binary sequence',13,'RS decode/Rising edge',13,'detector/Relational',13,'Operator']],...
		'position',[135,18,165,42])

add_block('built-in/Inport',[sys,'/',['Binary sequence',13,'RS decode/Rising edge',13,'detector/in_1']])
set_param([sys,'/',['Binary sequence',13,'RS decode/Rising edge',13,'detector/in_1']],...
		'position',[10,15,30,35])

add_block('built-in/Memory',[sys,'/',['Binary sequence',13,'RS decode/Rising edge',13,'detector/Memory1']])
set_param([sys,'/',['Binary sequence',13,'RS decode/Rising edge',13,'detector/Memory1']],...
		'position',[125,80,165,110])

add_block('built-in/Constant',[sys,'/',['Binary sequence',13,'RS decode/Rising edge',13,'detector/one1']])
set_param([sys,'/',['Binary sequence',13,'RS decode/Rising edge',13,'detector/one1']],...
		'Value','thld',...
		'position',[30,135,60,155])

add_block('built-in/Relational Operator',[sys,'/',['Binary sequence',13,'RS decode/Rising edge',13,'detector/Relational',13,'Operator2']])
set_param([sys,'/',['Binary sequence',13,'RS decode/Rising edge',13,'detector/Relational',13,'Operator2']],...
		'Operator','>',...
		'position',[265,23,295,47])

add_block('built-in/Relational Operator',[sys,'/',['Binary sequence',13,'RS decode/Rising edge',13,'detector/Relational',13,'Operator1']])
set_param([sys,'/',['Binary sequence',13,'RS decode/Rising edge',13,'detector/Relational',13,'Operator1']],...
		'position',[200,88,230,112])

add_block('built-in/Outport',[sys,'/',['Binary sequence',13,'RS decode/Rising edge',13,'detector/out_1']])
set_param([sys,'/',['Binary sequence',13,'RS decode/Rising edge',13,'detector/out_1']],...
		'position',[340,25,360,45])
add_line([sys,'/',['Binary sequence',13,'RS decode/Rising edge',13,'detector']],[300,35;335,35])
add_line([sys,'/',['Binary sequence',13,'RS decode/Rising edge',13,'detector']],[235,100;240,100;240,40;260,40])
add_line([sys,'/',['Binary sequence',13,'RS decode/Rising edge',13,'detector']],[170,30;260,30])
add_line([sys,'/',['Binary sequence',13,'RS decode/Rising edge',13,'detector']],[65,145;100,145;100,35;130,35])
add_line([sys,'/',['Binary sequence',13,'RS decode/Rising edge',13,'detector']],[100,145;180,145;180,105;195,105])
add_line([sys,'/',['Binary sequence',13,'RS decode/Rising edge',13,'detector']],[170,95;195,95])
add_line([sys,'/',['Binary sequence',13,'RS decode/Rising edge',13,'detector']],[35,25;130,25])
add_line([sys,'/',['Binary sequence',13,'RS decode/Rising edge',13,'detector']],[45,25;45,95;120,95])
set_param([sys,'/',['Binary sequence',13,'RS decode/Rising edge',13,'detector']],...
		'Mask Display','Edge\nDetector',...
		'Mask Type','Edge Detection',...
		'Mask Dialogue','Rising edge detection:|Threshold:',...
		'Mask Translate','thld=@1;')
set_param([sys,'/',['Binary sequence',13,'RS decode/Rising edge',13,'detector']],...
		'Mask Help','When the input signal rising edge passed the given shreshold, this block outputs one, otherwise outputs zero.',...
		'Mask Entries','0.1\/')


%     Finished composite block ['Binary sequence',13,'RS decode/Rising edge',13,'detector'].

set_param([sys,'/',['Binary sequence',13,'RS decode/Rising edge',13,'detector']],...
		'ForeGround',2,...
		'position',[280,51,330,79])


%     Subsystem  ['Binary sequence',13,'RS decode/Pulses deivide',13,'same sample time'].

new_system([sys,'/',['Binary sequence',13,'RS decode/Pulses deivide',13,'same sample time']])
set_param([sys,'/',['Binary sequence',13,'RS decode/Pulses deivide',13,'same sample time']],'Location',[55,192,283,286])

add_block('built-in/Outport',[sys,'/',['Binary sequence',13,'RS decode/Pulses deivide',13,'same sample time/Outport']])
set_param([sys,'/',['Binary sequence',13,'RS decode/Pulses deivide',13,'same sample time/Outport']],...
		'position',[160,35,180,55])

add_block('built-in/S-Function',[sys,'/',['Binary sequence',13,'RS decode/Pulses deivide',13,'same sample time/S-function1']])
set_param([sys,'/',['Binary sequence',13,'RS decode/Pulses deivide',13,'same sample time/S-function1']],...
		'function name','homopuls',...
		'parameters','samp_time, out_divid, offset',...
		'position',[40,31,105,59])
add_line([sys,'/',['Binary sequence',13,'RS decode/Pulses deivide',13,'same sample time']],[110,45;155,45])
set_param([sys,'/',['Binary sequence',13,'RS decode/Pulses deivide',13,'same sample time']],...
		'Mask Display','Vector\nPulse',...
		'Mask Type','Pulse generator')
set_param([sys,'/',['Binary sequence',13,'RS decode/Pulses deivide',13,'same sample time']],...
		'Mask Dialogue','Vector pulses with sample rate of each pulse as Sample_time ./Divider. |Sample time (scalar, sec):|Divider (integer vector):|Offset (have same dimension as Divider):')
set_param([sys,'/',['Binary sequence',13,'RS decode/Pulses deivide',13,'same sample time']],...
		'Mask Translate','samp_time=@1;out_divid=@2;offset=@3;')
set_param([sys,'/',['Binary sequence',13,'RS decode/Pulses deivide',13,'same sample time']],...
		'Mask Help','This block generator a vector pulses. The output length is the same as the vector defined in "Divider". The ith element of the output has the sample time Sample_time/Divider(i). All elements in "Divider" must be integers.')
set_param([sys,'/',['Binary sequence',13,'RS decode/Pulses deivide',13,'same sample time']],...
		'Mask Entries','samp*k\/[n k]\/[0 0]\/')


%     Finished composite block ['Binary sequence',13,'RS decode/Pulses deivide',13,'same sample time'].

set_param([sys,'/',['Binary sequence',13,'RS decode/Pulses deivide',13,'same sample time']],...
		'position',[15,86,75,124])

add_block('built-in/Demux',[sys,'/',['Binary sequence',13,'RS decode/Demux']])
set_param([sys,'/',['Binary sequence',13,'RS decode/Demux']],...
		'outputs','2',...
		'position',[95,85,140,120])


%     Subsystem  ['Binary sequence',13,'RS decode/Register-shift'].

new_system([sys,'/',['Binary sequence',13,'RS decode/Register-shift']])
set_param([sys,'/',['Binary sequence',13,'RS decode/Register-shift']],'Location',[237,587,592,717])

add_block('built-in/Inport',[sys,'/',['Binary sequence',13,'RS decode/Register-shift/in_2']])
set_param([sys,'/',['Binary sequence',13,'RS decode/Register-shift/in_2']],...
		'Port','2',...
		'position',[30,70,50,90])

add_block('built-in/Inport',[sys,'/',['Binary sequence',13,'RS decode/Register-shift/in_1']])
set_param([sys,'/',['Binary sequence',13,'RS decode/Register-shift/in_1']],...
		'position',[30,55,50,75])

add_block('built-in/Outport',[sys,'/',['Binary sequence',13,'RS decode/Register-shift/out_2']])
set_param([sys,'/',['Binary sequence',13,'RS decode/Register-shift/out_2']],...
		'Port','2',...
		'position',[330,70,350,90])

add_block('built-in/Outport',[sys,'/',['Binary sequence',13,'RS decode/Register-shift/out_1']])
set_param([sys,'/',['Binary sequence',13,'RS decode/Register-shift/out_1']],...
		'position',[330,55,350,75])

add_block('built-in/Mux',[sys,'/',['Binary sequence',13,'RS decode/Register-shift/Mux']])
set_param([sys,'/',['Binary sequence',13,'RS decode/Register-shift/Mux']],...
		'inputs','2',...
		'position',[80,56,115,89])

add_block('built-in/S-Function',[sys,'/',['Binary sequence',13,'RS decode/Register-shift/S-function']])
set_param([sys,'/',['Binary sequence',13,'RS decode/Register-shift/S-function']],...
		'function name','regshift',...