ooknyqmf.m

很多通信系统的matlab程序,便于修改调用

function [ret,x0,str,ts,xts]=ooknyqmf(t,x,u,flag);
%OOKNYQMF	is the M-file description of the SIMULINK system named OOKNYQMF.
%	The block-diagram can be displayed by typing: OOKNYQMF.
%
%	SYS=OOKNYQMF(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 OOKNYQMF 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 OOKNYQMF with a FLAG of zero:
%	[SIZES]=OOKNYQMF([],[],[],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',[194,42,1259,320])
     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.000001')
set_param(sys,'Max step size', '10')
set_param(sys,'Relative error','1e-6')
set_param(sys,'Return vars',   '')

add_block('built-in/Gain',[sys,'/','Gain1'])
set_param([sys,'/','Gain1'],...
		'Gain','1/20',...
		'position',[810,112,835,138])


%     Subsystem  'Nyquist-Filter'.

new_system([sys,'/','Nyquist-Filter'])
set_param([sys,'/','Nyquist-Filter'],'Location',[131,545,589,761])


%     Subsystem  ['Nyquist-Filter/Pulses deivide',13,'same sample time1'].

new_system([sys,'/',['Nyquist-Filter/Pulses deivide',13,'same sample time1']])
set_param([sys,'/',['Nyquist-Filter/Pulses deivide',13,'same sample time1']],'Location',[55,192,283,286])

add_block('built-in/S-Function',[sys,'/',['Nyquist-Filter/Pulses deivide',13,'same sample time1/S-function1']])
set_param([sys,'/',['Nyquist-Filter/Pulses deivide',13,'same sample time1/S-function1']],...
		'function name','homopuls',...
		'parameters','samp_time, out_divid, offset',...
		'position',[40,31,105,59])

add_block('built-in/Outport',[sys,'/',['Nyquist-Filter/Pulses deivide',13,'same sample time1/Outport']])
set_param([sys,'/',['Nyquist-Filter/Pulses deivide',13,'same sample time1/Outport']],...
		'position',[160,35,180,55])
add_line([sys,'/',['Nyquist-Filter/Pulses deivide',13,'same sample time1']],[110,45;155,45])
set_param([sys,'/',['Nyquist-Filter/Pulses deivide',13,'same sample time1']],...
		'Mask Display','Vector\nPulse',...
		'Mask Type','Pulse generator')
set_param([sys,'/',['Nyquist-Filter/Pulses deivide',13,'same sample time1']],...
		'Mask Dialogue','Vector pulses with sample rate of each pulse as Sample_time ./Divider. |Sample time (scalar, sec):|Divider (integer vector):|Offset (have same dimension as Divider:')
set_param([sys,'/',['Nyquist-Filter/Pulses deivide',13,'same sample time1']],...
		'Mask Translate','samp_time=@1;out_divid=@2;offset=@3;')
set_param([sys,'/',['Nyquist-Filter/Pulses deivide',13,'same sample time1']],...
		'Mask Help','This block generator a vector pulses. The output length is the same as the vector defined in "Divider". The ith element of the output has the sample time Sample_time/Divider(i). All elements in "Divider" must be integers.')
set_param([sys,'/',['Nyquist-Filter/Pulses deivide',13,'same sample time1']],...
		'Mask Entries','td\/1\/off\/')


%     Finished composite block ['Nyquist-Filter/Pulses deivide',13,'same sample time1'].

set_param([sys,'/',['Nyquist-Filter/Pulses deivide',13,'same sample time1']],...
		'position',[35,71,95,109])

add_block('built-in/Inport',[sys,'/','Nyquist-Filter/in_1'])
set_param([sys,'/','Nyquist-Filter/in_1'],...
		'position',[35,30,55,50])

add_block('built-in/Outport',[sys,'/','Nyquist-Filter/out_1'])
set_param([sys,'/','Nyquist-Filter/out_1'],...
		'position',[330,40,350,60])

add_block('built-in/Filter',[sys,'/','Nyquist-Filter/Filter'])
set_param([sys,'/','Nyquist-Filter/Filter'],...
		'Numerator','num',...
		'Denominator','den',...
		'Sample time','ts',...
		'position',[255,31,310,69])


%     Subsystem  'Nyquist-Filter/Memory Impulse'.

new_system([sys,'/','Nyquist-Filter/Memory Impulse'])
set_param([sys,'/','Nyquist-Filter/Memory Impulse'],'Location',[541,520,939,745])

add_block('built-in/Inport',[sys,'/','Nyquist-Filter/Memory Impulse/in_2'])
set_param([sys,'/','Nyquist-Filter/Memory Impulse/in_2'],...
		'Port','2',...
		'position',[25,30,45,50])


%     Subsystem  ['Nyquist-Filter/Memory Impulse/discrete',13,'Rising edge',13,'detector'].

new_system([sys,'/',['Nyquist-Filter/Memory Impulse/discrete',13,'Rising edge',13,'detector']])
set_param([sys,'/',['Nyquist-Filter/Memory Impulse/discrete',13,'Rising edge',13,'detector']],'Location',[419,521,808,713])

add_block('built-in/Unit Delay',[sys,'/',['Nyquist-Filter/Memory Impulse/discrete',13,'Rising edge',13,'detector/Unit Delay']])
set_param([sys,'/',['Nyquist-Filter/Memory Impulse/discrete',13,'Rising edge',13,'detector/Unit Delay']],...
		'Sample time','ts',...
		'position',[120,87,170,103])

add_block('built-in/Outport',[sys,'/',['Nyquist-Filter/Memory Impulse/discrete',13,'Rising edge',13,'detector/out_1']])
set_param([sys,'/',['Nyquist-Filter/Memory Impulse/discrete',13,'Rising edge',13,'detector/out_1']],...
		'position',[340,25,360,45])

add_block('built-in/Relational Operator',[sys,'/',['Nyquist-Filter/Memory Impulse/discrete',13,'Rising edge',13,'detector/Relational',13,'Operator1']])
set_param([sys,'/',['Nyquist-Filter/Memory Impulse/discrete',13,'Rising edge',13,'detector/Relational',13,'Operator1']],...
		'position',[200,88,230,112])

add_block('built-in/Relational Operator',[sys,'/',['Nyquist-Filter/Memory Impulse/discrete',13,'Rising edge',13,'detector/Relational',13,'Operator2']])
set_param([sys,'/',['Nyquist-Filter/Memory Impulse/discrete',13,'Rising edge',13,'detector/Relational',13,'Operator2']],...
		'Operator','>',...
		'position',[265,23,295,47])

add_block('built-in/Constant',[sys,'/',['Nyquist-Filter/Memory Impulse/discrete',13,'Rising edge',13,'detector/one1']])
set_param([sys,'/',['Nyquist-Filter/Memory Impulse/discrete',13,'Rising edge',13,'detector/one1']],...
		'Value','thld',...
		'position',[30,135,60,155])

add_block('built-in/Inport',[sys,'/',['Nyquist-Filter/Memory Impulse/discrete',13,'Rising edge',13,'detector/in_1']])
set_param([sys,'/',['Nyquist-Filter/Memory Impulse/discrete',13,'Rising edge',13,'detector/in_1']],...
		'position',[10,15,30,35])

add_block('built-in/Relational Operator',[sys,'/',['Nyquist-Filter/Memory Impulse/discrete',13,'Rising edge',13,'detector/Relational',13,'Operator']])
set_param([sys,'/',['Nyquist-Filter/Memory Impulse/discrete',13,'Rising edge',13,'detector/Relational',13,'Operator']],...
		'position',[135,18,165,42])
add_line([sys,'/',['Nyquist-Filter/Memory Impulse/discrete',13,'Rising edge',13,'detector']],[35,25;130,25])
add_line([sys,'/',['Nyquist-Filter/Memory Impulse/discrete',13,'Rising edge',13,'detector']],[45,25;45,95;115,95])
add_line([sys,'/',['Nyquist-Filter/Memory Impulse/discrete',13,'Rising edge',13,'detector']],[175,95;195,95])
add_line([sys,'/',['Nyquist-Filter/Memory Impulse/discrete',13,'Rising edge',13,'detector']],[65,145;100,145;100,35;130,35])
add_line([sys,'/',['Nyquist-Filter/Memory Impulse/discrete',13,'Rising edge',13,'detector']],[100,145;180,145;180,105;195,105])
add_line([sys,'/',['Nyquist-Filter/Memory Impulse/discrete',13,'Rising edge',13,'detector']],[170,30;260,30])
add_line([sys,'/',['Nyquist-Filter/Memory Impulse/discrete',13,'Rising edge',13,'detector']],[235,100;240,100;240,40;260,40])
add_line([sys,'/',['Nyquist-Filter/Memory Impulse/discrete',13,'Rising edge',13,'detector']],[300,35;335,35])
set_param([sys,'/',['Nyquist-Filter/Memory Impulse/discrete',13,'Rising edge',13,'detector']],...
		'Mask Display','Edge\nDetector',...
		'Mask Type','Edge Detection',...
		'Mask Dialogue','Rising edge detection:|Threshold:|Sampling time:')
set_param([sys,'/',['Nyquist-Filter/Memory Impulse/discrete',13,'Rising edge',13,'detector']],...
		'Mask Translate','thld=@1;ts=@2;')
set_param([sys,'/',['Nyquist-Filter/Memory Impulse/discrete',13,'Rising edge',13,'detector']],...
		'Mask Help','When the input signal rising edge passed the given shreshold, this block outputs one, otherwise outputs zero.')
set_param([sys,'/',['Nyquist-Filter/Memory Impulse/discrete',13,'Rising edge',13,'detector']],...
		'Mask Entries','thr\/ts\/')


%     Finished composite block ['Nyquist-Filter/Memory Impulse/discrete',13,'Rising edge',13,'detector'].

set_param([sys,'/',['Nyquist-Filter/Memory Impulse/discrete',13,'Rising edge',13,'detector']],...
		'position',[65,26,120,54])

add_block('built-in/Switch',[sys,'/','Nyquist-Filter/Memory Impulse/Switch'])
set_param([sys,'/','Nyquist-Filter/Memory Impulse/Switch'],...
		'Threshold','.5',...
		'position',[295,124,320,156])

add_block('built-in/Constant',[sys,'/','Nyquist-Filter/Memory Impulse/Constant'])
set_param([sys,'/','Nyquist-Filter/Memory Impulse/Constant'],...
		'Value','x',...
		'position',[210,155,230,175])

add_block('built-in/Outport',[sys,'/','Nyquist-Filter/Memory Impulse/out_1'])
set_param([sys,'/','Nyquist-Filter/Memory Impulse/out_1'],...
		'position',[355,130,375,150])

add_block('built-in/Product',[sys,'/','Nyquist-Filter/Memory Impulse/Product'])
set_param([sys,'/','Nyquist-Filter/Memory Impulse/Product'],...
		'position',[155,90,180,110])

add_block('built-in/Zero-Order Hold',[sys,'/',['Nyquist-Filter/Memory Impulse/Zero-Order',13,'Hold']])
set_param([sys,'/',['Nyquist-Filter/Memory Impulse/Zero-Order',13,'Hold']],...
		'Sample time','ts',...
		'position',[205,84,240,116])

add_block('built-in/Inport',[sys,'/','Nyquist-Filter/Memory Impulse/in_1'])
set_param([sys,'/','Nyquist-Filter/Memory Impulse/in_1'],...
		'position',[35,185,55,205])
add_line([sys,'/','Nyquist-Filter/Memory Impulse'],[50,40;60,40])
add_line([sys,'/','Nyquist-Filter/Memory Impulse'],[125,40;130,40;130,95;150,95])
add_line([sys,'/','Nyquist-Filter/Memory Impulse'],[235,165;250,165;250,140;290,140])
add_line([sys,'/','Nyquist-Filter/Memory Impulse'],[325,140;350,140])
add_line([sys,'/','Nyquist-Filter/Memory Impulse'],[60,195;280,195;290,150])
add_line([sys,'/','Nyquist-Filter/Memory Impulse'],[95,195;95,105;150,105])
add_line([sys,'/','Nyquist-Filter/Memory Impulse'],[185,100;200,100])
add_line([sys,'/','Nyquist-Filter/Memory Impulse'],[245,100;260,100;260,130;290,130])
set_param([sys,'/','Nyquist-Filter/Memory Impulse'],...
		'Mask Display','1/sinc',...
		'Mask Type','Impulse',...
		'Mask Dialogue','1/sinc filter:|Holding time:|Threshold for clock pulse: ',...
		'Mask Translate','ts=@1;thr=@2;x=thr>0;')
set_param([sys,'/','Nyquist-Filter/Memory Impulse'],...
		'Mask Help','When the threshold is greater than zero, this block outputs an impulse signal, which begins at the raising edge of the second input and lasts the holding time. The impulse holds a constant amplitude, which equals the amplitude of the first input at the begining of the impulse. A clock raising edge will trigger next impulse. When the threshold is less or equal to zero, this block faithfully outputs the first input.')
set_param([sys,'/','Nyquist-Filter/Memory Impulse'],...
		'Mask Entries','ts\/x\/')


%     Finished composite block 'Nyquist-Filter/Memory Impulse'.

set_param([sys,'/','Nyquist-Filter/Memory Impulse'],...
		'position',[155,30,225,70])
add_line([sys,'/','Nyquist-Filter'],[315,50;325,50])
add_line([sys,'/','Nyquist-Filter'],[60,40;150,40])
add_line([sys,'/','Nyquist-Filter'],[230,50;250,50])
add_line([sys,'/','Nyquist-Filter'],[100,90;115,90;115,60;150,60])
set_param([sys,'/','Nyquist-Filter'],...
		'Mask Display','R-C\nfilter',...
		'Mask Type','R-C filter')
set_param([sys,'/','Nyquist-Filter'],...
		'Mask Dialogue','Design a raised cosine filter and use it to filter the input signal.|Symbol interval (sec):|Computation sampling time (sec):|Filter rolloff factor:|Delay step (integer to multiply symbol interval):|Filter type (FIR or IIR) / (Normal or Sqrt):|With or without sinc filter (1 or 0, suggest 1):')
set_param([sys,'/','Nyquist-Filter'],...
		'Mask Translate','td=@1;if length(td)>1,off=td(2);td=td(1);else,off=0;end;ts=@2;R=@3;delay=@4;fil=@5;[num,den]=rcosine(1/td(1),1/ts(1),fil,R,delay);x=@6/2;')
set_param([sys,'/','Nyquist-Filter'],...
		'Mask Help','This block designs a raised cosine filter using the specified parameter. The filter is used to filter the input signal. You can optionally use the sinc filter in the filtering computation.')
set_param([sys,'/','Nyquist-Filter'],...
		'Mask Entries','1\/.05\/.5\/6\/''FIR''\/0\/')


%     Finished composite block 'Nyquist-Filter'.

set_param([sys,'/','Nyquist-Filter'],...
		'position',[710,103,790,147])

add_block('built-in/Scope',[sys,'/','x_e(t)'])