clc;
clear?all;
close?all;
[XFsbits]?=?wavread（‘C:\Users\xiaodan\Desktop\如果這就是愛情?00_00_00-00_00_06.wav‘）;
X1?=?X（:1）;????????%取出雙通道中其中一個通道
X1=X1‘;
%sound（X14410016）;?%改變一下播放的頻率或者比特數，會有不同的聲音效果

i?=?length（X1）;
t=（0:i-1）/Fs;
n?=?（1:i）‘;
v?=?0.8*randn（i1）;
v=v‘;

ar?=?[11/2];
v1?=?filter（1arv）;
x?=?X1?+?v1;
%sound（xFsbits）;
ma?=?[1-0.80.4-0.2];
v2?=?filter（ma1v）;
%sound（v2Fsbits）;

%?初始化RLS算法
Worder=32;???????????????????????????%濾波器階數
lambda=1?;???????????????????????????%設置遺忘因子
Delta=0.001?;????????????????????????%輸入協方差估計，當SNR高時取較小正常數
p=（1/Delta）?*?eye?（?WorderWorder?）;?%時間相關矩陣的逆矩陣初值???
w=zeros（Worder1）;???????????????????%權向量，列向量
[a1N]?=?size（x）;????????????????????%將矩陣的行數返回到第一個輸出變量a1，將矩陣的列數返回到第二個輸出變量N
frefpad?=?[zeros（1Worder-1）?v2];????%輸入（行向量），補零可以使輸出從第0個開始

%RLS算法公式
for?i?=?