窄带随机过程地产生及其性能测试

窄带随机过程地产生及其性能测试

ID:29450828

大小:119.50 KB

页数:7页

时间:2018-12-20

窄带随机过程地产生及其性能测试_第1页
窄带随机过程地产生及其性能测试_第2页
窄带随机过程地产生及其性能测试_第3页
窄带随机过程地产生及其性能测试_第4页
窄带随机过程地产生及其性能测试_第5页
资源描述:

《窄带随机过程地产生及其性能测试》由会员上传分享,免费在线阅读,更多相关内容在工程资料-天天文库

1、实用标准文案实验四窄带随机过程的产生及其性能测试一、实验目的1、基于随机过程的莱斯表达式产生窄带随机过程。2、掌握窄带随机过程的特性,包括均值(数学期望)、方差、相关函数及功率谱密度等。二、实验原理窄带随机过程的产生如下图所示:三、实验内容1、按照上面的结构框图,基于随机过程的莱斯表达式,用Matlab产生一满足条件的窄带随机过程。m.文件如下:%产生一个p个点的高斯窄带随机过程functionf=suiji1(p)n=1:p;w=-pi:2*pi/1000:pi;R=100;C=0.001;wn=1/2*

2、pi*R*C;[b,a]=butter(1,wn);g=randn(1,1001);y=filter(b,a,g);at=y.*cos(w.*n);bt=y.*sin(w.*n);ft=at-bt;subplot(211)plot(ft)subplot(212)ksdensity(ft)在command命令框里写入:suiji1(1000)即产生一个1000个点的高斯窄带随机过程窄带随机过程波形及其概率密度图分别如下所示:精彩文档实用标准文案2、画出该随机过程的若干次实现,观察其形状。该随机过程的四次实现,

3、代码如下:>>fori=1:1:4symsRC;n=1:1001;w=-pi:2*pi/1000:pi;R=100;C=0.001;wn=1/2*pi*R*C;[b,a]=butter(1,wn);g=randn(1,1001);y=filter(b,a,g);at=y.*cos(w.*n);bt=y.*sin(w.*n);ft=at-bt;subplot(4,2,2*i-1)plot(ft)subplot(4,2,2*i)ksdensity(ft)end形状如下:精彩文档实用标准文案3、编写Matlab程

4、序计算该随机过程的均值函数,自相关函数,功率谱,包络,包络平方及相位的一维概率密度画出相应的图形并给出解释。1)各个随机过程的均值:在command命令框中输入:mu=mean(u)%均值计算,u是标准正态mat=mean(at)%at,bt是u经过低通后的信号mbt=mean(bt)mzt=mean(zt)%计算窄带信号的均值运行结果:mu=-0.0171mat=-2.9152e-004mbt=4.1458e-004mzt=-7.0610e-004分析:由于u是标准正态的,所以均值趋近于零,而at,bt是

5、由u通过一个线性系统得到的,所以输出均值不变,仍为零,从程序运行结果可以看出,u,at,bt均值都趋近于零。精彩文档实用标准文案2)自相关函数:程序代码:Rx=xcorr(u);%分别计算出u,at,bt,zt的自相关Rat=xcorr(at);Rbt=xcorr(bt);Rzt=xcorr(zt);subplot(2,2,1);plot(Rx);title('Rx');%并绘制图形subplot(2,2,2)plot(Rat);title('Rat');subplot(2,2,3)plot(Rbt);ti

6、tle('Rbt');Subplot(2,2,4)plot(Rzt);title('Rz');运行结果:图1各个随机过程的自相关函数分析:各个过程都是实的,中心点上相关程度最高,而且观察到:at,bt,zt三个过程在中心点两边对称位置上各有一个峰值,其他位置上,自相关函数接近于零。3)功率谱密度程序代码:subplot(2,2,1);精彩文档实用标准文案periodogram(u);title('Sx');subplot(2,2,2)periodogram(at);title('Sa');subplot(2

7、,2,3)periodogram(bt);title('Sb');Subplot(2,2,4)periodogram(zt);title('Sz');运行结果:图2各个过程的功率谱密度分析:由u的功率谱可以看出,u是白噪声。4)包络、包络平方、相位的一维概率密度程序代码:A2=(at).^2+(bt).^2;A=A2.^(1/2);subplot(3,1,1);x=0:0.002:0.14;hist(A,x);axistight;精彩文档实用标准文案xlabel('包络概率密度')subplot(3,1,2

8、);x=0:0.0002:0.018;hist(A2,x);axistight;xlabel('包络平方概率密度')subplot(3,1,3)Q=bt./at;Qt=atan(Q);x=-2:0.025:2;hist(Qt,x);axistight;xlabel('相位概率密度')运行结果:图3概率密度图像分析:从上图可以看出,窄带过程包络服从瑞丽分布,包络平方服从指数分布。总体分析:由上图可看出,随机过程包络

当前文档最多预览五页,下载文档查看全文

此文档下载收益归作者所有

当前文档最多预览五页,下载文档查看全文
温馨提示:
1. 部分包含数学公式或PPT动画的文件,查看预览时可能会显示错乱或异常,文件下载后无此问题,请放心下载。
2. 本文档由用户上传,版权归属用户,天天文库负责整理代发布。如果您对本文档版权有争议请及时联系客服。
3. 下载前请仔细阅读文档内容,确认文档内容符合您的需求后进行下载,若出现内容与标题不符可向本站投诉处理。
4. 下载文档时可能由于网络波动等原因无法下载或下载错误,付费完成后未能成功下载的用户请联系客服处理。