窄带高斯随机信号功率谱的产生
摘要
基于Visual Studio的开发环境，本文首先通过线性同余法产生高斯白噪声。接着利用FIR数字滤波器的原理，产生低通滤波器。然后让高斯信号过滤波器，最后通过莱斯表达式的方式表示窄带高斯白噪声，求出其对应的自相关函数，利用维纳—辛钦公式，对自相关函数进行傅里叶变换便可以得出功率谱密度函数。
一、问
题
快递公司问题件快递公司问题件货款处理关于圆的周长面积重点题型关于解方程组的题及答案关于南海问题
重述、分析
画出窄带高斯白噪声的功率谱。
要画出窄带高斯白噪声的功率谱，最主要的问题是求出它的自相关函数，然后通过维纳—辛钦公式便可求出。而对于自相关函数的求解主要是得出窄带高斯信号，所以本文基于莱斯表达式，求出它的窄带高斯信号。具体流程如图（1）：
图(1)算法流图
二、问题的解决与算法分析
1、产生高斯白噪声
基于前一次大作业，我们知道产生白噪声的方法有线性同余法。
void gauss(float b[])
{
SYSTEMTIME sys;
GetLocalTime(&sys);
int seed;
float a[100];
seed=sys.wMilliseconds; //随机种子输入,产生范围为0—2^31-1
int i,j,w,k;
float x,y;
for(w=0;w<5000;w++)
for(k=0;k<20;k++) //5000*20个均匀分布的统计近似服从正态分布
{
x=0;y=0;
for(i=0;i<100;i++)
{
seed=16807*(seed%127773); //利用乘同余法产生其他随机数,共100个
if (seed<0)
seed += 2147483647; //处理溢出的情况
a[i]=seed/2147483647.0; //a即为0-1均匀分布数
}
for(j=20;j<32;j++)
{
x=a[j]+x; //取其中12的个数的和，方便计算
}
y=x-6; //产生标准正态分布随机数
for(j=0;j
y)
b[j]=b[j]+1; //统计各个范围内随机数个数
}
}
int gdriver=DETECT,gmode,errorcode;
initgraph(&gdriver,&gmode,"");
initgraph(1000, 600);
int X0=100,Y0=500;
setcolor(RED);
for(i=0;i21&&as<44)//选择汉宁窗
{
n=ceil(8*PI/bt);
if(n%2==0)
n=n+1;
m=n/2;
for(i=0;i44&&as<53)//选择汉明窗
{
n=ceil(8*PI/bt);
if(n%2==0)
n=n+1;
m=n/2;
for(i=0;i53)//选择布莱克曼窗
{
n=ceil(12*PI/bt);
if(n%2==0)
n=n+1;
m=n/2;
for(i=0;ij)
{
a=xr[i];
bj=xi[i];
xi[i]=xi[j];
xr[i]=xr[j];
xr[j]=a;
xi[j]=bj;
}
}
}
void conv(float g[], float h[], float a[])//卷积
{
int i,la;
for(la=0;la=0)&&(la-i
心得体会
决胜全面小康心得体会学党史心得下载党史学习心得下载军训心得免费下载党史学习心得下载
通过本次大作业，对知识点的理解，还有编程能力的提高都有很大帮助。比如对于随机信号产生，高斯信号的产生，还有如何去处理时域与频域的关系，如何将这个问题转换到另一个问题，通过这种联想的思想去解决生活中的问题。但是，从宏观角度来看，我最大的收获是如何用有限的时间创造出高效率还有对问题的转化能力。另外通过本次试验可以感受到课本的教学思路和系统路线，有助于我们更深入探讨与研究相关的知识。
浙公网安备 33021202002483