原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
查看更多
随机信号分析大作
班级：021252
学号姓名：谭红光
利用Matlab程序设计一正弦型信号加高斯白噪声的复合信号。
（1）分析复合信号的功率谱密度、幅度分布特性；
（2）分析复合信号通过RC积分电路后的功率谱密度和相应的幅度分布特性；
（3）分析复合信号通过理想低通系统后的功率谱密度和相应的幅度分布特性。
解：设定正选信号的频率为10HZ，抽样频率为100HZ
x=sin(2*pi*fc*t)
（1）正弦函数加上高斯白噪声：
y=awgn(x,10)
y的幅度分布特性可以通过傅里叶变换得到：
Y(jw)=fft(y)
y 的功率谱密度：
G(w)=Y(jw).*conj(Y(jw)/length(Y(jw)))
随机序列自相关函数的无偏估计公式为：
（2）复合信号 y通过RC积分电路后得到信号y2
通过卷积计算可以得到y2 即：y2= conv2(y,b*pi^-b*t)
y2的幅度分布特性可以通过傅里叶变换得到：
Y2(jw)=fft(y2)
y2的功率谱密度：
G2(w)=Y2(jw).*conj(Y2(jw)/length(Y2(jw)))
（3）复合信号 y通过理想滤波器电路后得到信号y3
通过卷积计算可以得到y3 即：y3=conv2(y,sin(10*t)/(pi*t))
y3的幅度分布特性可以通过傅里叶变换得到：
Y3(jw)=fft(y3)
y3的功率谱密度：
G3(w)=Y3(jw).*conj(Y3(jw)/length(Y3(jw)))
源程序如下：
clear all;
fs=100;
fc=10;
n=201;
t=0:1/fs:2;
x=sin(2*pi*fc*t); ? y=awgn(x,10);
10dBy1=a+b*randn(1,n);
m=50;
i=-0.49:1/fs:0.49;
for j=1:m
R(j)=sum(y(1:n-j-1).*y(j:199),2)/(n-j);
Ry(49+j)=R(j);
Ry(51-j)=R(j);
end
subplot(5,2,1);
plot(t,x,r);
title(x = sin(2*pi*fc*t));
ylabel(x);
xlabel(t/20pi);
grid;
subplot(5,2,2);
plot(t,y,r);
title(y=awgn(x,10));
ylabel(y);
xlabel(t/20pi);
grid;
subplot(5,2,3);
plot(i,Ry,r);
ylabel(Ry);
xlabel(i);
grid;
FY=fft(y);
FY1=fftshift(FY);
f=(0:200)*fs/n-fs/2;
subplot(5,2,5);
plot(f,abs(FY1),r);
ylabel(F(jw));
xlabel(w);
grid;
P=FY1.*conj(FY1)/length(FY1);
subplot(5,2,6);
plot(f,P,r);
ylabel(G(w));
xlabel(w);
grid;
b=10;
y2=conv2(y,b*pi^-b*t);
Fy2=fftshift(fft(y2));
f=(0:400)*fs/n-fs/2;
subplot(5,2,7);
plot(f,abs(Fy2),r);
ylabel(Fy2(jw));
xlabel(w);
grid;
P2=Fy2.*conj(Fy2)/length(Fy2);
subplot(5,2,8);
plot(f,P2,r);
ylabel(Gy2(w));
xlabel(w);
grid;
y3=conv2(y,sin(10*t)/(pi*t));
Fy3=fftshift(fft(y3));
f3=(0:200)*fs/n-fs/2;
subplot(5,2,9);
plot(f3,abs(Fy3),r);
ylabel(Fy3(jw));
xlabel(w);
grid;
P3=Fy3.*conj(Fy3)/length(Fy3);
subplot(5,2,10);
plot(f3,P3,r);
ylabel(Gy3(w));
xlabel(w);
grid;
图形如下：
文档评论（0）