tstcyclc.m

通信系统的matlab仿真程序

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


%     Subsystem  'Error rate1'.

new_system([sys,'/','Error rate1'])
set_param([sys,'/','Error rate1'],'Location',[255,368,509,510])

add_block('built-in/S-Function',[sys,'/','Error rate1/S-function'])
set_param([sys,'/','Error rate1/S-function'],...
		'function name','sbiterr',...
		'parameters','num_lin, K, TDelay, Sample',...
		'position',[130,52,195,78])

add_block('built-in/Mux',[sys,'/','Error rate1/Mux'])
set_param([sys,'/','Error rate1/Mux'],...
		'inputs','2',...
		'position',[65,46,100,79])

add_block('built-in/Inport',[sys,'/','Error rate1/in_1'])
set_param([sys,'/','Error rate1/in_1'],...
		'position',[15,45,35,65])

add_block('built-in/Inport',[sys,'/','Error rate1/in_2'])
set_param([sys,'/','Error rate1/in_2'],...
		'Port','2',...
		'position',[15,60,35,80])
add_line([sys,'/','Error rate1'],[105,65;125,65])
add_line([sys,'/','Error rate1'],[40,55;60,55])
add_line([sys,'/','Error rate1'],[40,70;60,70])
set_param([sys,'/','Error rate1'],...
		'Mask Display','Error\nmeter',...
		'Mask Type','Symbol/bit error counter')
set_param([sys,'/','Error rate1'],...
		'Mask Dialogue','Use the input data from the 1st port as the reference signal to detect the number of errors and error rate of the 2nd port inputted signal.|Bit per symbol:|Number of digits on display:|Delay between input (1st port) and output (2nd port):|Sampling time (sec):')
set_param([sys,'/','Error rate1'],...
		'Mask Translate','K=@1;num_lin=@2;TDelay=@3;Sample=@4;')
set_param([sys,'/','Error rate1'],...
		'Mask Help','The signal from the first port must be a scalar. The second input port can be a vector with a same-time delay to the input signal. The comparison happens only at the sampling point.')
set_param([sys,'/','Error rate1'],...
		'Mask Entries','4\/20\/0\/1\/')


%     Finished composite block 'Error rate1'.

set_param([sys,'/','Error rate1'],...
		'hide name',0,...
		'position',[820,8,900,52])

add_block('built-in/Mux',[sys,'/','Mux'])
set_param([sys,'/','Mux'],...
		'inputs','3',...
		'position',[750,103,785,137])


%     Subsystem  ['Binary vector',13,'cyclic decode'].

new_system([sys,'/',['Binary vector',13,'cyclic decode']])
set_param([sys,'/',['Binary vector',13,'cyclic decode']],'Location',[267,160,850,320])

add_block('built-in/Outport',[sys,'/',['Binary vector',13,'cyclic decode/out_1']])
set_param([sys,'/',['Binary vector',13,'cyclic decode/out_1']],...
		'position',[535,90,555,110])

add_block('built-in/S-Function',[sys,'/',['Binary vector',13,'cyclic decode/Modulus',13,'operation2']])
set_param([sys,'/',['Binary vector',13,'cyclic decode/Modulus',13,'operation2']],...
		'function name','arymodu',...
		'parameters','md',...
		'Mask Display','Modulo',...
		'Mask Type','Modulo')
set_param([sys,'/',['Binary vector',13,'cyclic decode/Modulus',13,'operation2']],...
		'Mask Dialogue','Modulus operation:|Modulus base:',...
		'Mask Translate','md = @1;')
set_param([sys,'/',['Binary vector',13,'cyclic decode/Modulus',13,'operation2']],...
		'Mask Help','This block output modulus operation of the input signal with module base as specified in block entry. The length of the output vector is the same as the input vector.')
set_param([sys,'/',['Binary vector',13,'cyclic decode/Modulus',13,'operation2']],...
		'Mask Entries','2\/',...
		'position',[455,81,515,119])

add_block('built-in/State-Space',[sys,'/',['Binary vector',13,'cyclic decode/Matrix',13,'Gain1']])
set_param([sys,'/',['Binary vector',13,'cyclic decode/Matrix',13,'Gain1']],...
		'A','[]',...
		'B','[]',...
		'C','[]',...
		'D','K',...
		'Mask Display','K',...
		'Mask Type','Matrix Gain',...
		'Mask Dialogue','Matrix Gain.|Gain matrix:')
set_param([sys,'/',['Binary vector',13,'cyclic decode/Matrix',13,'Gain1']],...
		'Mask Translate','K = @1;',...
		'Mask Help','Multiplies input vector by entered matrix to produce output vector (y=Au).',...
		'Mask Entries','G\/')
set_param([sys,'/',['Binary vector',13,'cyclic decode/Matrix',13,'Gain1']],...
		'position',[415,86,440,114])

add_block('built-in/S-Function',[sys,'/',['Binary vector',13,'cyclic decode/Modulus',13,'operation1']])
set_param([sys,'/',['Binary vector',13,'cyclic decode/Modulus',13,'operation1']],...
		'function name','arymodu',...
		'parameters','md',...
		'Mask Display','Modulo',...
		'Mask Type','Modulo')
set_param([sys,'/',['Binary vector',13,'cyclic decode/Modulus',13,'operation1']],...
		'Mask Dialogue','Modulus operation:|Modulus base:',...
		'Mask Translate','md = @1;')
set_param([sys,'/',['Binary vector',13,'cyclic decode/Modulus',13,'operation1']],...
		'Mask Help','This block output modulus operation of the input signal with module base as specified in block entry. The length of the output vector is the same as the input vector.')
set_param([sys,'/',['Binary vector',13,'cyclic decode/Modulus',13,'operation1']],...
		'Mask Entries','2\/',...
		'position',[340,81,400,119])

add_block('built-in/Sum',[sys,'/',['Binary vector',13,'cyclic decode/Sum']])
set_param([sys,'/',['Binary vector',13,'cyclic decode/Sum']],...
		'position',[305,90,325,110])

add_block('built-in/Combinatorial Logic',[sys,'/',['Binary vector',13,'cyclic decode/Combinatorial',13,'Logic']])
set_param([sys,'/',['Binary vector',13,'cyclic decode/Combinatorial',13,'Logic']],...
		'Truth Table','trt',...
		'position',[225,33,255,57])

add_block('built-in/S-Function',[sys,'/',['Binary vector',13,'cyclic decode/Modulus',13,'operation']])
set_param([sys,'/',['Binary vector',13,'cyclic decode/Modulus',13,'operation']],...
		'function name','arymodu',...
		'parameters','md',...
		'Mask Display','Modulo',...
		'Mask Type','Modulo')
set_param([sys,'/',['Binary vector',13,'cyclic decode/Modulus',13,'operation']],...
		'Mask Dialogue','Modulus operation:|Modulus base:',...
		'Mask Translate','md = @1;')
set_param([sys,'/',['Binary vector',13,'cyclic decode/Modulus',13,'operation']],...
		'Mask Help','This block output modulus operation of the input signal with module base as specified in block entry. The length of the output vector is the same as the input vector.')
set_param([sys,'/',['Binary vector',13,'cyclic decode/Modulus',13,'operation']],...
		'Mask Entries','2\/',...
		'position',[140,26,200,64])

add_block('built-in/State-Space',[sys,'/',['Binary vector',13,'cyclic decode/Matrix',13,'Gain']])
set_param([sys,'/',['Binary vector',13,'cyclic decode/Matrix',13,'Gain']],...
		'A','[]',...
		'B','[]',...
		'C','[]',...
		'D','K',...
		'Mask Display','K',...
		'Mask Type','Matrix Gain',...
		'Mask Dialogue','Matrix Gain.|Gain matrix:')
set_param([sys,'/',['Binary vector',13,'cyclic decode/Matrix',13,'Gain']],...
		'Mask Translate','K = @1;',...
		'Mask Help','Multiplies input vector by entered matrix to produce output vector (y=Au).',...
		'Mask Entries','H\/')
set_param([sys,'/',['Binary vector',13,'cyclic decode/Matrix',13,'Gain']],...
		'position',[100,31,125,59])

add_block('built-in/Inport',[sys,'/',['Binary vector',13,'cyclic decode/in_1']])
set_param([sys,'/',['Binary vector',13,'cyclic decode/in_1']],...
		'position',[30,35,50,55])
add_line([sys,'/',['Binary vector',13,'cyclic decode']],[520,100;530,100])
add_line([sys,'/',['Binary vector',13,'cyclic decode']],[445,100;450,100])
add_line([sys,'/',['Binary vector',13,'cyclic decode']],[405,100;410,100])
add_line([sys,'/',['Binary vector',13,'cyclic decode']],[330,100;335,100])
add_line([sys,'/',['Binary vector',13,'cyclic decode']],[55,45;95,45])
add_line([sys,'/',['Binary vector',13,'cyclic decode']],[80,45;80,105;300,105])
add_line([sys,'/',['Binary vector',13,'cyclic decode']],[260,45;285,45;285,95;300,95])
add_line([sys,'/',['Binary vector',13,'cyclic decode']],[205,45;220,45])
add_line([sys,'/',['Binary vector',13,'cyclic decode']],[130,45;135,45])
set_param([sys,'/',['Binary vector',13,'cyclic decode']],...
		'Mask Display','plot(0,0,100,100,x,y);Cyclic de',...
		'Mask Type','Vector I/O cyclic decode')
set_param([sys,'/',['Binary vector',13,'cyclic decode']],...
		'Mask Dialogue','Cyclic decode with codeword length N\nand message length K.|Codeword length N:|Message length K, or generator polynomial Pg:')
set_param([sys,'/',['Binary vector',13,'cyclic decode']],...
		'Mask Translate','n=@1;p=@2;if length(p)==1, p=cyclpoly(n,p,''min'');end;if isempty(p),error(''No such cyclic code exist'');end;[H,G]=cyclgen(n,p);trt=htruthtb(G);[n,m]=size(G);G=[zeros(n,m-n),eye(n)];[x,y]=codeicon(2);')
set_param([sys,'/',['Binary vector',13,'cyclic decode']],...
		'Mask Help','This block decodes length N binary vector codeword input into length K binary vector message output using cyclic code. When Pg is specified instead of K, the block uses the user-defined cyclic polynomial Pg. Use function CYCLPOLY for a list of cyclic generator polynomials.')
set_param([sys,'/',['Binary vector',13,'cyclic decode']],...
		'Mask Entries','7\/4\/')


%     Finished composite block ['Binary vector',13,'cyclic decode'].

set_param([sys,'/',['Binary vector',13,'cyclic decode']],...
		'hide name',0,...
		'position',[525,163,605,207])

add_block('built-in/S-Function',[sys,'/',['Integer vector',13,'to scalar2']])
set_param([sys,'/',['Integer vector',13,'to scalar2']],...
		'hide name',0,...
		'function name','simbi2de',...
		'parameters','p',...
		'Mask Display','Vec to sca\nconverter',...
		'Mask Type','Base M vector to decimal integer')
set_param([sys,'/',['Integer vector',13,'to scalar2']],...
		'Mask Dialogue','Convert base M integer with each\nposition in the input vector into scalar decimal output.|Input integer base M:',...