keisan生活や実務に役立つ高精度計算サイト

Tanh-Sinh数値積分の分点と重み

Tanh-Sinh数値積分の分点(nodes)と重み(weights)を計算します。

\(\normalsize {\large\int_{\small -1}^{\small 1}}f(x)dx\simeq{\large\int_{-t_a}^{t_a}}f(x(t))x'(t)dt\simeq{\large\displaystyle \sum_{\small i=1}^{n}}w_{i}f(x_i)\\ \ nodes\\ \hspace{20px} x_i=\tanh(\frac{\pi}{2}\sinh(t_i)),\hspace{10px}t_i=-t_a+(i-1)h\\ \ weights\\ \hspace{20px} w_i={\large\frac{\frac{\pi}{2}cosh(t_i)}{cosh^2(\frac{\pi}{2}\sinh(t_i))}}h,\hspace{10px}h={\large\frac{2t_a}{n-1}}\\ \)
nodes
order n
n=2,3,4,..,100
ta
\(\normalsize Tanh-Sinh\ integration\\\hspace{5px} (1)\hspace{1px} {\large\int_{\small -1}^{\small 1}}f(x)dx={\large\int_{\small -\infty}^{\small \infty}}f(x(t))x'(t)dt\\ \hspace{90px}\simeq{\large\int_{-t_a}^{t_a}}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))}}\\ (2)\ Trapezoid\\ \hspace{20px} {\large\int_{\small -1}^{\small 1}}f(x)dx\simeq{\large\displaystyle \sum_{\small i=1}^{n}}w_{i}f(x_i)\\ \hspace{10px}nodes\\ \hspace{20px} x_i=\tanh(\frac{\pi}{2}\sinh(t_i)),\hspace{10px}t_i=-t_a+(i-1)h\\ \hspace{10px}weights\\ \hspace{20px} w_i={\large\frac{\frac{\pi}{2}\cosh(t_i)}{\cosh^2(\frac{\pi}{2}\sinh(t_i))}}h,\hspace{10px}h={\large\frac{2t_a}{n-1}}\\ \)

この計算についてお客様の声はまだありません。