rec_frommag_complex.m

来自「Second and Higher-Order Statistics based」· M 代码 · 共 37 行

function g=rec_frommag_complex(M,I,N);
%------------------------------------------------------------------------------
%	rec_frommag_complex.m
%	This function generate the fourth order cross cumulant of 4 signals
%
% Usage:
%	 [g] = rec_frommag_complex(M,I,N);
% Where
%	g         : the reconstructed time-domain signal.
%	M         : the magnitude response.
%	I         : the index of the frequencies to be used in the reconstruction.
%	N         : the length of the time-domain signal to be reconstructed.
% 
%  Communications and Signal Processing Laboratory
%  ECE Department, Drexel University
%  Philadelphia, PA 19104, USA
%  http://www.ece.drexel.edu/CSPL
%-------------------------------------------------------------------------------

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);