微分方程及其解的定义

微分方程及其解的定义
微分方程及其解的定义

微分方程

什么是微分方程它是怎样产生的这是首先要回答的问题.

300多年前,由牛顿(Newton,1642-1727)和莱布尼兹(Leibniz,1646-1716)所创立的微积分学,是人类科学史上划时代的重大发现,而微积分的产生和发展,又与求解微分方程问题密切相关.这是因为,微积分产生的一个重要动因来自于人们探求物质世界运动规律的需求.一般地,运动规律很难全靠实验观测认识清楚,因为人们不太可能观察到运动的全过程.然而,运动物体(变量)与它的瞬时变化率(导数)之间,通常在运动过程中按照某种己知定律存在着联系,我们容易捕捉到这种联系,而这种联系,用数学语言表达出来,其结果往往形成一个微分方程.一旦求出这个方程的解,其运动规律将一目了然.下面的例子,将会使你看到微分方程是表达自然规律的一种最为自然的数学语言.

例1 物体下落问题

设质量为m的物体,在时间t=0时,在距地面高度为H处以初始速度v(0) = v0垂直地面下落,求此物体下落时距离与时间的关系.

解如图1-1建立坐标系,设为t时刻物体的位置坐标.于是物体下落的速度为

加速度为

质量为m的物体,在下落的任一时刻所受到的外力有重力mg和空气阻力,当速度不太大时,空气阻力可取为与速度成正比.于是根据牛顿第二定律

F = ma (力=质量×加速度)

可以列出方程

(·= )

其中k >0为阻尼系数,g是重力加速度.

式就是一个微分方程,这里t是自变量,x是未知函数,是未知函数对t导数.现在,我们还不会求解方程,但是,如果考虑k=0的情形,即自由落体运动,此时方程可化为

将上式对t积分两次得

其中和是两个独立的任意常数,它是方程的解.

一般说来,微分方程就是联系自变量、未知函数以及未知函数的某些导数之间的关系式.如果其中的未知函数只与一个自变量有关,则称为常微分方程;如果未知函数是两个或两个以上自变量的函数,并且在方程中出现偏导数,则称为偏微分方程.本书所介绍的都是常微分方程,有时就简称微分方程或方程.

例如下面的方程都是常微分方程

()

()

(·=)()

(′=)()

在一个常微分方程中,未知函数最高阶导数的阶数,称为方程的阶.这样,一阶常微分方程的一般形式可表为

()如果在中能将y′解出,则得到方程

()

或()

称为一阶隐式方程,称为一阶显式方程,称为微分形式的一阶方程.

n 阶隐式方程的一般形式为

()n 阶显式方程的一般形式为

在方程中,如果左端函数F对未知函数y和它的各阶导数y′,y″,…,y(n)的全体而言是一次的,则称为线性常微分方程,否则称它为非线性常微分方程.这样,一个以y为未知函数,以x为自变量的n阶线性微分方程具有如下形式:

()

显然,方程是一阶线性方程;方程是一阶非线性方程;方程是二阶线性方程;方程是二阶非线性方程.

通解与特解

微分方程的解就是满足方程的函数,可定义如下.

定义1.1设函数在区间I上连续,且有直到n阶的导数.如果把

代入方程,得到在区间I上关于x的恒等式,

则称为方程在区间I上的一个解.

这样,从定义可以直接验证:

1. 函数y = x2+C是方程在区间(-∞,+∞)上的解,其中C是任意的常数.

2. 函数是方程在区间(-1,+1)上的解,其中C是任意常数.又方程有两个明显的常数解y =±1,这两个解不包含在上述解中.

3. 函数是方程在区间(-∞,+∞)上的解,其中和是独立的任意常数.

4. 函数是方程(1.7)在区间(-∞,+∞)上的解,其中和是独立的任意常数.

这里,我们仅验证3,其余留给读者完成.事实上,在(-∞,+∞)上有

所以在(-∞,+∞)上有

从而该函数是方程的解.

从上面的讨论中,可以看到一个重要事实,那就是微分方程的解中可以包含任意常数,其中任意常数的个数可以多到与方程的阶数相等,也可以不含任意常数.我们把n

阶常微分方程的含有n个独立的任意常数C1,C2,…,Cn的解,称为该方程的通解,如果方程的解不包含任意常数,则称它为特解.由隐式表出的通解称为通积分,而由隐式表出的特解称为特积分.

由上面的定义,不难看出,函数和

分别是方程,和的通解,函数是方程的通积分,而函数

y =±1是方程的特解.通常方程的特解可对通解中的任意常数以定值确定,这种确定过程,需要下面介绍的初始值条件,或简称初值条件.

初值问题

例1中的函数显然是方程的通解,由于和是两个任意常数,这表明方程有无数个解,解的图像见下面的图a和图b所示.

图a(C1>固定,C2>0)图b(C1=0,C2>0)而实际经验表明,一个自由落体运动仅能有一条运动轨迹.产生这种多解性的原因是因为方程所表达的是任何一个自由落体,在任意瞬时t所满足的关系式,并未考虑运动的初始状态,因此,通过积分求得的其通解所描述的是任何一个自由落体的运动规律.显然,在同一初始时刻,从不同的高度或以不同初速度自由下落的物体,应有不同的运动轨迹.为了求解满足初值条件的解,我们可以把例1中给出的两个初始值条件,即

初始位置x(0)= H 初始速度

代入到通解中,推得

于是,得到满足上述初值条件的特解为

它描述了初始高度为H,初始速度为v0的自由落体运动规律.

求微分方程满足初值条件的解的问题称为初值问题.

于是我们称是初值问题

的解.

对于一个n 阶方程,初值条件的一般提法是

其中是自变量的某个取定值,而是相应的未知函数及导数的给定值.方程的初值问题常记为

初值问题也常称为柯希(Cauchy)问题.

对于一阶方程,若已求出通解,只要把初值条件

代入通解中,得到方程

从中解出C,设为,代入通解,即得满足初值条件的解.

对于n 阶方程,若已求出通解后,代入初值条件,得

到n个方程式

如果能从式中确定出,代回通解,即得所求初值问题的. 例2 求方程

的满足初值条件的解.

解方程通解为

求导数后得

将初值条件代入,得到方程组

解出和得

故所求特解为

积分曲线

为了便于研究方程解的性质,我们常常考虑解的图象.一阶方程的一个特解的图象是xoy平面上的一条曲线,称为方程的积分曲线,而通解的图象是平面上的一族曲线,称为积分曲线族.例如,方程的通解+C是xoy平面上

的一族抛物曲线.而是过点(0,0)的一条积分曲线.以后,为了叙述简便,我们对解和积分曲线这两个名词一般不加以区别.对于二阶和二阶以上的方程,也有积分曲线和积分曲线族的概念,只不过此时积分曲线所在的空间维数不同,我们将在第4章详细讨论.

最后,我们要指出,本书中按习惯用

分别代表,

分别代表

本节要点:

1.常微分程的定义,方程的阶,隐式方程,显式方程,线性方程,非线性方程.

2.常微分方程解的定义,通解,特解,通积分,特积分.

3.初值问题及初值问题解的求法.

4.解的几何意义,积分曲线.

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

第七章 常微分方程初值问题的数值解法 --------学习小结 一、本章学习体会 通过本章的学习,我了解了常微分方程初值问题的计算方法,对于解决那些很难求解出解析表达式的,甚至有解析表达式但是解不出具体的值的常微分方程非常有用。在这一章里求解常微分方程的基本思想是将初值问题进行离散化,然后进行迭代求解。在这里将初值问题离散化的方法有三种,分别是差商代替导数的方法、Taylor 级数法和数值积分法。常微分方程初值问题的数值解法的分类有显示方法和隐式方法,或者可以分为单步法和多步法。在这里单步法是指计算第n+1个y 的值时,只用到前一步的值,而多步法则是指计算第n+1个y 的值时,用到了前几步的值。通过对本章的学习,已经能熟练掌握如何用Taylor 级数法去求解单步法中各方法的公式和截断误差,但是对线性多步法的求解理解不怎么透切,特别是计算过程较复杂的推理。 在本章的学习过程中还遇到不少问题,比如本章知识点多,公式多,在做题时容易混淆,其次对几种R-K 公式的理解不够透彻,处理一个实际问题时,不知道选取哪一种公式,通过课本里面几种方法的计算比较得知其误差并不一样,,这个还需要自己在往后的实际应用中多多实践留意并总结。 二、本章知识梳理 常微分方程初值问题的数值解法一般概念 步长h ,取节点0,(0,1,...,)n t t nh n M =+=,且M t T ≤,则初值问题000 '(,),()y f t y t t T y t y =≤≤?? =?的数值解法的一般形式是 1(,,,...,,)0,(0,1,...,)n n n n k F t y y y h n M k ++==-

微分方程的基础知识及解析解

微分方程的基础知识及解析解

微分方程的基础知识与练习 (一)微分方程基本概念: 首先通过一个具体的问题来给出微分方程的基本概念。 (1)一条曲线通过点(1,2),且在该曲线上任一点M (x ,y )处的切线的斜率为2x ,求这条曲线的方程。 解 设曲线方程为)(x y y =.由导数的几何意义可知函数)(x y y =满足 x dx dy 2= (1) 同时还满足以下条件: 1=x 时,2=y (2) 把(1)式两端积分,得 ?=xdx y 2 即 C x y +=2 (3) 其中C 是任意常数。 把条件(2)代入(3)式,得 1=C , 由此解出C 并代入(3)式,得到所求曲线方程: 12+=x y (4) (2)列车在水平直线路上以20s m /的速度行驶;当制动时列车获得加速度2/4.0s m -.问开始制动后多少时间列车才能停住,以及列车在这段时间里行驶了多少路程? 解 设列车开始制动后t 秒时行驶了s 米。根据题意,反映制动阶段列车运动规律的函数)(t s s =满足: 4.02 2-=dt s d (5) 此外,还满足条件: 0=t 时,20,0== =dt ds v s (6) (5)式两端积分一次得: 14.0C t dt ds v +-== (7) 再积分一次得

2122.0C t C t s ++-= (8) 其中21,C C 都是任意常数。 把条件“0=t 时20=v ”和“0=t 时0=s ”分别代入(7)式和(8)式,得 0 ,2021==C C 把21,C C 的值代入(7)及(8)式得 ,204.0+-=t v (9) t t s 202.02+-= (10) 在(9)式中令0=v ,得到列车从开始制动到完全停止所需的时间: )(504 .020s t ==。 再把5=t 代入(10)式,得到列车在制动阶段行驶的路程 ).(5005020502.02m s =?+?-= 上述两个例子中的关系式(1)和(5),(6)都含有未知函数的导数,它们都是微分方程。 1.微分方程的概念 一般地,凡含有未知函数、未知函数的导数及自变量的方程,叫做微分方程。未知函数是一元函数的方程叫做常微分方程;未知函数是多元函数的方程,叫做偏微分方程。我们只研究常微分方程。微分方程中所出现的未知函数的最高阶导数的阶数,叫做微分方程的阶。 例如,方程(1)是一阶微分方程;方程(5)是二阶微分方程方程。又如,方程 ()x y y y y y 2sin 5'12''10'''44=+-+-是四阶微分方程。 一般地,n 阶微分方程的形式是 ()(,,',...,)0,n F x y y y = (11) 其中F 是个2+n 变量的函数。这里必须指出,在方程(11)中,)(n y 是必须出现的,而 )1(,...,',,-n y y y x 等变量则可以不出现。例如n 阶微分方程

常微分方程数值解

第四章常微分方程数值解 [课时安排]6学时 [教学课型]理论课 [教学目的和要求] 了解常微分方程初值问题数值解法的一些基本概念,如单步法和多步法,显式和隐式,方法的阶数,整体截断误差和局部截断误差的区别和关系等;掌握一阶常微分方程初值问题的一些常用的数值计算方法,例如欧拉(Euler)方法、改进的欧拉方法、龙贝-库塔(Runge-Kutta)方法、阿达姆斯(Adams)方法等,要注意各方法的特点及有关的理论分析;掌握构造常微分方程数值解的数值积分的构造方法和泰勒展开的构造方法的基本思想,并能具体应用它们导出一些常用的数值计算公式及评估截断误差;熟练掌握龙格-库塔(R-K)方法的基本思想,公式的推导,R-K公式中系数的确定,特别是能应用“标准四阶R-K公式”解题;掌握数值方法的收敛性和稳定性的概念,并能确定给定方法的绝对稳定性区域。[教学重点与难点] 重点:欧拉方法,改进的欧拉方法,龙贝-库塔方法。 难点:R—K方法,预估-校正公式。 [教学内容与过程] 4.1 引言 本章讨论常微分方程初值问题 (4.1.1) 的数值解法,这也是科学与工程计算经常遇到的问题,由于只有很特殊的方程能用解析方法求解,而用计算机求解常微分方程的初值问题都要采用数值方法.通常我们假定(4.1.1)中 f(x,y)对y满足Lipschitz条件,即存在常数L>0,使对,有 (4.1.2) 则初值问题(4.1.1)的解存在唯一. 假定(4.1.1)的精确解为,求它的数值解就是要在区间上的一组离散点 上求的近似.通常取 ,h称为步长,求(4.1.1)的数值解是按节点的顺序逐步 推进求得.首先,要对方程做离散逼近,求出数值解的公式,再研究公式的局部截

常微分方程基本概念习题附解答

§1.2 常微分方程基本概念习题及解答 1.dx dy =2xy,并满足初始条件:x=0,y=1的特解。 解:y dy =2xdx 两边积分有:ln|y|=x 2+c y=e 2x +e c =cex 2另外y=0也是原方程的解,c=0时,y=0 原方程的通解为y= cex 2,x=0 y=1时 c=1 特解为y= e 2 x . 2. y 2dx+(x+1)dy=0 并求满足初始条件:x=0,y=1的特解。 解:y 2dx=-(x+1)dy 2y dy dy=-11+x dx 两边积分: -y 1=-ln|x+1|+ln|c| y=|)1(|ln 1+x c 另外y=0,x=-1也是原方程的解 x=0,y=1时 c=e 特解:y=| )1(|ln 1+x c 3.dx dy =y x xy y 32 1++ 解:原方程为:dx dy =y y 21+31x x + y y 21+dy=3 1x x +dx 两边积分:x(1+x 2)(1+y 2)=cx 2 4. (1+x)ydx+(1-y)xdy=0 解:原方程为: y y -1dy=-x x 1+dx 两边积分:ln|xy|+x-y=c

另外 x=0,y=0也是原方程的解。 5.(y+x )dy+(x-y)dx=0 解:原方程为: dx dy =-y x y x +- 令 x y =u 则dx dy =u+x dx du 代入有: -112++u u du=x 1dx ln(u 2+1)x 2=c-2arctgu 即 ln(y 2+x 2)=c-2arctg 2x y . 6. x dx dy -y+22y x -=0 解:原方程为: dx dy =x y +x x ||-2)(1x y - 则令 x y =u dx dy =u+ x dx du 211 u - du=sgnx x 1dx arcsin x y =sgnx ln|x|+c 7. tgydx-ctgxdy=0 解:原方程为:tgy dy =ctgx dx 两边积分:ln|siny|=-ln|cosx|-ln|c| siny=x c cos 1=x c cos 另外y=0也是原方程的解,而c=0时,y=0. 所以原方程的通解为sinycosx=c. 8 dx dy +y e x y 32 +=0 解:原方程为:dx dy =y e y 2e x 3

常微分方程数值解法

第七章 常微分方程数值解法 常微分方程中只有一些典型方程能求出初等解(用初等函数表示的解),大部分的方程是求不出初等解的。另外,有些初值问题虽然有初等解,但由于形式太复杂不便于应用。因此,有必要探讨常微分方程初值问题的数值解法。本章主要介绍一阶常微分方程初值问题的欧拉法、龙格-库塔法、阿达姆斯方法,在此基础上推出一阶微分方程组与高阶方程初值问题的 数值解法;此外,还将简要介绍求解二阶常微分方程值问题的差分方法、试射法。 第一节 欧拉法 求解常微分方程初值问题 ?????==0 0)() ,(y x y y x f dx dy (1) 的数值解,就是寻求准确解)(x y 在一系列离散节点 <<<<

常微分方程基本知识点

常微分方程基本知识点 第一章 绪论 1. 微分方程的概念(常微分与偏微),什么是方程的阶数,线性与非线性,齐次与非齐次,解、特解、部分解和通解的概念及判断! (重要) 例:03)(22=-+y dx dy x dx dy (1阶非线性); x e dx y d y =+22sin 。 2.运用导数的几何意义建立简单的微分方程。(以书后练习题为主) (习题1,2,9题) 例:曲线簇cx x y -=3满足的微分方程是:__________. 第二章 一阶方程的初等解法 1.变量分离方程的解法(要能通过适当的变化化成变量分离方程);(重要) 2.齐次方程的解法(变量代换);(重要) 3.线性非齐次方程的常数变易法; 4.分式线性方程、贝努利方程、恰当方程的概念及判断(要能熟练的判断各种类型的一阶方程)(重要); 例题:(1).经变换_____y c u os =___________后, 方程1cos sin '+=+x y y y 可化为___线性_____方程; (2).经变换_____y x u 32-=____________后, 方程1 )32(1 '2+-=y x y 可化为____变量分离__方程; (3).方程0)1(222=+-dy e dx ye x x x 为:线性方程。

(4).方程221 'y x y -=为:线性方程。 5.积分因子的概念,会判断某个函数是不是方程的积分因子; 6.恰当方程的解法(分项组合方法)。(重要) 第三章 一阶方程的存在唯一性定理 1.存在唯一性定理的内容要熟记,并能准确确定其中的h ; 2.会构造皮卡逐步逼近函数序列来求第k 次近似解!(参见书上例题和习题 3.1的1,2,3题) 第四章 高阶微分方程 1.n 阶线性齐次(非齐次)微分方程的概念,解的概念,基本解组,解的线性相关与线性无关,齐次与非齐次方程解的性质; 2.n 阶线性方程解的Wronskey 行列式与解的线性相关与线性无关的关系; 3.n 阶线性齐次(非齐次)微分方程的通解结构定理!!(重要) 4.n 阶线性非齐次微分方程的常数变易法(了解); 5.n 阶常系数线性齐次与非齐次微分方程的解法(Eurler 待定指数函数法确定基本解组),特解的确定(比较系数法、复数法);(重要) 例题:t te x x 24=-'',确定特解类型? (习题4.2相关题目) 6.2阶线性方程已知一个特解的解法(作线性齐次变换)。(重要) 7.其他如Euler 方程、高阶方程降阶、拉普拉斯变换法等了解。

常微分方程的数值解

实验4 常微分方程的数值解 【实验目的】 1.掌握用MATLAB软件求微分方程初值问题数值解的方法; 2.通过实例用微分方程模型解决简化的实际问题; 3.了解欧拉方法和龙格-库塔方法的基本思想和计算公式,及稳定性等概念。 【实验内容】 题3 小型火箭初始重量为1400kg,其中包括1080kg燃料。火箭竖直向上发射时燃料燃烧率为18kg/s,由此产生32000N的推力,火箭引擎在燃料用尽时关闭。设火箭上升时空气阻力正比于速度的平方,比例系数为m,求引擎关闭瞬间火箭的高度、速度、加速度,及火箭到达最高点的时的高度和加速度,并画出高度、速度、加速度随时间变化的图形。 模型及其求解 火箭在上升的过程可分为两个阶段,在全过程中假设重力加速度始终保持不变,g=s2。 在第一个过程中,火箭通过燃烧燃料产生向上的推力,同时它还受到自身重力(包括自重和该时刻剩余燃料的重量)以及与速度平方成正比的空气阻力的作用,根据牛顿第二定律,三个力的合力产生加速度,方向竖直向上。因此有如下二式: a=dv/dt=/m=/(1400-18t) dh/dt=v 又知初始时刻t=0,v=0,h=0。记x(1)=h,x(2)=v,根据MATLAB 可以求出0到60秒内火箭的速度、高度、加速度随时间的变化情况。程序如下: function [ dx ] = rocket( t,x ) a=[*x(2)^2)/(1400-18*t)]; dx=[x(2);a]; end ts=0:1:60;

x0=[0,0]; [t,x]=ode45(@rocket,ts,x0); h=x(:,1); v=x(:,2); a=[*(v.^2))./(1400-18*t)]; [t,h,v,a]; 数据如下: t h v a 000

微分方程的基本概念

求函数关系是数学中的重要问题。然而,在实际中有时很难直接找出函数关系,我们所得到的仅是含有未知函数及其导数的关系式,称之为微分方程.我们的任务就是求解微分方程,找出未知函数。本章将介绍一些微分方程的基本概念和几种常用的微分方程的解法. 微分方程的基本概念 下面通过几个例题来说明微分方程的基本概念. 例1 一曲线通过)2,1(点,且在该曲线上任一点),(y x 处 的切线的斜率为x 2,求曲线的方程. 解 由导数的几何意义可得 x dx dy 2= ① 此外,未知函数)(x y y =还应满足条件 1=x 时,2=y (或写成21==x y ) ② 在式①两端积分,得 C x y +=2 , ③ 其中C 为任意常数.将条件②代入式③中,得1=C , 于是得所求曲线的方程为 ④ 12+=x y

我们知道式③表示一族曲线, 曲线族中的每一条曲线的函数 代入式①中都成为恒等式, 而式④仅表示是其中的一条,它是通过点()2,1的. 从以上例子中,可归纳出如下一些基本概念. (一)微分方程:含有自变量、未知函数以及未知函数导数或微分的方程叫微分方程(以下简称方程)。在方程中出现的未知函数导数的最高阶数成为微分方程的阶,n 阶微分方程的一般形式为 ()(,,,,,)0n F x y y y y '''=L ⑤ 如式①为一阶微分方程.

(二)解:一个函数代入微分方程后,使其成为恒等式,则该函数称为微分方程的解. 含有任意常数,且独立的任意常数的个数和微分方程的阶数相等的解,称为微分方程的通解或一般解.不含任意常数的解叫特解. 若I x x y ∈=),(?为方程⑤的解,则有 ()[,(),(),,()]0n F x x x x φφφ'≡L , I x ∈. 方程⑤的通解应含有n 个独立的任意常数, 其通解有时用隐函数表达式 12(,,,,,)0n x y C C C Φ=L 表示. ⑥ 例如:式③为方程①的通解.

常微分方程教材

第九章 微分方程 一、教学目标及基本要求 (1) 了解微分方程及其解、通解、初始条件和特解的概念。 (2) 掌握变量可分离的方程和一阶线性方程的解法,会解齐次方程。 (3) 会用降阶法解下列方程:),(),,(),()(y y f y y x f y x f y n '='''=''=。 (4) 理解二阶线性微分方程解的性质以及解的结构定理。 (5) 掌握二阶常系数齐次线性微分方程的解法,并会解某些高于二阶的常系数齐次线性微分方程。 (6) 会求自由项多项式、指数函数、正弦函数、余弦函数,以及它们的和与二阶常系数非齐次线性微分方程的 特解和通解。 (7) 会用微分方程解决一些简单的应用问题。 二、本章教学内容的重点和难点 1、理解和熟悉微分方程的一些基本概念; 2、掌握一阶和高阶微分方程的各种初等积分法; 3、熟悉线性方程的基础理论,掌握常系数二阶线性齐次与非齐次方程的解法; 4、会列微分方程及其始值问题去解决实际问题。 三、本章教学内容的深化和拓宽: 1、分离变量法的理论根据; 2、常用的变量代换; 3、怎样列微分方程解应用题; 4、黎卡提方程; 5、全微分方程的推广; 6、二阶齐次方程; 7、高阶微分方程的补充; 8、求线性齐次方程的另一个线性无关的解; 9、求线性非齐次方程的一个特解; 10、常数变易法。 本章的思考题和习题 解下列方程(第1-6题) 1、2)0(,)1(==+'+y x y y x 2、()[]f dx x f e e x f x x x ,)(02?+=可微 3、212 22sin 22sin 1X e y x y y x ++='?+ 4、0)3(24=+-xydx dy x y 5、21)0(,1)0(,022- ='=='+''y y y x y 6、2y y y x y '-'+'= 7、已知可微函数)(x f 满足 ?-=+x x f f x f x x f dx x f 12)()1(,1)()()(和求; 8、已知)(,,1)(2 1)(10x f f x f da ax f 求可微+= ?; 9、求与曲线族C y x =+2232相交成ο45角的曲线; 10、一容器的容积为100L ,盛满盐水,含10kg 的盐,现以每分钟3L 的速度向容器内注入淡水冲淡盐水,又以同样的速度将盐水抽入原先盛满淡水的同样大小的另一容器内,多余的水便从容器内流出,问经过多少时间,两容器内的含盐量相等?

微分方程(习题及解答)

第十二章 微分方程 § 微分方程基本概念、可分离变量的微分方程、齐次微分方程 一、单项选择题 1. 下列所给方程中,不是微分方程的是( ) . (A)2xy y '=; (B)222x y C +=; (C)0y y ''+=; (D)(76)d ()d 0x y x x y y -++=. 答(B). 2. 微分方程4(3)520y y xy y '''+-=的阶数是( ). (A)1; (B)2; (C)3; (D)4; 答(C). 3. 下列所给的函数,是微分方程0y y ''+=的通解的是( ). (A)1cos y C x =; (B)2sin y C x =; (C)cos sin y x C x =+; (D)12cos sin y C x C x =+ 答(D). 4. 下列微分方程中,可分离变量的方程是( ). (A)x y y e +'=; (B)xy y x '+=; (C)10y xy '--=; (D)()d ()d 0x y x x y y -++=. 答(A). 5. 下列微分方程中,是齐次方程是微分方程的是( ). (A)x y y e +'=; 2(B)xy y x '+=; (C)0y xy x '--=; (D)()d ()d 0x y x x y y -++=. 答(D). 二、填空题 1.函数25y x =是否是微分方程2xy y '=的解 . 答:是 . 2.微分方程 3d d 0,4x x y y y x =+==的解是 . 答:2225x y +=. 3.微分方程2 3550x x y '+-=的通解是 . 答:32 52 x x y C =++. 4.微分方程ln 0xy y y '-=的通解是 . 答: Cx y e =. 5'=的通解是 . 答:arcsin arcsin y x C =+. 6.微分方程 (ln ln )xy y y y x '-=-的通解是. 答: Cx y e x =. 三、解答题 1.求下列微分方程的通解. (1) 22sec tan d sec tan d 0x y x y x y +=; (2) 2()y xy a y y '''-=+; 解: 解: (3) d 10d x y y x +=; (4) 23d (1)0.d y y x x ++=

《常微分方程》课程大纲

《常微分方程》课程大纲 一、课程简介 课程名称:常微分方程学时/学分:3/54 先修课程:数学分析,高等代数,空间解析几何,或线性代数(行列式,矩阵与线性方程组,线性空间F n,欧氏空间R n,特征值与矩阵的对角化), 高等数学(多元微积分,无穷级数)。 面向对象:本科二年级或以上学生 教学目标:围绕基本概念与基本理论、具体求解和实际应用三条主线开展教学活动,通过该课程的教学,希望学生正确理解常微分方程的基本概念,掌握基本理论和主要方法,具有一定的解题能力和处理相关应用问题的思维方式,如定性分析解的性态和定量近似求解等思想,并希望学生初步了解常微分方程的近代发展,为学习动力系统学科的近代内容和后续课程打下基础。 二、教学内容和要求 常微分方程的教学内容分为七部分,对不同的内容提出不同的教学要求。(数字表示供参考的相应的学时数,第一个数为课堂教学时数,第二个数为习题课时数) 第一章基本概念(2,0) (一)本章教学目的与要求: 要求学生正确掌握微分方程,通解,线性与非线性,积分曲线,线素场(方

向场),定解问题等基本概念。本章教学重点解释常微分方程解的几何意义。 (二)教学内容: 1.由实际问题:质点运动即距离与时间关系(牛顿第二运动定律),放射性元素衰变过程,人口总数发展趋势估计等,通过建立数学模型,导出微分方程。 2.基本概念(常微分方程,偏微分方程,阶,线性,非线性,解,定解问题,特解,通解等)。 3.一阶微分方程组的几何定义,线素场(方向场),积分曲线。 4.常微分方程所讨论的基本问题。 第二章初等积分法(4,2) (一)本章教学目的与要求: 要求学生熟练掌握分离变量法,常数变易法,初等变换法,积分因子法等初等解法。 本章教学重点对经典的几类方程介绍基本解法,勾通初等积分法与微积分学基本定理的关系。并通过习题课进行初步解题训练,提高解题技巧。 (二)教学内容: 1. 恰当方程(积分因子法); 2. 分离变量法 3. 一阶线性微分方程(常数变易法) 4. 初等变换法(齐次方程,伯努利方程,黎卡提方程)

微分方程中的几个基础概念

微分方程中的几个基础概念 微分方程—基础 微分方程(Differential equation, DFQ)是一种用来描述函数与其导数之间关系的数学方程。与之前所接触初等数学代数方程的解不同,它的解不是数,而是符合方程关系的函数。 微分方程的起源约在十七世纪末,为了解决自然科学发展中遇到物理及天文学问题而产生,随着微积分的诞生与在各个科学领域中的广泛应用,很多问题被归化为某类微分方程的问题。 在微分方程分支中,存在很多各种各样已知类型的微分方程。实事上,提高对微分方程的理解的最好的方法之一是首先处理基本的分类系统。为什么?因为你可能永远不会遇到完全陌生的微分方程。大多数微分方程已经被解决了,因此,普遍适用的解决方法很可能已经存在。 除了描述方程本身的性质外,对微分方程进行分类和识别的真正附加值来自于为跳转点提供一张导图。求解微分方程的诀窍不是创造原始解法,而是对已证明的解法进行分类和应用;有时,可能需要几步把一类方程转换为另一类等效方程,以获得可实现的广义解。 最常用于描述微分方程的四个属性是: ?常微分与偏微分 ?线性与非线性 ?齐次与非齐次

?微分阶数 虽然这个列表并非详尽无遗,但是它是我们学习首先要掌握的知识,通常在微分方程学期课程的前几周会进行回顾;通过快速回顾每一个类别,我们将会配备基本的入门工具包来处理常见的微分方程问题。 常微分与偏微分 首先,我们在自然中所发现的微分方程最常见的分类来源于从我们手边的问题中所发现的导数类型;简单地说,方程是否包含偏导数? 如果不包含,那么它是一个常微分方程(, Ordinary differential equation)。如果包含,那么它是一个偏微分方程(, Partial differential equation)。 常微分方程是未知函数只含有一个自变量的微分方程,其微分基于该单一的自变量,通常是时间。一个常微分方程有一组离散的(有限的)变量;它们通常是一维动力系统的模型,例如:钟摆随时间的摆动。 另一方面,偏微分方程相当复杂,因为它们通常涉及多个自变量,其多种多样的偏微分方程可能基于也可能并不基于一个已知的自变量。偏微分方程常被用来描述自然界中各种各样的现象,例如:热,空间中的流体速度,或电动力学。这些似乎完全不同的物理现象被化为偏微分方程;它们在随机偏微分方程中得到推广。 下面的这些例子有助于我们分辨微分方程的导数类型包括:

常微分方程数值解法

i.常微分方程初值问题数值解法 常微分方程初值问题的真解可以看成是从给定初始点出发的一条连续曲线。差分法是常微分方程初值问题的主要数值解法,其目的是得到若干个离散点来逼近这条解曲线。有两个基本途径。一个是用离散点上的差商近似替代微商。另一个是先对微分方程积分得到积分方程,再利用离散点作数值积分。 i.1 常微分方程差分法 考虑常微分方程初值问题:求函数()u t 满足 (,), 0du f t u t T dt =<≤ (i.1a ) 0(0)u u = (i.1b) 其中(,)f t u 是定义在区域G : 0t T ≤≤, u <∞上的连续函数,0u 和T 是给定的常数。我们假设(,)f t u 对u 满足Lipschitz 条件,即存在常数L 使得 121212(,)(,), [0,]; ,(,)f t u f t u L u u t T u u -≤-?∈∈-∞∞ (i.2) 这一条件保证了(i.1)的解是适定的,即存在,唯一,而且连续依赖于初值0u 。 通常情况下,(i.1)的精确解不可能用简单的解析表达式给出,只能求近似解。本章讨论常微分方程最常用的近似数值解法-差分方法。先来讨论最简单的Euler 法。为此,首先将求解区域[0,]T 离散化为若干个离散点: 0110N N t t t t T -=<< <<= (i.3) 其中n t hn =,0h >称为步长。 在微积分课程中我们熟知,微商(即导数)是差商的极限。反过来,差商就是微商的近似。在0t t =处,在(i.1a )中用向前差商 10()()u t u t h -代替微商du dt ,便得 10000()()(,())u t u t hf t u t ε=++ 如果忽略误差项0ε,再换个记号,用i u 代替()i u t 便得到 1000(,)u u hf t u -= 一般地,我们有 1Euler (,), 0,1, ,1n n n n u u hf t u n N +=+=-方法: (i.4) 从(i.1b) 给出的初始值0u 出发,由上式可以依次算出1,,N t t 上的差分解1,,N u u 。

微分方程的基础知识与练习

微分方程的基础知识与练习 (一)微分方程基本概念: 首先通过一个具体的问题来给出微分方程的基本概念。 (1)一条曲线通过点(1,2),且在该曲线上任一点M (x ,y )处的切线的斜率为2x ,求这条曲线的方程。 解 设曲线方程为)(x y y =.由导数的几何意义可知函数)(x y y =满足 x dx dy 2= (1) 同时还满足以下条件: 1=x 时,2=y (2) 把(1)式两端积分,得 ?=xdx y 2 即 C x y +=2 (3) 其中C 是任意常数。 把条件(2)代入(3)式,得 1=C , 由此解出C 并代入(3)式,得到所求曲线方程: 12+=x y (4) (2)列车在水平直线路上以20s m /的速度行驶;当制动时列车获得加速度 2/4.0s m -.问开始制动后多少时间列车才能停住,以及列车在这段时间里行驶了 多少路程? 解 设列车开始制动后t 秒时行驶了s 米。根据题意,反映制动阶段列车运 动规律的函数)(t s s =满足: 4.02 2-=dt s d (5) 此外,还满足条件: 0=t 时,20,0== =dt ds v s (6) (5)式两端积分一次得: 14.0C t dt ds v +-== (7) 再积分一次得

2122.0C t C t s ++-= (8) 其中21,C C 都是任意常数。 把条件“0=t 时20=v ”和“0=t 时0=s ”分别代入(7)式和(8)式,得 0 ,2021==C C 把21,C C 的值代入(7)及(8)式得 ,204.0+-=t v (9) t t s 202.02+-= (10) 在(9)式中令0=v ,得到列车从开始制动到完全停止所需的时间: )(504 .020 s t == 。 再把5=t 代入(10)式,得到列车在制动阶段行驶的路程 ).(5005020502.02m s =?+?-= 上述两个例子中的关系式(1)和(5),(6)都含有未知函数的导数,它们 都是微分方程。 1.微分方程的概念 一般地,凡含有未知函数、未知函数的导数及自变量的方程,叫做微分方程。未知函数是一元函数的方程叫做常微分方程;未知函数是多元函数的方程,叫做偏微分方程。我们只研究常微分方程。微分方程中所出现的未知函数的最高阶导数的阶数,叫做微分方程的阶。 例如,方程(1)是一阶微分方程;方程(5)是二阶微分方程方程。又如,方程 ()x y y y y y 2sin 5'12''10'''44=+-+-是四阶微分方程。 一般地,n 阶微分方程的形式是 ()(,,',...,)0,n F x y y y = (11) 其中F 是个2+n 变量的函数。这里必须指出,在方程(11)中,)(n y 是必须出现的,而 )1(,...,',,-n y y y x 等变量则可以不出现。例如n 阶微分方程

§1常微分方程的基本概念

第十三章 常微分方程简介 本章介绍微分方程的有关概念及某些简单微分方程的解法。 微分方程是包含未知函数及其导数的方程。由微分方程能够求出未知函数的解析表达式,从而掌握所研究的客观现象的变化规律和发展趋势。因此,掌握这方面的知识,用之分析解决问题是非常重要的。 由于在大多数情况下,微分方程很难求出初等解(即解的形式是初等函数)。那么,就需要研究解的存在理论,借助计算机求出微分方程的数值解。 本章的内容,仅仅包含常微分方程的一些最初步的知识,特殊的一阶和部分二阶微分方程的初等解法;最后一节讨论微分方程的简单应用。 §1 常微分方程的基本概念 像过去我们研究其他许多问题一样,首先通过具体实际例子来引入微分方程的概念。 两个实例 例1.1 设某一平面曲线上任意一点),(y x 处的切线斜率等于该点处横坐标x 的2倍,且曲线通过点)2,1(,求该曲线的方程。 解 平面上的曲线可由一元函数来表示 设所求的曲线方程为)(x f y =,根据导数的几何意义,由题意得 x dx dy 2=(这是一个含未知函数)(x f y =的导数的方程)。 另外,由题意,曲线通过点)2,1(,所以,所求函数)(x f y =还满足2|1==x y 。 从而得到 12 (1.1)|2(1.2) x dy x dx y =ì??=?í??=??,。 为了解出)(x f y =,我们只要将的两端积分,得 ?+=+==C x C x xdx y 22 2 22, 我们说 C x y +=2对于任意常数C 都满足方程。 再由条件,将2|1==x y 代入C x y +=2,即

C +=2121=?C 。 故所求曲线的方程为12+=x y 。 再看一个例子: 例1.2 设质点以匀加速度a 作直线运动,且0=t 时0,0v v s ==。求质点运 动的位移与时间t 的关系。 解 这是一个物理上的运动问题。 设质点运动的位移与时间的关系为 )(t s s =。 则由二阶导数的物理意义,知a t d s d =22,这是一个含有二阶导数的方程。 再由题意000 |0 |t t s v v ==ì=??í ?=??,因此,)(t S S =应满足问题 22 000 (1.3)|0|(1.4)t t d s a dt s v v ==ì??=?í??==???,,。 要解这个问题,我们可以将两边连续积分两次,即 1C at dt ds +=, ??++=21C dt C tdt a s ,即 2122 C t C t a s ++=, 其中21,C C 为任意常数。 由条件,因为0|0==t s ,代入,得02=C ; 再由00|v v t ==,代入,得01v C =。 故得 t v t a s 02 2 += 为所求。 下面我们将通过分析这两个具体的例子,给出微分方程的一些基本概念。 微分方程的基本概念 总结所给出的两个具体的例子,我们看到: (1) 例的)1(式和例 的)1(式都是含有未知函数的导数的等式(例1含一阶导数,例2含二阶导数); (2) 通过积分可以解出满足这等式的函数;

常微分方程数值解法

第八章 常微分方程数值解法 考核知识点: 欧拉法,改进欧拉法,龙格-库塔法,单步法的收敛性与稳定性。 考核要求: 1. 解欧拉法,改进欧拉法的基本思想;熟练掌握用欧拉法,改进欧拉法、求微 分方程近似解的方法。 2. 了解龙格-库塔法的基本思想;掌握用龙格-库塔法求微分方程近似解的方 法。 3. 了解单步法的收敛性、稳定性与绝对稳定性。 例1 用欧拉法,预估——校正法求一阶微分方程初值问题 ? ??=-='1)0(y y x y ,在0=x (0,1)0.2近似解 解 (1)用1.0=h 欧拉法计算公式 n n n n n n x y y x y y 1.09.0)(1.01+=-+=+,1.0=n 计算得 9.01=y 82.01.01.09.09.02=?+?=y (2)用预估——校正法计算公式 1,0)(05.01.09.0)0(111)0(1=???-+-+=+=++++n y x y x y y x y y n n n n n n n n n 计算得 91.01=y ,83805.02=y 例2 已知一阶初值问题 ???=-='1 )0(5y y y 求使欧拉法绝对稳定的步长n 值。 解 由欧拉法公式 n n n n y h y h y y )51(51-=-=+ n n y h y ~)51(~1-=+

相减得01)51()51(e h e h e n n n -==-=-Λ 当 151≤-h 时,4.00≤

微分方程及其解的定义

微分方程 什么是微分方程它是怎样产生的这是首先要回答的问题. 300多年前,由牛顿(Newton,1642-1727)和莱布尼兹(Leibniz,1646-1716)所创立的微积分学,是人类科学史上划时代的重大发现,而微积分的产生和发展,又与求解微分方程问题密切相关.这是因为,微积分产生的一个重要动因来自于人们探求物质世界运动规律的需求.一般地,运动规律很难全靠实验观测认识清楚,因为人们不太可能观察到运动的全过程.然而,运动物体(变量)与它的瞬时变化率(导数)之间,通常在运动过程中按照某种己知定律存在着联系,我们容易捕捉到这种联系,而这种联系,用数学语言表达出来,其结果往往形成一个微分方程.一旦求出这个方程的解,其运动规律将一目了然.下面的例子,将会使你看到微分方程是表达自然规律的一种最为自然的数学语言. 例1 物体下落问题 设质量为m的物体,在时间t=0时,在距地面高度为H处以初始速度v(0) = v0垂直地面下落,求此物体下落时距离与时间的关系. 解如图1-1建立坐标系,设为t时刻物体的位置坐标.于是物体下落的速度为 加速度为 质量为m的物体,在下落的任一时刻所受到的外力有重力mg和空气阻力,当速度不太大时,空气阻力可取为与速度成正比.于是根据牛顿第二定律 F = ma (力=质量×加速度) 可以列出方程 (·= ) 其中k >0为阻尼系数,g是重力加速度. 式就是一个微分方程,这里t是自变量,x是未知函数,是未知函数对t导数.现在,我们还不会求解方程,但是,如果考虑k=0的情形,即自由落体运动,此时方程可化为 将上式对t积分两次得 其中和是两个独立的任意常数,它是方程的解. 一般说来,微分方程就是联系自变量、未知函数以及未知函数的某些导数之间的关系式.如果其中的未知函数只与一个自变量有关,则称为常微分方程;如果未知函数是两个或两个以上自变量的函数,并且在方程中出现偏导数,则称为偏微分方程.本书所介绍的都是常微分方程,有时就简称微分方程或方程.

常微分方程的基本概念

考点:常微分方程的基本概念【☆☆☆☆☆】 1.微分方程:含有未知函数的导数或微分的方程称为微分方程. 若未知函数是一元函数,则称为常微分方程; 若未知函数是多元函数,则称为偏微分方程. 考题链接: 例:*320y x y x y xdy ydx ''=++=+=,, 2.阶:未知函数的最高阶导数的阶数. 考题链接: 例:微分方程()2 420x y y x y '''+-=的阶数是( ) A.1 B.2 C.3 D.4 3.性微分方程: ()()()()()()*012n n f x y f x y f x y f x y f x '?+?+?+ +?= 考题链接: 例:判断下列函数是否为线性方程. (1)2y x y '=+ (2)2sin y x y x '=++ (3)sin 0y x y '-+= (4)2y yy x '''-= (5)()2 3y x y '=+ 4.解:若()y x ?=代入方程成为恒等式,则称()y x ?=为方程的一个解. (1)通解:含有相互独立(不能合并,212y C x C x =+与12y C x C x =+)的任意常数,且任意常数的个数与方程的阶数相同的微分方程的解. (2)特解:不含任意常数的解. 例1:某二阶常微分方程的下列解中为通解的是( ) A.sin y C x = B.12sin cos y C x C x =+ C.sin cos y x x =+ D.()12cos y C C x =+

例2:函数sin y C x =(其中C 为任意常数)是微分方程0y y ''+=的( ) A.通解 B.特解 C.解 D.不是解 例3:已知微分方程x y ay e '+=的一个特解为x y xe =,则a =________. 考点:可分离变量的微分方程【☆☆☆☆☆】 (1)标准形式:()()f y dy g x dx = (2)解法:①分离变量,化为标准形式;②两边同时积分. 例1:微分方程0dx dy y x +=的通解是( ) A.2225x y += B.34x y C += C.22x y C += D.227y x -= 例2:方程22sec tan sec tan 0x ydx y xdy +=的通解为________. 例3:微分方程220dy xy dx -=满足条件()11y =-的特解是( ) A.21 y x = B.21y x =- C.2y x = D.2y x =- 考点:齐次方程【☆☆☆☆☆】 (1)标准形式:y y f x ?? = ??? 考题链接: 例:22x y x y '=+不是 222x y x y '=+是 (2)解法:①化为标准形式; ②令y u x = ,代入方程消去y ; ③化为x 与u 的可分离变量的微分方程,求解. 例:求sin 0y xy x y x '--=的通解. 考点:一阶线性微分方程【☆☆☆☆☆】 (1)标准形式:()()y P x y Q x '+=

常微分方程数值解法

第八章 常微分方程的数值解法 一.内容要点 考虑一阶常微分方程初值问题:?????==0 0)() ,(y x y y x f dx dy 微分方程的数值解:设微分方程的解y (x )的存在区间是[a,b ],在[a,b ]内取一系列节 点a= x 0< x 1<…< x n =b ,其中h k =x k+1-x k ;(一般采用等距节点,h=(b-a)/n 称为步长)。在每个节点x k 求解函数y(x)的近似值:y k ≈y(x k ),这样y 0 , y 1 ,...,y n 称为微分方程的数值解。 用数值方法,求得f(x k )的近似值y k ,再用插值或拟合方法就求得y(x)的近似函数。 (一)常微分方程处置问题解得存在唯一性定理 对于常微分方程初值问题:?????==0 0)() ,(y x y y x f dx dy 如果: (1) 在B y y A x x 00≤-≤≤,的矩形内),(y x f 是一个二元连续函数。 (2) ),(y x f 对于y 满足利普希茨条件,即 2121y y L y x f y x f -≤-),(),(则在C x x 0≤≤上方程?????==0 0)() ,(y x y y x f dx dy 的解存在且唯一,这里C=min((A-x 0),x 0+B/L),L 是利普希茨常数。 定义:任何一个一步方法可以写为),,(h y x h y y k k k 1k Φ+=+,其中),,(h y x k k Φ称为算法的增量函数。 收敛性定理:若一步方法满足: (1)是p 解的. (2) 增量函数),,(h y x k k Φ对于y 满足利普希茨条件. (3) 初始值y 0是精确的。则),()()(p h O x y kh y =-kh =x -x 0,也就是有 0x y y lim k x x kh 0h 0 =--=→)( (一)、主要算法 1.局部截断误差 局部截断误差:当y(x k )是精确解时,由y(x k )按照数值方法计算出来的1~ +k y 的误差y (x k+1)- 1~ +k y 称为局部截断误差。 注意:y k+1和1~ +k y 的区别。因而局部截断误差与误差e k +1=y (x k +1) -y k +1不同。 如果局部截断误差是O (h p+1),我们就说该数值方法具有p 阶精度。

相关文档
最新文档