tstrsbin.m

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

function [ret,x0,str,ts,xts]=tstrsbin(t,x,u,flag);
%TSTRSBIN     is the M-file description of the SIMULINK system named TSTRSBIN.
%       The block-diagram can be displayed by typing: TSTRSBIN.
%
%       SYS=TSTRSBIN(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 TSTRSBIN 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 TSTRSBIN with a FLAG of zero:
%  [SIZES]=TSTRSBIN([],[],[],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',[137,230,951,530])
     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_rs_demo=7;K_rs_demo=5;M_rs_demo=3;');

add_block('built-in/Sum',[sys,'/','Sum3'])
set_param([sys,'/','Sum3'],...
           'BackGround',7,...
              'inputs','+',...
                'position',[385,160,405,180])


%     Subsystem  ['Random-integer',13,'generator'].

new_system([sys,'/',['Random-integer',13,'generator']])
set_param([sys,'/',['Random-integer',13,'generator']],'Location',[373,639,729,808])

add_block('built-in/Constant',[sys,'/',['Random-integer',13,'generator/Constant']])
set_param([sys,'/',['Random-integer',13,'generator/Constant']],...
         'Value','x',...
         'position',[65,105,95,125])

add_block('built-in/State-Space',[sys,'/',['Random-integer',13,'generator/Matrix',13,'Gain']])
set_param([sys,'/',['Random-integer',13,'generator/Matrix',13,'Gain']],...
             'A','[]',...
            'B','[]',...
            'C','[]',...
            'D','K',...
             'Mask Display','K',...
          'Mask Type','Matrix Gain',...
           'Mask Dialogue','Matrix Gain.|Gain matrix:')
set_param([sys,'/',['Random-integer',13,'generator/Matrix',13,'Gain']],...
          'Mask Translate','K = @1;',...
          'Mask Help','Multiplies input vector by entered matrix to produce output vector (y=Au).',...
            'Mask Entries','ones(1,N-1)\/')
set_param([sys,'/',['Random-integer',13,'generator/Matrix',13,'Gain']],...
               'position',[230,35,275,65])

add_block('built-in/Outport',[sys,'/',['Random-integer',13,'generator/out_1']])
set_param([sys,'/',['Random-integer',13,'generator/out_1']],...
               'position',[310,40,330,60])

add_block('built-in/Relational Operator',[sys,'/',['Random-integer',13,'generator/Relational',13,'Operator']])
set_param([sys,'/',['Random-integer',13,'generator/Relational',13,'Operator']],...
             'position',[155,38,185,62])

add_block('built-in/Signal Generator',[sys,'/',['Random-integer',13,'generator/Uniform',13,'noise gen.']])
set_param([sys,'/',['Random-integer',13,'generator/Uniform',13,'noise gen.']],...
          'Peak','1.000000',...
           'Peak Range','5.000000',...
             'Freq','1234.000000',...
                'Freq Range','6170.000000',...
          'Wave','Rnd',...
                'Units','Rads')
set_param([sys,'/',['Random-integer',13,'generator/Uniform',13,'noise gen.']],...
                'position',[35,28,80,62])
add_line([sys,'/',['Random-integer',13,'generator']],[85,45;150,45])
add_line([sys,'/',['Random-integer',13,'generator']],[100,115;115,115;115,55;150,55])
add_line([sys,'/',['Random-integer',13,'generator']],[190,50;225,50])
add_line([sys,'/',['Random-integer',13,'generator']],[280,50;305,50])
set_param([sys,'/',['Random-integer',13,'generator']],...
           'Mask Display','Random\nNumber',...
             'Mask Type','Random number generator.')
set_param([sys,'/',['Random-integer',13,'generator']],...
                'Mask Dialogue','Randomly generate integers in the range of [0,1,...,N-1]. N is the M-ary number defined below:|M-ary number:',...
                'Mask Translate','N=@1;x=[(2-N)/N:2/N:(N-2)/N];')
set_param([sys,'/',['Random-integer',13,'generator']],...
              'Mask Help','This block generate integer number 0 to M-1, where M is a integer. The probability of the apperance for each number is not guaranteed to be the same.')
set_param([sys,'/',['Random-integer',13,'generator']],...
           'Mask Entries','2^(K_rs_demo*M_rs_demo)\/')


%     Finished composite block ['Random-integer',13,'generator'].

set_param([sys,'/',['Random-integer',13,'generator']],...
           'BackGround',5,...
              'position',[25,35,85,75])

add_block('built-in/Mux',[sys,'/','Mux1'])
set_param([sys,'/','Mux1'],...
               'inputs','2',...
                'position',[665,145,700,180])


%     Subsystem  ['Graph',13,'err'].

new_system([sys,'/',['Graph',13,'err']])
set_param([sys,'/',['Graph',13,'err']],'Location',[0,59,274,252])

add_block('built-in/Inport',[sys,'/',['Graph',13,'err/x']])
set_param([sys,'/',['Graph',13,'err/x']],...
               'position',[65,55,85,75])

add_block('built-in/S-Function',[sys,'/',['Graph',13,'err/S-function',13,'M-file which plots',13,'lines',13,'']])
set_param([sys,'/',['Graph',13,'err/S-function',13,'M-file which plots',13,'lines',13,'']],...
                'function name','sfuny',...
             'parameters','ax, color,dt',...
         'position',[130,55,180,75])
add_line([sys,'/',['Graph',13,'err']],[90,65;125,65])
set_param([sys,'/',['Graph',13,'err']],...
              'Mask Display','plot(0,0,100,100,[90,10,10,10,90,90,10],[65,65,90,40,40,90,90],[90,78,69,54,40,31,25,10],[77,60,48,46,56,75,81,84])',...
                'Mask Type','Graph scope.')
set_param([sys,'/',['Graph',13,'err']],...
           'Mask Dialogue','Graph scope using MATLAB graph window.\nEnter plotting ranges and line type.|Time range:|y-min:|y-max:|Line type (rgbw-:*). Seperate each plot by ''/'':')
set_param([sys,'/',['Graph',13,'err']],...
           'Mask Translate','color = @4; ax = [0, @1, @2, @3]; dt = -1;')
set_param([sys,'/',['Graph',13,'err']],...
                'Mask Help','This block plots to the MATLAB graph window and can be used as an improved version of the Scope block. Look at the m-file sfuny.m to see how it works. This block can take scalar or vector input signal.')
set_param([sys,'/',['Graph',13,'err']],...
              'Mask Entries','10\/-1.2\/.2+floor((N_rs_demo-K_rs_demo)/2)\/''mx/c+/m-/w-/b+''\/')


%     Finished composite block ['Graph',13,'err'].

set_param([sys,'/',['Graph',13,'err']],...
         'position',[755,146,785,184])

add_block('built-in/Note',[sys,'/','In this demo, you may chooset the variable code word length N, message length K, and M.'])
set_param([sys,'/','In this demo, you may chooset the variable code word length N, message length K, and M.'],...
            'position',[315,220,320,225])

add_block('built-in/Note',[sys,'/','The default set of this demo is N = N_rs_demo=7; K =K_rs_demo= 5; M =M_rs_demo = 3. One burst error-correction.'])
set_param([sys,'/','The default set of this demo is N = N_rs_demo=7; K =K_rs_demo= 5; M =M_rs_demo = 3. One burst error-correction.'],...
            'position',[390,240,395,245])


%     Subsystem  ['Binary vector',13,'noise generator'].

new_system([sys,'/',['Binary vector',13,'noise generator']])
set_param([sys,'/',['Binary vector',13,'noise generator']],'Location',[334,411,752,588])

add_block('built-in/S-Function',[sys,'/',['Binary vector',13,'noise generator/S-function']])
set_param([sys,'/',['Binary vector',13,'noise generator/S-function']],...
          'function name','simrdbit',...
          'parameters','n, prob, seed',...
                'position',[90,35,140,55])

add_block('built-in/Inport',[sys,'/',['Binary vector',13,'noise generator/Inport']])
set_param([sys,'/',['Binary vector',13,'noise generator/Inport']],...
             'position',[30,65,50,85])

add_block('built-in/Switch',[sys,'/',['Binary vector',13,'noise generator/Switch']])
set_param([sys,'/',['Binary vector',13,'noise generator/Switch']],...
              'Threshold','0.5',...
           'position',[215,59,245,91])

add_block('built-in/Memory',[sys,'/',['Binary vector',13,'noise generator/Memory']])
set_param([sys,'/',['Binary vector',13,'noise generator/Memory']],...
            'orientation',2,...
             'x0','zeros(1,n)',...
           'position',[215,110,255,140])

add_block('built-in/Outport',[sys,'/',['Binary vector',13,'noise generator/out_1']])
set_param([sys,'/',['Binary vector',13,'noise generator/out_1']],...
           'position',[340,65,360,85])
add_line([sys,'/',['Binary vector',13,'noise generator']],[250,75;335,75])
add_line([sys,'/',['Binary vector',13,'noise generator']],[145,45;175,45;175,65;210,65])
add_line([sys,'/',['Binary vector',13,'noise generator']],[55,75;210,75])
add_line([sys,'/',['Binary vector',13,'noise generator']],[70,75;70,45;85,45])
add_line([sys,'/',['Binary vector',13,'noise generator']],[210,125;175,125;175,85;210,85])
add_line([sys,'/',['Binary vector',13,'noise generator']],[280,75;280,125;260,125])
set_param([sys,'/',['Binary vector',13,'noise generator']],...
         'Mask Display','Random\nB-vector',...
           'Mask Type','Binary vector noise generator.')
set_param([sys,'/',['Binary vector',13,'noise generator']],...
             'Mask Dialogue','Generatea binary vector with at most given number of ones randomly distributed in the vector.|Vector length(can''t be changed during simulation):|Probability of a "one", two "ones",three "ones"...:|Seed:')
set_param([sys,'/',['Binary vector',13,'noise generator']],...
            'Mask Translate','n=@1;prob=@2;seed=@3;')
set_param([sys,'/',['Binary vector',13,'noise generator']],...
         'Mask Help','This block outputs a zero vector with ''1''s randomly distributed in any bit of the vector. The maximum number of ''1''s appears in the vector is same as the length of the vector "probability". The value of ith element of the ''probability'' is the probability of i number of ''1''s in the vector, e.g. [1/4 1/4 1/4] means probability of single, double, three ''1''s or no ''1''s in the vector is 1/4. The block refresh its output whenever the input is no-zero.')
set_param([sys,'/',['Binary vector',13,'noise generator']],...
              'Mask Entries','N_rs_demo*M_rs_demo\/1/2\/12345\/')


%     Finished composite block ['Binary vector',13,'noise generator'].

set_param([sys,'/',['Binary vector',13,'noise generator']],...
         'position',[245,153,300,187])

add_block('built-in/Constant',[sys,'/','Const'])
set_param([sys,'/','Const'],...
            'position',[190,160,210,180])

add_block('built-in/Constant',[sys,'/','Const1'])
set_param([sys,'/','Const1'],...
          'position',[195,65,215,85])


%     Subsystem  ['Binary-vector',13,'R-S encode'].

new_system([sys,'/',['Binary-vector',13,'R-S encode']])
set_param([sys,'/',['Binary-vector',13,'R-S encode']],'Location',[160,111,633,282])

add_block('built-in/Outport',[sys,'/',['Binary-vector',13,'R-S encode/out_1']])
set_param([sys,'/',['Binary-vector',13,'R-S encode/out_1']],...
          'position',[425,60,445,80])

add_block('built-in/Inport',[sys,'/',['Binary-vector',13,'R-S encode/in_2']])
set_param([sys,'/',['Binary-vector',13,'R-S encode/in_2']],...
          'Port','2',...
          'position',[20,110,40,130])

add_block('built-in/Inport',[sys,'/',['Binary-vector',13,'R-S encode/in_1']])
set_param([sys,'/',['Binary-vector',13,'R-S encode/in_1']],...
          'position',[30,25,50,45])

add_block('built-in/Mux',[sys,'/',['Binary-vector',13,'R-S encode/Mux']])