随机信号大作业

《随机信号大作业》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、随机信号大作业021151班02115091薛雪歌第一章上机题：设有随机初相信号X(t)=5cos（t+），其中相位是在区间（0,2）上均匀分布的随机变量。试用Matlab编程产生其三个样本函数。解：由Matlab产生的三个样本函数如下图所示：程序源代码：clcclearm=unifrnd(0,2*pi,1,10);fork=1:3t=1:0.1:10;X=5*cos(t+m(k));plot(t,X);holdonendxlabel('t');ylabel('X(t)');gridon;axistight;第二章上机题：利用Matlab程序设

2、计一正弦型信号加高斯白噪声的复合信号。（1）分析复合信号的功率谱密度，幅度分布的特性；（2）分析复合信号通过RC积分电路后的功率谱密度和相应的幅度分布特性；（3）分析复合信号通过理想低通系统后的功率谱密度和相应的幅度分布特性。解：设正弦信号的频率为10HZ，抽样频率为100HZx=sin(2*pi*fc*t)正弦曲线图：程序块代码：clearall;fs=100;fc=10;n=201;t=0:1/fs:2;x=sin(2*pi*fc*t);y=awgn(x,10);m=50;i=-0.49:1/fs:0.49;forj=1:mR(j)=sum

3、(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('正弦信号曲线');ylabel('x');xlabel('t/20pi');grid;（1）正弦信号加上高斯白噪声产生复合信号y：y=awgn(x,10)对复合信号进行傅里叶变换得到傅里叶变换：Y(jw)=fft(y)复合信号的功率谱密度函数为：G(w)=Y(jw).*conj(Y(jw)/length(Y(jw)))复合信号的曲线图，频谱图和功率谱图

4、：程序块代码：plot(t,y,'r');title('复合信号曲线');ylabel('y');xlabel('t/20pi');grid;程序块代码：FY=fft(y);FY1=fftshift(FY);f=(0:200)*fs/n-fs/2;plot(f,abs(FY1),'r');title('复合信号频谱图');ylabel('F(jw)');xlabel('w');grid;程序块代码：P=FY1.*conj(FY1)/length(FY1);plot(f,P,'r');title('复合信号功率谱密度图');ylabel('G(w

5、)');xlabel('w');grid;（2）正弦曲线的复合信号通过RC积分电路后得到信号为：通过卷积计算可以得到y2即：y2=conv2(y,b*pi^-b*t)y2的幅度分布特性可以通过傅里叶变换得到Y2(jw)=fft(y2)y2的功率谱密度G2(w)=Y2(jw).*conj(Y2(jw)/length(Y2(jw)))复合信号通过RC积分电路后的曲线频谱图和功率谱图：程序块代码：b=10;y2=conv2(y,b*pi^-b*t);Fy2=fftshift(fft(y2));f=(0:400)*fs/n-fs/2;plot(f,ab

6、s(Fy2),'r');title('复合信号通过RC系统后频谱图');ylabel('Fy2(jw)');xlabel('w');grid;程序代码：P2=Fy2.*conj(Fy2)/length(Fy2);plot(f,P2,'r');title('复合信号通过ＲＣ系统后功率密度图');ylabel('Gy2(w)');xlabel('w');grid;（3）复合信号y通过理想滤波器电路后得到信号y3通过卷积计算可以得到y3即：y3=conv2(y,sin(10*t)/(pi*t))y3的幅度分布特性可以通过傅里叶变换得到Y3(jw)=f

7、ft(y3)，y3的功率谱密度G3(w)=Y3(jw).*conj(Y3(jw)/length(Y3(jw)))复合信号通过理想滤波器后的频谱图和功率密度图：程序块代码：y3=conv2(y,sin(10*t)/(pi*t));Fy3=fftshift(fft(y3));f3=(0:200)*fs/n-fs/2;plot(f3,abs(Fy3),'r');title('复合信号通过理想滤波器频谱图');ylabel('Fy3(jw)');xlabel('w');grid;程序块代码：P3=Fy3.*conj(Fy3)/length(Fy3);p

8、lot(f3,P3,'r');title('理想信号通过理想滤波器功率密度图');ylabel('Gy3(w)');xlabel('w');grid;