您正在使用新版界面(Beta)。返回原版界面

本页面由机器翻译自动翻译而成,翻译内容可能存在不准确之处,敬请谅解。

keisan服务于日常生活和工作的高精度计算网站

龙格-库塔法(4次,2阶常微分方程)

管理书签
2阶常微分方程 y’’=F(x,y,y’)的解 y=f(x)将通过4阶的龙格-库塔法求解。在初始条件 y0=f(x0), y’0=f’(x0) 下,求解 x 在 x0≦x≦xn 的范围内。
\(\normalsize y''=F(x,y,y')\hspace{30px} y_0=f(x_0),\ y'_0=f'(x_0) \rightarrow\ y=f(x)\\\)
F(x,y,p(=y'))
x0
初始条件
y0
=f(x0) 初始条件
y'0=p0
=f'(x0) 初始条件
xn
x0≦x≦xn
[ n]分割数
\(\normalsize Runge-Kutta\ method\\ (1)\ y''=F(x,y,y')\hspace{30px} p=y',\ y_0=f(x_0),\ y'_0=f'(x_0) \rightarrow\ y=f(x)\\ (2)\ p_{n+1}=p_n+{\large\frac{1}{6}}(j_1+2j_2+2j_3+j_4)+{\small O}(h^5)\\\\ \hspace{25px} y_{n+1}=y_n+{\large\frac{1}{6}}(k_1+2k_2+2k_3+k_4)+{\small O}(h^5)\\\\ \hspace{25px} j_1=hF(x_n,\ y_n,\ p_n)\\ \hspace{25px} k_1=h \cdot p_n \\ \hspace{25px} j_2=hF(x_n+{\large\frac{h}{2}},\ y_n+{\large\frac{k_1}{2}},\ p_n+{\large\frac{j_1}{2}})\\ \hspace{25px} k_2=h ( p_n + {\large\frac{j_1}{2}} )\\ \hspace{25px} j_3=hF(x_n+{\large\frac{h}{2}},\ y_n+{\large\frac{k_2}{2}},\ p_n+{\large\frac{j_2}{2}})\\ \hspace{25px} k_3=h ( p_n + {\large\frac{j_2}{2}} )\\ \hspace{25px} j_4=hF(x_n+h,\ y_n+k_3,\ p_n+j_3)\\ \hspace{25px} k_4=h ( p_n + j_3)\\ \)

目前还没有关于此计算的用户评价。