y=randn(1,2500);
y=y/std(y);
y=y-mean(y);
a=0;
b=sqrt(5);
y=a+b*y;
就得到了 N ( 0,5 ) 的高斯分布序列.
MATLAB中產(chǎn)生高斯白噪聲的兩個函數(shù)
MATLAB中產(chǎn)生高斯白噪聲非常方便,可以直接應用兩個函數(shù),一個是WGN,另一個是AWGN.WGN用于產(chǎn)生高斯白噪聲,AWGN則用于在某一信號中加入高斯白噪聲.
1.WGN：產(chǎn)生高斯白噪聲
y = wgn(m,n,p) 產(chǎn)生一個m行n列的高斯白噪聲的矩陣,p以dBW為單位指定輸出噪聲的強度.
y = wgn(m,n,p,imp) 以歐姆(Ohm)為單位指定負載阻抗.
y = wgn(m,n,p,imp,state) 重置RANDN的狀態(tài).
在數(shù)值變量后還可附加一些標志性參數(shù)：
y = wgn(…,POWERTYPE) 指定p的單位.POWERTYPE可以是'dBW','dBm'或'linear'.線性強度(linear power)以瓦特(Watt)為單位.
y = wgn(…,OUTPUTTYPE) 指定輸出類型.OUTPUTTYPE可以是'real'或'complex'.
2.AWGN：在某一信號中加入高斯白噪聲
y = awgn(x,SNR) 在信號x中加入高斯白噪聲.信噪比SNR以dB為單位.x的強度假定為0dBW.如果x是復數(shù),就加入復噪聲.
y = awgn(x,SNR,SIGPOWER) 如果SIGPOWER是數(shù)值,則其代表以dBW為單位的信號強度；如果SIGPOWER為'measured',則函數(shù)將在加入噪聲之前測定信號強度.
y = awgn(x,SNR,SIGPOWER,STATE) 重置RANDN的狀態(tài).
y = awgn(…,POWERTYPE) 指定SNR和SIGPOWER的單位.POWERTYPE可以是'dB'或'linear'.如果POWERTYPE是'dB',那么SNR以dB為單位,而SIGPOWER以dBW為單位.如果POWERTYPE是'linear',那么SNR作為比值來度量,而SIGPOWER以瓦特為單位.