自适应滤波的MATLAB实现.doc

自适应滤波的MATLAB实现

发表于:2006-12-0608:56来源:中国振动联盟

1.LMS算法的仿真程序：

%lms算法

clearall

closeall

holdoff%系统信道权数

sysorder=5;%抽头数

N=1000;%总采样次数

inp=randn(N,1);%产生高斯随机系列

n=randn(N,1);

[b,a]=butter(2,0.25);

Gz=tf(b,a,-1);%逆变换函数

h=[0.0976;0.2873;0.3360;0.2210;0.0964;];%信道特性向量

y=lsim(Gz,inp);%加入噪声

n=n__std(y)/(10__std(n));%噪声信号

d=y+n;%期望输出信号

totallength=size(d,1);%步长

N=60;%60节点作为训练序列

%算法的开始

w=zeros(sysorder,1);%初始化

forn=sysorder:N

u=inp(n:-1:n-sysorder+1);%u的矩阵

y(n)=w__u;%系统输出

e(n)=d(n)-y(n);%误差

ifn20

mu=0.32;

else

mu=0.15;

end

w=w+mu__u__e(n);%