matlab - 在 MATLAB 中计算 FFT 图下的面积

Question

目前我对一组数据进行了 FFT，得到了一个 x 轴为频率、y 轴为振幅的图。我想计算图表下方的面积来给我能量。

我不确定如何确定面积，因为我没有方程，而且我只想要地 block 的某个区域而不是地 block 下的整个区域。有什么办法可以做到吗？

Answer 1

有很多方法可以用 Matlab 进行数值积分。这是一个例子:

%# 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')