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