自适应信号处理 LMS算法

%lms 算法
clear all
close all
hold off%系统信道权数
sysorder = 5 ;%抽头数
N=1000;%总采样次数
inp = randn（N1）;%产生高斯随机系列
n = randn（N1）;
[ba] = butter（20.25）;
Gz = tf（ba-1）;%逆变换函数
h= [0.0976;0.2873;0.3360;0.2210;0.0964;];%信道特性向量
y = lsim（Gzinp）;%加入噪声
n = n * std（y）/（10*std（n））;%噪声信号
d = y + n;%期望输出信号
totallength=size（d1）;%步长
N=60 ; %60节点作为训练序列
%算法的开始
w = zeros （ sysorder  1 ） ;%初始化
for n = sysorder : N 
u = inp（n:-1:n-sysorder+1） ;% u的矩阵
y（n）= w‘ * u;%系统输出
e（n） = d（n） -