原创力文档- 1、本文档共11页，可阅读全部内容。
- 2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
- 3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
随机信号分析大作业
學院：
電子工程學院
指导教師：
李琦
12月11曰
1.设有随机初相信号Xt=5cos?(
代码：%sjf均匀分布的随机相位
sjf=unifrnd(0,2*pi,4,1);
fork=1:3
t=1:0.1:10;
x=5*cos(t+sjf(k));
plot(t,x);
holdon
end
xlabel(t);
ylabel(X(t));
gridon;
axistight;
仿真成果：
随机初相信号的三個样本函数
2．运用matlab程序设计设计
（1）分析复合信号的功率谱密度，幅度分布的特性；
（2）分析复合信号通過RC积分電路後的功率谱密度和對应的幅度分布特性；
（3）分析复合信号通過理想低通系统後的功率谱密度和對应的幅度分布特性。
答：（1）
代码：%初始条件:正弦信号频率為10Hz，采样频率為100Hz,加入高斯白噪声的信噪比為10
fc=20;
fs=100;
XZB=10;
T=2;
t=0:0.01:2;
n=2*fs+1;
x=sin(2*pi*fc*t);
y=awgn(x,XZB);
figure(1);
hist(y,100);
FF1=fft(y);
FF2=fftshift(FF1);
f=(0:n-1)*fs/n-fs/2;
[Cy,a]=xcorr(y,unbiased);
a=a*T/n;
figure(2);
Cyk=abs(fftshift(fft(Cy,n)));
plot(f,Cyk);
xlabel(频率);
ylabel(G(f));
gridon;
仿真成果：
复合信号的功率谱幅度分布直方图
复合信号的双边功率谱密度曲线
（2）代码：clearall;
fc=20;
fs=100;
XZB=10;
T=2;
t=0:0.01:2;
n=2*fs+1;
x=sin(2*pi*fc*t);
y=awgn(x,XZB);
R=2;C=0.5;%设置积分電路的参数
a=1/(R*C);
ht=a*exp(-a*t);%RC积分電路的冲激响应
yo=conv2(y,ht);
yo=yo(2*fs+1:4*fs+1);
figure(1);
plot(t,y);
xlabel(t);ylabel(幅度);
figure(2);
hist(yo,100);
ylabel(频率/组距);xlabel(组距);
gridon;
holdon;
FY1=fft(yo);
FY2=fftshift(FY1);
f=(0:n-1)*fs/n-fs/2;
[Cy,a]=xcorr(yo,unbiased);
a=a*T/n;
figure(3);
Cyk=abs(fftshift(fft(Cy,n)));
plot(f,Cyk);
ylabel(G(f));xlabel(频率f/Hz);
gridon;
仿真成果：
复合信号的時域波形
复合信号RC积分後的幅度分布直方图
复合信号RC积分後的双边功率谱密度曲线
（3）代码：clear;clc;
fc=10;
fs=100;
XZB=10;
T=2;
t=0:1/fs:2;
n=2*fs+1;
x=sin(2*pi*fc*t);%输入信号
yi=awgn(x,XZB);%输入的复合信号
figure(1);
plot(t,yi);
xlabel(t);ylabel(幅度);
figure(2);
N=32;
Wn=0.24;
b=fir1(N,Wn);
yo=filter(b,1,yi);
plot(t,yo);
xlabel(t);ylabel(幅度);
figure(3);
hist(yo,100);
ylabel(频率/组距);xlabel(组距);
gridon;
holdon;
FF_Y1=fft(yo);
FF_Y2=fftshift(FF_Y1);
f=(0:n-1)*fs/n-fs/2;
figure(4);
plot(f,abs(FF_Y2));
ylabel(Yo(f));xlabel(频率f/Hz);gridon;
[Cy,a]=xcorr(yo,unbiased);
a=a*T/n;
figure(5);
Cyk=abs(fftshift(fft(Cy,n)));
plot(f,Cyk);
ylabel(G(f));xlabel(频率f/Hz);
gridon;
仿真成果：
复合信号的時域波形
复合信号低通滤波後的時域波形
复合信号低通滤波後的幅度分布直方图
复合信号低通滤波後的频谱特性
复合信号低通滤波後的双边功率谱密度曲线
文档评论（0）