af.c

自适应滤波器设计的仿真程序

void main()
{
	int n,k;
	float x[64],y1[64],y2[64],h[64],d[64],e[64];
	
	for(n=0;n<64;n++)
	{
		x[n] = n;
		h[n] = 0.1;
		d[n] = 202.0/64*n;
		y1[n] = 0;
		y2[n] = 0;
		e[n] = 0;
	}
	
	for(n=0;n<64;n++)
	{
		for(k=0;k<=n;k++)
		{
			y1[n] += h[k]*x[n-k];
		}
	}
	
	e[0]=0;
	for(n=0;n<64;n++)
	{
		for(k=0;k<=n;k++)
		{
			y2[n] += h[k]*x[n-k];
		}
		e[n] = d[n] - y2[n];
		for(k=0;k<=n;k++)
		{
			h[k] += 0.005*e[n]*x[n-k];
		}
	}
	
	e[0]=0;
}