龙格-库塔法(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)\\\) |
| \(\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)\\ \) |
目前还没有关于此计算的用户评价。