第八章 常微分方程初步

常微分方程和偏微分方程的数值解法教学大纲

上海交通大学致远学院 《常微分方程和偏微分方程的数值解法》教学大纲 一、课程基本信息 课程名称（中文）：常微分方程和偏微分方程的数值解法 课程名称（英文）：Numerical Methods for Ordinary and Partial Differential Equations 课程代码：MA300 学分 / 学时：4学分 / 68学时 适用专业：致远学院与数学系相关专业 先修课程：偏微分方程，数值分析 后续课程：相关课程 开课单位：理学院数学系计算与运筹教研室 Office hours: 每周二19：00—21：00，地点：数学楼1204 二、课程性质和任务 本课程是致远学院和数学系应用数学和计算数学方向的一门重要专业基础课程，其主要任务是通过数学建模、算法设计、理论分析和上机实算“四位一体”的教学方法，使学生掌握常微分方程与偏微分方程数值解的基本方法、基本原理和基本理论，进一步提升同学们利用计算机解决实际问题的能力。在常微分方程部分，将着重介绍常微分方程初值问题的单步法，含各类Euler方法和Runge-Kutta方法，以及线性多步法。将简介常微分方程组和高阶常微分方程的数值方法。在偏微分方程部分，将系统介绍求解椭圆、双曲、抛物型方程的差分方法的构造方法和理论分析技巧，对于椭圆型方程的边值问题将介绍相应变分原理与有限元方法。将在课堂上实时演示讲授的核心算法的计算效果，以强调其直观效果与应用性。本课程重视实践环节建设，学生要做一定数量的大作业。 三、教学内容和基本要求 第一部分：常微分方程数值解法 1 引论 1.1回顾：一阶常微分方程初值问题及解的存在唯一性定理

考研数学三(常微分方程与差分方程)-试卷4

考研数学三（常微分方程与差分方程）-试卷4 (总分：58.00，做题时间：90分钟) 一、选择题(总题数：3，分数：6.00) 1.选择题下列每题给出的四个选项中，只有一个选项符合题目要求。（分数： 2.00） __________________________________________________________________________________________ 解析： 2.设函数y 1 (x)，y 2 (x)，y 3 (x)线性无关，而且都是非齐次线性方程(6．2)的解，C 1，C 2为任意常数，则该非齐次方程的通解是 （分数：2.00） A.C 1 y 1 +C 2 y 2 +y 3． B.C 1 y 1 +C 2 y 2 -(C 1 +C 2 )y 3． C.C 1 y 1 +C 2 y 2 -(1-C 1 -C 2 )y 3． D.C 1 y 1 +C 2 y 2 +(1-C 1 -C 2 )y 3．√ 解析：解析：对于选项(D)来说，其表达式可改写为 y 3 +C 1 (y 1 -y 3 )+C 2 (y 2 -y 3 )，而且y 3是非齐次方程(6．2)的一个特解，y 1 -y 3与y 2 -y 3是(6．4)的两个线性无关的解，由通解的结构可知它就是(6．2)的通解．故应选(D)． 3.已知sin 2 x，cos 2 x是方程y""+P(x)y"+Q(x)y=0的解，C 1，C 2为任意常数，则该方程的通解不是（分数：2.00） A.C 1 sin 2 x+C 2 cos 2 x． B.C 1 +C 2 cos2x． C.C 1 sin 2 2x+C 2 tan 2 x．√ D.C 1 +C 2 cos 2 x． 解析：解析：容易验证sin 2 x与cos 2 x是线性无关的两个函数，从而依题设sin 2 x，cos 2 x为该方程的两个线性无关的解，故C 1 sin 2 x+C 2 cos 2 x为方程的通解．而(B)，(D)中的解析式均可由C 1 sin 2 x+C 2 cos 2 x恒等变换得到，因此，由排除法，仅C 1 sin 2 2x+C 2 tan 2 x不能构成该方程的通解．事实上，sin 2 2x，tan 2 x都未必是方程的解，故选(C)． 二、填空题(总题数：1，分数：2.00) 4.当y＞0时的通解是y= 1． （分数：2.00） 填空项1:__________________ （正确答案：正确答案：[*]） 解析：解析：将原方程改写成，然后令y=ux，则y"=u+xu"．代入后将会发现该变形计算量较大．于 是可转换思维方式，将原方程改写成分离变量，然后积分得 三、解答题(总题数：25，分数：50.00) 5.解答题解答应写出文字说明、证明过程或演算步骤。（分数：2.00） __________________________________________________________________________________________ 解析： 6.求微分方程x(y 2 -1)dx+y(x 2 -1)dy=0的通解． （分数：2.00） __________________________________________________________________________________________ 正确答案：(正确答案：用(x 2 -1)(y 2 -1)除方程的两端，则原方程化为由此可见这是一个变量可

常微分方程边值问题的数值解法

第8章 常微分方程边值问题的数值解法 引 言 第7章介绍了求解常微分方程初值问题的常用的数值方法；本章将介绍常微分方程的边值问题的数值方法。 只含边界条件(boundary-value condition)作为定解条件的常微分方程求解问题称为常微分方程的边值问题(boundary-value problem). 为简明起见，我们以二阶边值问题为 则边值问题(8.1.1)有唯一解。 推论 若线性边值问题 ()()()()()(),, (),()y x p x y x q x y x f x a x b y a y b αβ'''=++≤≤?? ==? (8.1.2) 满足 (1) (),()p x q x 和()f x 在[,]a b 上连续； (2) 在[,]a b 上, ()0q x >, 则边值问题(8.1.1)有唯一解。 求边值问题的近似解，有三类基本方法： (1) 差分法(difference method)，也就是用差商代替微分方程及边界条件中的导数，最终化为代数方程求解； (2) 有限元法(finite element method)；

(3) 把边值问题转化为初值问题，然后用求初值问题的方法求解。 差分法 8.2.1 一类特殊类型二阶线性常微分方程的边值问题的差分法 设二阶线性常微分方程的边值问题为 (8.2.1)(8.2.2) ()()()(),,(),(), y x q x y x f x a x b y a y b αβ''-=<

常微分方程的初等解法与求解技巧

师大学本科毕业论文（设计） 常微分方程的初等解法与求解技巧 姓名娟 院系数学与计算机科学学院 专业信息与计算科学 班级12510201 学号1251020126 指导教师王晓锋 答辩日期 成绩

常微分方程的初等解法与求解技巧 容摘要 常微分方程在数学中发挥着举足轻重的作用，同时它的应用在日常生活里随处可见，因此掌握常微分方程的初等解法与求解技巧是非常必要的.本论文主要论述了其发展、初等解法与求解技巧，前者主要有变量分离、积分因子、一阶隐式微分方程的参数表示，通过举例从中总结出其求解技巧，目的是掌握其求解技巧. 【关键词】变量分离一阶隐式微分方程积分因子求解技巧

Elementary Solution and Solving Skills of Ordinary Differential Equation Abstract Ordinary differential equations take up significant position in mathematics, and at the same time, the application of it can be seen everywhere in our daily life, therefore, it’s necessary to grasp the elementary solution of ordinary differential equations and solving skills. This paper mainly introduced the definition of ordinary differential equations, elementary solution method and solving skills, the former mainly included the separation of variables, integral factor, a parameter-order differential equations implicit representation, by way of examples to sum up their solving skills, the purpose is to master the skills to solve. 【Key Words】the separation of variables the first order implicit differential equation integrating factor solution techniques

微积分(B)常微分方程与差分方程 练习题

For personal use only in study and research; not for commercial use 2013-2014(2) 大学数学(B) 练习题 第六章 For personal use only in study and research; not for commercial use 一、选择题 1. 微分方程xy y 2='的通解为 （ ） A. C e y x +=2 ； B. 2 x Ce y =； For personal use only in study and research; not for commercial use C. 2 C x y e =； D. x Ce y =. 2. 函数221x c y c e +=是微分方程20y y y '''--=的 （ ） A. 通解； B. 特解； C. 不是解； D. 是解, 但既不是通解, 也不是特解. 3. 设线性无关的函数321,,y y y 都是二阶非齐次线性微分方程)()()(x f y x q y x p y =+'+''的解， 21,C C 是任意常数，则该方程的通解是 （ ） A. 32211y y C y C ++； B. 3212211)(y C C y C y C +-+； C. 3212211)1(y C C y C y C ---+； D. 3212211)1(y C C y C y C --++. 4. 微分方程22y x y y x += +'是 （ ） A. 可分离变量的微分方程； B. 齐次微分方程； C. 一阶线性齐次微分方程； D. 一阶线性非齐次微分方程. 二、填空题 1. 微分方程y y y x ln ='的通解是 . 2. 方程x y y sin 2='的奇解为_______________.

常微分方程组(边值)

常微分方程组边值问题解法 打靶法Shooting Method （shooting.m ） %打靶法求常微分方程的边值问题 function [x,a,b,n]=shooting(fun,x0,xn,eps) if nargin<3 eps=1e-3; end x1=x0+rand; [a,b]=ode45(fun,[0,10],[0,x0]'); c0=b(length(b),1); [a,b]=ode45(fun,[0,10],[0,x1]'); c1=b(length(b),1); x2=x1-(c1-xn)*(x1-x0)/(c1-c0); n=1; while (norm(c1-xn)>=eps & norm(x2-x1)>=eps) x0=x1;x1=x2; [a,b]=ode45(fun,[0,10],[0,x0]'); c0=b(length(b),1); [a,b]=ode45(fun,[0,10],[0,x1]'); c1=b(length(b),1) x2=x1-(c1-xn)*(x1-x0)/(c1-c0); n=n+1; end x=x2; 应用打靶法求解下列边值问题： ()()??? ????==- =010004822y y y dx y d 解：将其转化为常微分方程组的初值问题

()????? ? ?????==-==t dx dy y y y dx dy y dx dy x 0011221 048 命令： x0=[0:0.1:10]; y0=32*((cos(5)-1)/sin(5)*sin(x0/2)-cos(x0/2)+1); 真实解 plot(x0,y0,'r') hold on [x,y]=ode45('odebvp',[0,10],[0,2]'); plot(x,y(:,1)) [x,y]=ode45('odebvp',[0,10],[0,5]'); plot(x,y(:,1)) [x,y]=ode45('odebvp',[0,10],[0,8]'); plot(x,y(:,1)) [x,y]=ode45('odebvp',[0,10],[0,10]'); plot(x,y(:,1))

第八章常微分方程数值解

151 第八章 常微分方程数值解 在工程和科学技术的实际问题中，常常需求解常微分方程。但由常微分方程理论可 知，常微分方程中往往只有少数较简单和典型的方程可求出其解析解。在大多数情况下，常微分方程只能用近似法求解。这种近似解法可分为两大类：一类是近似解析法，如级数解法、逐次逼近法等；另一类则是数值解法，它给出方程在一些离散点上的近似解。 本章主要讨论一阶常微分方程的初值问题： ()()?????==0 ,y a y y x f dx dy b x a ≤≤ （8.1） 从理论上讲，只要方程中的()y x f ,连续且关于y 满足李普希兹（Lipschitz ）条件，即存在常数L ，使 ()()2121,,y y L y x f y x f -≤- 则常微分方程存在唯一解)(x y y =。 所谓微分方程数值解，就是求微分方程的解()x y 在一系列离散节点 b x x x x a n n =<<<<=-110 处()i x y 的近似值i y ),,1,0(n i =. 相邻的两个节点之间的距离i i i x x h -=+1称为由i x 到1+i x 的步长，通常取为常数h 。 求数值解，首先应将微分方程离散化，常用的方法有： （1） 用差商代替微商 若用向前差商代替微商，即 ()() ()()()i i i i i x y x f x y h x y x y ,1='≈-+ )1,,1,0(-=n i 代入（8.1）中的微分方程，则得 ()1+i x y ()()()i i i x y x hf x y ,+≈

152 记)(i x y 的近似值i y ，则由上式右端可计算出)(1+i x y 的近似值，即 ()i i i i y x hf y y ,1+=+ )1,,1,0(-=n i （8.2） （2） 数值积分法 利用数值积分法左矩形公式 ()()i i x y x y -+1=()()()i i x x y x hf dx x y x f i i ,,1 ≈? + 可得同样算法 ()i i i i y x hf y y ,1+=+ （3） 用泰勒（Taylor ）公式 将函数)(x y 在i x 处展开，取一次Taylor 多项式近似，则得 ()()h x y x y i i +=+1()()i i x y h x y '+≈()()()i i i x y x hf x y ,+= 从而也得到离散化得计算公式 ()i i i i y x hf y y ,1+=+ §1 欧拉（Euler ）方法 1.1欧拉方法 对一阶微分方程（8.1），把区间[]b a ,作n 等分：b x x x x a n n =<<<<=-110 ， 则分点为 ih a x i +=， n a b h -= ),2,1(n i = 由以上讨论可知，无论用一阶向前差商，还是用数值积分法左矩形公式，或者用泰勒公式取前两项都可得到同样的离散化计算公式 ()i i i i y x hf y y ,1+=+ 并将初值条件代入，则得到数值算法： () ()? ? ?=+=+a y y y x hf y y i i i i 01, ),2,1(n i = （8.3） 称其为欧拉方法。 几何上欧拉方法就是用一条折线近似表示曲线()x y y =（如图8-1）。因此欧拉方法又称为欧拉折线方法。

常微分方程和差分方程解法归纳

常微分方程解法归纳 1. 一阶微分方程部分 ① 可分离变量方程（分离变量法） 如果一阶微分方程),(y x f dx dy =中的二元函数),(y x f 可表示为) ()(),(y h x g y x f =的形式，我们称)()(y h x g dx dy =为可分离变量的方程。 对于这类方程的求解我们首先将其分离变量为 dx x g y h dy )() (=的形式，再对此式两边积分得到 C dx x g y h dy +=??)()(从而解出)()(y h x g dx dy =的解，其中C 为任意常数。 具体例子可参考书本P10—P11的例题。 ②一阶线性齐次、非齐次方程（常数变易法） 如果一阶微分方程),(y x f dx dy =中的二元函数),(y x f 可表示为 y x P x Q y x f )()(),(-=的形式，我们称由此形成的微分方程)()(x Q y x P dx dy =+为一阶线 性微分方程，特别地，当0)(≡x Q 时我们称其为一阶线性齐次微分方程，否则为一阶线性非齐次微分方程。 对于这类方程的解法，我们首先考虑一阶线性齐次微分方程 0)(=+y x P dx dy ，这是可分离变量的方程，两边积分即可得到?=-dx x P Ce y )(，其中C 为任意常数。这也是一阶线性 非齐次微分方程的特殊情况，两者的解存在着对应关系，设)(x C 来替换C ，于是一阶线性 非齐次微分方程存在着形如?=-dx x P e x C y )()(的解。将其代入)()(x Q y x P dx dy =+我们就可 得到)()()()()()()()()(x Q e x C x P e x C x P e x C dx x P dx x P dx x P =?+?-?'---这其实也就是 ? ='dx x P e x Q x C )()()(，再对其两边积分得C dx e x Q x C dx x P +? =? )()()(，于是将其回代入 ? =-dx x P e x C y )()(即得一阶线性微分方程)()(x Q y x P dx dy =+的通解? ? ? ??+??=?-C dx e x Q e y dx x P dx x P )()()(。 具体例子可参照书本P16—P17的例题。

第十一章 常微分方程边值问题的数值解法汇总

第十一章 常微分方程边值问题的数值解法 工程技术与科学实验中提出的大量问题是常微分方程边值问题．本章将研究常微分方程边值问题的数值求解方法.主要介绍三种边界条件下的定解问题和两大类求解边值问题的数值方法,打靶法算法和有限差分方法. 11.1 引言 在很多实际问题中都会遇到求解常微分方程边值问题. 考虑如下形式的二阶常微分方程 ),,(y y x f y '='', b x a <<, (11.1.1) 在如下三种边界条件下的定解问题: 第一种边界条件: α=)(a y , β=)(b y (11.1.2) 第二种边界条件: α=')(a y , β=')(b y (11.1.2) 第三种边界条件: ? ? ?=-'=-'101 0)()()()(b b y b y a a y a y βα, (11.1.13) 其中0 0, ,00000>+≥≥b a b a . 常微分方程边值问题有很多不同解法, 本书仅介绍打靶方法和有限差分方法. 11.2 打靶法 对于二阶非线性边值问题 ()()().,,βα==≤≤'=''b y a y b x a y y x f y ，，， (11.2.1) 打靶法近似于使用初值求解的情况． 我们需要利用一个如下形式问题初值解的序列： ()()v a w a w b x a w w x f w ='=≤≤'='')(,,，，，α， （11.2.2） 引进参数v 以近似原边界值问题的解．选择参数k v v =，以使： ()()β==∞ →b y v b w k k ,lim ， （11.2.3）

其中),(k v x w 定义为初值问题（11.2.2）在k v v =时的解，同时()x y 定义为边值问题（11.2.1）的解． 首先定义参数0v ，沿着如下初值问题解的曲线，可以求出点),(αa 对应的初始正视图 ()()v a w a w b x a w w x f w ='=≤≤'='')(,,，，，α． （11.2.4） 如果),(0v b w 不严格收敛于β，那么我们选择1v 等值以修正近似值，直到),(0v b w 严格逼近β． 为了取得合适的参数k v ，现在假定边值问题（11.2.1）有唯一解，如果),(v x w 定义为初始问题（11.2.2）的解，那么v 可由下式确定： 0),(=-βv b w ． （11.2.5） 由于这是一个非线性方程，我们可以利用Newton 法求解．首先选择初始值0v ，然后由下式生成序列 ),)(()),((111----- =k k k k v b dv dw v b w v v β，此处),(),)(( 11--=k k v b dv dw v b dv dw ， （11.2.6） 同时要求求得),)(( 1-k v b dv dw ，因为),(v b w 的表达式未知，所以求解这个有一点难度；我们只能得到这么一系列的值。 ，，，),(),(),(),(1210-??k v b w v b w v b w v b w 假如我们如下改写初值问题（11.2.2），使其强调解对x 和v 的依赖性 ()()v v a w v a w b x a v x w v x w x f w ='=≤≤'=''),(,),(),,(,，，，α，（11.2.7） 保留初始记号以显式与x 的微分相关．既然要求当k v v =时),)((v b dv dw 的值，那么我们需要求出表达式（11.2.7）关于v 的偏导数．过程如下： )),(),,(,(),(v x w v x w x v f v x v w '??=?''? ),()),(),,(,()),(),,(,(v x v w v x w v x w x w f v x v x w v x w x x f ??'??+??'??= ) ,()),(),,(,(v x v w v x w v x w x w f ?'?''??+ 又因为x 跟v 相互独立，所以当b x a ≤≤上式如下；

常微分方程

中国海洋大学本科生课程大纲 课程属性：公共基础/通识教育/学科基础/专业知识/工作技能，课程性质：必修、选修 一、课程介绍 1.课程描述： 常微分方程是数学学科的一门基础理论课程,是一门专业必修课，是数学分析，高等代数和解析几何的综合应用和发展。通过学习不仅可加强先修课程中已学过的概念和方法，且为后续课程的学习准备解决问题的方法和工具。常微分方程与微积分同时诞生，以解决天文学、力学等实际问题而闻名于世，是研究事物、物体和现象运动、演化和变化规律的最基本的数学理论和方法。现实生活中许多原理和规律都可以描述成适当的常微分方程，如物体运动、生物群体竞争、疾病的传播等。对这些规律的描述、认识和分析，往往可以归结为用常微分方程描述的数学模型的分析和研究。由此可知，它是数学理论联系实际的重要渠道之一，它有着深刻而生动的实际背景，它从生产实践与科学技术中产生，而又成为现代科学技术中分析问题与解决问题的一种强有力的工具。 2.设计思路： 本课程适用于数学与应用数学专业、信息与计算科学专业二年级本科生。本课程主要包括六个部分内容：初等积分法；基本定理；一阶线性微分方程组；n 阶线性微分方程；定性理论与稳定性理论简介；一阶偏微分方程初步。 - 1 -

初等积分法主要讲解几类能用初等（积分）解法求解的方程类型及其求解方法。要求学生掌握各种类型的解法，具有判断一个给定方程的类型和正确求解的能力。重点是求解方法，难点是识别方程的类型以及熟练掌握求解方法。 基本定理包括解的存在唯一性定理，解的延展定理，解对初值的连续依赖性定理和解的可微性定理，构成了常微分方程主要理论部分。解的存在唯一性定理表明，若右端函数满足连续和利布希兹条件，则保证方程的解存在性与唯一性。它是常微分方程理论中最基本的定理，有其重大的理论意义。另一方面，由于能求得精确解的方程不多，所以该定理给出的求近似解法就具有重要的实际意义。解的延拓定理及解对初值的连续依赖性与可微性定理揭示了微分方程的重要性质。要求学生必需理解本章定理的条件和结论，掌握证明方法，能运用定理证明有关问题。重点是证明的思路和方法，特别是逐次逼近法，难点是贯穿定理证过程的利布希兹条件运用和证明过程中不等式技巧的把握； 一阶线性微分方程组主要讲线性微分方程组的理论。线性微分方程组理论是微分方程理论中的重要部分，无论从实用的角度或从理论的角度来说，本章所提供的方法和结果都是非常重要的，它是进一步学习常微分方程理论和其它有关课程必不可少的基本知识，基本要求：（1）理解线性微分方程组解的存在与唯一性定理，熟练掌握逐步逼近法；（2）掌握线性微分方程组的一般理论，把握解空间的代数结构；（3）基解矩阵求法。一般齐次线性微分方程组的基解矩阵是难以通过积分求得，但当系数矩阵是常系数矩阵时，可以通过代数方法（Jordan标准型、矩阵指数）求出基解矩阵。重点掌握一阶线性微分方程组的解空间结构和常系数线性微分方程组的解法，难点是证明一阶齐次常微分方程组的解空间是n 维线性空间和一阶常系数齐次或非齐次微分方程组的求解。 - 1 -

郑州大学研究生课程数值分析复习---第八章 常微分方程数值解法

郑州大学研究生课程（2012-2013学年第一学期）数值分析 Numerical Analysis 习题课 第八章常微分方程数值解法

待求解的问题：一阶常微分方程的初值问题/* Initial-Value Problem */: ?????=∈=0 )(] ,[),(y a y b a x y x f dx dy 解的存在唯一性（“常微分方程”理论）：只要f (x , y ) 在[a , b ] ×R 1 上连续，且关于y 满足Lipschitz 条件，即存在与x , y 无关的常数L 使 对任意定义在[a , b ] 上的y 1(x ) 和y 2(x ) 都成立，则上述IVP 存在唯一解。 1212|(,)(,)||| f x y f x y L y y ?≤?一、要点回顾

§8.2 欧拉(Euler)法 通常取(常数),则Euler 法的计算格式 h h x x i i i ==?+1?? ?=+=+) (),(001x y y y x hf y y i i i i i =0,1,…,n ( 8.2 )

§8.2 欧拉(Euler)法(1) 用差商近似导数 )) (,()()()()(1n n n n n n x y x hf x y x y h x y x y +=′+≈+?? ?=+=+) (),(01a y y y x hf y y n n n n 差分方程初值问题向前Euler 方法h x y x y x y n n n ) ()()(1?≈ ′+)) (,() ()(1n n n n x y x f h x y x y ≈?+))(,()(n n n x y x f x y =′

常微分方程的初等解法与求解技巧

山西师范大学本科毕业论文（设计） 常微分方程的初等解法与求解技巧 姓名张娟 院系数学与计算机科学学院 专业信息与计算科学 班级12510201 学号1251020126 指导教师王晓锋 答辩日期 成绩

常微分方程的初等解法与求解技巧 内容摘要 常微分方程在数学中发挥着举足轻重的作用，同时它的应用在日常生活里随处可见，因此掌握常微分方程的初等解法与求解技巧是非常必要的.本论文主要论述了其发展、初等解法与求解技巧，前者主要有变量分离、积分因子、一阶隐式微分方程的参数表示，通过举例从中总结出其求解技巧，目的是掌握其求解技巧. 【关键词】变量分离一阶隐式微分方程积分因子求解技巧

Elementary Solution and Solving Skills of Ordinary Differential Equation Abstract Ordinary differential equations take up significant position in mathematics, and at the same time, the application of it can be seen everywhere in our daily life, therefore, it’s necessary to grasp the elementary solution of ordinary differential equations and solving skills. This paper mainly introduced the definition of ordinary differential equations, elementary solution method and solving skills, the former mainly included the separation of variables, integral factor, a parameter-order differential equations implicit representation, by way of examples to sum up their solving skills, the purpose is to master the skills to solve. 【Key Words】the separation of variables the first order implicit differential equation integrating factor solution techniques

常微分方程组(边值)

常微分方程组边值问题解法 打靶法Shooti ng Method (shoot in g.m ) % 丁靶法求常微分方程的边值问题 function [x,a,b ,n]=shooti ng(fu n, xO,x n, eps) if nargin<3 eps=1e-3; end x1=x0+ra nd; [a,b]=ode45(fu n, [0,10],[0,x0]'); c0=b(le ngth(b),1); [a,b]=ode45(fu n, [0,10],[0,x1]'); c1=b(le ngth(b),1); x2=x1-(c1-x n)*(x1-x0)/(c1-c0); n=1; while (no rm(c1-x n)>=eps & no rm(x2-x1)>=eps) x0=x1;x 仁x2; [a,b]=ode45(fu n,[ 0,10],[0,x0]'); cO=b(le ngth(b),1); [a,b]=ode45(fu n,[ 0,10],[0,x1]'); c1= b(le ngth(b),1) x2=x1-(c1-x n)*(x1-x0)/(c1-c0); n=n+1; end x=x2; 应用打靶法求解下列边值问题: y 10 0 解：将其转化为常微分方程组的初值问题

命令: xO=[O:O.1:1O]; y0=32*((cos(5)-1)/si n( 5)*si n(x0/2)-cos(x0/2)+1); plot(xO,yO,'r') hold on [x,y]=ode45('odebvp',[0,10],[0,2]'); plot(x,y(:,1)) [x,y]=ode45('odebvp',[0,10],[0,5]'); plot(x,y(:,1)) [x,y]=ode45('odebvp',[0,10],[0,8]'); plot(x,y(:,1)) [x,y]=ode45('odebvp',[0,10],[0,10]'); plot(x,y(:,1)) dy i dx y 2 dy 2 dx y i 0 y 4 y o dy dx X0 真实解 30 ' 12^4567^9 10

第8章 常微分方程数值解法 本章主要内容： 1.欧拉法、改进欧拉法 2

第8章 常微分方程数值解法 本章主要内容： 1．欧拉法、改进欧拉法. 2．龙格-库塔法。 3．单步法的收敛性与稳定性。 重点、难点 一、微分方程的数值解法 在工程技术或自然科学中，我们会遇到的许多微分方程的问题，而我们只能对其中具有较简单形式的微分方程才能够求出它们的精确解。对于大量的微分方程问题我们需要考虑求它们的满足一定精度要求的近似解的方法，称为微分方程的数值解法。本章我们主要 讨论常微分方程初值问题?????==00 )() ,(y x y y x f dx dy 的数值解法。 数值解法的基本思想是：在常微分方程初值问题解的存在区间[a,b]内，取n+1个节点a=x 0＜x 1＜…＜x N =b （其中差h n = x n –x n-1称为步长，一般取h 为常数，即等步长），在这些节点上把常微分方程的初值问题离散化为差分方程的相应问题，再求出这些点的上的差分方程值作为相应的微分方程的近似值（满足精度要求）。 二、欧拉法与改进欧拉法 欧拉法与改进欧拉法是用数值积分方法对微分方程进行离散化的一种方法。 将常微分方程),(y x f y ='变为() *+=?++1 1))(,()()(n x n x n n dt t y t f x y x y 1．欧拉法（欧拉折线法） 欧拉法是求解常微分方程初值问题的一种最简单的数值解法。 欧拉法的基本思想：用左矩阵公式计算（＊）式右端积分，则得欧拉法的计算公式为：N a b h N n y x hf y y n n n n -= -=+=+)1,...,1,0(),(1 欧拉法局部截断误差 11121 )(2 ++++≤≤''=n n n n n x x y h R ξξ或简记为O （h 2）。

微分方程与差分方程详细讲解与例题

第七章 常微分方程与差分方程 常微分方程是高等数学中理论性和应用性都较强的一部分，是描述客观规律的一种重要方法，是处理物理、力学、几何等应用问题的一个重要工具，微分和积分的知识是研究微分方程的基础。微分方程作为考试的重点容，每年研究生考试均会考到。特别是微分方程的应用问题，既是重点，也是难点，在复习时必须有所突破。 【数学一大纲容】常微分方程的基本概念；变量可分离的方程；齐次方程；一阶线性方程；伯努利（Bernoulli ）方程；全微分方程；可用简单的变量代换求解的某些微分方程；可降阶的高阶微分方程；线性微分方程解的性质及解的结构定理；二阶常系数齐次线性微分方程；高于二阶的某些常系数齐次线性微分方程；简单的二阶常系数非齐次线性微分方程；欧拉（Euler ）方程；微分方程的简单应用。 【数学二大纲容】常微分方程的基本概念；变量可分离的方程；齐次方程；一阶线性微分方程；可降阶的高阶微分方程；线性微分方程解的性质及解的结构定理；二阶常数齐次线性微分方程；高于二阶的某些常系数齐次线性微分方程；简单的二阶常系数非齐次线性微分方程；微分方程的一些简单应用。 【大纲要求】要理解微分方程的有关概念，如阶、解、通解、特解、定解条件等，掌握几类方程的解法：如变量可分离方程，齐次方程，一阶线性微分方程，伯努利方程，可降阶方程等。理解线性微分方程解的性质和解的结构，掌握求解常系数齐次线性方程的方法，掌握求解某些自由项的常系数非齐次线性方程的待定系数法。了解欧拉方程的概念，会求简单的欧拉方程。会用微分方程处理物理、力学、几何中的简单问题。 【考点分析】本章包括三个重点容： 1．常见的一阶、二阶微分方程求通解或特解。求解常微分方程重要的是判断方程为哪种类型，并记住解法的推导过程。 2．微分方程的应用问题，这是一个难点，也是重点。利用微分方程解决实际问题时，若是几何问题，要根据问题的几何特性建立微分方程。若是物理问题，要根据某些物理定律建立微分方程，也有些问题要利用微元法建立微分方程。 3．数学三要求掌握一阶常系数线性差分方程的求解方法，了解差分与差分方程及其通解与特解等概念，会用差分方程求解简单的经济应用问题。 【考点八十三】形如()()y f x g y '=的一阶微分方程称为变量可分离微分方程。可分离变量的微分方程的解题程序： 当()0,()()()() dy g y y f x g y f x dx g y '≠=? =时，然后左、右两端积分 (),()dy f x dx C g y =+?? 上式即为变量可分离微分方程的通解。其中，C 为任意常数，1 ()() dy g y g y ? 表示函数的一个原函数，()f x dx ?表示函数()f x 的一个原函数. 【例7.1】微分方程1+++='y x xy y 的通解为____________。

第八章 常微分方程答案(2012[1].6)

例1 微分方程221y x y xy '=-+-满足1)0(=y 的特解为 ． 解：2 2 2(1)(1)(1)(1)11dy dy y x y x dx x dx y y '=-+? =-?=-++?? 解得 2 arctan 2 x y x C =-+，由0 14 x y C π ==?= 则方程的特解为 2arctan 24 x y x π=-+ 或 2tan()24x y x π =-+ 例2 解微分方程3 23 x xy y y -='． 解：323x xy y y -='即为3 2 1y x y y x ?? ???'=?? - ??? ，为齐次微分方程.令y u y xu y u xu x ''=?=?=+， 由已知321 u y u '=-，整理得211 u du dx u x -=， 两边积分得 2 22ln ln ln ln 2ln 22u u y u x C Cy Cy x ?? -=+?=?= ??? 则方程的通解为 2 2ln y Cy x ?? = ??? . 例3 微分方程x y y x ln =+'满足1)1(=y 的特解为 ． 解：原方程整理得1ln x y y x x '+ =，为一阶线性非齐次微分方程. 由通解公式得 11 ln 1ln ln 1dx dx x x x C y e e dx C xdx C x x x x - ??????=+=+=-+???????? 由1)1(=y 解得2C =，所以微分方程x y y x ln =+'满足1)1(=y 的特解为2 ln 1.y x x =-+

例4 微分方程3 1 y xy y += '的通解为 ． 解： 3 3dx dx xy y yx y dy dy =+? -=， 通解为 2 22 32 22232y y y ydy ydy e y e dy C Ce y x e y e dy C --???? +????? ???? =-?-? ?=+=?? 例5 解微分方程y x y y x 24=-'． ……① 解 原方程可化为y x y x y =?- '4 （2 1 =α的贝努里方程），即 x y x y y =?-'4 1 ……② 作换元y u = ，则 y y dx du 2' = ，②可化为 22x u x dx du =-（一阶线性非齐次方程） ……③ 由常数变易法可得③的通解为： )2ln (2x C x u + =， 故原方程通解为 )2 ln (2x C x y + =． 例6 已知函数(),()f x g x 满足x e x g x f x f x g x g x f 2)()(),()(),()(=+='='，且()00f =， 求)()()(x g x f x F =所满足的一阶微分方程，并求)(x F 的表达式． 解：(1) 由)()()()()(x g x f x g x f x F '+'='=)()(2 2x f x g + =)()(2)]()([2 x g x f x g x f -+)(242x F e x -=, 可见，)(x F 所满足的一阶微分方程为 2()2()4(0)0 x F x F x e F '?+=? =?．