窄带高斯白噪声.pdfVIP

窄带高斯随机信号功率谱的产生

摘要

基于VisualStudio的开发环境，本文首先通过线性同余法产生高斯白噪声。接着利用FIR

数字滤波器的原理，产生低通滤波器。然后让高斯信号过滤波器，最后通过莱斯表达式的方

式表示窄带高斯白噪声，求出其对应的自相关函数，利用维纳—辛钦公式，对自相关函数进

行傅里叶变换便可以得出功率谱密度函数。

一、问题重述、分析

画出窄带高斯白噪声的功率谱。

要画出窄带高斯白噪声的功率谱，最主要的问题是求出它的自相关函数，然后通过维纳

—辛钦公式便可求出。而对于自相关函数的求解主要是得出窄带高斯信号，所以本文基于莱

斯表达式，求出它的窄带高斯信号。具体流程如图（1）：

sinwt

0

低通滤

Nt

a波器

减

窄带高斯随

高斯白噪声加法

器机信号

加

低通滤

Nt

b波器

coswt

0

图(1)算法流图

二、问题的解决与算法分析

1、产生高斯白噪声

基于前一次大作业，我们知道产生白噪声的方法有线性同余法。

voidgauss(floatb[])

{

SYSTEMTIMEsys;

GetLocalTime(sys);

intseed;

floata[100];

seed=sys.wMilliseconds;//随机种子输入,产生范围为0—2^31-1

inti,j,w,k;

floatx,y;

for(w=0;w5000;w++)

for(k=0;k20;k++)//5000*20个均匀分布的统计近似服从正态分布

{

x=0;y=0;

for(i=0;i100;i++)

{

seed=16807*(seed%127773);//利用乘同余法产生其他随机数,共100个

if(seed0)

seed+=2147483647;//处理溢出的情况

a[i]=seed/2147483647.0;//a即为0-1均匀分布数

}

for(j=20;j32;j++)

{

x=a[j]+x;//取其中12的个数的和，方便计算

}

y=x-6;//产生标准正态分布随机数

for(j=0;