数学建模社区-数学中国
标题:
常微分方程的解法 (三): 龙格—库塔(Runge—Kutta)方法 、线性多步法
[打印本页]
作者:
浅夏110
时间:
2020-6-9 14:53
标题:
常微分方程的解法 (三): 龙格—库塔(Runge—Kutta)方法 、线性多步法
§ 4 龙格—库塔(Runge—Kutta)方法
z# _: F/ i v: G0 F/ f1 ?
1 龙格—库塔方法的基本思想
* K3 A6 H. t, P. h9 E0 j2 x
: i- Z/ A9 ]: Q3 ?9 r8 V1 z
4 T8 R7 X8 f( j8 F4 Z) ~
! ^) S" W5 v- Q( L
确定系数以提高精度
, S$ ]) x: W& o4 U+ D
3 _" M3 s+ ?* P0 o3 [' O
" A% Z H- e0 i, j- M+ [
7 l0 F! G6 w9 v+ R; e) E& k( Z
; y5 ?/ s2 o0 R! r2 P7 v
8 x9 I J0 Y9 n
3 }' a$ l& i: F: R
2 . RK方法 : 4 阶龙格—库塔公式
d, G) x3 U0 U7 I
要进一步提高精度,必须取更多的点,如取 4 点构造如下形式的公式:
5 O7 i$ J5 q% Z5 d
7 o3 q1 ]) n, `3 p/ ^6 B, ^$ K
2 i: U! j0 v* D8 w
8 x4 M( S" r& U2 m1 v
8 u5 z/ T' s7 G2 T# Q t( k
& E7 r1 N7 n% e( [9 r# h9 `! E( y6 E. [7 Z
这就是常用的 4 阶龙格—库塔方法(简称 RK 方法).
6 l! _! G3 H8 x, k
' y3 R0 T+ ]+ L9 b9 |. \
" L! U( Y1 R5 p. p
B6 @; q) d; F9 h4 `
5 线性多步法
6 k8 K5 f! _ ?
多步法的基本思想 、增量函数
4 ?+ Y7 X% n9 X7 Z
9 y/ @: o* z) ?+ J* X% m/ n" s: d
F: x3 l) G! E) M* r
% c: t* a/ L6 `) \8 w5 K- t
' j$ B, W; E( {$ V2 j5 G
0 h* U( ]1 J" U+ c1 S t7 k& U
. a2 Y- y( D5 b8 x
2 Z }* K9 y2 |( p6 f. m$ G# J
§6 一阶微分方程组与高阶微分方程的数值解法
. r* A3 p3 }% \% M7 ]" ~
6.1 一阶微分方程组的数值解法
5 c8 S0 M/ N" R/ S; ]: F/ i% K
L+ P8 Q) [: u
4 K' e3 Z) o. S
1 E) V9 x# k' C' _% @
么问题(25)在[a,b] 上存在唯一解 y = y(x) 。 问题(25)与(1)形式上完全相同,故对初值问题(1)所建立的各种数值解法可 全部用于求解问题(25)。
2 @& q$ U# \$ i
/ _% Y! E% @: v8 W
6.2 高阶微分方程的数值解法
$ B- n# y$ D N9 u# p! ]" {" e. w
高阶微分方程的初值问题可以通过变量代换化为一阶微分方程组初值问题。
$ J: [' ^/ B' V- a
& w, t9 m* x5 b# B' Y
' A! a7 Q- F* A! G' U6 |
7 b+ X& {( A: |/ W: o+ y3 u+ i
! x; Z+ J s, f! C! w6 G1 c
% a" v l' f% N2 r' R3 S1 J$ \
刚性方程组、Stiff 方程组
! i. T' h0 G3 g+ b1 Z
最后需要指出的是,在化学工程及自动控制等领域中,所涉及的常微分方程组初值 问题常常是所谓的“刚性”问题。具体地说,对一阶线性微分方程组
/ |% ]' C' ?5 c+ }5 }3 ~4 H
. s$ D$ g, S* C e# ^, y3 u
5 h* t& K& P0 b) E
————————————————
) A: e' S! i' D9 ` A% s; y
版权声明:本文为CSDN博主「wamg潇潇」的原创文章,遵循CC 4.0 BY-SA版权协议,转载请附上原文出处链接及本声明。
' t" ]9 @. w, A# E# A1 f: L
原文链接:https://blog.csdn.net/qq_29831163/article/details/89703598
1 P4 p# @# B6 F5 D- `9 z0 { H
8 u: Y) F& k Y+ t
4 @/ d- K8 L: U+ V) Y
欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)
Powered by Discuz! X2.5