How to take the engineering FFT (magnitude spectrum scaled correctly) of Signal in Matlab?

17 次查看（过去 30 天）

farzad 2020-11-29

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/668058-how-to-take-the-engineering-fft-magnitude-spectrum-scaled-correctly-of-signal-in-matlab

评论： farzad 2020-11-30

Hi All

the fft function in matlab gives a theoretical mathematical fft of the signal, and is not scaled to units and is real immaginary and not amplitude-frequency

I think in better terms, I want to plot the amplitude spectrum of the signal,, with the X axis being the frequencies of the input signal.

is there an alternative way of fft to do this ? and is the following link correct?

for example if this is my signal

using the code in the link I get

which its extremums do not meet the extremums of the time signal

8 个评论
显示 6更早的评论隐藏 6更早的评论

Walter Roberson 2020-11-30

Data would help, along with the code you are using.

farzad 2020-11-30

在 MATLAB Online 中打开

It was not possible to attach the signal file here I dont know why. It is an xlsx file of 2.8 MB. but does not get attached. I will add the code to the question as well. I have tried both of these. you can generate a random signal and check with it. my problem is not the signal

function [frq, amp] = simpleFFT(signal)
            
                n = length(signal); 
                 %generate the vector of frequencies
                
                halfn = floor(n / 2)+1;
                
                deltaf = 1 / ( n / ScanRate);
                
                frq = (0:(halfn-1)) * deltaf;
                z = fft(signal, n); %do the actual work
                x=signal(:,1);
                xFFT = abs(fft(x))/length(x);
                amp = abs(fft(x).^2);
                
end

function [frq, amp] = simpleFFT(signal,ScanRate)
            
                N=length(signal);
                N1=2^nextpow2(N);
                X=fft(signal,N1);
                X=X(1:N1/2);%Discard Half of Points
                X_mag=abs(X); %Magnitude Spectrum
                X_phase=phase(X); %Phase Spectrum
                frq=ScanRate*(0:N1/2-1)/N1; %Frequency axis
                amp= X_mag/N1;
end