tstqask.m

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

function [ret,x0,str,ts,xts]=tstqask(t,x,u,flag);
%TSTQASK       is the M-file description of the SIMULINK system named TSTQASK.
%        The block-diagram can be displayed by typing: TSTQASK.
%
%        SYS=TSTQASK(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 TSTQASK 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 TSTQASK with a FLAG of zero:
%   [SIZES]=TSTQASK([],[],[],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',[2,429,830,856])
     open_system(sys)
end;
set_param(sys,'algorithm',     'RK-45')
set_param(sys,'Start time',    '0.0')
set_param(sys,'Stop time',     '10')
set_param(sys,'Min step size', '0.0001')
set_param(sys,'Max step size', '.003')
set_param(sys,'Relative error','1e-3')
set_param(sys,'Return vars',   '')
set_param(sys,'Load callback','N=16;T=.1;Fc=100;');


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

new_system([sys,'/',['Random-number',13,'generator']])
set_param([sys,'/',['Random-number',13,'generator']],'Location',[115,400,471,569])

add_block('built-in/White Noise',[sys,'/',['Random-number',13,'generator/Random',13,'Number']])
set_param([sys,'/',['Random-number',13,'generator/Random',13,'Number']],...
              'Seed','seed',...
               'position',[40,27,85,63])

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

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

add_block('built-in/State-Space',[sys,'/',['Random-number',13,'generator/Matrix',13,'Gain']])
set_param([sys,'/',['Random-number',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-number',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-number',13,'generator/Matrix',13,'Gain']],...
                'position',[230,35,275,65])

add_block('built-in/Constant',[sys,'/',['Random-number',13,'generator/Constant']])
set_param([sys,'/',['Random-number',13,'generator/Constant']],...
          'Value','x',...
         'position',[65,105,95,125])
add_line([sys,'/',['Random-number',13,'generator']],[280,50;305,50])
add_line([sys,'/',['Random-number',13,'generator']],[190,50;225,50])
add_line([sys,'/',['Random-number',13,'generator']],[100,115;115,115;115,55;150,55])
add_line([sys,'/',['Random-number',13,'generator']],[90,45;150,45])
set_param([sys,'/',['Random-number',13,'generator']],...
              'Mask Display','Random\nNumber',...
             'Mask Type','Random number generator',...
               'Mask Dialogue','Random number generator:|M-ary number:|Seed:')
set_param([sys,'/',['Random-number',13,'generator']],...
              'Mask Translate','N=@1;x=[(2-N)/N:2/N:(N-2)/N];seed=@2;')
set_param([sys,'/',['Random-number',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-number',13,'generator']],...
            'Mask Entries','N\/4432\/')


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

set_param([sys,'/',['Random-number',13,'generator']],...
             'BackGround',5,...
              'position',[50,30,110,70])


%     Subsystem  ['Clock',13,'pulse'].

new_system([sys,'/',['Clock',13,'pulse']])
set_param([sys,'/',['Clock',13,'pulse']],'Location',[272,165,909,406])

add_block('built-in/Unit Delay',[sys,'/',['Clock',13,'pulse/Unit Delay1']])
set_param([sys,'/',['Clock',13,'pulse/Unit Delay1']],...
             'orientation',2,...
             'Sample time','[Ts/2,st2]',...
          'x0','ini',...
          'position',[140,165,190,185])

add_block('built-in/Unit Delay',[sys,'/',['Clock',13,'pulse/Unit Delay']])
set_param([sys,'/',['Clock',13,'pulse/Unit Delay']],...
          'orientation',2,...
             'Sample time','[Ts/2,st1]',...
          'position',[145,75,195,95])

add_block('built-in/Gain',[sys,'/',['Clock',13,'pulse/Gain1']])
set_param([sys,'/',['Clock',13,'pulse/Gain1']],...
            'Gain','g',...
          'position',[235,122,255,138])

add_block('built-in/Constant',[sys,'/',['Clock',13,'pulse/Constant']])
set_param([sys,'/',['Clock',13,'pulse/Constant']],...
                'position',[25,30,45,50])

add_block('built-in/Saturation',[sys,'/',['Clock',13,'pulse/Saturation']])
set_param([sys,'/',['Clock',13,'pulse/Saturation']],...
              'Lower Limit','0',...
           'Upper Limit','1',...
           'position',[365,40,390,60])

add_block('built-in/Sum',[sys,'/',['Clock',13,'pulse/Sum2']])
set_param([sys,'/',['Clock',13,'pulse/Sum2']],...
               'hide name',0,...
               'inputs','+-',...
               'position',[300,40,320,60])

add_block('built-in/Sum',[sys,'/',['Clock',13,'pulse/Sum']])
set_param([sys,'/',['Clock',13,'pulse/Sum']],...
         'hide name',0,...
               'inputs','+-',...
               'position',[135,35,155,55])

add_block('built-in/Sum',[sys,'/',['Clock',13,'pulse/Sum1']])
set_param([sys,'/',['Clock',13,'pulse/Sum1']],...
               'hide name',0,...
               'inputs','+-',...
               'position',[130,120,150,140])

add_block('built-in/Relational Operator',[sys,'/',['Clock',13,'pulse/Relational',13,'operator']])
set_param([sys,'/',['Clock',13,'pulse/Relational',13,'operator']],...
             'hide name',0,...
               'position',[370,102,400,133])

add_block('built-in/Gain',[sys,'/',['Clock',13,'pulse/Gain']])
set_param([sys,'/',['Clock',13,'pulse/Gain']],...
            'hide name',0,...
               'Gain','ht',...
         'position',[540,102,565,128])

add_block('built-in/Product',[sys,'/',['Clock',13,'pulse/Product']])
set_param([sys,'/',['Clock',13,'pulse/Product']],...
           'hide name',0,...
               'position',[470,105,495,125])

add_block('built-in/Clock',[sys,'/',['Clock',13,'pulse/Clock1']])
set_param([sys,'/',['Clock',13,'pulse/Clock1']],...
               'hide name',0,...
               'position',[295,100,315,120])

add_block('built-in/Constant',[sys,'/',['Clock',13,'pulse/Constant1']])
set_param([sys,'/',['Clock',13,'pulse/Constant1']],...
              'Value','stt',...
               'position',[295,149,320,171])

add_block('built-in/Outport',[sys,'/',['Clock',13,'pulse/out_1']])
set_param([sys,'/',['Clock',13,'pulse/out_1']],...
               'hide name',0,...
               'position',[590,105,610,125])
add_line([sys,'/',['Clock',13,'pulse']],[260,130;275,130;275,55;295,55])
add_line([sys,'/',['Clock',13,'pulse']],[155,130;230,130])
add_line([sys,'/',['Clock',13,'pulse']],[325,50;360,50])
add_line([sys,'/',['Clock',13,'pulse']],[160,45;295,45])
add_line([sys,'/',['Clock',13,'pulse']],[50,40;70,40;70,125;125,125])
add_line([sys,'/',['Clock',13,'pulse']],[135,175;100,175;100,135;125,135])
add_line([sys,'/',['Clock',13,'pulse']],[155,130;210,130;210,175;195,175])
add_line([sys,'/',['Clock',13,'pulse']],[50,40;130,40])
add_line([sys,'/',['Clock',13,'pulse']],[160,45;230,45;230,85;200,85])
add_line([sys,'/',['Clock',13,'pulse']],[140,85;105,85;105,50;130,50])
add_line([sys,'/',['Clock',13,'pulse']],[320,110;365,110])
add_line([sys,'/',['Clock',13,'pulse']],[405,120;465,120])
add_line([sys,'/',['Clock',13,'pulse']],[500,115;535,115])
add_line([sys,'/',['Clock',13,'pulse']],[570,115;585,115])
add_line([sys,'/',['Clock',13,'pulse']],[325,160;345,160;345,125;365,125])
add_line([sys,'/',['Clock',13,'pulse']],[395,50;430,50;430,110;465,110])
set_param([sys,'/',['Clock',13,'pulse']],...
               'Mask Display','plot(0,0,100,100,[90,75,75,60,60,35,35,20,20,10],[20,20,80,80,20,20,80,80,20,20])',...
          'Mask Type','Pulse Generator')
set_param([sys,'/',['Clock',13,'pulse']],...
              'Mask Dialogue','Pulse Generator.|Pulse period (secs):|Pulse width:|Pulse height:|Pulse start time:')
set_param([sys,'/',['Clock',13,'pulse']],...
               'Mask Translate','Ts=@1; du=@2; ht=@3; stt=@4; g=1; ini=ones(length(@3),1); if (Ts<2*du), g=-1; du=du-Ts/2; ini=zeros(length(@3),1); end; st1=rem(stt,Ts); st2=rem(stt+du,Ts);')
set_param([sys,'/',['Clock',13,'pulse']],...
            'Mask Help','Pulse generator which ensures pulse\ntransitions are hit. Provides a vector of pulses when the height is entered as a vector.\nUnmask to see how it works.')
set_param([sys,'/',['Clock',13,'pulse']],...
           'Mask Entries','T\/.05\/1\/0\/')


%     Finished composite block ['Clock',13,'pulse'].

set_param([sys,'/',['Clock',13,'pulse']],...
                'BackGround',4,...
              'position',[65,267,100,303])


%     Subsystem  ['Encoding only',13,'no-modu'].

new_system([sys,'/',['Encoding only',13,'no-modu']])