matlab产生窄带高斯随机过程,窄带随机过程的产生及其性能测试

实验四 窄带随机过程的产生及其性能测试

一、实验目的

1、基于随机过程的莱斯表达式产生窄带随机过程。

2、掌握窄带随机过程的特性，包括均值(数学期望)、方差、相关函数及功率谱密度等。 二、实验原理

窄带随机过程的产生如下图所示：

00)()sin(2)

f t b t f t p p -

三、实验内容

1、按照上面的结构框图，基于随机过程的莱斯表达式t t b t t a t X 00sin )(cos )()(ωω-=，用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)

在command 命令框里写入： suiji1(1000)

即产生一个1000个点的高斯窄带随机过程

窄带随机过程波形及其概率密度图分别如下所示：