增加四個節點中的節點數量

Question

這可能是一個簡單的問題，但我試圖增加四個集成中的節點數量。我找不到這方面的相關文件。

我有一個一維函數（Dirac delta-like），它大多爲零，但在非常窄的時間間隔內非零。所以自適應四元組無法捕獲該間隔並返回零。我試圖通過更多點抽樣來捕捉峯值。你怎麼做到這一點？ SciPy文檔並沒有多大意義。

這個窄峯的位置發生變化。我不知道高峯的確切位置。

Answer 1

quadpy（我的一個項目）用Gauss-Kronrod在1D中進行自適應正交運算。您可以指定Kronrod規則的程度，並通過增加域中的點數來啓動它。

安裝與

pip install quadpy 

，並嘗試

from numpy import sin 
import quadpy 

val, error_estimate = quadpy.line_segment.adaptive_integrate(
     lambda x: sin(5*x), 
     [0.0, 1.0], 
     1.0e-10, 
     kronrod_degree=10 
     ) 

print(val) 

殺青kronrod_degree域中獲得更多的積分。