2022年西电随机信号分析大作业参考 .pdf
随机信号分析大作业学院：电子工程学院班级： 021151 学号： 02115037 姓名：隋伟哲名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 1 页，共 9 页 - - - - - - - - - 第一题 ：设有随机信号X（t）=5cos(t+a),其中相位a 是在区间（ 0,2）上均匀分布的随机变量，使用Matlab 编程产生其三个样本函数。解：源程序如下：clc;clear; C=2*pi*rand(1,3);% 在0,2产生均匀分布的相位角t=1:.1:80; y1=5*cos(t+C(1); %将产生的随机相位角逐一代入随机过程中y2=5*cos(t+C(2); %将产生的随机相位角逐一代入随机过程中y3=5*cos(t+C(3); %将产生的随机相位角逐一代入随机过程中plot(t,y1,r-); hold on; plot(t,y2,g-); hold on; plot(t,y3,k-); xlabel(t);ylabel(X(t); grid on;axis(0 30 -8 8); title(随机相位的三条样本曲线); 产生的三条样本曲线：名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 2 页，共 9 页 - - - - - - - - - 第二题 ：利用 Matlab 程序设计一正弦型信号加高斯白噪声的复合信号。 （1）分析复合信号的功率谱密度、幅度分布特性；（2）分析复合信号通过RC 积分电路后的功率谱密度和相应的幅度分布特性；（3） 分析复合信号通过理想低通系统后的功率谱密度和相应的幅度分布特性。解：设定正选信号的频率为10HZ，抽样频率为 100HZ x=sin(2*pi*fc*t) 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 3 页，共 9 页 - - - - - - - - - （1）正弦函数加上高斯白噪声：y=awgn(x,10) y 的幅度分布特性可以通过傅里叶变换得到：Y(jw)=fft(y) y 的功率谱密度：G(w)=Y(jw).*conj(Y(jw)/length(Y(jw) 随机序列自相关函数的无偏估计公式为：101()( ) ()NmxxnRmx n x nmNm01mN（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) 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 4 页，共 9 页 - - - - - - - - - 源程序如下：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); 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 5 页，共 9 页 - - - - - - - - - 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); 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 6 页，共 9 页 - - - - - - - - - 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); 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 7 页，共 9 页 - - - - - - - - - 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; 图形如下：名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 8 页，共 9 页 - - - - - - - - - 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 9 页，共 9 页 - - - - - - - - -
