python中的2变量数值积分积分

Question

在下面的测试中，我想在每个频率步骤的a和b之间进行积分。我需要提取数组中的数据，这样我就可以显示每个freq步骤的积分结果。最终，我想为一个复杂的可积函数重现这一点，在这里我将绘制真实和想象的部分。由于我没有相同的变量维度，我得到了一个广播错误。我不知道如何用numpy写这篇文章。

f = np.sin(x)*freq  # for each freq calculate the integrale and store the result 
ValueError: operands could not be broadcast together with shapes (11,) (5,)

import numpy as np
a = 0
b = 10
n = 11
h = (b - a) / (n - 1)
x = np.linspace(a, b, n)
freq = np.linspace(0.001, 100, 5)
f = np.exp(x*freq)  # for each freq calculate the integrale and store the result of whatever function of 2 variables
#integration
for f in freq:
I_simp = ((h/3) * (f[x[0],freq[0]] + 2*sum(f[x[:n-2:2],freq[f]]) 
+ 4*sum(f[x[1:n-1:2],freq[f]]) + f[x[n-1],freq[-1]]))
print(I_simp)    #print the array , in case of complex i will then extract real and imag

Answer 1

我建议您在集成之前绘制函数。

Wolfram的这个图，如果正确的话，表明函数是奇异的。此集成会出现问题。