tstcycsq.m

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

function [ret,x0,str,ts,xts]=tstcycsq(t,x,u,flag);
%TSTCYCSQ	is the M-file description of the SIMULINK system named TSTCYCSQ.
%	The block-diagram can be displayed by typing: TSTCYCSQ.
%
%	SYS=TSTCYCSQ(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 TSTCYCSQ 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 TSTCYCSQ with a FLAG of zero:
%	[SIZES]=TSTCYCSQ([],[],[],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',[54,72,691,289])
     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', '1')
set_param(sys,'Relative error','1e-3')
set_param(sys,'Return vars',   '')
set_param(sys,'Load callback','N_cyc_demo=7;K_cyc_demo=4;GWN=.01;');

add_block('built-in/Note',[sys,'/',['This demo uses the MATLAB workspace variable to set the code world length N_cyc_demo=7,',13,'and message length K_cyc_demo= 4. You may assign the variables to be different values.      ']])
set_param([sys,'/',['This demo uses the MATLAB workspace variable to set the code world length N_cyc_demo=7,',13,'and message length K_cyc_demo= 4. You may assign the variables to be different values.      ']],...
		'position',[315,165,320,170])


%     Subsystem  ['Binary sequence',13,'cyclic encode'].

new_system([sys,'/',['Binary sequence',13,'cyclic encode']])
set_param([sys,'/',['Binary sequence',13,'cyclic encode']],'Location',[149,109,788,279])


%     Subsystem  ['Binary sequence',13,'cyclic encode/Binary-vector',13,'Cyclic encode'].

new_system([sys,'/',['Binary sequence',13,'cyclic encode/Binary-vector',13,'Cyclic encode']])
set_param([sys,'/',['Binary sequence',13,'cyclic encode/Binary-vector',13,'Cyclic encode']],'Location',[235,592,599,711])

add_block('built-in/Inport',[sys,'/',['Binary sequence',13,'cyclic encode/Binary-vector',13,'Cyclic encode/in_1']])
set_param([sys,'/',['Binary sequence',13,'cyclic encode/Binary-vector',13,'Cyclic encode/in_1']],...
		'position',[30,35,50,55])

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

add_block('built-in/S-Function',[sys,'/',['Binary sequence',13,'cyclic encode/Binary-vector',13,'Cyclic encode/Modulus',13,'operation']])
set_param([sys,'/',['Binary sequence',13,'cyclic encode/Binary-vector',13,'Cyclic encode/Modulus',13,'operation']],...
		'function name','arymodu',...
		'parameters','md',...
		'Mask Display','Modulo',...
		'Mask Type','Modulo')
set_param([sys,'/',['Binary sequence',13,'cyclic encode/Binary-vector',13,'Cyclic encode/Modulus',13,'operation']],...
		'Mask Dialogue','Modulus operation:|Modulus base:',...
		'Mask Translate','md = @1;')
set_param([sys,'/',['Binary sequence',13,'cyclic encode/Binary-vector',13,'Cyclic encode/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 sequence',13,'cyclic encode/Binary-vector',13,'Cyclic encode/Modulus',13,'operation']],...
		'Mask Entries','2\/',...
		'position',[180,26,240,64])

add_block('built-in/Outport',[sys,'/',['Binary sequence',13,'cyclic encode/Binary-vector',13,'Cyclic encode/out_1']])
set_param([sys,'/',['Binary sequence',13,'cyclic encode/Binary-vector',13,'Cyclic encode/out_1']],...
		'position',[290,35,310,55])
add_line([sys,'/',['Binary sequence',13,'cyclic encode/Binary-vector',13,'Cyclic encode']],[55,45;105,45])
add_line([sys,'/',['Binary sequence',13,'cyclic encode/Binary-vector',13,'Cyclic encode']],[140,45;175,45])
add_line([sys,'/',['Binary sequence',13,'cyclic encode/Binary-vector',13,'Cyclic encode']],[245,45;285,45])
set_param([sys,'/',['Binary sequence',13,'cyclic encode/Binary-vector',13,'Cyclic encode']],...
		'Mask Display','Cyclic\nencode',...
		'Mask Type','Cyclic encoding')
set_param([sys,'/',['Binary sequence',13,'cyclic encode/Binary-vector',13,'Cyclic encode']],...
		'Mask Dialogue','Cyclic encoding:|Code word length:|Generator polynomial or message length:')
set_param([sys,'/',['Binary sequence',13,'cyclic encode/Binary-vector',13,'Cyclic encode']],...
		'Mask Translate','n=@1;p=@2;if length(p)==1,p=cyclpoly(n,p,''min'');if isempty(p),error(''No such cyclic code exist'');end;end;[h,G]=cyclgen(n, p);G=G'';')
set_param([sys,'/',['Binary sequence',13,'cyclic encode/Binary-vector',13,'Cyclic encode']],...
		'Mask Help','This block encodes the input binary signal using cyclic coding method. Let k=code_word_length - degree_of_generator_poly. The input of this block is a k-array binary message vector. The output of this block is an n-array binary code word vector. If the entry for "Generator polynomial" is a number instead of polynomial, the number is taken as the message length (input vector length).')
set_param([sys,'/',['Binary sequence',13,'cyclic encode/Binary-vector',13,'Cyclic encode']],...
		'Mask Entries','n\/p\/')


%     Finished composite block ['Binary sequence',13,'cyclic encode/Binary-vector',13,'Cyclic encode'].

set_param([sys,'/',['Binary sequence',13,'cyclic encode/Binary-vector',13,'Cyclic encode']],...
		'position',[320,15,385,55])

add_block('built-in/Outport',[sys,'/',['Binary sequence',13,'cyclic encode/out_2']])
set_param([sys,'/',['Binary sequence',13,'cyclic encode/out_2']],...
		'Port','2',...
		'position',[525,45,545,65])


%     Subsystem  ['Binary sequence',13,'cyclic encode/Triggered',13,'buffer down'].

new_system([sys,'/',['Binary sequence',13,'cyclic encode/Triggered',13,'buffer down']])
set_param([sys,'/',['Binary sequence',13,'cyclic encode/Triggered',13,'buffer down']],'Location',[306,291,744,461])

add_block('built-in/Outport',[sys,'/',['Binary sequence',13,'cyclic encode/Triggered',13,'buffer down/out_1']])
set_param([sys,'/',['Binary sequence',13,'cyclic encode/Triggered',13,'buffer down/out_1']],...
		'position',[390,60,410,80])

add_block('built-in/Outport',[sys,'/',['Binary sequence',13,'cyclic encode/Triggered',13,'buffer down/out_2']])
set_param([sys,'/',['Binary sequence',13,'cyclic encode/Triggered',13,'buffer down/out_2']],...
		'Port','2',...
		'position',[355,80,375,100])

add_block('built-in/Mux',[sys,'/',['Binary sequence',13,'cyclic encode/Triggered',13,'buffer down/Mux']])
set_param([sys,'/',['Binary sequence',13,'cyclic encode/Triggered',13,'buffer down/Mux']],...
		'inputs','[in_size 1 1]',...
		'position',[110,59,140,101])

add_block('built-in/Demux',[sys,'/',['Binary sequence',13,'cyclic encode/Triggered',13,'buffer down/Demux']])
set_param([sys,'/',['Binary sequence',13,'cyclic encode/Triggered',13,'buffer down/Demux']],...
		'outputs','[out_size 1]',...
		'position',[290,62,330,98])

add_block('built-in/S-Function',[sys,'/',['Binary sequence',13,'cyclic encode/Triggered',13,'buffer down/S-function']])
set_param([sys,'/',['Binary sequence',13,'cyclic encode/Triggered',13,'buffer down/S-function']],...
		'function name','regdown',...
		'parameters','ini_idx, idx_inc, thrshld, cycl_flag',...
		'position',[185,70,235,90])

add_block('built-in/Inport',[sys,'/',['Binary sequence',13,'cyclic encode/Triggered',13,'buffer down/in_1']])
set_param([sys,'/',['Binary sequence',13,'cyclic encode/Triggered',13,'buffer down/in_1']],...
		'position',[15,55,35,75])

add_block('built-in/Inport',[sys,'/',['Binary sequence',13,'cyclic encode/Triggered',13,'buffer down/in_2']])
set_param([sys,'/',['Binary sequence',13,'cyclic encode/Triggered',13,'buffer down/in_2']],...
		'Port','2',...
		'position',[40,70,60,90])

add_block('built-in/Inport',[sys,'/',['Binary sequence',13,'cyclic encode/Triggered',13,'buffer down/in_3']])
set_param([sys,'/',['Binary sequence',13,'cyclic encode/Triggered',13,'buffer down/in_3']],...
		'Port','3',...
		'position',[65,85,85,105])
add_line([sys,'/',['Binary sequence',13,'cyclic encode/Triggered',13,'buffer down']],[240,80;285,80])
add_line([sys,'/',['Binary sequence',13,'cyclic encode/Triggered',13,'buffer down']],[145,80;180,80])
add_line([sys,'/',['Binary sequence',13,'cyclic encode/Triggered',13,'buffer down']],[40,65;105,65])
add_line([sys,'/',['Binary sequence',13,'cyclic encode/Triggered',13,'buffer down']],[335,70;385,70])
add_line([sys,'/',['Binary sequence',13,'cyclic encode/Triggered',13,'buffer down']],[65,80;105,80])
add_line([sys,'/',['Binary sequence',13,'cyclic encode/Triggered',13,'buffer down']],[335,90;350,90])
add_line([sys,'/',['Binary sequence',13,'cyclic encode/Triggered',13,'buffer down']],[90,95;105,95])
set_param([sys,'/',['Binary sequence',13,'cyclic encode/Triggered',13,'buffer down']],...
		'Mask Display','Triggered\nbuffer dn',...
		'Mask Type','Triggered buffer down.')
set_param([sys,'/',['Binary sequence',13,'cyclic encode/Triggered',13,'buffer down']],...
		'Mask Dialogue','Triggered buffer down.|Input signal vector size:|Initial index (integer vector):|Increament for each index (integer vector):|Trigger threshold:|When index exceed register size, 1 for cyclic ouput, 0 for zero output:')
set_param([sys,'/',['Binary sequence',13,'cyclic encode/Triggered',13,'buffer down']],...
		'Mask Translate','in_size=@1; ini_idx=@2; idx_inc=@3; thrshld=@4; cycl_flag=@5; out_size=length(ini_idx);')
set_param([sys,'/',['Binary sequence',13,'cyclic encode/Triggered',13,'buffer down']],...
		'Mask Help','The register inside this block take the signal vector at the 1st inport by the raising edge of the second inport. The fisrt outport outputs a signal vector at the raising edge of the third inport.')
set_param([sys,'/',['Binary sequence',13,'cyclic encode/Triggered',13,'buffer down']],...
		'Mask Entries','n\/0\/1\/.1\/0\/')


%     Finished composite block ['Binary sequence',13,'cyclic encode/Triggered',13,'buffer down'].

set_param([sys,'/',['Binary sequence',13,'cyclic encode/Triggered',13,'buffer down']],...
		'BackGround',3,...
		'position',[435,28,495,62])

add_block('built-in/Outport',[sys,'/',['Binary sequence',13,'cyclic encode/out_1']])
set_param([sys,'/',['Binary sequence',13,'cyclic encode/out_1']],...
		'position',[565,30,585,50])

add_block('built-in/Demux',[sys,'/',['Binary sequence',13,'cyclic encode/Demux']])
set_param([sys,'/',['Binary sequence',13,'cyclic encode/Demux']],...
		'outputs','2',...
		'position',[130,80,170,115])


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

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

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

add_block('built-in/S-Function',[sys,'/',['Binary sequence',13,'cyclic encode/Pulses deivide',13,'same sample time/S-function1']])
set_param([sys,'/',['Binary sequence',13,'cyclic encode/Pulses deivide',13,'same sample time/S-function1']],...
		'function name','homopuls',...
		'parameters','samp_time, out_divid, offset',...