In MATLAB 7.11.2(R2010b) Service Pack 2, in the Signal Processing Toolbox file pmtm.m (Multitaper), line 232, the current code uses abs(Xx).^2, where XX is the set of raw Fourier coefficients, to eventually calculate the PSD. Shouldn't that be Xx.conj(Xx) since power is defined as the square of the norm of the Fourier coefficients?
% Compute DFT using FFT or Goertzel
[Xx,w] = computeDFT(E(:,1:k).*x(:,ones(1,k)),nfft,Fs); %Sk = abs(Xx).^2; % OLD CODE HERE
Sk = Xx.*conj(Xx); % AK CHANGE HERE
Thanks, A K
Best Answer