N = 128;    %抽頭數
Dlen = 1024;  %輸出長度
Step = 0.005; %步長
%注意：需要給出輸入data和參考信號s(n)才能運行本程序
w=zeros(N,1);           %N為抽頭個數
p=(1/delta)*eye(N);       %delta為正則化參數
for n= N:Dlen
    u(1:N)=data(n:-1:n-N+1);  %濾波器輸入
    v=p*u;
    k=(1/step)*v/(1+(1/step)*u'*v);
    y(n) = w'*u;            %濾波器輸出
    e(n)=s(n)-y(n);
    w=w+k*conj(e(n));
    p=(1/step)*(eye(N)-k*u')*p;
end