rec_frommag_complex.m

来自「用MATLAB实现MIMO系统盲辨识」· M 代码 · 共 20 行

function g=rec_frommag_complex(M,I,N);

NF=length(M);

theta=zeros(size(M));
tmp=M;

for i=1:600
   	G=tmp.*exp(j*theta);
      g=(ifft(G));
      g(N+1:NF)=zeros(1,NF-N);	
      g=real(g);	
      G=fft(g,NF);
      theta=phase(G);
      
      tmp=abs(G);
		tmp(I)=M(I);
		%tmp=[tmp(1:NF/2+1) fliplr(tmp(2:NF/2))];
end
g=g(1:N);