com_sour.m

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

set_param([sys,'/',['Rayleigh noise',13,'generator']],...
		'Mask Help','This block generates Rayleigh noise with the given sigma value. The sigma value is the standard deviations of the zero mean Gaussian noise components in generating the Rayleigh noise. The output vector size of this block is the same as the vector size of the seed.')
set_param([sys,'/',['Rayleigh noise',13,'generator']],...
		'Mask Entries','1\/17658\/')


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

set_param([sys,'/',['Rayleigh noise',13,'generator']],...
		'position',[455,413,535,457])


%     Subsystem  ['Gaussian noise',13,'generator'].

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

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

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

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

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

add_block('built-in/State-Space',[sys,'/',['Gaussian noise',13,'generator/Matrix',13,'Gain']])
set_param([sys,'/',['Gaussian noise',13,'generator/Matrix',13,'Gain']],...
		'A','[]',...
		'B','[]',...
		'C','[]',...
		'D','K',...
		'Mask Display','K',...
		'Mask Type','Matrix Gain',...
		'Mask Dialogue','Matrix Gain.|Gain matrix:')
set_param([sys,'/',['Gaussian noise',13,'generator/Matrix',13,'Gain']],...
		'Mask Translate','K = @1;',...
		'Mask Help','Multiplies input vector by entered matrix to produce output vector (y=Au).',...
		'Mask Entries','d\/')
set_param([sys,'/',['Gaussian noise',13,'generator/Matrix',13,'Gain']],...
		'position',[140,61,165,89])
add_line([sys,'/',['Gaussian noise',13,'generator']],[85,75;135,75])
add_line([sys,'/',['Gaussian noise',13,'generator']],[170,75;215,75])
add_line([sys,'/',['Gaussian noise',13,'generator']],[245,70;280,70])
add_line([sys,'/',['Gaussian noise',13,'generator']],[165,30;190,30;190,65;215,65])
set_param([sys,'/',['Gaussian noise',13,'generator']],...
		'Mask Display','plot(0,0,100,100,x,y);Gaussian',...
		'Mask Type','Gaussian noise generator')
set_param([sys,'/',['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,'/',['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,'/',['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,'/',['Gaussian noise',13,'generator']],...
		'Mask Entries','[0 1]\/[1 2]\/[123454 4321]\/')


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

set_param([sys,'/',['Gaussian noise',13,'generator']],...
		'position',[455,113,535,157])


%     Subsystem  ['Rician noise',13,'generator K'].

new_system([sys,'/',['Rician noise',13,'generator K']])
set_param([sys,'/',['Rician noise',13,'generator K']],'Location',[52,446,521,649])


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

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

add_block('built-in/State-Space',[sys,'/',['Rician noise',13,'generator K/Gaussian noise',13,'generator1/Matrix',13,'Gain']])
set_param([sys,'/',['Rician noise',13,'generator K/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 K/Gaussian noise',13,'generator1/Matrix',13,'Gain']],...
		'Mask Dialogue','Matrix Gain.|Gain matrix:',...
		'Mask Translate','K = @1;')
set_param([sys,'/',['Rician noise',13,'generator K/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 K/Gaussian noise',13,'generator1/Matrix',13,'Gain']],...
		'position',[140,61,165,89])

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

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

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

add_block('built-in/Outport',[sys,'/',['Rician noise',13,'generator K/Gaussian noise',13,'generator1/out_1']])
set_param([sys,'/',['Rician noise',13,'generator K/Gaussian noise',13,'generator1/out_1']],...
		'position',[285,60,305,80])
add_line([sys,'/',['Rician noise',13,'generator K/Gaussian noise',13,'generator1']],[165,30;190,30;190,65;215,65])
add_line([sys,'/',['Rician noise',13,'generator K/Gaussian noise',13,'generator1']],[245,70;280,70])
add_line([sys,'/',['Rician noise',13,'generator K/Gaussian noise',13,'generator1']],[170,75;215,75])
add_line([sys,'/',['Rician noise',13,'generator K/Gaussian noise',13,'generator1']],[85,75;135,75])
set_param([sys,'/',['Rician noise',13,'generator K/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 K/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 K/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 K/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 K/Gaussian noise',13,'generator1']],...
		'Mask Entries','0\/s\/seed+max(seed)\/')


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

set_param([sys,'/',['Rician noise',13,'generator K/Gaussian noise',13,'generator1']],...
		'position',[40,108,120,152])


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

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

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

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

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

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

add_block('built-in/State-Space',[sys,'/',['Rician noise',13,'generator K/Gaussian noise',13,'generator/Matrix',13,'Gain']])
set_param([sys,'/',['Rician noise',13,'generator K/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,'generator K/Gaussian noise',13,'generator/Matrix',13,'Gain']],...
		'Mask Dialogue','Matrix Gain.|Gain matrix:',...
		'Mask Translate','K = @1;')
set_param([sys,'/',['Rician noise',13,'generator K/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,'generator K/Gaussian noise',13,'generator/Matrix',13,'Gain']],...
		'position',[140,61,165,89])
add_line([sys,'/',['Rician noise',13,'generator K/Gaussian noise',13,'generator']],[85,75;135,75])
add_line([sys,'/',['Rician noise',13,'generator K/Gaussian noise',13,'generator']],[170,75;215,75])
add_line([sys,'/',['Rician noise',13,'generator K/Gaussian noise',13,'generator']],[245,70;280,70])
add_line([sys,'/',['Rician noise',13,'generator K/Gaussian noise',13,'generator']],[165,30;190,30;190,65;215,65])
set_param([sys,'/',['Rician noise',13,'generator K/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,'generator K/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 K/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 K/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 K/Gaussian noise',13,'generator']],...
		'Mask Entries','m\/s\/seed\/')


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

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

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

add_block('built-in/S-Function',[sys,'/',['Rician noise',13,'generator K/Array',13,'function']])
set_param([sys,'/',['Rician noise',13,'generator K/Array',13,'function']],...
		'function name','arysin',...
		'parameters','opr',...
		'Mask Display','Array\nFcn',...
		'Mask Type','Array function')
set_param([sys,'/',['Rician noise',13,'generator K/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,'generator K/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,'generator K/Sum']])
set_param([sys,'/',['Rician noise',13,'generator K/Sum']],...
		'position',[245,55,265,75])

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

add_block('built-in/Product',[sys,'/',['Rician noise',13,'generator K/Product']])
set_param([sys,'/',['Rician noise',13,'generator K/Product']],...
		'position',[160,48,190,72])
add_line([sys,'/',['Rician noise',13,'generator K']],[375,65;405,65])
add_line([sys,'/',['Rician noise',13,'generator K']],[270,65;290,65])
add_line([sys,'/',['Rician noise',13,'generator K']],[200,130;210,130;210,70;240,70])
add_line([sys,'/',['Rician noise',13,'generator K']],[125,130;135,130;135,125;160,125])
add_line([sys,'/',['Rician noise',13,'generator K']],[135,130;135,135;160,135])
add_line([sys,'/',['Rician noise',13,'generator K']],[195,60;240,60])
add_line([sys,'/',['Rician noise',13,'generator K']],[125,55;155,55])
add_line([sys,'/',['Rician noise',13,'generator K']],[135,55;135,65;155,65])
set_param([sys,'/',['Rician noise',13,'generator K']],...
		'Mask Display','plot(0,0,100,100,x,y);Rician (K)',...
		'Mask Type','Rician noise generator')
set_param([sys,'/',['Rician noise',13,'generator K']],...
		'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 K']],...
		'Mask Translate','K=@1;s=@2;seed=@3;m=sqrt(2*K(:)).*s;[x,y]=srcsicon(3);')
set_param([sys,'/',['Rician noise',13,'generator K']],...
		'Mask Help','The Rician noise is composed with two independant Gaussian noises. These two Gau