资源描述：
《西电随机信号分析大作业》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、随机信号分析大作业学院：电子工程学院班级：021151学号：02115037姓名：隋伟哲第一题：设有随机信号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));%将产生的随机相位角
2、逐一代入随机过程中plot(t,y1,'r-');holdon;plot(t,y2,'g--');holdon;plot(t,y3,'k-');xlabel('t');ylabel('X(t)');gridon;axis([030-88]);title('随机相位的三条样本曲线');产生的三条样本曲线：第二题：利用Matlab程序设计一正弦型信号加高斯白噪声的复合信号。（1）分析复合信号的功率谱密度、幅度分布特性；（2）分析复合信号通过RC积分电路后的功率谱密度和相应的幅度分布特性；（3）分析复合信号通过理想低通系统后的功率谱密度和相应的幅度分布特
3、性。解：设定正选信号的频率为10HZ，抽样频率为100HZx=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
4、(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)))源程序如下：clearall;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(
5、1,n);m=50;i=-0.49:1/fs:0.49;forj=1:mR(j)=sum(y(1:n-j-1).*y(j:199),2)/(n-j);Ry(49+j)=R(j);Ry(51-j)=R(j);endsubplot(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');gr
6、id;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=
7、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));
8、f3=(0:200)*fs/n-fs/2;subplot(5,2,9);plot(f3,abs(Fy3),'r');yla
