untitled7.asv

来自「基于MMSE的波束形成器」· ASV 代码 · 共 23 行

I=eye(m);
P=1*I;
lamda=0.05;%  遗忘因子
w(1:m,1)=0;
for k=2:N
    L=P*S(:,k)/(lamda+(S(:,k))'*P*S(:,k));
    d=s(1,k);% 期待响应信号
    e(k)=d-(w(:,k-1))'*S(:,k);
    w(:,k)=w(:,k-1)+L*conj(e(k));
    P=(1-L*(S(:,k))')*P/lamda;
end



%%%LMS

   u=0.001;
   w(1:m,1)=0;
  for k=1:N
     d=s(1,k);
     e(k)=d-w(:,k)'*S(:,k);
     w(:,k+1)=w(:,k)+2*u*conj(e(k))*S(:,k);
  end