untitled7.m

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

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