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

N
Nalani · Mar 4, 2021 · 2K views
Question
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
Expert Answer
Profile picture of Kshitij Singh
Kshitij Singh PhD Expert
Answered Aug 29, 2026
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
Z = 1+1i;
abs(Z)^2  
Z*conj(Z)
100% Run Guarantee 3-Hour Fast-Track Delivery

Need a Custom Version or Complete Simulation for This Problem?

Our 500+ PhD engineers build, debug, and optimize working MATLAB scripts and Simulink (.slx) models tailored to your exact assignment rubrics with zero plagiarism.

Tested on MATLAB R2024b / R2026a
Turnitin 0% Plagiarism Report
Free 7-Day Revisions Guarantee
Have a different question? Ask here

Get a Free Consultation or a Sample Assignment Review!