楂樻柉鐧藉櫔澹扮殑matlab瀹炵幇 - 鐧惧害鏂囧簱

通信系统建模与仿真

实验一、高斯白噪声的matlab实现

要求：

样本点：100 1000 标准差：0.2 2 10 均值： 0 0.2

白噪声

如果噪声的功率谱密度在所有的频率上均为一常数，即

Pn(f)?n0,(???f???),(W/Hz)式中：n0为常数，责成该噪声为白噪声，用n(t)表示。

高斯白噪声的matlab实现

1．样本点为1000、均值为0、标准差为0.2时，高斯白噪声分布为下图所示：

程序如下所示：

% White background nois clear all

f = 1:1:1000;

for i = 1:length(f)

K = (0.2) * randn(1,1) - 0;

P(i) = 10.^(K - 3.95*(10^-5)*f(i)); A(i) = sqrt(2*P(i)); end

xifft = ifft(A);

realx = real(xifft);

ti = [1:length(xifft)-1]/1000;

realx2(1:length(xifft)-1) = realx(2:length(xifft)); plot(ti,realx2)

2．样本点为1000、均值为0、标准差为2时，高斯白噪声分布为下图所示：

程序如下所示：

% White background nois clear all

f = 1:1:1000;

for i = 1:length(f)

K = (2) * randn(1,1) - 0;

P(i) = 10.^(K - 3.95*(10^-5)*f(i)); A(i) = sqrt(2*P(i)); end

xifft = ifft(A);

realx = real(xifft);

ti = [1:length(xifft)-1]/1000;

realx2(1:length(xifft)-1) = realx(2:length(xifft)); plot(ti,realx2)

3．样本点为1000、均值为0、标准差为10时，高斯白噪声分布为下图所示：

程序如下所示：

% White background nois clear all

f = 1:1:1000;

for i = 1:length(f)

K = (10) * randn(1,1) - 0;

P(i) = 10.^(K - 3.95*(10^-5)*f(i)); A(i) = sqrt(2*P(i)); end

xifft = ifft(A);

realx = real(xifft);

ti = [1:length(xifft)-1]/1000;

realx2(1:length(xifft)-1) = realx(2:length(xifft)); plot(ti,realx2)