pmtm (Multitaper) calculates PSD from Fourier coefficients as abs(Xx)^2; shouldn't that be Xx.conj(Xx)?

Question

Ajoy 2012-12-13

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/56531-pmtm-multitaper-calculates-psd-from-fourier-coefficients-as-abs-xx-2-shouldn-t-that-be-xx-conj-x

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

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Wayne King 2012-12-13

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/56531-pmtm-multitaper-calculates-psd-from-fourier-coefficients-as-abs-xx-2-shouldn-t-that-be-xx-conj-x#answer_68451

编辑：Wayne King 2012-12-13

在 MATLAB Online 中打开

There is no difference between

abs(Z).^2

and

Z.*conj(Z)

for a complex vector. They are equivalent. In fact, for a complex number the modulus (norm) is defined to be the square root of the product of Z and its conjugate.

Just convince yourself for a single complex number