#include <stdio.h>
#define LEN 2000//序列长度
#define M 6//滤波器阶数
#define var_v 1
void uniform(double a,double b,long int *seed,double w[],double variance);
int main()
{
double x[LEN],h[M]={0},y[LEN]={0},e,d[LEN],w[LEN];
double u=0.00001;
int i,j,k,m;
long int seed=9;
uniform(0,1,&seed,w,var_v); //产生噪声v
for(i=0;i<LEN;i++)
{
d[i]=5;
x[i]=d[i]+w[i];
}
for(i=0;i<LEN;i++)
{
if(i<M)
{
for(j=0;j<i+1;j++)
y[i]+=h[j]*x[i-j];
e=d[i]-y[i];
for(k=0;k<i+1;k++)
h[k]+=2*u*e*x[i-k];
//////////////////////
for(m=0;m<i;m++)
printf("%f\t",h[m]);
printf("\n");
}
else
{
for(j=0;j<M;j++)
y[i]+=h[j]*x[i-j];
e=d[i]-y[i];
for(k=0;k<M;k++)
h[k]+=2*u*e*x[i-k];
//////////////////////
for(m=0;m<M;m++)
printf("%f\t",h[m]);
printf("\n");
}
}
//////////////////////
return 0;
}
LMS.rar_lms 收敛_matlab 均方误差_最小均方误差_计算误差
版权申诉
45 浏览量
2022-09-24
15:39:31
上传
评论
收藏 182KB RAR 举报
LMS.rar （17个子文件） 
uniform.cpp 562B
lms lms.ncb 33KB lms.dsw 529B lms.plg 240B uniform.cpp 563B main.cpp 871B lms.opt 53KB Debug vc60.idb 33KB main.obj 4KB uniform.obj 3KB lms.pch 209KB lms.ilk 207KB vc60.pdb 44KB lms.pdb 441KB
lms.exe 196KB lms.dsp 4KB main.cpp 846B
评论0