tceps2c2.m

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

set_param([sys,'/',['C-map QASK',13,'mod baseband/QASK',13,'circle-map/Imag']],...
		'Port','2',...
		'position',[275,100,295,120])
add_line([sys,'/',['C-map QASK',13,'mod baseband/QASK',13,'circle-map']],[50,30;75,30])
add_line([sys,'/',['C-map QASK',13,'mod baseband/QASK',13,'circle-map']],[120,30;175,30])
add_line([sys,'/',['C-map QASK',13,'mod baseband/QASK',13,'circle-map']],[145,30;145,110;170,110])
add_line([sys,'/',['C-map QASK',13,'mod baseband/QASK',13,'circle-map']],[220,30;270,30])
add_line([sys,'/',['C-map QASK',13,'mod baseband/QASK',13,'circle-map']],[215,110;270,110])
set_param([sys,'/',['C-map QASK',13,'mod baseband/QASK',13,'circle-map']],...
		'Mask Display','QASK\nC-map',...
		'Mask Type','QASK C-map')
set_param([sys,'/',['C-map QASK',13,'mod baseband/QASK',13,'circle-map']],...
		'Mask Dialogue','QASK mapping with circle constellation. The constellation is defined in NIC, RIC and PIC.|Number of symbols in each circle NIC:|Radii in each circle RIC:|Phase shift in each circle PIC (rad):|Symbol interval and offset (sec):')
set_param([sys,'/',['C-map QASK',13,'mod baseband/QASK',13,'circle-map']],...
		'Mask Translate','N=@1;A=@2;P=@3;td=@4;y1=apkconst(N,A,P);y1=[y1;y1];y1=y1(:);y2=imag(y1);y1=real(y1);x=[-sum(N) 1:sum(N)-1;1:sum(N)-1 10*sum(N)]-.5;x=x(:);')
set_param([sys,'/',['C-map QASK',13,'mod baseband/QASK',13,'circle-map']],...
		'Mask Help','This block maps a digital signal onto in-phase and quadratic components, which are outputs in the first and second output port respectively. The circle constellation is defined in NIC, RIC, and PIC. 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,'/',['C-map QASK',13,'mod baseband/QASK',13,'circle-map']],...
		'Mask Entries','NIR\/AIR\/PIR\/td\/')


%     Finished composite block ['C-map QASK',13,'mod baseband/QASK',13,'circle-map'].

set_param([sys,'/',['C-map QASK',13,'mod baseband/QASK',13,'circle-map']],...
		'hide name',0,...
		'position',[65,23,145,67])

add_block('built-in/Inport',[sys,'/',['C-map QASK',13,'mod baseband/in_1']])
set_param([sys,'/',['C-map QASK',13,'mod baseband/in_1']],...
		'hide name',0,...
		'position',[25,35,40,55])

add_block('built-in/Outport',[sys,'/',['C-map QASK',13,'mod baseband/out_1']])
set_param([sys,'/',['C-map QASK',13,'mod baseband/out_1']],...
		'hide name',0,...
		'position',[290,35,310,55])
add_line([sys,'/',['C-map QASK',13,'mod baseband']],[150,35;175,35])
add_line([sys,'/',['C-map QASK',13,'mod baseband']],[150,55;175,55])
add_line([sys,'/',['C-map QASK',13,'mod baseband']],[265,45;285,45])
add_line([sys,'/',['C-map QASK',13,'mod baseband']],[45,45;60,45])
set_param([sys,'/',['C-map QASK',13,'mod baseband']],...
		'Mask Display','plot(0,0,100,100,x,y,u,v);C-QASK',...
		'Mask Type','Baseband C-QASK modulation')
set_param([sys,'/',['C-map QASK',13,'mod baseband']],...
		'Mask Dialogue','Output the complex envelope of a QASK\nmodulated signal. The QASK circle conste-\nllation is defined by the NIC, RIC, and PIC.|Number of symbols in each circle (NIC):|Radii in each circle (RIC):|Phase shift in each circle (PIC, rad):|Initial phase (rad):|Symbol interval (sec):|Sampling time (sec):')
set_param([sys,'/',['C-map QASK',13,'mod baseband']],...
		'Mask Translate','NIR=@1;AIR=@2;PIR=@3;Ph=@4;td=@5;ts=@6;[x,y]=moduicon(5);[u,v]=cmplxicn(1,50,get_param(gcb,''orientation''));')
set_param([sys,'/',['C-map QASK',13,'mod baseband']],...
		'Mask Help','This block outputs the complex envelope of a QASK modulated signal. This block uses circle constellation. The vector lengths of numbers on each circle, radii on each circle, and phases on each circle should be the same. Use function DMODCE to view the constellation of the signal distribution.')
set_param([sys,'/',['C-map QASK',13,'mod baseband']],...
		'Mask Entries','NIR\/AIR\/PIR\/0\/td\/ts\/')


%     Finished composite block ['C-map QASK',13,'mod baseband'].

set_param([sys,'/',['C-map QASK',13,'mod baseband']],...
		'hide name',0,...
		'position',[140,38,220,82])


%     Subsystem  ['Sampled read',13,'from wksp'].

new_system([sys,'/',['Sampled read',13,'from wksp']])
set_param([sys,'/',['Sampled read',13,'from wksp']],'Location',[355,528,737,675])

add_block('built-in/S-Function',[sys,'/',['Sampled read',13,'from wksp/Message source']])
set_param([sys,'/',['Sampled read',13,'from wksp/Message source']],...
		'BackGround',5,...
		'function name','trigwksp',...
		'parameters','var,thd,cyclFlag,ini',...
		'Mask Display','Triggered\nread',...
		'Mask Type','Read from workspace')
set_param([sys,'/',['Sampled read',13,'from wksp/Message source']],...
		'Mask Dialogue','The raising edge of input pulse trigger the block output the next row from the workspace variable.|Read the variable:|Threshold for the input pulse:|Cyclic read variable (1) or amend zeros (0).|Output before first trigger:')
set_param([sys,'/',['Sampled read',13,'from wksp/Message source']],...
		'Mask Translate','var=@1;thd=@2;cyclFlag=@3;ini=@4;')
set_param([sys,'/',['Sampled read',13,'from wksp/Message source']],...
		'Mask Help','This block reads the data from a workspace variable. The output is refreshed at the raising edge of the input pulse.',...
		'Mask Entries','var\/.1\/cyc\/ini\/')
set_param([sys,'/',['Sampled read',13,'from wksp/Message source']],...
		'position',[180,53,255,97])

add_block('built-in/Outport',[sys,'/',['Sampled read',13,'from wksp/out_1']])
set_param([sys,'/',['Sampled read',13,'from wksp/out_1']],...
		'position',[295,65,315,85])


%     Subsystem  ['Sampled read',13,'from wksp/Pulses deivide',13,'same sample time1'].

new_system([sys,'/',['Sampled read',13,'from wksp/Pulses deivide',13,'same sample time1']])
set_param([sys,'/',['Sampled read',13,'from wksp/Pulses deivide',13,'same sample time1']],'Location',[55,192,283,286])

add_block('built-in/Outport',[sys,'/',['Sampled read',13,'from wksp/Pulses deivide',13,'same sample time1/Outport']])
set_param([sys,'/',['Sampled read',13,'from wksp/Pulses deivide',13,'same sample time1/Outport']],...
		'position',[160,35,180,55])

add_block('built-in/S-Function',[sys,'/',['Sampled read',13,'from wksp/Pulses deivide',13,'same sample time1/S-function1']])
set_param([sys,'/',['Sampled read',13,'from wksp/Pulses deivide',13,'same sample time1/S-function1']],...
		'function name','homopuls',...
		'parameters','samp_time, out_divid, offset',...
		'position',[40,31,105,59])
add_line([sys,'/',['Sampled read',13,'from wksp/Pulses deivide',13,'same sample time1']],[110,45;155,45])
set_param([sys,'/',['Sampled read',13,'from wksp/Pulses deivide',13,'same sample time1']],...
		'Mask Display','Vector\nPulse',...
		'Mask Type','Pulse generator')
set_param([sys,'/',['Sampled read',13,'from wksp/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,'/',['Sampled read',13,'from wksp/Pulses deivide',13,'same sample time1']],...
		'Mask Translate','samp_time=@1;out_divid=@2;offset=@3;')
set_param([sys,'/',['Sampled read',13,'from wksp/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,'/',['Sampled read',13,'from wksp/Pulses deivide',13,'same sample time1']],...
		'Mask Entries','ts\/1\/off\/')


%     Finished composite block ['Sampled read',13,'from wksp/Pulses deivide',13,'same sample time1'].

set_param([sys,'/',['Sampled read',13,'from wksp/Pulses deivide',13,'same sample time1']],...
		'position',[50,56,110,94])
add_line([sys,'/',['Sampled read',13,'from wksp']],[115,75;175,75])
add_line([sys,'/',['Sampled read',13,'from wksp']],[260,75;290,75])
set_param([sys,'/',['Sampled read',13,'from wksp']],...
		'Mask Display','plot(0,0,100,100,x,y);Rd wksp',...
		'Mask Type','Read from workspace')
set_param([sys,'/',['Sampled read',13,'from wksp']],...
		'Mask Dialogue','Read from a workspace variable at\nsampling time point.|Variable (column_number=block_output_length):|Data output sample time (sec):|Cyclic control (1: cyclic read, 0: output zeros after eof).|Initial output:')
set_param([sys,'/',['Sampled read',13,'from wksp']],...
		'Mask Translate','var=@1;td=@2;ts=td(1);if (length(td)>=2),off=td(2);else off=0;end;cyc=@3;ini=@4;[x,y]=srcsicon(7);')
set_param([sys,'/',['Sampled read',13,'from wksp']],...
		'Mask Help','This block reads from a workspace variable at a given sample time. The column number is the block output vector length. When the simulation reached the end of the variable (eof) row, the block will cyclically output from the first row of the variable or output zeros based on the cyclic control entree.')
set_param([sys,'/',['Sampled read',13,'from wksp']],...
		'Mask Entries','randint(1000,1,M)\/td\/1\/0\/')


%     Finished composite block ['Sampled read',13,'from wksp'].

set_param([sys,'/',['Sampled read',13,'from wksp']],...
		'hide name',0,...
		'position',[10,38,85,82])


%     Subsystem  ['Rician noise channel',13,'w vary para CE'].

new_system([sys,'/',['Rician noise channel',13,'w vary para CE']])
set_param([sys,'/',['Rician noise channel',13,'w vary para CE']],'Location',[107,289,478,470])

add_block('built-in/Outport',[sys,'/',['Rician noise channel',13,'w vary para CE/out_1']])
set_param([sys,'/',['Rician noise channel',13,'w vary para CE/out_1']],...
		'position',[325,60,345,80])

add_block('built-in/Sum',[sys,'/',['Rician noise channel',13,'w vary para CE/Sum']])
set_param([sys,'/',['Rician noise channel',13,'w vary para CE/Sum']],...
		'position',[275,60,295,80])


%     Subsystem  ['Rician noise channel',13,'w vary para CE/Gaussian',13,'noise'].

new_system([sys,'/',['Rician noise channel',13,'w vary para CE/Gaussian',13,'noise']])
set_param([sys,'/',['Rician noise channel',13,'w vary para CE/Gaussian',13,'noise']],'Location',[482,386,971,557])

add_block('built-in/White Noise',[sys,'/',['Rician noise channel',13,'w vary para CE/Gaussian',13,'noise/Random',13,'Number']])
set_param([sys,'/',['Rician noise channel',13,'w vary para CE/Gaussian',13,'noise/Random',13,'Number']],...
		'Seed','s',...
		'position',[40,82,85,118])

add_block('built-in/Inport',[sys,'/',['Rician noise channel',13,'w vary para CE/Gaussian',13,'noise/in_1']])
set_param([sys,'/',['Rician noise channel',13,'w vary para CE/Gaussian',13,'noise/in_1']],...
		'position',[185,40,205,60])

add_block('built-in/Inport',[sys,'/',['Rician noise channel',13,'w vary para CE/Gaussian',13,'noise/in_2']])
set_param([sys,'/',['Rician noise channel',13,'w vary para CE/Gaussian',13,'noise/in_2']],...
		'Port','2',...
		'position',[105,45,125,65])

add_block('built-in/Outport',[sys,'/',['Rician noise channel',13,'w vary para CE/Gaussian',13,'noise/out_1']])
set_param([sys,'/',['Rician noise channel',13,'w vary para CE/Gaussian',13,'noise/out_1']],...
		'position',[350,80,370,100])

add_block('built-in/Sum',[sys,'/',['Rician noise channel',13,'w vary para CE/Gaussian',13,'noise/Sum']])
set_param([sys,'/',['Rician noise channel',13,'w vary para CE/Gaussian',13,'noise/Sum']],...
		'position',[265,80,285,100])

add_block('built-in/Product',[sys,'/',['Rician noise channel',13,'w vary para CE/Gaussian',13,'noise/Product']])
set_param([sys,'/',['Rician noise channel',13,'w vary para CE/Gaussian',13,'noise/Product']],...
		'position',[175,85,200,105])
add_line([sys,'/',['Rician noise channel',13,'w vary para CE/Gaussian',13,'noise']],[210,50;230,50;230,85;260,85])
add_line([sys,'/',['Rician noise channel',13,'w vary para CE/Gaussian',13,'noise']],[130,55;145,55;145,90;170,90])
add_line([sys,'/',['Rician noise channel',13,'w vary para CE/Gaussian',13,'noise']],[90,100;170,100])
add_line([sys,'/',['Rician noise channel',13,'w vary para CE/Gaussian',13,'noise']],[205,95;260,95])
add_line([sys,'/',['Rician noise channel',13,'w vary para CE/Gaussian',13,'noise']],[290,90;345,90])
set_param([sys,'/',['Rician noise channel',13,'w vary para CE/Gaussian',13,'noise']],...
		'Mask Display','Gaussian\nnoise',...
		'Mask Type','Gaussian noise.',...
		'Mask Dialogue','Gaussian distribution noise.|Initial Seed(s):')
set_param([sys,'/',['Rician noise channel',13,'w vary para CE/Gaussian',13,'noise']],...
		'Mask Translate','s=@1;')
set_param([sys,'/',['Rician noise channel',13,'w vary para CE/Gaussian',13,'noise']],...
		'Mask Help','The meanvalue is determined by value from the frist input. The standard deviation is determined by the second input. ')
set_param([sys,'/',['Rician noise channel',13,'w vary para CE/Gaussian',13,'noise']],...
		'Mask Entries','s\/')


%     Finished composite block ['Rician noise channel',13,'w vary para CE/Gaussian',13,'noise'].

set_param([sys,'/',['Rician noise channel',13,'w vary para CE/Gaussian',13,'noise']],...
		'position',[155,55,230,95])

add_block('built-in/Inport',[sys,'/',['Rician noise channel',13,'w vary para CE/in_2']])
set_param([sys,'/',['Rician noise channel',13,'w vary para CE/in_2']],...
		'Port','2',...
		'position',[15,45,35,65])

add_block('built-in/Mux',[sys,'/',['Rician noise channel',13,'w vary para CE/Mux']])
set_param([sys,'/',['Rician noise channel',13,'w vary para CE/Mux']],...