随机过程实验4.docx
- 文档编号:5915842
- 上传时间:2023-01-02
- 格式:DOCX
- 页数:12
- 大小:232.85KB
随机过程实验4.docx
《随机过程实验4.docx》由会员分享,可在线阅读,更多相关《随机过程实验4.docx(12页珍藏版)》请在冰豆网上搜索。
随机过程实验4
课程实验报告
课程名称:
随机过程实验
实验项目名称:
线性系统对随机过程的响应
专业班级:
通信工程1301班
姓名:
王少丹
学号:
201308030104
指导教师:
何松华
1.
实验目的
通过本仿真实验了解正态白色噪声随机过程通过线性系统后相关函数以及功率谱的变化;培养计算机编程能力。
2.实验要求
采用MATLAB或VB语言进行编程
(1)运用正态分布随机数产生函数产生均值为零、根方差=1
的白色噪声样本序列[或可参考实验1的正态分布产生方法]
{u(n)|n=1,2,…,2000};画出噪声u(n)的波形图
(2)设离散时间线性系统的差分方程为
画出x(n)的波形图。
(3)随机过程x(n)的理论上的功率谱密度函数为
在[0,]范围内对进行采样,采样间隔0.001,计算S(i×0.001)(i=1,2,…,1000);画出波形图。
(4)根据步骤
(2)产生的数据序列x(n)计算相关函数的估计值
与理论值1.1296、-0.666、0.85、0、0、0的差异。
(5)根据相关函数的估计值对随机过程的功率谱密度函数进行估计
在[0,]范围内对进行采样,采样间隔0.001,计算S1(i×0.001)(i=1,2,…,1000);画出波形图;比较其与理论上的功率谱密度函数S()的差异。
(6)仿照实验1的方法统计数据x(n)在不同区间出现的概率,计算其理论概率,观察二者是否基本一致。
3.程序代码
(见后)
4.实验结果
可以看出噪声的随机性
均值为零、根方差=1的白色噪声样本序列即上图所示
的离散时间线性系统即如上图所示。
随机过程的功率谱密度函数图像均匀分布,函数曲线光滑。
与理论值相较,有一定差距,但是基本一样。
在[0,]范围内对进行采样,采样间隔0.001,计算S1(i×0.001)
与上次图片相较,不如上次的曲线光滑,但是大致走向相同,说明取样和理论有差距但是略有不同。
仿照实验1的方法统计数据x(n)在不同区间出现的概率,计算其理论概率,观察到二者基本一致。
5.实验体会
通过本仿真实验了解了正态白色噪声随机过程通过线性系统后相关函数以及功率谱的变化;培养了计算机编程能力。
6.完整代码附录:
3_1
n=1:
2000;
u1(n)=rand(1,2000);
u2(n)=rand(1,2000);
u(n)=sqrt(-2*log(u1(n))).*cos(2*pi*u2(n));
stem(u,'.');
title('u(n)');
3_2
n=3:
2000;
x(n)=u(n)-0.36*u(n-1)+0.85*u(n-2);
stem(x,'.');
title('x(n)');
3_3
i=1:
1000;
w=0.001*pi.*i;
s=(abs(1-0.36.*exp((-1j).*w)+0.85.*exp((-2j).*w))).*(abs(1-0.36.*exp((-1j).*w)+0.85.*exp((-2j).*w)));
stem(s,'.');
title('s(i*0.001*pi)');
3_4
N=2000;
Rx=rand(1,6);
form=1:
1:
6
sum=0;
forn=(2+m):
1:
N
sum=sum+x(n)*x(n-m+1);
end
Rx(m)=sum/(1999-m);
end
Rx
3_5
S1=rand(1,1000);
fori=1:
1:
1000
S1(i)=Rx
(1)+2*Rx
(2)*cos(i*0.001*pi)+2*Rx(3)*cos(2*i*0.001*pi);
end
figure
stem(S1)
title('S1(w),w=[0,pi]')
3_6
N=2000;
P1=0;P2=0;P3=0;P4=0;
forn=3:
1:
N
if(x(n)<-1)
P1=P1+1;
elseif(x(n)>=-1&x(n)<=0)
P2=P2+1;
elseif(x(n)>0&x(n)<=1)
P3=P3+1;
else
P4=P4+1;
end
end
end
end
p1=P1/N
p2=P2/N
p3=P3/N
p4=P4/N
p=p1+p2+p3+p4
figure
hist(x,1000)
return
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 随机 过程 实验