tstrand.m

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

set_param([sys,'/','Error-Meter'],...
		'position',[345,83,420,127])


%     Subsystem  'Mean&std'.

new_system([sys,'/','Mean&std'])
set_param([sys,'/','Mean&std'],'Location',[16,148,871,412])

add_block('built-in/Inport',[sys,'/','Mean&std/Actual'])
set_param([sys,'/','Mean&std/Actual'],...
		'position',[25,35,45,55])

add_block('built-in/Digital Clock',[sys,'/','Mean&std/Digital Clock'])
set_param([sys,'/','Mean&std/Digital Clock'],...
		'Sample time','Ts',...
		'position',[45,80,80,110])

add_block('built-in/Fcn',[sys,'/','Mean&std/Fcn2'])
set_param([sys,'/','Mean&std/Fcn2'],...
		'Expr','u[1]/Ts +1',...
		'position',[105,83,165,107])

add_block('built-in/Unit Delay',[sys,'/','Mean&std/Unit Delay1'])
set_param([sys,'/','Mean&std/Unit Delay1'],...
		'orientation',2,...
		'Sample time','Ts',...
		'position',[350,87,400,103])

add_block('built-in/Outport',[sys,'/','Mean&std/mean value'])
set_param([sys,'/','Mean&std/mean value'],...
		'position',[545,45,565,65])

add_block('built-in/Mux',[sys,'/','Mean&std/Mux1'])
set_param([sys,'/','Mean&std/Mux1'],...
		'inputs','3',...
		'position',[235,39,265,71])

add_block('built-in/Fcn',[sys,'/','Mean&std/Fcn3'])
set_param([sys,'/','Mean&std/Fcn3'],...
		'Expr','(u[3]*(u[2]-1) + u[1])/u[2]',...
		'position',[290,41,465,69])

add_block('built-in/Mux',[sys,'/','Mean&std/Mux'])
set_param([sys,'/','Mean&std/Mux'],...
		'inputs','3',...
		'position',[235,134,265,166])

add_block('built-in/Unit Delay',[sys,'/','Mean&std/Unit Delay'])
set_param([sys,'/','Mean&std/Unit Delay'],...
		'orientation',2,...
		'Sample time','Ts',...
		'position',[345,192,395,208])

add_block('built-in/Fcn',[sys,'/','Mean&std/Fcn'])
set_param([sys,'/','Mean&std/Fcn'],...
		'Expr','(u[3]*(u[2]-1) + u[1]*u[1])/u[2]',...
		'position',[290,136,465,164])

add_block('built-in/Sum',[sys,'/','Mean&std/Sum'])
set_param([sys,'/','Mean&std/Sum'],...
		'inputs','-+',...
		'position',[645,135,665,155])

add_block('built-in/Outport',[sys,'/',['Mean&std/Mean',13,'Squared',13,'Error1']])
set_param([sys,'/',['Mean&std/Mean',13,'Squared',13,'Error1']],...
		'Port','2',...
		'position',[790,135,810,155])

add_block('built-in/Fcn',[sys,'/','Mean&std/Fcn1'])
set_param([sys,'/','Mean&std/Fcn1'],...
		'Expr','sqrt(u[1])',...
		'position',[690,132,755,158])

add_block('built-in/Mux',[sys,'/','Mean&std/Mux2'])
set_param([sys,'/','Mean&std/Mux2'],...
		'inputs','2',...
		'position',[515,101,545,134])

add_block('built-in/Fcn',[sys,'/','Mean&std/Fcn4'])
set_param([sys,'/','Mean&std/Fcn4'],...
		'Expr','u[1]*u[1]',...
		'position',[555,107,620,133])
add_line([sys,'/','Mean&std'],[85,95;100,95])
add_line([sys,'/','Mean&std'],[340,200;220,200;230,160])
add_line([sys,'/','Mean&std'],[170,95;190,95;190,150;230,150])
add_line([sys,'/','Mean&std'],[270,150;285,150])
add_line([sys,'/','Mean&std'],[270,55;285,55])
add_line([sys,'/','Mean&std'],[470,55;540,55])
add_line([sys,'/','Mean&std'],[760,145;785,145])
add_line([sys,'/','Mean&std'],[50,45;210,45;230,45])
add_line([sys,'/','Mean&std'],[190,95;190,55;230,55])
add_line([sys,'/','Mean&std'],[495,55;495,95;405,95])
add_line([sys,'/','Mean&std'],[345,95;210,95;210,65;230,65])
add_line([sys,'/','Mean&std'],[200,45;200,140;230,140])
add_line([sys,'/','Mean&std'],[470,150;480,150;480,200;400,200])
add_line([sys,'/','Mean&std'],[670,145;685,145])
add_line([sys,'/','Mean&std'],[480,150;640,150])
add_line([sys,'/','Mean&std'],[625,120;625,140;640,140])
add_line([sys,'/','Mean&std'],[495,95;495,110;510,110])
add_line([sys,'/','Mean&std'],[190,125;500,125;510,125])
add_line([sys,'/','Mean&std'],[550,120;550,120])
set_param([sys,'/','Mean&std'],...
		'Mask Display','Mean\n& std',...
		'Mask Type','Mean Squared Error',...
		'Mask Dialogue','Calculates the mean squared error of the two inputs sampling at the delay time.|Delay time:')
set_param([sys,'/','Mean&std'],...
		'Mask Translate','Ts=@1;',...
		'Mask Help','For the inputs x and xHat and delay time T, the output, e, at the n''th time step will be\ne[nT] = sqrt( sum( (x-xHat)^2 )/n ).',...
		'Mask Entries','td\/')


%     Finished composite block 'Mean&std'.

set_param([sys,'/','Mean&std'],...
		'position',[245,83,300,127])


%     Subsystem  'Mean&var'.

new_system([sys,'/','Mean&var'])
set_param([sys,'/','Mean&var'],'Location',[16,148,774,399])

add_block('built-in/Outport',[sys,'/',['Mean&var/Mean',13,'Squared',13,'Error1']])
set_param([sys,'/',['Mean&var/Mean',13,'Squared',13,'Error1']],...
		'Port','2',...
		'position',[700,135,720,155])

add_block('built-in/Fcn',[sys,'/','Mean&var/Fcn4'])
set_param([sys,'/','Mean&var/Fcn4'],...
		'Expr','u[1]*u[1]',...
		'position',[555,107,620,133])

add_block('built-in/Mux',[sys,'/','Mean&var/Mux2'])
set_param([sys,'/','Mean&var/Mux2'],...
		'inputs','2',...
		'position',[515,101,545,134])

add_block('built-in/Sum',[sys,'/','Mean&var/Sum'])
set_param([sys,'/','Mean&var/Sum'],...
		'inputs','-+',...
		'position',[645,135,665,155])

add_block('built-in/Fcn',[sys,'/','Mean&var/Fcn'])
set_param([sys,'/','Mean&var/Fcn'],...
		'Expr','(u[3]*(u[2]-1) + u[1]*u[1])/u[2]',...
		'position',[290,136,465,164])

add_block('built-in/Unit Delay',[sys,'/','Mean&var/Unit Delay'])
set_param([sys,'/','Mean&var/Unit Delay'],...
		'orientation',2,...
		'Sample time','Ts',...
		'position',[345,192,395,208])

add_block('built-in/Mux',[sys,'/','Mean&var/Mux'])
set_param([sys,'/','Mean&var/Mux'],...
		'inputs','3',...
		'position',[235,134,265,166])

add_block('built-in/Fcn',[sys,'/','Mean&var/Fcn3'])
set_param([sys,'/','Mean&var/Fcn3'],...
		'Expr','(u[3]*(u[2]-1) + u[1])/u[2]',...
		'position',[290,41,465,69])

add_block('built-in/Mux',[sys,'/','Mean&var/Mux1'])
set_param([sys,'/','Mean&var/Mux1'],...
		'inputs','3',...
		'position',[235,39,265,71])

add_block('built-in/Outport',[sys,'/','Mean&var/mean value'])
set_param([sys,'/','Mean&var/mean value'],...
		'position',[545,45,565,65])

add_block('built-in/Unit Delay',[sys,'/','Mean&var/Unit Delay1'])
set_param([sys,'/','Mean&var/Unit Delay1'],...
		'orientation',2,...
		'Sample time','Ts',...
		'position',[350,87,400,103])

add_block('built-in/Fcn',[sys,'/','Mean&var/Fcn2'])
set_param([sys,'/','Mean&var/Fcn2'],...
		'Expr','u[1]/Ts +1',...
		'position',[105,83,165,107])

add_block('built-in/Digital Clock',[sys,'/','Mean&var/Digital Clock'])
set_param([sys,'/','Mean&var/Digital Clock'],...
		'Sample time','Ts',...
		'position',[45,80,80,110])

add_block('built-in/Inport',[sys,'/','Mean&var/Actual'])
set_param([sys,'/','Mean&var/Actual'],...
		'position',[25,35,45,55])
add_line([sys,'/','Mean&var'],[670,145;695,145])
add_line([sys,'/','Mean&var'],[550,120;550,120])
add_line([sys,'/','Mean&var'],[170,95;190,95;190,150;230,150])
add_line([sys,'/','Mean&var'],[190,125;510,125])
add_line([sys,'/','Mean&var'],[470,55;540,55])
add_line([sys,'/','Mean&var'],[495,55;495,95;405,95])
add_line([sys,'/','Mean&var'],[495,95;495,110;510,110])
add_line([sys,'/','Mean&var'],[625,120;625,140;640,140])
add_line([sys,'/','Mean&var'],[470,150;480,150;480,200;400,200])
add_line([sys,'/','Mean&var'],[480,150;640,150])
add_line([sys,'/','Mean&var'],[50,45;230,45])
add_line([sys,'/','Mean&var'],[200,45;200,140;230,140])
add_line([sys,'/','Mean&var'],[345,95;210,95;210,65;230,65])
add_line([sys,'/','Mean&var'],[190,95;190,55;230,55])
add_line([sys,'/','Mean&var'],[270,55;285,55])
add_line([sys,'/','Mean&var'],[270,150;285,150])
add_line([sys,'/','Mean&var'],[340,200;220,200;230,160])
add_line([sys,'/','Mean&var'],[85,95;100,95])
set_param([sys,'/','Mean&var'],...
		'Mask Display','Mean\n& var',...
		'Mask Type','Mean & Variant',...
		'Mask Dialogue','Calculates the mean squared error of the two inputs sampling at the delay time.|Delay time:')
set_param([sys,'/','Mean&var'],...
		'Mask Translate','Ts=@1;',...
		'Mask Help','For the inputs x and xHat and delay time T, the output, e, at the n''th time step will be\ne[nT] = sqrt( sum( (x-xHat)^2 )/n ).',...
		'Mask Entries','td\/')


%     Finished composite block 'Mean&var'.

set_param([sys,'/','Mean&var'],...
		'position',[175,160,245,200])
add_line(sys,[305,95;340,95])
add_line(sys,[305,115;340,115])
add_line(sys,[135,240;185,240;185,105;240,105])

drawnow

% Return any arguments.
if (nargin | nargout)
	% Must use feval here to access system in memory
	if (nargin > 3)
		if (flag == 0)
			eval(['[ret,x0,str,ts,xts]=',sys,'(t,x,u,flag);'])
		else
			eval(['ret =', sys,'(t,x,u,flag);'])
		end
	else
		[ret,x0,str,ts,xts] = feval(sys);
	end
else
	drawnow % Flash up the model and execute load callback
end