Question

目前我对一组数据进行了FFT，这给出了一个频率为x轴，幅度为y轴的图。我想计算图下的面积给我能量。

我不确定如何确定该区域因为我没有方程式而且我只想要一个区域的情节而不是整个区域的情节。我有办法吗？

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

Answer 2

FFT是离散的，而不是连续的 - 您只需要对所有bin值求和。如果您正在查看功率谱（幅度平方），则bin值以W / Hz为单位，因此您需要将每个值（或者仅仅是总和）乘以以Hz为单位的bin宽度来获得功率（和因此输入样本中的总能量。）

在MATLAB中计算FFT图下的面积

2 个答案: