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

窄带随机过程地产生及其性能测试 - - - - - -精品可编辑word学习资料 gF3A9Q8X4G8 — — hA2Q9W10Y1N5 — — lQ7H3U2W4D4 试验四 窄带随机过程的产生与其性能测试 一、试验目的 1、基于随机过程的莱斯表达式产生窄带随机过程； 2、把握窄带随机过程的特性，包括均值〔数学期望〕 、方差、相关函数与功率谱密度等；二、试验原理 窄带随机过程的产生如如下图所示： |精. |品. |可. |编. |辑. |学. |习. |资. |料. * | * | * | * | |欢. |迎. |下. |载.  高斯白噪声 高斯白噪声 三、试验内容  低通滤波器 低通滤波器  a(t ) b(t ) cos(2 sin(2 f0t ) f0t )  a(t )cos(2  f0t)  b(t )sin(2  f0t) 1、根据上面的结构框图， 基于随机过程的莱斯表达式 X (t) a(t ) cos 0t b(t) sin 0t ，用 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,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) 在 mand命令框里写入： suiji1(1000) 即产生一个 1000 个点的高斯窄带随机过程 窄带随机过程波形与其概率密度图分别如下所示： |精. |品. |可. |编. |辑. |学. |习. |资. |料. * | * | * | * | |欢. |迎. |下. |载.  1 21 2 1 0 -1 -2 0 200 400 600 800 1000 1200 .5 1 .5 0 -2.5 -2 -1.5 -1 -0.5 0 0.5 1 1.5 2 2、画出该随机过程的假如干次实现，观看其外形；该随机过程的四次实现，代码如下： >> for i=1:1:4 syms 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,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 外形如下： 21 2 1 0 500 1000 1500 0 -4 2 -2 0 2 4 1 0 500 1000 1500 0 -2 2 0 2 4 0 -5 5 0 -5 2 |精. |品. |可. |编. 0 1 |辑. |学. |习. |资. |料. * -2 0 500 1000 1500 0 -4 -2 0 2 4 | * 2 2 | * | * 0 1 | |欢. |迎. |下. |载. -2 0 500 1000 1500 0 -2 0 2 4 3、编写 Matlab 程序运算该随机过程的均值函数，自相关函数，功率谱，包络，包络平方与相位的一维概率密度画出相应的图形并给出说明； 1〕各个随机过程的均值： 在 mand命令框中输入： mu=mean(u) % 均值运算 ,u 是标准正态mat=mean(at) %at ， bt 是u经过低通后的信号mbt=mean(bt) mzt=mean(zt) % 运算窄带信号的均值 运行结果： mu = ma t= mbt = mzt = 分析： 由于 u 是标准正态的， 所以均值趋近于零， 而 at ，bt 是由 u 通过一个线性系统得到的， 所以输出均值不变，仍为零，从程序运行结果可以看出， u，at ， bt 均值都趋近于零； 2〕自相关函数： 2〕自相关函数： 程序代码： Rx=xcorr(u); Rat=xcorr(at); Rbt=xcorr(bt); Rzt=xcorr(zt); % 分别运算出 u， at ， bt ， zt 的自相关 subplot(2,2,1); plot(Rx);title( subplot(2,2,2) plot(Rat);title( subplot(2,2,3) plot(Rbt);ti