realifft.m

来自「TSTOOL应用软件,内有说明文件,可以适当修改,应用比较方便.」· M 代码 · 共 15 行

function rs = realifft(data, N)

% inverse fft for fourier coefficents from real valued  original data
% needs the length of the original time series from which the fft was computed
% see also : realfft

data = data(:);
if mod(N,2) 	% odd length
	data = [data ; conj(data(end:-1:2))];  
else			% even length
	data = [data ; conj(data(end-1:-1:2))];  
end 

rs = real(ifft(data)); % remove small imaginary components introduced by rounding errors