实验四窄带随机过程的产生及其性能测试一、实验目的、方差、相关函数及功率谱密度等。1、基于随机过程的莱斯表达式产生窄带随机过程。2、掌握窄带随机过程的特性，包括均值(数学期望)、实验原理 窄带随机过程的产生如下图所示:-b(t)sin(2pft)三、实验内容1、按照上面的结构框图，基于随机过程的莱斯表达式X(t) a(t)cos ot b(t)sin t，用Matlab产生一满足条件的窄带随机过程。m.文件如下: %产生一个p个点的高斯窄带随机过程 function f=suiji1(p) n=1:p;w=-pi:2*pi/1000:pi;R=100;C=0.001;wn=1/2*pi*R*C;b,a=butter(1,w n);g=ra ndn (1,1001);y=filter(b,a,g);at=y.*cos(w.* n);bt=y.*si n( w.* n);ft=at-bt;subplot(211)plot(ft)subplot(212)ksde nsity(ft)在comma nd命令框里写入：suiji1(1000)即产生一个1000个点的高斯窄带随机过程窄带随机过程波形及其概率密度图分别如下所示:210-12、画出该随机过程的若干次实现，观察其形状。 该随机过程的四次实现，代码如下： for i=1:1:4syms R C;n=1:1001;w=-pi:2*pi/1000:pi;R=100;C=0.001;wn=1/2*pi*R*C;b,a=butter(1,w n);g=ra ndn (1,1001);y=filter(b,a,g); at=y.*cos(w.* n); bt=y.*si n( w.* n); ft=at-bt; subplot(4,2,2*i-1) plot(ft) subplot(4,2,2*i) ksde nsity(ft) end形状如下：1000500150020-2500100015003、编写Matlab程序计算该随机过程的均值函数，自相关函数，功率谱，包络，包络平方及相位的一维概率密度画出相应的图形并给出解释。1 ）各个随机过程的均值：在 comma nd命令框中输入：mu=mea n(u)%均值计算,u是标准正态mat=mea n(at)%at， bt是u经过低通后的信号mbt=mea n(bt)mzt=mea n(zt)%计算窄带信号的均值运行结果：mu =-0.0171ma t=-2.9152e-004mbt =4.1458e-004mzt =-7.0610e-004分析：由于u是标准正态的，所以均值趋近于零，而at, bt是由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);title( Rbt);Subplot(2,2,4)plot(Rzt);title( Rz);运行结果：Rat6000400020000-200050001000015000图1各个随机过程的自相关函数分析：各个过程都是实的，中心点上相关程度最高，而且观察到：at，bt，zt三个过程在中心点两边对称位置上各有一个峰值，其他位置上，自相关函数接近于零。3) 功率谱密度 程序代码： subplot(2,2,1);periodogram(u); title( Sx);subplot(2,2,2) periodogram(at);title( Sa); subplot(2,2,3) periodogram(bt);title( Sb);Subplot(2,2,4) periodogram(zt); title( Sz);运行结果：益-dEO 里P&VCDP)斗*uanba 上/JUAAOd t-dLURfpetcnpj 乂 2LJ2nbal0vQdw:分20-20-40Sa-20 Ir-一叫口ornnalized FrequencySb-200Q匸nbp/Jo50 亡益UJezpe、Bs販u匚 mnbalnlda/sd/40so1 I0051 I【ormalized Firequen匚y (xn rad/sampleP-Normalized Frequency (xi rad/sample)Sz-80L0 0 6 1normalized Frequency (xi rad/samplej图2 各个过程的功率谱密度由u的功率谱可以看出，u是白噪声。4) 包络、包络平方、相位的一维概率密度 程序代码：A2=(at)A2+(btL2;A=A2.a(1/2);subplot(3,1,1);x=0:0.002:0.14;hist(A,x);axis tight ;xlabel(包络概率密度)subplot(3,1,2);x=0:0.0002:0.018;hist(A2,x);axis tight ;xlabel(包络平方概率密度)subplot(3,1,3)Q=bt./at;Qt=ata n(Q);x=-2:0.025:2;hist(Qt,x);axis tight ;xlabel(相位概率密度)运行结果:20010000 020.040.060 OS 010.120 14包络概率密度10&相位概率密度图3概率密度图像分析：从上图可以看出，窄带过程包络服从瑞丽分布，包络平方服从指数分布。总体分析：由上图可看出，随机过程包络的一维概率密度服从瑞利分布，相位的一维概率密度服从均匀分布，而包络平方的一维概率密度近似指数分布，与书本上的结论相符合。四、实验心得（1）能够基于随机过程的莱斯表达式产生窄带随机过程;(2) 通过做图，掌握了窄带随机过程的特性，包括数学期望、方差、相关函数及功率谱 密度等。(3) 通过观察和分析所得图形，验证了书本上学到的有关窄带随机过程的包络、包络平 方和相位的一维概率密度的知识。