[xFs]=audioread（‘rf.MP3‘）;
x=x（1:220501）;
x1=fft（x）;
p2=abs（x1/22050）;
p1=p2（1:11025+1）;
f=Fs*（0:11025）/22050;
%b=[0.10.20.40.20.1];
%d=filter（b1x）;
k=x.*（-1）
mu=0.05;
h=adaptfilt.lms（5mu）;
[ye]=filter（hxk）;
m=x+y;
figure（1）;
subplot（311）;
plot（1:22050x‘-b‘）;
xlabel（‘時間指數‘）;ylabel（‘幅值‘）;
legend（‘采樣信號‘）;
subplot（312）;
plot（1:22050k‘--r‘）;
xlabel（‘時間指數‘）;ylabel（‘幅值‘）;
legend（‘降噪信號‘）;
subplot（313）;
plot（1:22050y‘-.c‘）;
xlabel（‘時間指數‘）;ylabel（‘幅值‘）;
legend（‘迭代信號‘）;
figure（2）;
subplo