com_sour.m

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

function [ret,x0,str,ts,xts]=com_sour(t,x,u,flag);
%COM_SOUR	is the M-file description of the SIMULINK system named COM_SOUR.
%	The block-diagram can be displayed by typing: COM_SOUR.
%
%	SYS=COM_SOUR(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 COM_SOUR 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 COM_SOUR with a FLAG of zero:
%	[SIZES]=COM_SOUR([],[],[],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',[271,85,1003,719])
     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,'/','Source//Sink Library'])
set_param([sys,'/','Source//Sink Library'],...
		'ForeGround',3,...
		'Font Weight','bold',...
		'Font Size',12,...
		'position',[325,5,330,10])


%     Subsystem  ['Triggered write',13,'to file'].

new_system([sys,'/',['Triggered write',13,'to file']])
set_param([sys,'/',['Triggered write',13,'to file']],'Location',[105,608,345,754])

add_block('built-in/Inport',[sys,'/',['Triggered write',13,'to file/in_2']])
set_param([sys,'/',['Triggered write',13,'to file/in_2']],...
		'Port','2',...
		'position',[15,75,35,95])

add_block('built-in/Inport',[sys,'/',['Triggered write',13,'to file/in_1']])
set_param([sys,'/',['Triggered write',13,'to file/in_1']],...
		'position',[15,60,35,80])

add_block('built-in/Mux',[sys,'/',['Triggered write',13,'to file/Mux']])
set_param([sys,'/',['Triggered write',13,'to file/Mux']],...
		'inputs','2',...
		'position',[65,61,100,94])

add_block('built-in/S-Function',[sys,'/',['Triggered write',13,'to file/S-function']])
set_param([sys,'/',['Triggered write',13,'to file/S-function']],...
		'function name','swritfil',...
		'parameters','filename, datatype, nump, threshold',...
		'position',[135,70,185,90])
add_line([sys,'/',['Triggered write',13,'to file']],[40,85;60,85])
add_line([sys,'/',['Triggered write',13,'to file']],[40,70;60,70])
add_line([sys,'/',['Triggered write',13,'to file']],[105,80;130,80])
set_param([sys,'/',['Triggered write',13,'to file']],...
		'Mask Display','plot(0,0,100,100,x,y,v,w);Wt file',...
		'Mask Type','Write to file')
set_param([sys,'/',['Triggered write',13,'to file']],...
		'Mask Dialogue','The rising edge of the signal from the\n2nd port triggers the block to write\nanother row into the specified file.|File name (string, include extension):|Data type (''ascii'', ''float'', ''integer''):|Number of trigger pulses between saved data:|Threshold in detecting trigger signal:')
set_param([sys,'/',['Triggered write',13,'to file']],...
		'Mask Translate','filename=@1;datatype=@2;nump=@3;threshold=@4;[x,y]=srcsicon(8);[v,w]=trigicon(0,25,get_param(gcb,''orientation''),2);')
set_param([sys,'/',['Triggered write',13,'to file']],...
		'Mask Help','This block writes the data inputted from the first port to the next row of the specified file at the rising edge of the second input trigger pulse. The first pulse triggers to write the first record. You can specify the number of pulses between writting records. When you specify ASCII format, there is no return added between two records.')
set_param([sys,'/',['Triggered write',13,'to file']],...
		'Mask Entries','''junk.dat''\/''ascii''\/0\/0.5\/')


%     Finished composite block ['Triggered write',13,'to file'].

set_param([sys,'/',['Triggered write',13,'to file']],...
		'position',[175,113,255,157])


%     Subsystem  ['QASK Mo//Dem',13,'demo4'].

new_system([sys,'/',['QASK Mo//Dem',13,'demo4']])
set_param([sys,'/',['QASK Mo//Dem',13,'demo4']],'Location',[175,552,335,697])
set_param([sys,'/',['QASK Mo//Dem',13,'demo4']],...
		'Mask Display','Triggered\nfile I/O\ndemo',...
		'Mask Dialogue','eval(''junk=magic(5);ttrigfil;save(''''junk.tmp'''',''''junk'''',''''-ascii'''')'')')


%     Finished composite block ['QASK Mo//Dem',13,'demo4'].

set_param([sys,'/',['QASK Mo//Dem',13,'demo4']],...
		'BackGround',7,...
		'hide name',0,...
		'position',[315,113,395,157])


%     Subsystem  ['Eye-diagram',13,'scatter plot'].

new_system([sys,'/',['Eye-diagram',13,'scatter plot']])
set_param([sys,'/',['Eye-diagram',13,'scatter plot']],'Location',[165,242,460,428])

add_block('built-in/Inport',[sys,'/',['Eye-diagram',13,'scatter plot/in_1']])
set_param([sys,'/',['Eye-diagram',13,'scatter plot/in_1']],...
		'position',[15,55,35,75])

add_block('built-in/Inport',[sys,'/',['Eye-diagram',13,'scatter plot/in_2']])
set_param([sys,'/',['Eye-diagram',13,'scatter plot/in_2']],...
		'Port','2',...
		'position',[20,115,40,135])

add_block('built-in/Mux',[sys,'/',['Eye-diagram',13,'scatter plot/Mux']])
set_param([sys,'/',['Eye-diagram',13,'scatter plot/Mux']],...
		'inputs','2',...
		'position',[65,70,100,105])

add_block('built-in/S-Function',[sys,'/',['Eye-diagram',13,'scatter plot/Plot1']])
set_param([sys,'/',['Eye-diagram',13,'scatter plot/Plot1']],...
		'function name','eyediasi',...
		'parameters','time_range, boundary, kept_length, eye_line, scatter_line, two_d_line',...
		'position',[130,70,200,110])
add_line([sys,'/',['Eye-diagram',13,'scatter plot']],[40,65;45,65;45,80;60,80])
add_line([sys,'/',['Eye-diagram',13,'scatter plot']],[45,125;50,125;60,95])
add_line([sys,'/',['Eye-diagram',13,'scatter plot']],[105,90;125,90])
set_param([sys,'/',['Eye-diagram',13,'scatter plot']],...
		'Mask Display','plot(0,0,100,100,x,y,v,w);',...
		'Mask Type','Eye-diagram/scatter plot')
set_param([sys,'/',['Eye-diagram',13,'scatter plot']],...
		'Mask Dialogue','Produce eye-pattern diagram/scatter\nplot/x-y plot with the rising edge of the second input as the decision line.|Eye-pattern plot time frame:|Lower and upper bound of in-coming signal (1x2 vector):|Number of saved traces:|Line-type for eye-pattern diagram (0 for no plot):|Line-type for scatter plot (0 for no plot):|Line-type for x-y plot (0 for no plot):')
set_param([sys,'/',['Eye-diagram',13,'scatter plot']],...
		'Mask Translate','time_range=@1;boundary=@2;kept_length=@3;eye_line=@4;scatter_line=@5;two_d_line=@6;[x,y]=srcsicon(10);[v,w]=trigicon(0,25,get_param(gcb,''orientation''),2);')
set_param([sys,'/',['Eye-diagram',13,'scatter plot']],...
		'Mask Help','This plot opens one, two, or three plots by assigning the line type parameters. The scatter plot point is selected by the rising edge of a trigger signal. A decision line is drawn at the rising edge of the trigger signal.')
set_param([sys,'/',['Eye-diagram',13,'scatter plot']],...
		'Mask Entries','pi\/[-1.1 1.1]\/10\/''r-/g--/b-.''\/''r+/gx/b*''\/''r-/g--/b-.''\/')


%     Finished composite block ['Eye-diagram',13,'scatter plot'].

set_param([sys,'/',['Eye-diagram',13,'scatter plot']],...
		'position',[175,188,255,232])


%     Subsystem  ['QASK Mo//Dem',13,'demo3'].

new_system([sys,'/',['QASK Mo//Dem',13,'demo3']])
set_param([sys,'/',['QASK Mo//Dem',13,'demo3']],'Location',[175,552,335,697])
set_param([sys,'/',['QASK Mo//Dem',13,'demo3']],...
		'Mask Display','Eye-pattern\ndemo',...
		'Mask Dialogue','eval(''ts=.01;Fc=100/pi/2;[num,den]=butter(2,Fc*2*ts);td=.1;M=16;tsteyesc'')')


%     Finished composite block ['QASK Mo//Dem',13,'demo3'].

set_param([sys,'/',['QASK Mo//Dem',13,'demo3']],...
		'BackGround',7,...
		'hide name',0,...
		'position',[315,188,395,232])


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

new_system([sys,'/',['Random-integer',13,'generator']])
set_param([sys,'/',['Random-integer',13,'generator']],'Location',[98,518,454,687])

add_block('built-in/S-Function',[sys,'/',['Random-integer',13,'generator/Uniform',13,'number']])
set_param([sys,'/',['Random-integer',13,'generator/Uniform',13,'number']],...
		'function name','srandint',...
		'parameters','seed,mul',...
		'position',[65,45,135,85])

add_block('built-in/Outport',[sys,'/',['Random-integer',13,'generator/out_1']])
set_param([sys,'/',['Random-integer',13,'generator/out_1']],...
		'position',[165,55,185,75])
add_line([sys,'/',['Random-integer',13,'generator']],[140,65;160,65])
set_param([sys,'/',['Random-integer',13,'generator']],...
		'Mask Display','plot(0,0,100,100,x,y);Random int',...
		'Mask Type','Random integer generator')
set_param([sys,'/',['Random-integer',13,'generator']],...
		'Mask Dialogue','Generate integers randomly distributed in range [0, M-1], where M is the M-ary number.|M-ary number:|Seed:',...
		'Mask Translate','seed=@2;mul=@1;[x,y]=srcsicon(4);')
set_param([sys,'/',['Random-integer',13,'generator']],...
		'Mask Help','This block generates an integer number in the range [0 M-1], where M is the M-ary number. The probability of the appearance for each number is equal.')
set_param([sys,'/',['Random-integer',13,'generator']],...
		'Mask Entries','8\/[1 2]\/')


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

set_param([sys,'/',['Random-integer',13,'generator']],...
		'position',[455,188,535,232])


%     Subsystem  ['QASK Mo//Dem',13,'demo8'].

new_system([sys,'/',['QASK Mo//Dem',13,'demo8']])
set_param([sys,'/',['QASK Mo//Dem',13,'demo8']],'Location',[175,552,335,697])
set_param([sys,'/',['QASK Mo//Dem',13,'demo8']],...
		'Mask Display','Gaussian\nrandom\ndemo',...
		'Mask Dialogue','eval(''td=.1;tgrand'')')


%     Finished composite block ['QASK Mo//Dem',13,'demo8'].

set_param([sys,'/',['QASK Mo//Dem',13,'demo8']],...
		'BackGround',7,...
		'hide name',0,...
		'position',[595,113,675,157])


%     Subsystem  ['QASK Mo//Dem',13,'demo9'].

new_system([sys,'/',['QASK Mo//Dem',13,'demo9']])
set_param([sys,'/',['QASK Mo//Dem',13,'demo9']],'Location',[175,552,335,697])
set_param([sys,'/',['QASK Mo//Dem',13,'demo9']],...
		'Mask Display','Poisson\nrandom\ndemo',...
		'Mask Dialogue','eval(''td=.1;tprand'')')


%     Finished composite block ['QASK Mo//Dem',13,'demo9'].

set_param([sys,'/',['QASK Mo//Dem',13,'demo9']],...
		'BackGround',7,...
		'hide name',0,...
		'position',[595,263,675,307])