Как лучше всего интегрировать произведение функции с функцией косинуса? - PullRequest
Новая
       26

Как лучше всего интегрировать произведение функции с функцией косинуса?

0 голосов
/ 12 ноября 2018

Я хочу интегрировать функцию этой формы $ \int_0^{s_f} V_{pot}(s)\cos(\frac{k\pi}{s_f}s) \mathrm{d}s$, где $V_{pot}$ - массив и имеет эту форму Vpot и k равно $k=0,1,2 ... 2*N_f$. Итак, это мой код 

import numpy as np
from scipy.integrate import simps
Nf = 200 #número de terminos de la integral
s = np.zeros(NP)
for i in range(0,NP):
    s[i] = (i-1)*H
def f(k):
    return [(1/SF)*np.cos(np.pi*k*s/SF)*Vpot]
IVp = np.zeros(2*Nf)
for k in range(0,2*Nf):
       Func = f(k)
       y1 = np.array(Func,dtype=float)
       I = simps(y1,s)
       IVp[k] = I

Затем я получил некоторые данные, но я не знаю, правильно ли это. Как я могу продолжить и какой метод интеграции лучше в этом случае?

...