matlab 白噪声实现

clc

fs=input('请输入信号的频率:fs=');

fc=5;    %抽样平率5Hz

x=(0:1/fs:2);

n=(fs*2)+1;

y1=sin(2*pi*fc*x);   %原正弦信号，频率为10

a=0;b=0.5;    %均值为a，方差为b^2

subplot(2,2,1);

plot(x,y1,'m');

title('y=sin(20pi*x)');

ylabel('y轴');

xlabel('x轴');

grid;

y2=a+b*randn(1,n);  %高斯白噪声

subplot(2,2,2);

plot(x,y2,'g');

title('N(0，0.25)的高斯白噪声');

ylabel('y轴');

xlabel('x轴');

grid;

y=y1+y2;  %加入噪声之后的信号

subplot(2,2,3);

plot(x,y,'r');

title('叠加了高斯白噪声的sinx');

ylabel('y轴');

xlabel('x轴');

grid;

FY=fft(y);  %傅里叶变换得出频谱函数

FY1=fftshift(FY);  %频谱校正

f=(1:n)*fs/n-fs/2;

subplot(2,2,4);

plot(f,abs(FY1),'b');

title('函数频谱图');

ylabel('F(jw)');

xlabel('w');

grid;d;