function [e,wop]=rls(m,N,X,d)
I=eye(m);
P=0.1*I;
lamda=0.1;%  遺忘因子
w(1:m,1)=0;
for k=2:N
    L=P*X(:,k)/(lamda+X(:,k)'*P*X(:,k));
    e(k)=d(k)-w(:,k-1)'*X(:,k);
    w(:,k)=w(:,k-1)+L*conj(e(k));
    P=(1-L*X(:,k)')*P/lamda;
end
wop=w(:,N);