tstricia.m

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

set_param([sys,'/',['Rician noise',13,'generator/Gaussian noise',13,'generator']],...
		'Mask Dialogue','Generate Gaussian distributed noise with\ngiven mean and variance value.|Mean value:|Variance (vector or matrix):|Initial seed:')
set_param([sys,'/',['Rician noise',13,'generator/Gaussian noise',13,'generator']],...
		'Mask Translate','m=@1;d=@2;s=@3;[m,d,s]=gwnoise(m,d,s);[x,y]=srcsicon(3);')
set_param([sys,'/',['Rician noise',13,'generator/Gaussian noise',13,'generator']],...
		'Mask Help','When the mean value and seed are vectors, the output of this block is a vector. The length of the vector is the same as the length of the mean value vector or the length of the seed. When the variance entry is a vector, the variance is a diagonal matrix.')
set_param([sys,'/',['Rician noise',13,'generator/Gaussian noise',13,'generator']],...
		'Mask Entries','m\/s\/seed\/')


%     Finished composite block ['Rician noise',13,'generator/Gaussian noise',13,'generator'].

set_param([sys,'/',['Rician noise',13,'generator/Gaussian noise',13,'generator']],...
		'position',[40,33,120,77])


%     Subsystem  ['Rician noise',13,'generator/Gaussian noise',13,'generator1'].

new_system([sys,'/',['Rician noise',13,'generator/Gaussian noise',13,'generator1']])
set_param([sys,'/',['Rician noise',13,'generator/Gaussian noise',13,'generator1']],'Location',[22,119,370,243])

add_block('built-in/Outport',[sys,'/',['Rician noise',13,'generator/Gaussian noise',13,'generator1/out_1']])
set_param([sys,'/',['Rician noise',13,'generator/Gaussian noise',13,'generator1/out_1']],...
		'position',[285,60,305,80])

add_block('built-in/Sum',[sys,'/',['Rician noise',13,'generator/Gaussian noise',13,'generator1/Sum']])
set_param([sys,'/',['Rician noise',13,'generator/Gaussian noise',13,'generator1/Sum']],...
		'position',[220,60,240,80])

add_block('built-in/Constant',[sys,'/',['Rician noise',13,'generator/Gaussian noise',13,'generator1/Constant']])
set_param([sys,'/',['Rician noise',13,'generator/Gaussian noise',13,'generator1/Constant']],...
		'Value','m',...
		'position',[140,20,160,40])

add_block('built-in/White Noise',[sys,'/',['Rician noise',13,'generator/Gaussian noise',13,'generator1/Random',13,'Number']])
set_param([sys,'/',['Rician noise',13,'generator/Gaussian noise',13,'generator1/Random',13,'Number']],...
		'Seed','s',...
		'position',[35,57,80,93])

add_block('built-in/State-Space',[sys,'/',['Rician noise',13,'generator/Gaussian noise',13,'generator1/Matrix',13,'Gain']])
set_param([sys,'/',['Rician noise',13,'generator/Gaussian noise',13,'generator1/Matrix',13,'Gain']],...
		'A','[]',...
		'B','[]',...
		'C','[]',...
		'D','K',...
		'Mask Display','K',...
		'Mask Type','Matrix Gain')
set_param([sys,'/',['Rician noise',13,'generator/Gaussian noise',13,'generator1/Matrix',13,'Gain']],...
		'Mask Dialogue','Matrix Gain.|Gain matrix:',...
		'Mask Translate','K = @1;')
set_param([sys,'/',['Rician noise',13,'generator/Gaussian noise',13,'generator1/Matrix',13,'Gain']],...
		'Mask Help','Multiplies input vector by entered matrix to produce output vector (y=Au).',...
		'Mask Entries','d\/')
set_param([sys,'/',['Rician noise',13,'generator/Gaussian noise',13,'generator1/Matrix',13,'Gain']],...
		'position',[140,61,165,89])
add_line([sys,'/',['Rician noise',13,'generator/Gaussian noise',13,'generator1']],[85,75;135,75])
add_line([sys,'/',['Rician noise',13,'generator/Gaussian noise',13,'generator1']],[170,75;215,75])
add_line([sys,'/',['Rician noise',13,'generator/Gaussian noise',13,'generator1']],[245,70;280,70])
add_line([sys,'/',['Rician noise',13,'generator/Gaussian noise',13,'generator1']],[165,30;190,30;190,65;215,65])
set_param([sys,'/',['Rician noise',13,'generator/Gaussian noise',13,'generator1']],...
		'Mask Display','plot(0,0,100,100,x,y);Gaussian',...
		'Mask Type','Gaussian noise generator')
set_param([sys,'/',['Rician noise',13,'generator/Gaussian noise',13,'generator1']],...
		'Mask Dialogue','Generate Gaussian distributed noise with\ngiven mean and variance value.|Mean value:|Variance (vector or matrix):|Initial seed:')
set_param([sys,'/',['Rician noise',13,'generator/Gaussian noise',13,'generator1']],...
		'Mask Translate','m=@1;d=@2;s=@3;[m,d,s]=gwnoise(m,d,s);[x,y]=srcsicon(3);')
set_param([sys,'/',['Rician noise',13,'generator/Gaussian noise',13,'generator1']],...
		'Mask Help','When the mean value and seed are vectors, the output of this block is a vector. The length of the vector is the same as the length of the mean value vector or the length of the seed. When the variance entry is a vector, the variance is a diagonal matrix.')
set_param([sys,'/',['Rician noise',13,'generator/Gaussian noise',13,'generator1']],...
		'Mask Entries','0\/s\/seed+max(seed)\/')


%     Finished composite block ['Rician noise',13,'generator/Gaussian noise',13,'generator1'].

set_param([sys,'/',['Rician noise',13,'generator/Gaussian noise',13,'generator1']],...
		'position',[40,108,120,152])
add_line([sys,'/',['Rician noise',13,'generator']],[125,55;155,55])
add_line([sys,'/',['Rician noise',13,'generator']],[135,55;135,65;155,65])
add_line([sys,'/',['Rician noise',13,'generator']],[195,60;240,60])
add_line([sys,'/',['Rician noise',13,'generator']],[125,130;135,130;135,125;160,125])
add_line([sys,'/',['Rician noise',13,'generator']],[135,130;135,135;160,135])
add_line([sys,'/',['Rician noise',13,'generator']],[200,130;210,130;210,70;240,70])
add_line([sys,'/',['Rician noise',13,'generator']],[270,65;290,65])
add_line([sys,'/',['Rician noise',13,'generator']],[375,65;405,65])
set_param([sys,'/',['Rician noise',13,'generator']],...
		'Mask Display','plot(0,0,100,100,x,y);Rician (K)',...
		'Mask Type','Rician noise generator')
set_param([sys,'/',['Rician noise',13,'generator']],...
		'Mask Dialogue','Generate Rician distributed noise.\nThe output vector size of this block is the same as the vector size of the seed.|Rician K-factor:|sigma:|Seed:')
set_param([sys,'/',['Rician noise',13,'generator']],...
		'Mask Translate','K=@1;s=@2;seed=@3;m=sqrt(2*K(:)).*s;[x,y]=srcsicon(3);')
set_param([sys,'/',['Rician noise',13,'generator']],...
		'Mask Help','The Rician noise is composed with two independant Gaussian noises. These two Gaussain noises have a same std, sigma. Their mean values are m1 and m2 respectively. The Rician K-factor=(m1^2+m2^2)/2/sigma^2. The output vector size of this block is the same as the vector size of the seed.')
set_param([sys,'/',['Rician noise',13,'generator']],...
		'Mask Entries','2\/1\/12345\/')


%     Finished composite block ['Rician noise',13,'generator'].

set_param([sys,'/',['Rician noise',13,'generator']],...
		'hide name',0,...
		'position',[30,33,110,77])

add_block('built-in/Constant',[sys,'/','Expected value'])
set_param([sys,'/','Expected value'],...
		'Value','[6 20]',...
		'position',[230,134,275,156])


%     Subsystem  ['Rician noise',13,'generator2'].

new_system([sys,'/',['Rician noise',13,'generator2']])
set_param([sys,'/',['Rician noise',13,'generator2']],'Location',[22,451,491,654])

add_block('built-in/Outport',[sys,'/',['Rician noise',13,'generator2/out_1']])
set_param([sys,'/',['Rician noise',13,'generator2/out_1']],...
		'position',[410,55,430,75])

add_block('built-in/S-Function',[sys,'/',['Rician noise',13,'generator2/Array',13,'function']])
set_param([sys,'/',['Rician noise',13,'generator2/Array',13,'function']],...
		'function name','arysin',...
		'parameters','opr',...
		'Mask Display','Array\nFcn',...
		'Mask Type','Array function')
set_param([sys,'/',['Rician noise',13,'generator2/Array',13,'function']],...
		'Mask Dialogue','Multi-input/multi-output functions, such as sin, cos, tan,exp,asin, acos,ln,log10,sgn.|Operation (string):',...
		'Mask Translate','opr=@1;')
set_param([sys,'/',['Rician noise',13,'generator2/Array',13,'function']],...
		'Mask Help','The vector length of the output is the same as the input.',...
		'Mask Entries','''sqrt''\/',...
		'position',[295,43,370,87])

add_block('built-in/Sum',[sys,'/',['Rician noise',13,'generator2/Sum']])
set_param([sys,'/',['Rician noise',13,'generator2/Sum']],...
		'position',[245,55,265,75])

add_block('built-in/Product',[sys,'/',['Rician noise',13,'generator2/Product1']])
set_param([sys,'/',['Rician noise',13,'generator2/Product1']],...
		'position',[165,118,195,142])

add_block('built-in/Product',[sys,'/',['Rician noise',13,'generator2/Product']])
set_param([sys,'/',['Rician noise',13,'generator2/Product']],...
		'position',[160,48,190,72])


%     Subsystem  ['Rician noise',13,'generator2/Gaussian noise',13,'generator'].

new_system([sys,'/',['Rician noise',13,'generator2/Gaussian noise',13,'generator']])
set_param([sys,'/',['Rician noise',13,'generator2/Gaussian noise',13,'generator']],'Location',[22,119,370,243])

add_block('built-in/State-Space',[sys,'/',['Rician noise',13,'generator2/Gaussian noise',13,'generator/Matrix',13,'Gain']])
set_param([sys,'/',['Rician noise',13,'generator2/Gaussian noise',13,'generator/Matrix',13,'Gain']],...
		'A','[]',...
		'B','[]',...
		'C','[]',...
		'D','K',...
		'Mask Display','K',...
		'Mask Type','Matrix Gain')
set_param([sys,'/',['Rician noise',13,'generator2/Gaussian noise',13,'generator/Matrix',13,'Gain']],...
		'Mask Dialogue','Matrix Gain.|Gain matrix:',...
		'Mask Translate','K = @1;')
set_param([sys,'/',['Rician noise',13,'generator2/Gaussian noise',13,'generator/Matrix',13,'Gain']],...
		'Mask Help','Multiplies input vector by entered matrix to produce output vector (y=Au).',...
		'Mask Entries','d\/')
set_param([sys,'/',['Rician noise',13,'generator2/Gaussian noise',13,'generator/Matrix',13,'Gain']],...
		'position',[140,61,165,89])

add_block('built-in/White Noise',[sys,'/',['Rician noise',13,'generator2/Gaussian noise',13,'generator/Random',13,'Number']])
set_param([sys,'/',['Rician noise',13,'generator2/Gaussian noise',13,'generator/Random',13,'Number']],...
		'Seed','s',...
		'position',[35,57,80,93])

add_block('built-in/Constant',[sys,'/',['Rician noise',13,'generator2/Gaussian noise',13,'generator/Constant']])
set_param([sys,'/',['Rician noise',13,'generator2/Gaussian noise',13,'generator/Constant']],...
		'Value','m',...
		'position',[140,20,160,40])

add_block('built-in/Sum',[sys,'/',['Rician noise',13,'generator2/Gaussian noise',13,'generator/Sum']])
set_param([sys,'/',['Rician noise',13,'generator2/Gaussian noise',13,'generator/Sum']],...
		'position',[220,60,240,80])

add_block('built-in/Outport',[sys,'/',['Rician noise',13,'generator2/Gaussian noise',13,'generator/out_1']])
set_param([sys,'/',['Rician noise',13,'generator2/Gaussian noise',13,'generator/out_1']],...
		'position',[285,60,305,80])
add_line([sys,'/',['Rician noise',13,'generator2/Gaussian noise',13,'generator']],[165,30;190,30;190,65;215,65])
add_line([sys,'/',['Rician noise',13,'generator2/Gaussian noise',13,'generator']],[245,70;280,70])
add_line([sys,'/',['Rician noise',13,'generator2/Gaussian noise',13,'generator']],[170,75;215,75])
add_line([sys,'/',['Rician noise',13,'generator2/Gaussian noise',13,'generator']],[85,75;135,75])
set_param([sys,'/',['Rician noise',13,'generator2/Gaussian noise',13,'generator']],...
		'Mask Display','plot(0,0,100,100,x,y);Gaussian',...
		'Mask Type','Gaussian noise generator')
set_param([sys,'/',['Rician noise',13,'generator2/Gaussian noise',13,'generator']],...
		'Mask Dialogue','Generate Gaussian distributed noise with\ngiven mean and variance value.|Mean value:|Variance (vector or matrix):|Initial seed:')
set_param([sys,'/',['Rician noise',13,'generator2/Gaussian noise',13,'generator']],...
		'Mask Translate','m=@1;d=@2;s=@3;[m,d,s]=gwnoise(m,d,s);[x,y]=srcsicon(3);')
set_param([sys,'/',['Rician noise',13,'generator2/Gaussian noise',13,'generator']],...
		'Mask Help','When the mean value and seed are vectors, the output of this block is a vector. The length of the vector is the same as the length of the mean value vector or the length of the seed. When the variance entry is a vector, the variance is a diagonal matrix.')
set_param([sys,'/',['Rician noise',13,'generator2/Gaussian noise',13,'generator']],...
		'Mask Entries','m1\/s\/seed\/')


%     Finished composite block ['Rician noise',13,'generator2/Gaussian noise',13,'generator'].

set_param([sys,'/',['Rician noise',13,'generator2/Gaussian noise',13,'generator']],...
		'position',[40,38,120,82])


%     Subsystem  ['Rician noise',13,'generator2/Gaussian noise',13,'generator1'].

new_system([sys,'/',['Rician noise',13,'generator2/Gaussian noise',13,'generator1']])
set_param([sys,'/',['Rician noise',13,'generator2/Gaussian noise',13,'generator1']],'Location',[22,119,370,243])

add_block('built-in/Outport',[sys,'/',['Rician noise',13,'generator2/Gaussian noise',13,'generator1/out_1']])
set_param([sys,'/',['Rician noise',13,'generator2/Gaussian noise',13,'generator1/out_1']],...
		'position',[285,60,305,80])

add_block('built-in/Sum',[sys,'/',['Rician noise',13,'generator2/Gaussian noise',13,'generator1/Sum']])
set_param([sys,'/',['Rician noise',13,'generator2/Gaussian noise',13,'generator1/Sum']],...
		'position',[220,60,240,80])

add_block('built-in/Constant',[sys,'/',['Rician noise',13,'generator2/Gaussian noise',13,'generator1/Constant']])
set_param([sys,'/',['Rician noise',13,'generator2/Gaussian noise',13,'generator1/Constant']],...
		'Value','m',...
		'position',[140,20,160,40])

add_block('built-in/White Noise',[sys,'/',['Rician noise',13,'generator2/Gaussian noise',13,'generator1/Random',13,'Number']])
set_param([sys,'/',['Rician noise',13,'generator2/Gaussian noise',13,'generator1/Random',13,'Number']],...
		'Seed','s',...
		'position',[35,57,80,93])

add_block('built-in/State-Space',[sys,'/',['Rician noise',13,'generator2/Gaussian noise',13,'generator1/Matrix',13,'Gain']])
set_param([sys,'/',['Rician noise',13,'generator2/Gaussian noise',13,'generator1/Matrix',13,'Gain']],...
		'A','[]',...
		'B','[]',...
		'C','[]',...
		'D','K',...
		'Mask Display','K',...
		'Mask Type','Matrix Gain')
set_param([sys,'/',['Rician noise',13,'generator2/Gaussian noise',13,'generator1/Matrix',13,'Gain']],...
		'Mask Dialogue','Matrix Gain.|Gain matrix:',...
		'Mask Translate','K = @1;')
set_param([sys,'/',['Rician noise',13,'generator2/Gaussian noise',13,'generator1/Matrix',13,'Gain']],...
		'Mask Help','Multiplies input vector by entered matrix to produce output vector (y=Au).',...
		'Mask Entries','d\/')
set_param([sys,'/',['Rician noise',13,'generator2/Gaussian noise',13,'generator1/Matrix',13,'Gain']],...