有限区间(a,b)的数值积分
使用双指数型(DE)数值积分法计算有限区间(a,b)的积分。也适用于端点为奇异点(±∞)的函数的积分。
被积函数 f(x) 需满足以下前提条件:除端点外为解析函数,且非周期函数。

| \( \normalsize Double\ exponential\ integration\\ (1)\ x\rightarrow \frac{b-a}{2}y+\frac{b+a}{2}\\ \hspace{30px} {\large\int_{\small a}^{\small b}}f(x)dx= {\large\int_{\small -1}^{\small 1}}f(\frac{b-a}{2}x+\frac{b+a}{2})\frac{b-a}{2}dx\\ (2) x\rightarrow \tanh(\frac{\pi}{2}\sinh(t))\\ \hspace{30px}{\large\int_{\small -1}^{\small 1}}f(x)dx={\large\int_{\small -\infty}^{\small \infty}}f(x(t))x'(t)dt\\ \hspace{140px}x(t)= \tanh(\frac{\pi}{2}\sinh(t))\\ \hspace{140px}x'(t)={\large\frac{\frac{\pi}{2}\cosh(t)}{\cosh^2(\frac{\pi}{2}\sinh(t))}}\\ (3)\ Trapezoid\\ \hspace{20px}S\simeq{\large\int_{-t_a}^{t_a}}f(x(t))x'(t)dt=h{\large\displaystyle \sum_{{\small j=-}\frac{N}{2}}^{\frac{N}{2}}}f(x(jh))x'(jh)\\ \hspace{240px}h={\large\frac{2t_a}{N}}\\ \) |
目前还没有关于此计算的用户评价。