【文件属性】:
文件名称:rls滤波器,用于去除噪声
文件大小:7KB
文件格式:C
更新时间:2013-06-09 14:07:00
rls
用于滤波,rls自适应滤波
#include
#include
#include
#define N 110000
float x[N];
float x1[110050];
float xout[N];
float w[50];
float p[50][50];
float u[50];
float kn[50];
float kn1[50];
float dn;
float e;
float outmiddle;
//阶数50//
//****************************************初始化程序用于计算w[][].p[]**********************************//
float initial ()
{
char i1,j1;
for(i1=0;i1<50;i1++)
{
w[i1]=0;
}
for(i1=0;i1<50;i1++)
{
for(j1=0;j1<50;j1++)
{
p[i1][j1]=0;
}
}
for(i1=0;i1<50;i1++)
{
p[i1][i1]=9999999.9;
}
return(0);
}
//****************************************************************************************************************//
//****************************************e[n]的求解**************************************************************//
float en(unsigned long int datanumber)
{
/*float result=0;
int i2,k2;
int j2;
j2=datanumber;
k2=datanumber+5;
e=0;
for(i2=0;i2<50;i2++)//计算w()*u()
{
e=e+w[i2]*u[i2];
}
dn=0;
for(i2=j2;i2