transformada - Calcular el área bajo el gráfico FFT en MATLAB

transformada de fourier matlab pdf

Question

Actualmente hice una FFT de un conjunto de datos que me da una gráfica con la frecuencia en el eje xy la amplitud en el eje y. Me gustaría calcular el área debajo del gráfico para darme la energía.

No estoy seguro de cómo determinar el área porque estoy sin la ecuación y también solo quiero un área determinada de la trama en lugar de toda el área debajo de la trama. ¿Hay alguna manera de hacerlo?

Answer 1

Hay muchas maneras de hacer integración numérica con Matlab. Aquí hay un ejemplo:

%# create some data
x = linspace(0,pi/2,100); %# 100 equally spaced points between 0 and pi/2
y = sin(x);

%# integrate using trapz, which calculates the area in the trapezoid defined by 
%# x(k),x(k+1),y(k),y(k+1) for k=1:length(x)
integral = trapz(x,y);

%# if you only want to integrate part of the data, do
partialIntegral = trapz(x(10:20),y(10:20));

%# show the integrated area
figure, 
area(x,y); 
hold on, 
area(x(10:20),y(10:20),''FaceColor'',''red'')

Answer 2

La FFT es discreta, no continua, solo debe sumar todos los valores de la bandeja. Si observa el espectro de potencia (magnitud al cuadrado), entonces los valores de bin están en W / Hz, por lo que necesitaría multiplicar cada valor (o alternativamente solo la suma), por el ancho del compartimiento en Hz para obtener energía (y de ahí la energía total en su muestra de entrada).