当前位置:文档之家› 随机信号分析大作业:利用matlab程序设计设计一个正旋信号加高斯白噪声的复合信号

随机信号分析大作业:利用matlab程序设计设计一个正旋信号加高斯白噪声的复合信号

班级:021012班学号:0210111X姓名:李X
随机信号大作业
利用matlab程序设计设计一个正旋信号加高斯白噪声的复合信号。

源代码:
正旋sinx信号
x=(0:0.01:2);
y1=sin(10*pi*x);
plot(x,y1,'r');
title('y=sin(10*pi*x)');
ylabel('y');
xlabel('x/10pi');
grid;正旋信号如下图:
高斯白噪声,当白噪声的方差为10的-4次方时
y2=0.01*randn(1,201);
plot(x,y2,'r');
title('高斯白噪声');
ylabel('y');
xlabel('x/10pi');
grid;
1:加入噪声之后的信号。

因为白噪声方差过小,变化过快,叠加信号大致显示
出的是正旋波形
y=y1+y2;
plot(x,y,'r');
title('叠加了高斯白噪声的sinx');
ylabel('y');
xlabel('x/10pi');
grid;
2:当增加白噪声的方差到0.01时,观察复合信号,可以发现,复合信号波形没有第一次平滑,较第一个复合波形而言更显无序性
3:方差增加到1时:
取方差为0.01时的白噪声作函数频谱图和白噪声自相关函数图:FY=fft(y);
FY1=fftshift(FY);
f=(0:200)*100/201-50;
subplot(1,2,1);
plot(f,abs(FY1),'r');
ylabel('F(jw)');
xlabel('w');
grid;
i=-0.49:1/100:0.49;
for j=1:50
R(j)=sum(y2(1:201-j-1).*y2(j:199),2)/(201-j);
Rx(49+j)=R(j);
Rx(51-j)=R(j);
end
subplot(1,2,2);
plot(i,Rx,'r');
ylabel('Rx');
xlabel('x');
grid;
左图为函数频谱图;右图为白噪声自相关函数图。

相关主题