应用统计与随机过程实验报告

实验三线性系统对随机过程的响应一、实验目的通过本仿真实验了解正态白色噪声随机过程通过线性系统后相关函数以及功率谱的变化；培养计算机编程能力。二、实验要求采用MATLAB或VB语言进行编程1)运用正态分布随机数产生函数产生均值为零、根方差=1的白色噪声样本序列[或可参考实验1的正态分布产生方法]{u(n)|n=1,2,…,2000}；画出噪声u(n)的波形图。2)设离散时间线性系统的差分方程为画出x(n)的波形图。3)随机过程x(n)的理论上的功率谱密度函数为在[0,π]范围内对w进行采样，采样间隔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,π]范围内对w进行采样，采样间隔0.001π，计算S1(i×()()-0.36(-1)0.85(-2)(3,4,...,2000)xnunununn22()|10.360.85|jjSee200031ˆ()()()(0,1,2,3,4,5)1998XnmRmxnxnmmm1ˆˆˆ()(0)2(1)cos()2(2)cos(2)XXXSRRR0.001π)(i=1,2,…,1000);画出波形图；比较其与理论上的功率谱密度函数S(w)的差异。6)仿照实验1的方法统计数据x(n)在不同区间出现的概率,计算其理论概率，观察二者是否基本一致。三、实验代码及结果1.运用正态分布随机数产生函数产生均值为零、根方差=1的白色噪声样本序列[或可参考实验1的正态分布产生方法]{u(n)|n=1,2,…,2000}；画出噪声u(n)的波形图。代码：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)');波形图：分析：运用正态分布随机数产生函数产生均值为零、根方差=1的白色噪声样本序列。2.设离散时间线性系统的差分方程为画出x(n)的波形图。代码：n=3:2000;x(n)=u(n)-0.36*u(n-1)+0.85*u(n-2);stem(x,'.');title('x(n)');波形图：分析：正态随机序列通过线性离散系统生成的还是正态随机序列。()()-0.36(-1)0.85(-2)(3,4,...,2000)xnunununn3.随机过程x(n)的理论上的功率谱密度函数为在[0,π]范围内对w进行采样，采样间隔0.001π，计算S(i×0.001π)(i=1,2,…,1000);画出波形图。代码：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)');波形图：22()|10.360.85|jjSee4.根据步骤(2)产生的数据序列x(n)计算相关函数的估计值与理论值1.1296、-0.666、0.85、0、0、0的差异。代码：Rx=rand(1,6);form=1:1:6sum=0;forn=(3+m):1:2000sum=sum+x(n)*x(n-m+1);endRx(m)=sum/(1999-m);endS1=rand(1,1000);fori=1:1:1000S1(i)=Rx(1)+2*Rx(2)*cos(i*0.001*pi)+2*Rx(3)*cos(2*i*0.001*pi);endfigurestem(S1)运行结果：分析：所得的数据与理论值1.1296、-0.666、0.85、0、0、0存在一定的差异。200031ˆ()()()(0,1,2,3,4,5)1998XnmRmxnxnmmm5.根据相关函数的估计值对随机过程的功率谱密度函数进行估计在[0,π]范围内对w进行采样，采样间隔0.001π，计算S1(i×0.001π)(i=1,2,…,1000);画出波形图；比较其与理论上的功率谱密度函数S(w)的差异。代码：N=1000;P1=0;P2=0;P3=0;P4=0;forn=3:1:NIf(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;elseP4=P4+1;endendendendp1=P1/Np2=P2/Np3=P3/Np4=P4/Np=p1+p2+p3+p4figurehist(x,1000)return1ˆˆˆ()(0)2(1)cos()2(2)cos(2)XXXSRRR运行结果：分析：采样计算得到的功率谱密度函数比较其与理论上的功率谱密度函数相比，没有完全成偶对称。数据的概率分布没有理论那样均匀。6.分析：理论概率Rx=1.8315-0.64300.8528-0.0473-0.0096-0.0102。所以二者基本一致。