hdslsimu.m

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

function [ret,x0,str,ts,xts]=hdslsimu(t,x,u,flag);
%HDSLSIMU	is the M-file description of the SIMULINK system named HDSLSIMU.
%	The block-diagram can be displayed by typing: HDSLSIMU.
%
%	SYS=HDSLSIMU(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 HDSLSIMU 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 HDSLSIMU with a FLAG of zero:
%	[SIZES]=HDSLSIMU([],[],[],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',[95,335,595,590])
     open_system(sys)
end;
set_param(sys,'algorithm',     'RK-45')
set_param(sys,'Start time',    '0.0')
set_param(sys,'Stop time',     '200')
set_param(sys,'Min step size', '0.002')
set_param(sys,'Max step size', '10')
set_param(sys,'Relative error','1e-3')
set_param(sys,'Return vars',   '')
set_param(sys,'Load callback','hdsldata');

add_block('built-in/Note',[sys,'/','High-Bit-Rate Digital Subscriber Line'])
set_param([sys,'/','High-Bit-Rate Digital Subscriber Line'],...
		'ForeGround',2,...
		'Font Weight','bold',...
		'Font Size',12,...
		'position',[245,20,250,25])


%     Subsystem  ['K-step',13,'delay'].

new_system([sys,'/',['K-step',13,'delay']])
set_param([sys,'/',['K-step',13,'delay']],'Location',[12,432,357,702])

add_block('built-in/Digital Clock',[sys,'/',['K-step',13,'delay/Digital Clock']])
set_param([sys,'/',['K-step',13,'delay/Digital Clock']],...
		'Sample time','ts',...
		'position',[60,140,95,170])

add_block('built-in/Inport',[sys,'/',['K-step',13,'delay/in_1']])
set_param([sys,'/',['K-step',13,'delay/in_1']],...
		'position',[20,80,40,100])

add_block('built-in/Outport',[sys,'/',['K-step',13,'delay/out_1']])
set_param([sys,'/',['K-step',13,'delay/out_1']],...
		'position',[325,120,345,140])

add_block('built-in/Filter',[sys,'/',['K-step',13,'delay/K-step',13,'delay']])
set_param([sys,'/',['K-step',13,'delay/K-step',13,'delay']],...
		'Numerator','[zeros(1,n), 1]',...
		'Denominator','1',...
		'Sample time','ts',...
		'position',[60,70,185,110])

add_block('built-in/Switch',[sys,'/',['K-step',13,'delay/Switch']])
set_param([sys,'/',['K-step',13,'delay/Switch']],...
		'Threshold','.1',...
		'position',[270,114,295,146])

add_block('built-in/Constant',[sys,'/',['K-step',13,'delay/Constant']])
set_param([sys,'/',['K-step',13,'delay/Constant']],...
		'Value','ini',...
		'position',[120,230,140,250])

add_block('built-in/Constant',[sys,'/',['K-step',13,'delay/Constant1']])
set_param([sys,'/',['K-step',13,'delay/Constant1']],...
		'Value','n*ts(1)+ts(2)',...
		'position',[15,194,100,216])

add_block('built-in/Relational Operator',[sys,'/',['K-step',13,'delay/Relational',13,'Operator']])
set_param([sys,'/',['K-step',13,'delay/Relational',13,'Operator']],...
		'position',[150,168,180,192])
add_line([sys,'/',['K-step',13,'delay']],[300,130;320,130])
add_line([sys,'/',['K-step',13,'delay']],[45,90;55,90])
add_line([sys,'/',['K-step',13,'delay']],[190,90;195,90;195,120;265,120])
add_line([sys,'/',['K-step',13,'delay']],[145,240;240,240;240,140;265,140])
add_line([sys,'/',['K-step',13,'delay']],[100,155;115,155;115,175;145,175])
add_line([sys,'/',['K-step',13,'delay']],[105,205;115,205;115,185;145,185])
add_line([sys,'/',['K-step',13,'delay']],[185,180;220,180;220,130;265,130])
set_param([sys,'/',['K-step',13,'delay']],...
		'Mask Display','    -k\nz   ',...
		'Mask Type','Multi-step delay',...
		'Mask Dialogue','Output is K-step delay of the input signal.|Delay steps (integer):|Sample time (sec):|Initial output:')
set_param([sys,'/',['K-step',13,'delay']],...
		'Mask Translate','n=@1;ts=@2;if length(ts)<2, ts=[ts 0];end;ini=@3;')
set_param([sys,'/',['K-step',13,'delay']],...
		'Mask Help','This block holds the input signal for K steps and then outputs the signal. The delay step cannot be changed during the simulation.',...
		'Mask Entries','43\/.0025\/0\/')


%     Finished composite block ['K-step',13,'delay'].

set_param([sys,'/',['K-step',13,'delay']],...
		'ForeGround',2,...
		'hide name',0,...
		'position',[145,69,190,101])


%     Subsystem  ['K-step',13,'delay1'].

new_system([sys,'/',['K-step',13,'delay1']])
set_param([sys,'/',['K-step',13,'delay1']],'Location',[12,432,357,702])

add_block('built-in/Relational Operator',[sys,'/',['K-step',13,'delay1/Relational',13,'Operator']])
set_param([sys,'/',['K-step',13,'delay1/Relational',13,'Operator']],...
		'position',[150,168,180,192])

add_block('built-in/Constant',[sys,'/',['K-step',13,'delay1/Constant1']])
set_param([sys,'/',['K-step',13,'delay1/Constant1']],...
		'Value','n*ts(1)+ts(2)',...
		'position',[15,194,100,216])

add_block('built-in/Constant',[sys,'/',['K-step',13,'delay1/Constant']])
set_param([sys,'/',['K-step',13,'delay1/Constant']],...
		'Value','ini',...
		'position',[120,230,140,250])

add_block('built-in/Switch',[sys,'/',['K-step',13,'delay1/Switch']])
set_param([sys,'/',['K-step',13,'delay1/Switch']],...
		'Threshold','.1',...
		'position',[270,114,295,146])

add_block('built-in/Filter',[sys,'/',['K-step',13,'delay1/K-step',13,'delay']])
set_param([sys,'/',['K-step',13,'delay1/K-step',13,'delay']],...
		'Numerator','[zeros(1,n), 1]',...
		'Denominator','1',...
		'Sample time','ts',...
		'position',[60,70,185,110])

add_block('built-in/Outport',[sys,'/',['K-step',13,'delay1/out_1']])
set_param([sys,'/',['K-step',13,'delay1/out_1']],...
		'position',[325,120,345,140])

add_block('built-in/Inport',[sys,'/',['K-step',13,'delay1/in_1']])
set_param([sys,'/',['K-step',13,'delay1/in_1']],...
		'position',[20,80,40,100])

add_block('built-in/Digital Clock',[sys,'/',['K-step',13,'delay1/Digital Clock']])
set_param([sys,'/',['K-step',13,'delay1/Digital Clock']],...
		'Sample time','ts',...
		'position',[60,140,95,170])
add_line([sys,'/',['K-step',13,'delay1']],[185,180;220,180;220,130;265,130])
add_line([sys,'/',['K-step',13,'delay1']],[105,205;115,205;115,185;145,185])
add_line([sys,'/',['K-step',13,'delay1']],[100,155;115,155;115,175;145,175])
add_line([sys,'/',['K-step',13,'delay1']],[145,240;240,240;240,140;265,140])
add_line([sys,'/',['K-step',13,'delay1']],[190,90;195,90;195,120;265,120])
add_line([sys,'/',['K-step',13,'delay1']],[45,90;55,90])
add_line([sys,'/',['K-step',13,'delay1']],[300,130;320,130])
set_param([sys,'/',['K-step',13,'delay1']],...
		'Mask Display','    -k\nz   ',...
		'Mask Type','Multi-step delay',...
		'Mask Dialogue','Output is K-step delay of the input signal.|Delay steps (integer):|Sample time (sec):|Initial output:')
set_param([sys,'/',['K-step',13,'delay1']],...
		'Mask Translate','n=@1;ts=@2;if length(ts)<2, ts=[ts 0];end;ini=@3;')
set_param([sys,'/',['K-step',13,'delay1']],...
		'Mask Help','This block holds the input signal for K steps and then outputs the signal. The delay step cannot be changed during the simulation.',...
		'Mask Entries','43\/.0025\/0\/')


%     Finished composite block ['K-step',13,'delay1'].

set_param([sys,'/',['K-step',13,'delay1']],...
		'ForeGround',2,...
		'hide name',0,...
		'position',[160,194,205,226])


%     Subsystem  'PAM Signal'.

new_system([sys,'/','PAM Signal'])
set_param([sys,'/','PAM Signal'],'Location',[113,193,448,317])


%     Subsystem  'PAM Signal/MASK map'.

new_system([sys,'/','PAM Signal/MASK map'])
set_param([sys,'/','PAM Signal/MASK map'],'Location',[520,162,902,335])

add_block('built-in/Constant',[sys,'/','PAM Signal/MASK map/Gain diff'])
set_param([sys,'/','PAM Signal/MASK map/Gain diff'],...
		'Value','Ac',...
		'position',[180,96,215,124])

add_block('built-in/Sum',[sys,'/','PAM Signal/MASK map/Sum1'])
set_param([sys,'/','PAM Signal/MASK map/Sum1'],...
		'inputs','+-',...
		'position',[250,63,270,87])

add_block('built-in/Gain',[sys,'/','PAM Signal/MASK map/Gain'])
set_param([sys,'/','PAM Signal/MASK map/Gain'],...
		'Gain','Kc',...
		'position',[180,25,220,55])

add_block('built-in/Inport',[sys,'/','PAM Signal/MASK map/in_1'])
set_param([sys,'/','PAM Signal/MASK map/in_1'],...
		'position',[40,30,60,50])

add_block('built-in/Zero-Order Hold',[sys,'/',['PAM Signal/MASK map/Zero-Order',13,'Hold']])
set_param([sys,'/',['PAM Signal/MASK map/Zero-Order',13,'Hold']],...
		'Sample time','td',...
		'position',[100,24,135,56])

add_block('built-in/Outport',[sys,'/','PAM Signal/MASK map/out_1'])
set_param([sys,'/','PAM Signal/MASK map/out_1'],...
		'position',[315,65,335,85])
add_line([sys,'/','PAM Signal/MASK map'],[140,40;175,40])
add_line([sys,'/','PAM Signal/MASK map'],[65,40;95,40])
add_line([sys,'/','PAM Signal/MASK map'],[275,75;310,75])
add_line([sys,'/','PAM Signal/MASK map'],[220,110;230,110;230,80;245,80])
add_line([sys,'/','PAM Signal/MASK map'],[225,40;230,40;230,70;245,70])
set_param([sys,'/','PAM Signal/MASK map'],...
		'Mask Display','MASK\nmap',...
		'Mask Type','MASK map')
set_param([sys,'/','PAM Signal/MASK map'],...
		'Mask Dialogue','Map a digital symbol onto an analog\nsignal in the range [-1, 1] using M-ary\namplitude shift keying modulation mapping.|M-ary number (input integers in range [0, M-1]):|Input symbol interval and offset (sec):')
set_param([sys,'/','PAM Signal/MASK map'],...
		'Mask Translate','N=@1;Ac=1;Kc=Ac*2/(N-1);td=@2;')
set_param([sys,'/','PAM Signal/MASK map'],...
		'Mask Help','This block maps an input digit signal onto a signal in the range [-1, 1]. The coded signal can be modulated using analog\nmodulation blocks. Input digits are in the range [0, M-1], where M is the M-ary number. The input symbol interval can be a two-element vector with the second element being offset. The default value offset is 0. Use the MODMAP function to see the constellation.')
set_param([sys,'/','PAM Signal/MASK map'],...
		'Mask Entries','2\/0.0025\/')


%     Finished composite block 'PAM Signal/MASK map'.