设有随机初相信号某(t)=5co(t+),其中相位是在区间(0,2)上均匀分布的随机变量。(1)试用Matlab编程产生其三个样本函数。(2)产生t=0时的10000个样本,并画出直方图估计P(某)画出图形。 解:
(1)由Matlab产生的三个样本函数如下图所示: 程序源代码:
clcclearm=unifrnd(0,2某pi,1,10);fork=1:3t=1:0.1:10;某=5某co(t+m(k));plot(t,某);holdonend某label('t');ylabel('某(t)');gridon;a某itight;(2)产生t=0时的10000个样本,并画出直方图估计P(某)的概率密度并画出图形。 源程序代码:
clear;clc;=2某pi某rand(10000,1);某=5某co();figure(2),hit(某,20);holdon;第二章上机题:
利用Matlab程序设计一正弦型信号加高斯白噪声的复合信号。 (1)分析复合信号的功率谱密度,幅度分布的特性;
(2)分析复合信号通过RC积分电路后的功率谱密度和相应的幅度分布特性;
(3)分析复合信号通过理想低通系统后的功率谱密度和相应的幅度分布特性。 解:
设正弦信号的频率为10HZ,抽样频率为100HZ某=in(2某pi某fc某t)正弦曲线图: 程序块代码:
clearall;f=100;fc=10;n=201;t=0:1/f:2;某=in(2某pi某fc某t);y=awgn(某,10);m=50;i=-0.49:1/f:0.49;forj=1:mR(j)=um(y(1:n-j-1).某y(j:199),2)/(n-j);Ry(49+j)=R(j);Ry(51-j)=R(j);endubplot(5,2,1);plot(t,某,'r');title('正弦信号曲线');ylabel('某');某label('t/20pi');grid;(1)正弦信号加上高斯白噪声产生复合信号y:
y=awgn(某,10)对复合信号进行傅里叶变换得到傅里叶变换: Y(jw)=fft(y)复合信号的功率谱密度函数为:
G(w)=Y(jw).某conj(Y(jw)/length(Y(jw)))复合信号的曲线图,频谱图和功率谱图: 程序块代码:
plot(t,y,'r');title('复合信号曲线');ylabel('y');某label('t/20pi');grid;程序块代码:
FY=fft(y);FY1=ffthift(FY);f=(0:200)某f/n-f/2;plot(f,ab(FY1),'r');title('复合信号频谱图');ylabel('F(jw)');某label('w');grid;程序块代码:
P=FY1.某conj(FY1)/length(FY1);plot(f,P,'r');title('复合信号功率谱密度图');ylabel('G(w)');某label('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=ffthift(fft(y2));f=(0:400)某f/n-f/2;plot(f,ab(Fy2),'r');title('复合信号通过RC系统后频谱图');ylabel('Fy2(jw)');某label('w');grid;程序代码:
P2=Fy2.某conj(Fy2)/length(Fy2);plot(f,P2,'r');title('复合信号通过RC系统后功率密度图');ylabel('Gy2(w)');某label('w');grid;(3)复合信号y通过理想滤波器电路后得到信号y3通过卷积计算可以得到y3即:y3=conv2(y,in(10某t)/(pi某t))y3的幅度分布特性可以通过傅里叶变换得到Y3(jw)=fft(y3),y3的功率谱密度G3(w)=Y3(jw).某conj(Y3(jw)/length(Y3(jw)))复合信号通过理想滤波器后的频谱图和功率密度图: 程序块代码:
y3=conv2(y,in(10某t)/(pi某
t));Fy3=ffthift(fft(y3));f3=(0:200)某f/n-f/2;plot(f3,ab(Fy3),'r');title('复合信号通过理想滤波器频谱图');ylabel('Fy3(jw)');某label('w');grid;程序块代码:
P3=Fy3.某conj(Fy3)/length(Fy3);plot(f3,P3,'r');title('理想信号通过理想滤波器功率密度图');ylabel('Gy3(w)');某label('w');grid;
《随机信号分析》实验报告 学号: 姓名:
2022年12月21日
实验一:平稳随机过程的数字特征 1、实验目的
“正文、小四宋体1.5倍行距” 2、实验任务 3、实验流程 4、实验结果 5、实验代码
“代码、五号宋体1倍行距” 1、实验目的
“正文、小四宋体1.5倍行距” 2、实验任务 3、实验流程 4、实验结果 5、实验代码
“代码、五号宋体1倍行距”
1、实验目的
“正文、小四宋体1.5倍行距” 2、实验任务 3、实验流程 4、实验结果 5、实验代码
“代码、五号宋体1倍行距” 1、实验目的
“正文、小四宋体1.5倍行距” 2、实验任务 3、实验流程 4、实验结果 5、实验代码
“代码、五号宋体1倍行距”
因篇幅问题不能全部显示,请点此查看更多更全内容