rls_main.c

对高斯白噪声进行rls滤波

#include "stdio.h"
#include "stdio.h"
#include "math.h"

void rls(x,d,n,w,m,r);
double gauss(double mean,double sigma,int *s);

main()
{int i,m,n;
long seed;
double r,w[10],d[100],x[100];
void rls();
seed=13579l;
n=100;
for(i=0;i<n;i++)
{x[i]=gauss(0.0,1.0,&seed);}
d[0]=2.0*x[0];
d[1]=2.0*x[1]-0.5*x[0];
d[2]=2.0*x[2]-0.5*x[1]+1.4*x[0];
for (i=3;i<n;i++)
{d[i]=2.0*x[i]-0.5*x[i-1]+1.4*x[i-2]+0.1*x[i-3];}
m=4;
for(i=0;i<m;i++)
{w[i]=0.0;};
r=1.0;
rls(x,d,n,w,m,r);
printf("\n The Coefficient of FIR Digita System\n");
for (i=0;i<m;i++)
{printf("   w(%d)=%10.7llf\n",i,w[i]);}
}