comcdcyc.m

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

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

add_block('built-in/Note',[sys,'/','Cyclic Coding//Decoding Library'])
set_param([sys,'/','Cyclic Coding//Decoding Library'],...
		'ForeGround',2,...
		'Font Weight','bold',...
		'Font Size',12,...
		'position',[200,10,205,15])


%     Subsystem  ['Vector cyclic',13,'co//dec demo'].

new_system([sys,'/',['Vector cyclic',13,'co//dec demo']])
set_param([sys,'/',['Vector cyclic',13,'co//dec demo']],'Location',[175,552,335,697])
set_param([sys,'/',['Vector cyclic',13,'co//dec demo']],...
		'Mask Display','Cyclic\nco/dec\ndemo',...
		'Mask Dialogue','eval(''tstcyclc'')')


%     Finished composite block ['Vector cyclic',13,'co//dec demo'].

set_param([sys,'/',['Vector cyclic',13,'co//dec demo']],...
		'BackGround',7,...
		'hide name',0,...
		'position',[290,58,370,102])


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

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

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

add_block('built-in/State-Space',[sys,'/',['Binary vector',13,'cyclic encode/Matrix',13,'Gain']])
set_param([sys,'/',['Binary vector',13,'cyclic encode/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 encode/Matrix',13,'Gain']],...
		'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 encode/Matrix',13,'Gain']],...
		'position',[110,31,135,59])

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

add_block('built-in/Outport',[sys,'/',['Binary vector',13,'cyclic encode/out_1']])
set_param([sys,'/',['Binary vector',13,'cyclic encode/out_1']],...
		'position',[290,35,310,55])
add_line([sys,'/',['Binary vector',13,'cyclic encode']],[55,45;105,45])
add_line([sys,'/',['Binary vector',13,'cyclic encode']],[140,45;175,45])
add_line([sys,'/',['Binary vector',13,'cyclic encode']],[245,45;285,45])
set_param([sys,'/',['Binary vector',13,'cyclic encode']],...
		'Mask Display','plot(0,0,100,100,x,y);Cyclic en',...
		'Mask Type','Vector I/O cyclic encode')
set_param([sys,'/',['Binary vector',13,'cyclic encode']],...
		'Mask Dialogue','Cyclic encode 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 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'';[x,y]=codeicon(1);')
set_param([sys,'/',['Binary vector',13,'cyclic encode']],...
		'Mask Help','This block encodes length K binary vector message input into length N binary vector codeword 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 encode']],...
		'Mask Entries','7\/4\/')


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

set_param([sys,'/',['Binary vector',13,'cyclic encode']],...
		'position',[50,58,130,102])


%     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/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_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/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/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/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/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/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,'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])