QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 2871|回复: 0
打印 上一主题 下一主题

[建模教程] 常微分方程的解法 (三): 龙格—库塔(Runge—Kutta)方法 、线性多步法

[复制链接]
字体大小: 正常 放大
浅夏110 实名认证       

542

主题

15

听众

1万

积分

  • TA的每日心情
    开心
    2020-11-14 17:15
  • 签到天数: 74 天

    [LV.6]常住居民II

    邮箱绑定达人

    群组2019美赛冲刺课程

    群组站长地区赛培训

    群组2019考研数学 桃子老师

    群组2018教师培训(呼伦贝

    群组2019考研数学 站长系列

    跳转到指定楼层
    1#
    发表于 2020-6-9 14:53 |只看该作者 |倒序浏览
    |招呼Ta 关注Ta |邮箱已经成功绑定
    § 4 龙格—库塔(Runge—Kutta)方法
    # M5 M6 ~0 n7 f1 龙格—库塔方法的基本思想
    4 k( j' t* `  r) ^) a9 K& ^% `+ y. ?8 p" X) |) o* N. n6 T

    4 Z% Q* Y0 F% Y$ M% b, m
    1 w# f( G: u. e确定系数以提高精度
      ?8 i4 k& W0 i$ _1 j0 s. l
    ( U* Z8 r3 t- R4 L/ Y& u9 u, A7 f4 D, \  v2 m1 W  R/ I3 T* |: C- k

    0 r& y* i! d% r  j! U2 H
    8 J2 U6 j0 `. u
    ' Z& b( g5 p& W# Q$ V! _% e! v' g6 ]8 W$ \5 x
    2 .  RK方法 :  4 阶龙格—库塔公式( N0 u( A3 p5 n) l5 d
    要进一步提高精度,必须取更多的点,如取 4 点构造如下形式的公式:
    ' s/ n  l( l0 M  S( q% @, ~3 C
    7 ~5 G! N2 P6 B! ]1 U4 [  W: g* e) _" D( o0 x, N

    8 h4 H* |+ J. M0 h# e: @% N& p+ u1 m* M4 m7 Z" V& r
    2 r. u7 m) T' A: i. k- m4 L$ h
    这就是常用的 4 阶龙格—库塔方法(简称 RK 方法).
    + h( S  _& o% V9 U7 ]& z- n. m1 }. I$ H
    % z, w1 r6 {) _7 J. u

      N8 V9 r  e5 r4 G) F5  线性多步法7 a* H- ?1 _) n* @% c; @& G
    多步法的基本思想 、增量函数: b0 f. i# Y+ r2 Y0 y" q6 Z7 l

    ( N( v- b! N( F, A; F8 b+ d, x$ ]7 E; N. b/ N" o- E9 E

    ) N+ X5 u# a4 H/ |& v: e+ Z0 m- p+ ~

    : Q# N& e3 v9 ^, q: g4 Q8 V
    + D% d1 ~$ g( ?  \! _1 b4 E4 q& A2 D' J6 ~3 i
    §6 一阶微分方程组与高阶微分方程的数值解法2 q! w! H$ m7 p) h% w5 W; u
    6.1 一阶微分方程组的数值解法
    - b! ]3 q) n" J$ o  J3 [2 t9 ]& z$ E# I4 F2 G$ q( k8 W
    7 A$ l. g% [1 v9 `6 D( Z

    ) p0 q4 ^! M" A- |; r么问题(25)在[a,b] 上存在唯一解 y = y(x) 。 问题(25)与(1)形式上完全相同,故对初值问题(1)所建立的各种数值解法可 全部用于求解问题(25)。" ], N5 H  i$ z) P' V

    - F' s, q. {/ \0 {: c6.2 高阶微分方程的数值解法
    : G; p$ s/ u% c9 U高阶微分方程的初值问题可以通过变量代换化为一阶微分方程组初值问题。
    3 |7 v, ?8 V: ~8 V
    6 D. \/ b6 g% W" u$ C' u. v" f! s  {' Y* W9 o7 D2 s- q
    ; T) B9 y6 V2 W/ @
    - [- A$ {0 X& H3 E. S! P/ h, j

    % O. q, I* u7 l7 O1 x刚性方程组、Stiff 方程组. R7 G/ f9 D+ A8 I0 E6 P
    最后需要指出的是,在化学工程及自动控制等领域中,所涉及的常微分方程组初值 问题常常是所谓的“刚性”问题。具体地说,对一阶线性微分方程组
    3 D' Y% x3 [# r9 E, a; \" Z& p) I* W8 j: b5 R0 l8 b6 g; b" Z6 v' s

    3 T8 P4 K% w: t9 r- r, [( X( r————————————————+ ^& Z4 G8 Y3 E9 X# |3 |1 U7 f
    版权声明:本文为CSDN博主「wamg潇潇」的原创文章,遵循CC 4.0 BY-SA版权协议,转载请附上原文出处链接及本声明。! K) r2 H& z& Z6 U! x
    原文链接:https://blog.csdn.net/qq_29831163/article/details/897035984 ~! N" y# M# Z% L
    & d7 |# x- ]' v+ z8 s
    ! g0 I1 P* W; ~0 E% R
    zan
    转播转播0 分享淘帖0 分享分享0 收藏收藏0 支持支持0 反对反对0 微信微信
    您需要登录后才可以回帖 登录 | 注册地址

    qq
    收缩
    • 电话咨询

    • 04714969085
    fastpost

    关于我们| 联系我们| 诚征英才| 对外合作| 产品服务| QQ

    手机版|Archiver| |繁體中文 手机客户端  

    蒙公网安备 15010502000194号

    Powered by Discuz! X2.5   © 2001-2013 数学建模网-数学中国 ( 蒙ICP备14002410号-3 蒙BBS备-0002号 )     论坛法律顾问:王兆丰

    GMT+8, 2026-6-11 04:07 , Processed in 0.361003 second(s), 50 queries .

    回顶部