LMS程序算法

clear;%清除工作空间
close all%...
hold off%...
sysorder=5;%抽头数...
samples_per_period=50; %周期采样频率
N=10*samples_per_period-1;%总采样次数
signal=cos（pi*0.02*[0:N-1]）+sin（2*pi*0.02*[0:N-1]）+cos（3*pi*0.02*[0:N-1]）+sin（4*pi*0.02*[0:N-1]）+cos（5*pi*0.02*[0:N-1]）...
    +sin（6*pi*0.02*[0:N-1]）+cos（7*pi*0.02*[0:N-1]）+sin（8*pi*0.02*[0:N-1]）+cos（9*pi*0.02*[0:N-1]）;%初始输入信号/期望信号
figure;%作图1 的第一子图
subplot（211）;
plot（signal）;
grid;
title（‘自适应滤波器的理想输入‘）;
nvar=0.5;%噪声方差
noise=nvar*randn（1N ）;%噪声信号
X=signal+noise;%带加性噪声的输入信号
delayX=[0 X];% 输入信号延时
subplot（212）;%作图1的第二子图
plot（delayX）;
grid;
title（‘自适应滤波器的噪声输入‘）;
SIGNAL=[signal 0];
%M=32;%滤波器长度
M=length（signal）;%M为接收数据长度
mu=0.002;%迭代步长
t