数值分析——第一章:常微分方程初值问题数值方法
第一章MATLAB与数值分析之常微分 方程初值问题的数值方法。 本章内容主要参考的是李荣华老师的书 (参见[1]) 中内容第一章,感兴趣的读者可以参考该书。前言
随着计算机不断地更新迭代和计算机技术的日新月异的蓬勃发展,科学计算迎来发展的重要机遇。自计算数学(Computational Mathematics)引入我国以来,计算数学已发展成为数学领域里中的一个重要分支,往往与基础数学、应用数学、概率论与数理统计并列为数学中的“四大高峰”。计算机的计算能力的提升为设计高效数值算法、大规模科学计算、并行计算、云计算等研究课题注入了新的能量。近年来,随着的人工智能的汹涌发展,越来越多的应用数学学者以及活跃在人工智能前沿的研究者意识到了计算数学对该领域的重要性和学习计算数学的必要性。高校计算数学专业的学生所授课程中,最为基础的课程是:数值分析,数值线性代数,以及优化与控制。本文作者将集中于数值分析这门课出发,分享系列Matlab代码,进而夯实数值分析的功底,二来进一步地加深数值分析的理解,同时也为初入计算数学以及对计算数学感兴趣的同仁、学者和学生描绘一幅了解计算数学和科学计算的概念画。
由于作者水平所限,因此难免有这样那样的错误,敬请诸位批评校正,谢谢!
提示:以下是本篇文章正文内容,下面案例可供参考
一、一阶常微分初值问题
设函数 f ( t , x ) f(t,x) f(t,x)在区域 G : 0 ≤ t ≤ T , ∣ u < ∞ G:0\leq t\leq T , | u<\infty G:0≤t≤T,∣u<∞ 上连续,求 u = u ( t ) u=u(t) u=u(t), 使得其满足
d u d t = f ( t , u ) \frac{du}{dt}=f(t,u) dtdu=f(t,u), 0 < t ≤ T , \quad\quad\quad\quad\quad\quad\quad\quad 0<t\leq T,~ ~ 0<t≤T, ( 1. 1 a 1.1_a 1.1a)
u ( 0 ) = u 0 u(0)=u_0 u(0)=u0, \quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad~ ( 1. 1 b 1.1_b 1.1b)
其中 u 0 u_0 u0是给定的值,一般称之为初值。系统 (1.1) 称为一阶常微分方程的初值问为了能使得问题 (1.1) 的解存在、唯一且连续依赖于初值,也就是说问题 (1.1) 适定,则需要对右端 f ( t , u ) f(t,u) f(t,u)加适当的要求,通常要求 f f f关于 u u u满足Lipschits条件,即,存在常数 L L L, 使得
∣ f ( t , u 1 ) − f ( t , u 2 ) ∣ ≤ L ∣ u 1 − U 2 ∣ |f(t,u_1)-f(t,u_2)|\leq L|u_1-U_2| ∣f(t,u1)−f(t,u2)∣≤L∣u