资源描述：
《随机信号实验-高斯过程的产生.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、随机信号实验--高斯过程的产生赵骏一、实验目的掌握产生具有一定相关性的高斯过程的方法。二、实验内容（1）产生协方差函数为,的高斯过程。（2）分析其概率密度函数Error!Nobookmarknamegiven.、均值、方差、相关函数和功率谱。三、实验仪器和设备（1）计算机一台。（2）Matlab软件。四、实验原理一般情况下，认为所要进行仿真的随机过程都满足平稳性和各态历经性。严格的说，产生的随机序列并非是随机过程的样本函数，但是可以近似地认为是一个样本函数。仿真协方差函数为的高斯过程的方法如下：1．产生个均值为0，方差为1并且相互独立的高斯分布随机数。2．
2、根据递推公式计算出一组随机数，其中初值为，采样间隔是。3．如果要仿真的随机信号数学期望不为0，将数学期望加到随机数上，就可以得到对应均值和方差的随机过程或随机序列的一个样本，生成随机序列样本函数。4．估计随机信号的期望和方差分别为，。5．计算产生的随机过程的自相关函数为。6．计算产生的随机过程的功率谱密度为。7．利用自相关函数傅立叶变换计算随机信号的理论功率谱密度和估计值进行比较。理论方法计算功率谱密度为。五、实验步骤与结果记录（1）高斯过程的产生按照实验原理中所描述的产生高斯过程的步骤，编写实验程序如下：clc,clearN=250;alp=1;sig=
3、1;delt=1;a=1;mm=zeros(1,N);x1=rand(1,N);x2=rand(1,N);x3=a.*sqrt(-2*log(x1)).*cos(2*pi*x2)+mm;%产生随高斯分布的随机数y(1)=sig*x3(1);forn=2:N%产生高斯过程一个样本y(n)=exp(-alp)*y(n-1)+sig*sqrt(1-exp(-2*alp*delt))*x3(n);endi=1:N;plot(i,y);holdon;plot(i,mm,'-');title('高斯过程样本函数')运行程序后得到如图1所示的高斯过程的样本函数。图1高斯
4、过程的样本函数（2）高斯过程统计特性的分析以上已经产生了高斯过程的一个样本函数，设其具有各态历经性，则可以利用这个样本函数的统计特性来近似代替整个高斯过程的统计特性。首先，根据公式计算样本的均值。程序如下：（接着上面的程序）M=0;fori=1:NM=M+y(i);endM=M/N计算得出样本的均值为0.0102，进而可以近似地认为该高斯过程的均值也为0.0102。由高斯信号的概率密度函数形式，以及均值为0，方差为1可以得出该高斯过程的实际概率密度函数为。其均值为0，与实验得到的0.0102基本相同。然后，根据公式计算样本的方差。程序如下：（接着上面的程序
5、）D=0;fori=1:ND=D+(y(i)-M)^2;endD=D/N计算得出样本的均值为0.9886，进而可以近似地认为该高斯过程的均值也为0.9886。可以看出，它与实际的方差1还是十分接近的。接着，根据公式计算样本的自相关函数。程序如下：（接着上面的程序）form=1:N%计算自相关函数正半轴forn=1:N-m+1rr(n)=y(n)*y(n+m-1);endr2(m)=sum(rr)/N;endform=1:N-1%由对称性计算负半轴r1(N-m)=r2(m+1);endj=-N+1:N-1;plot(j,[r1,r2]);title('自相关
6、函数')运行程序后得到如图2所示的自相关函数。图2自相关函数通过观察以上得出的自相关函数可以发现它的几个特征：a．偶函数；b．在零点取得最大值。最后，根据公式计算样本的功率谱密度函数。程序如下：（接着上面的程序）ss=fft(y,512);s=(ss.*conj(ss))/512;%计算功率谱密度估值f=1000*(0:255)/512;plot(f,s(1:256))title('功率谱密度')运行程序后得到如图3所示的功率谱密度函数。图3功率谱密度函数（3）实际自相关函数和功率谱密度函数实际自相关函数为。编写程序如下：clc,cleart=-250:2
7、50;r=exp(-abs(t));plot(t,r)title('自相关函数')得出其图像如图4所示。图4实际自相关函数可以看出它与实验得出的自相关函数在形状和数值上大致相同。通过编写程序如下：clc,clearsymstf=exp(-abs(t));g=fourier(f)可以得到实际功率谱密度函数为。编写程序如下：clc,clearw=0:0.01:20;h=2./(w.^2+1);plot(w,h)title('实际功率谱密度')运行程序后得到实际功率谱密度如图5所示。图5实际功率谱密度函数图与样本函数得到的功率谱密度函数相比，可以看出，实验得到的
8、函数在总体的变化趋势上与实际的函数大致相同。但是，实验得到的函数，
