擅长:python、mysql、java
<p>如果使用简单的和(或梯形)积分:</p>
<pre><code>ft = np.fft.fft(y)
integral = sum(ft[105:121])
</code></pre>
<p>或者</p>
^{pr2}$
<p>看来应该行得通。在</p>
<pre><code>>>> import numpy as np
>>> x = x=np.arange(0,15,0.01)
>>> from numpy import exp,pi
>>> y=exp(-0.3*x)*exp(x*pi*20j)+exp(-0.9*x)*exp(x*pi*25j)+exp(-0.9*x)*exp(x*pi*15j)
>>> ft = np.fft.fft(y)
>>> np.trapz(ft[105:121])
(642.14009362811771+142.9776425340925j)
>>> sum(ft[105:121])
(652.29308789751224+152.70583448308713j)
</code></pre>