差分方程真题

差分方程真题
差分方程真题

第九章 微分方程与差分方程(历年考研真题)

1、设函数()y y x =满足条件440(0)2,(0)4

y y y y y '''++=??'==? ,求广义积分0()d y x x +∞? 。 2、已知连续函数()f x 满足条件320()()d e 3

x

x t f x f t =+?,求()f x 。 3

、求微分方程d d y x =的通解。 4、设函数()f t 在[0,)+∞上连续,且满足方程

22242

4()d e x y t t f t f x y π+≤=+??,求()f t 。 5、设函数()f x 在[1,)+∞上连续。若由曲线()y f x =,直线1,(1)x x t t ==>与x 轴所围成的平面图形绕x 轴旋转一周所形成的旋转体的体积为 2()[()(1)]3V t t f t f π=

-。试求()y f x =所满足的微分方程,并求该微分方程满足条件229

x y ==的解。 6、设有微分方程2()y y x ?'-=,其中2,1;()0, 1.x x x ??

试求在(,)-∞+∞内的连续函数()y y x =,使之在(,1)-∞和(1,)+∞内都满足所给方程,且满足条件(0)0y =。

7、求微分方程22e 0x y y '''--=满足条件(0)1,(0)1y y '==的解。

8、已知()n f x 满足 1()()e n x n n f x f x x -'=+(n 为正整数),且e (1)n f n =,求函数项级数1

()n n f x ∞=∑ 之和。

9、(1)验证函数 3693()1()3!6!9!(3)!

n

x x x x y x x n =++++++-∞<<+∞L L 满足微分方程e x y y y '''++=;

(2)利用(1)的结果求幂级数30(3)!n

n x n ∞

=∑的和函数。 10、设()()()F x f x g x =,其中函数(),()f x g x 在(,)-∞+∞内满足以下条件:

()(),()(),(0)0,()()2e x

f x

g x g x f x f f x g x ''===+=且。

(1) 求()F x 所满足的一阶微分方程; (2)求出()F x 的表达式。

11、设级数 468

()242462468

x x x x +++-∞<<+∞L g g g g g g 的和函数为()S x 。求: (1)()S x 所满足的一阶微分方程; (2)()S x 的表达式。

12、设()y f x =是第一象限内连接点(0,1),(1,0)A B 的一段连续曲线,(,)M x y 为该曲线上任意一点,点C 为M 在x 轴上的投影,O 为坐标原点,若梯形OCMA 的面积与曲边三角形CBM 的面积之和为3163

x +,求()f x 的表达式。 13、设(,)f u v 具有连续偏导数,且满足(,)(,)u v f u v f u v uv ''+=,求2()(,)x y x e

f x x -=所满足的一

阶微分方程,并求其通解。

14、差分方程12t t t y y t +-=的通解为 。 15、差分方程121050t t y y t ++-=的通解为 。

16、某公司每年的工资总额在比上一年增加20%的基础上再追加2百万元,若以t W 表示第t 年的工资总额(单位:百万元),则t W 满足的差分方程是 。

附答案:

1、1

2、32()3e

2e x x f x =- 3、(0)y C x +=> 4、224()(41)e t f t t ππ=+ 5、2232x y y xy '=-,31x y x =+ 6、222e 1,1;()(1e )e , 1.

x x x y x x -?-≤=?->? 7、211(32)e 44x y x =++

8、e ln(1)x

x -- 9、(2)221()e e ()33x x y x x x -=+-∞<<+∞ 10、(1)2()2()4e x F x F x '+=; (2)22()e e x x F x -=-

11、(1)3

,(0)02

x y xy y '=+=; (2)222()e 12x x S x =-+- 12、2()(1)f x x =- 13、222e x y y x -'+=;3

2()e 3

x x y C -=+ 14、(2)2,0,1,2,t t y c t t =+-=L 15、51(5)(),0,1,2,126

t t y c t t =-+-=L 16、11.22t t W W -=+

习题详解-第10章微分方程与差分方程初步

习题10-1 1. 指出下列方程的阶数: (1)4620x y y x y '''''-+=. (2)2 2 d d 0d d Q Q Q L R t c t ++=. (3)2d cos d ρ ρθθ +=. (4)2()d 2d 0y x y x x y -+=. 解:(1)三阶(2)二阶(3)一阶(4)一阶 2. 验证下列给出的函数是否为相应方程的解: (1)2x y y '=, 2y Cx =. (2)2(+1)d d x y y x =, +1y x =. (3)20y y y '''++=, x y x e -=. (4)22d 0.4d s t =-, 2120.2s t c t c =-++. 解:(1)是,代入即可. (2)是,代入即可; (3)是,因为 ,2x x x x y e xe y e xe ----'''=-=-+,满足20y y y '''++=; (4)是,代入,2 12d d 0.4,0.4d d s s t C t t =-+=-,显然满足. 3. 验证:函数x =C 1cos kt +C 2sin kt (k ≠0)是微分方程 222d 0d x k x t += 的通解. 解:221212()sin cos ,()cos sin ,x t C k kt C k kt x t C k kt C k kt '''=-+=--满足2 22 d 0d x k x t +=,所以是解,又因为含有两个任意常数12,C C ,且方程是二阶的,故是通解. 4. 已知函数x =C 1cos kt +C 2sin kt (k ≠0)是微分方程222d 0d x k x t +=的通解,求满足初始条件 x | t 2 x | t 的特解. 解:上题可知是微分方程通解,且12()sin cos ,x t C k kt C k kt '=-+代入初值条件0|02,|0t t x x ='===,得122,0C C ==,所以特解为2cos (0).x kt k =≠ 习题10-2 1. 求下列微分方程的通解: (1)()2 310y y x '++=; (2) 2 +'=x y y ; (3) d d sin xcos y y sin y cos x x =; (4) 2 d d d d x xy y y x y y +=+; (5) 22 d d d d y y y x xy x x +=; (6) d d y x y x x y -= +; (7) 22 d d y y x xy x =+; (8) )2(tan 21 2y x y +='. 解:(1)这是可分离变量方程,分离变量得 () 2 31d =d y y x x +- 两端分别积分:

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

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

差分方程真题

第九章 微分方程与差分方程(历年考研真题) 1、设函数()y y x =满足条件440(0)2,(0)4 y y y y y '''++=??'==? ,求广义积分0()d y x x +∞? 。 2、已知连续函数()f x 满足条件320()()d e 3 x x t f x f t =+?,求()f x 。 3 、求微分方程d d y x =的通解。 4、设函数()f t 在[0,)+∞上连续,且满足方程 22242 4()d e x y t t f t f x y π+≤=+??,求()f t 。 5、设函数()f x 在[1,)+∞上连续。若由曲线()y f x =,直线1,(1)x x t t ==>与x 轴所围成的平面图形绕x 轴旋转一周所形成的旋转体的体积为 2()[()(1)]3V t t f t f π= -。试求()y f x =所满足的微分方程,并求该微分方程满足条件229 x y ==的解。 6、设有微分方程2()y y x ?'-=,其中2,1;()0, 1.x x x ?? 试求在(,)-∞+∞内的连续函数()y y x =,使之在(,1)-∞和(1,)+∞内都满足所给方程,且满足条件(0)0y =。 7、求微分方程22e 0x y y '''--=满足条件(0)1,(0)1y y '==的解。 8、已知()n f x 满足 1()()e n x n n f x f x x -'=+(n 为正整数),且e (1)n f n =,求函数项级数1 ()n n f x ∞=∑ 之和。 9、(1)验证函数 3693()1()3!6!9!(3)! n x x x x y x x n =++++++-∞<<+∞L L 满足微分方程e x y y y '''++=; (2)利用(1)的结果求幂级数30(3)!n n x n ∞ =∑的和函数。 10、设()()()F x f x g x =,其中函数(),()f x g x 在(,)-∞+∞内满足以下条件: ()(),()(),(0)0,()()2e x f x g x g x f x f f x g x ''===+=且。 (1) 求()F x 所满足的一阶微分方程; (2)求出()F x 的表达式。

差分方程的解法

1、常系数线性差分方程的解 方程( 8)其中为常数,称方程(8)为常系数线性方程。 又称方程(9) 为方程(8)对应的齐次方程。 如果(9)有形如的解,带入方程中可得: (10) 称方程(10)为方程(8)、(9)的特征方程。 显然,如果能求出(10)的根,则可以得到(9)的解。 基本结果如下: (1)若(10)有k个不同的实根,则(9)有通解: , (2)若(10)有m重根,则通解中有构成项: (3)若(10)有一对单复根,令:,,则(9)的通解中有构成项: (4)若有m 重复根:,,则(9)的通项中有成项:

综上所述,由于方程(10)恰有k 个根,从而构成方程 (9)的通解中必有k个独立的任意常数。通解可记为: 如果能得到方程(8)的一个特解:,则(8)必有通解: + (11) (1)的特解可通过待定系数法来确定。 例如:如果为n 的多项式,则当b不是特征根时,可设成形如形式的特解,其中为m次多项式;如果b是r重根时,可设特解:,将其代入(8)中确定出系数即可。 2、差分方程的z变换解法 对差分方程两边关于取Z变换,利用的Z 变换F(z)来表示出的Z变换,然后通过解代数方程求出F(z),并把F(z)在z=0的解析圆环域中展开成洛朗级数,其系数就是所要求的 例1设差分方程,求 解:解法1:特征方程为,有根: 故:为方程的解。 由条件得: 解法2:设F(z)=Z(),方程两边取变换可得:

由条件得 由F(z)在中解析,有 所以, 3、二阶线性差分方程组 设,,形成向量方程组 (12)则 (13)(13)即为(12)的解。 为了具体求出解(13),需要求出,这可以用高等代数的方法计算。常用的方法有: (1)如果A为正规矩阵,则A必可相似于对角矩阵,对角线上的元素就是A的特征值,相似变换矩阵由A的特征向量构成:。 (2)将A 分解成为列向量,则有 从而,

差分方程模型的稳定性分析分析解析

分类号 学号密题 目 (中、英文) 作者姓名 指导教师 学科门类 提交论文日期专业名称 成绩评定 数学与应用数学 理 学

咸阳师范学院2016届本科毕业设计(论文) 摘要 微分方程是研究数学的一个重要分支,是本科期间我们必须掌握的基本知识,而本文我们研究的是一个递推关系式,也称差分方程。它是一种离散化的微分方程,是利用描述客观事物的数量关系的一种重要的数学思想来建立模型的。而利用差分方程建立模型解决问题的方法在生活中随处可见,比如在自由竞争市场经济中的蛛网模型是利用差分方程分析经济何时趋于稳定,又如金融问题中的养老保险也是利用差分方程来分析保险品种的实际投资价值。而差分方程模型是描述客观世界中随离散时间变量演化规律的有力建模工具。本文首先给出差分方程的定义以及求解过程并给出判断差分方程稳定性的判断方法,随后以同一环境下的羊群和草群的相互作用为模型分析其种群的数量变化过程,进而研究线性差分方程的稳定性,最后用一个实际模型来更好的说明差分方程的稳定性对解决实际问题有非常大的帮助。 关键字:差分方程;差分方程模型;平衡点;稳定性

差分方程模型的稳定性分析 Abstract Difference equation is also called recursive equation, it is to describe the relationship between the number of objective things of a kind of important mathematical model. And the use of the differential equation model of the solution can be found everywhere in life. Such as cobweb model in the free market economy is to use the difference equation analysis when the economic stability, and as the financial problem of pension insurance breed difference equation is used to analysis the actual investment value. This paper gives the judge the stability of difference equation to judge method, then in the same group of sheep and grass under the environment of interaction analysis for the model a process, the number of the population change, in turn, study the stability of the linear difference equation. In the end, one practical model to better explain the stability of difference equation. Key words:Difference equation;Difference equation model ; Balance point; Stability

(完整版)差分方程模型(讲义)

差分方程模型 一. 引言 数学模型按照离散的方法和连续的方法,可以分为离散模型和连续模型。 1. 确定性连续模型 1) 微分法建模(静态优化模型),如森林救火模型、血管分支模型、最优价格模型。 2) 微分方程建模(动态模型),如传染病模型、人口控制与预测模型、经济增长模型。 3) 稳定性方法建模(平衡与稳定状态模型),如军备竞赛模型、种群的互相竞争模型、种群的互相依存模型、种群弱肉强食模型。 4) 变分法建模(动态优化模型),如生产计划的制定模型、国民收入的增长模型、渔业资源的开发模型。 2. 确定性离散模型 1) 逻辑方法建模,如效益的合理分配模型、价格的指数模型。 2) 层次分析法建模,如旅游景点的选择模型、科研成果的综合评价模型。 3)图的方法建模,如循环比赛的名次模型、红绿灯的调节模型、化学制品的存放模型。 4)差分方程建模,如市场经济中的蛛网模型、交通网络控制模型、借贷模型、养老基金设置模型、人口的预测与控制模型、生物种群的数量模型。 随着科学技术的发展,人们将愈来愈多的遇到离散动态系统的问题,差分方程就是建立离散动态系统数学模型的有效方法。 在一般情况下,动态连续模型用微分方程方法建立,与此相适应,当时间变量离散化以后,可以用差分方程建立动态离散模型。有些实际问题既可以建立连续模型,又可建立离散模型,究竟采用那种模型应视建模的目的而定。例如,人口模型既可建立连续模型(其中有马尔萨斯模型Malthus、洛杰斯蒂克Logistic模型),又可建立人口差分方程模型。这里讲讲差分方程在建立离散动态系统数学模型的的具体应用。

二. 差分方程简介 在实际中,许多问题所研究的变量都是离散的形式,所建立的数学模型也是离散的,譬如,像政治、经济和社会等领域中的实际问题。有些时候,即使所建立的数学模型是连续形式,例如像常见的微分方程模型、积分方程模型等。但是,往往都需要用计算机求数值解。这就需要将连续变量在一定的条件下进行离散化,从而将连续型模型转化为离散型模型。因此,最后都归结为求解离散形式的差分方程解的问题。关于差分方程理论和求解方法在数学建模和解决实际问题的过程中起着重要作用。 1. 差分方程的定义 给定一个数列{}n x , 把数列中的前1+n 项i x ),,2,1,0(n i Λ=关联起来得到的方程,则称这个方程为差分方程。 2. 常系数线性齐次差分方程 常系数线性齐次差分方程的一般形式为 02211=++++---k n k n n n x a x a x a x Λ, (1) 或者表示为 0),,,,(1=++k n n n x x x n F Λ (1’) 其中k 为差分方程的阶数,其中k a a a ,,,21Λ为差分方程的系数,且0≠k a )(n k ≤。 对应的代数方程 02211=++++--k k k k a a a Λλλλ (2) 称为差分方程(1)的对应的特征方程。(2)式中的k 个根k λλλ,,,21Λ称为(1)式的特征根。 2.1 差分方程的解 常系数线性齐次差分方程的解主要是由相应的特征根的不同情况有不同的形式。下面分别就特征根为单根、重根和复根的情况给出方程解的形式。 2.1.1 特征根为单根(互不相同的根) 设差分方程(1)有k 个单特征根(互不相同的根)k λλλ,,,21Λ,则

差分方程与概率计算

第4期 随着科学技术的发展,差分方程在各个领域得到越来越多的应用,本文将介绍差分方程的一个简单的应用,即如何利用差分方程来求概率问题,虽然差分方程及其解法在很多方面类似于微分方程,但由于很少书籍介绍差分方程的内容,现在先了解一下差分方程的基本概念。 1.差分的概念[1] 定义1 设y(t)为定义在整数集上的函数,则称△y(t)=y(t+1)-y(t)为函数y(t)的一阶差分, △(△y(t))=△2y(t)称为y(t)的二阶差分,△ny(t)=△(△n-1y(t))称为y(t)的n阶差分。 对于连续函数y(t),可以在区间[a,b]内插入n-1个分点:a<t0<t1<…<tn=b(为方便计算,可取等距离点),得函数值y(t0),y(t1),…,y(tn),同样定义y(t)的各阶差分。 上述定义也可以称为向前差分,还可以用不同的形式定义向后差分与中心差分,三者实质是相同的,可以互相转换。差分具有线性运算及类似微分的运算性质。 2.差分方程的概念[1] 定义2差分方程的一般形式为:F(y(t);△y(t),…,△ny(t))=0,方程中的最大足标i+n与最小足标 i之差为n时,称之为n阶的差分方程,其一般形式为:a0(t)y(t)+a1(t)y(t)+…+an(t)y(t)=b(t),当b(t)=0 时,称为其次的,否则称为非其次的。 在求概率中应用到的一类差分方程,是一类简单的特殊形式,常用到的只有一阶常系数线性差分方程和二阶常系数线性差分方程,其一般形式为:xn+1=axn+b (1);xn+2=axn+1+bxn(2) 3.一阶、二阶常系数线性差分方程的解[2]引理1 对于一阶常系数线性差分方程xn+1=axn+bxn,a,b为常数,若已知x1=c(c为常数), 则xn+1=an c+(1-an ) 1-a 引理2[3] 对于二阶常系数线性差分方程xn+2=axn+1+bx,a,b为常数,若x1=m1,x2=m2(m1,m2为常数), 则xn+1=λ1n (m1λ2-m2)λ2-λ1+λ2n (m2-m1λ2)λ2-λ1 ,其中λ1、λ2是方程λ 2 -aλ-b=0的两根。证明令xn=Aλn代入(2)得:Aλn(λ2-aλ-b)=0,称方程λ2 -aλ-b=0为差分方程(2)的特征方程,且(2) 的解与特征方程的解有关系式:xn+1=c1λ1n+1 +c2λ2n+2 , 因给定初值x1=m1x2=m2" , 代入上式得:m1=c1λ1+c2λ2 m2=c1λ12+c2λ2 2 " 差分方程与概率计算 唐燕玉 (安庆师范学院学报编辑部,安徽安庆246011) 摘要:全文介绍了差分方程的概念,并给出了一阶差分方程xn+1=axn+b的通解与给定初始条件x1=c的特解,同时又给出了二阶差分方程xn+2=axn+1+bxn的通解与给定初始条件x1=m1,x2=m2的特解,并详细讨论了这两种差分方程在概率论中的应用。 关键词:概率;差分;差分方程;试验;全概公式中图分类号:O211 文献标识码:A 文章编号:1007-4260(2006)04-0091-03 收稿日期:2006-01-28 作者简介:唐燕玉(1951-),女,安徽枞阳人,安庆师范学院学报(自然科学版)主编。 安庆师范学院学报(自然科学版) JournalofAnqingTeachersCollege(NaturalScienceEdition) 2006年11月 Nov.2006第12卷第4期 Vol.12No.4

差分方程的解法分析及MATLAB实现(程序)

差分方程的解法分析及MATLAB 实现(程序) 摘自:张登奇,彭仕玉.差分方程的解法分析及其MATLAB 实现[J]. 湖南理工学院学报.2014(03) 引言 线性常系数差分方程是描述线性时不变离散时间系统的数学模型,求解差分方程是分析离散时间系统的重要内容.在《信号与系统》课程中介绍的求解方法主要有迭代法、时域经典法、双零法和变换域 法[1]. 1 迭代法 例1 已知离散系统的差分方程为)1(3 1)()2(81)1(43)(-+=-+--n x n x n y n y n y ,激励信号为)()4 3()(n u n x n =,初始状态为21)2(4)1(=-=-y y ,.求系统响应. 根据激励信号和初始状态,手工依次迭代可算出24 59)1(,25)0(==y y . 利用MATLAB 中的filter 函数实现迭代过程的m 程序如下: clc;clear;format compact; a=[1,-3/4,1/8],b=[1,1/3,0], %输入差分方程系数向量,不足补0对齐 n=0:10;xn=(3/4).^n, %输入激励信号 zx=[0,0],zy=[4,12], %输入初始状态 zi=filtic(b,a,zy,zx),%计算等效初始条件 [yn,zf]=filter(b,a,xn,zi),%迭代计算输出和后段等效初始条件 2 时域经典法 用时域经典法求解差分方程:先求齐次解;再将激励信号代入方程右端化简得自由项,根据自由项形 式求特解;然后根据边界条件求完全解[3].用时域经典法求解例1的基本步骤如下. (1)求齐次解.特征方程为081432=+-αα,可算出4 1 , 2121==αα.高阶特征根可用MATLAB 的roots 函数计算.齐次解为. 0 , )4 1()21()(21≥+=n C C n y n n h (2)求方程的特解.将)()4 3()(n u n x n =代入差分方程右端得自由项为 ?????≥?==-?+-1,)4 3(9130 ,1)1()43(31)()43(1n n n u n u n n n 当1≥n 时,特解可设为n p D n y )4 3()(=,代入差分方程求得213=D . (3)利用边界条件求完全解.当n =0时迭代求出25)0(=y ,当n ≥1时,完全解的形式为 ,)4 3(213 )41()21()(21n n n C C n y ?++=选择求完全解系数的边界条件可参考文[4]选)1(),0(-y y .根据边界条件求得35,31721=-=C C .注意完全解的表达式只适于特解成立的n 取值范围,其他点要用 )(n δ及其延迟表示,如果其值符合表达式则可合并处理.差分方程的完全解为

差分方程的基本知识(3)

差分方程模型的理论和方法 1、差分方程:差分方程反映的是关于离散变量的取值与变化规律。通过建立一个或几个离散变量取值所满足的平衡关系,从而建立差分方程。 差分方程就是针对要解决的目标,引入系统或过程中的离散变量,根据实际背景的规律、性质、平衡关系,建立离散变量所满足的平衡关系等式,从而建立差分方程。通过求出和分析方程的解,或者分析得到方程解的特别性质(平衡性、稳定性、渐近性、振动性、周期性等),从而把握这个离散变量的变化过程的规律,进一步再结合其他分析,得到原问题的解。 2、应用:差分方程模型有着广泛的应用。实际上,连续变量可以用离散变量来近似和逼近,从而微分方程模型就可以近似于某个差分方程模型。差分方程模型有着非常广泛的实际背景。在经济金融保险领域、生物种群的数量结构规律分析、疾病和病虫害的控制与防治、遗传规律的研究等许许多多的方面都有着非常重要的作用。可以这样讲,只要牵涉到关于变量的规律、性质,就可以适当地用差分方程模型来表现与分析求解。 3、差分方程建模:在实际建立差分方程模型时,往往要将变化过程进行划分,划分成若干时段,根据要解决问题的目标,对每个时段引入相应的变量或向量,然后通过适当假设,根据事物系统的实际变化规律和数量相互关系,建立每两个相邻时段或几个相邻时段或者相隔某几个时段的量之间的变化规律和运算关系(即用相应设定的变量进行四则运算或基本初等函数运算或取最运算等)等式(可以多个并且应当充分全面反映所有可能的关系),从而建立起差分方程。或者对事物系统进行划分,划分成若干子系统,在每个子系统中引入恰当的变量或向量,然后分析建立起子过程间的这种量的关系等式,从而建立起差分方程。在这里,过程时段或子系统的划分方式是非常非常重要的,应当结合已有的信息和分析条件,从多种可选方式中挑选易于分析、针对性强的划分,同时,对划分后的时段或子过程,引入哪些变量或向量都是至关重要的,要仔细分析、选择,尽量扩大对过程或系统的数量感知范围,包括对已有的、已知的若干量进行结合运算、取最运算等处理方式,目的是建立起简洁、深刻、易于求解分析的差分方程。在后面我们所举的实际例子中,这方面的内容应当重点体会。

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

第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 αβ''-=<

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

第七章 常微分方程与差分方程 常微分方程是高等数学中理论性和应用性都较强的一部分,是描述客观规律的一种重要方法,是处理物理、力学、几何等应用问题的一个重要工具,微分和积分的知识是研究微分方程的基础。微分方程作为考试的重点容,每年研究生考试均会考到。特别是微分方程的应用问题,既是重点,也是难点,在复习时必须有所突破。 【数学一大纲容】常微分方程的基本概念;变量可分离的方程;齐次方程;一阶线性方程;伯努利(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 的通解为____________。

时间序列分析讲义 第01章 差分方程

第一章 差分方程 差分方程是连续时间情形下微分方程的特例。差分方程及其求解是时间序列方法的基础,也是分析时间序列动态属性的基本方法。经济时间序列或者金融时间序列方法主要处理具有随机项的差分方程的求解问题,因此,确定性差分方程理论是我们首先需要了解的重要内容。 §1.1 一阶差分方程 假设利用变量t y 表示随着时间变量t 变化的某种事件的属性或者结构,则t y 便是在时间t 可以观测到的数据。假设t y 受到前期取值1-t y 和其他外生变量t w 的影响,并满足下述方程: t t t w y y ++=-110φφ (1.1) 在上述方程当中,由于t y 仅线性地依赖前一个时间间隔自身的取值1-t y ,因此称具有这种结构的方程为一阶线性差分方程。如果变量t w 是确定性变量,则此方程是确定性差分方程;如果变量t w 是随机变量,则此方程是随机差分方程。在下面的分析中,我们假设t w 是确定性变量。 例1.1 货币需求函数 假设实际货币余额、实际收入、银行储蓄利率和商业票据利率的对数变量分别表示为t m 、t I 、bt r 和ct r ,则可以估计出美国货币需求函数为: ct bt t t t r r I m m 019.0045.019.072.027.01--++=- 上述方程便是关于t m 的一阶线性差分方程。可以通过此方程的求解和结构分析,判断其他外生变量变化对货币需求的动态影响。 1.1.1 差分方程求解:递归替代法 差分方程求解就是将方程变量表示为外生变量及其初值的函数形式,可以通过以前的数据计算出方程变量的当前值。 由于方程结构对于每一个时间点都是成立的,因此可以将(1.1)表示为多个方程: 0=t :01100w y y ++=-φφ 1=t :10101w y y ++=φφ t t =:t t t w y y ++=-110φφ 依次进行叠代可以得到: 1011211010110101)()1()(w w y w w y y ++++=++++=--φφφφφφφφ 0111122113121102)1(w w w y y φφφφφφφ++++++=- i t i i t t i i t w y y ∑∑=-=++=0 111 1 0φφφφ (1.2) 上述表达式(1.2)便是差分方程(1.1)的解,可以通过代入方程进行验证。上述通过叠代将 t y 表示为前期变量和初始值的形式,从中可以看出t y 对这些变量取值的依赖性和动态变化 过程。 1.1. 2. 差分方程的动态分析:动态乘子(dynamic multiplier) 在差分方程的解当中,可以分析外生变量,例如0w 的变化对t 阶段以后的t y 的影响。假设初始值1-y 和t w w ,,1 不受到影响,则有:

差分方程求解

例题:已知差分方程51 (2)(1)()(+1)+0.5()66 x k x k x k r k r k +-++=,其中r (k )=1,k ≥0,x (0)=1, x (1)=2。 (1) 试由迭代法求其全解的前5项; (2) 分别由古典法求其零输入解、零状态解,以及全解; (3) 用Z 变换法求解差分方程。 解:注:解题过程中出现的下标“zi ”和“zs ”分别表示零输入条件和零状态条件。 1. 迭代法 题目中给出的条件仅仅是零输入初始条件,进行迭代求解时的初始条件应该是全解初始条件。 (1) 零输入初始条件 本题已给出零输入时的两个初始条件x zi (0)=1,x zi (1)=2。 (2) 零状态初始条件 取k =-2时,则51 (0)(1)(2)(1)0.5(2)66x x x r r --+-=-+-,得x zs (0)=0; 取k =-1 时,则51 (1)(0)(1)(0)0.5(1)66 x x x r r -+-=+-,求得x zs (1)=1。 (3) 全解初始条件 x (0)= x zi (0)+ x zs (0)=1; x (1)= x zi (1)+ x zs (1)=3。 (4) 根据求出的全解x (0)和x (1),利用迭代法求解 取k =0时,则51(2)(1)(0)(1)0.5(0)66x x x r r -+=+,求得23(2)6x =; 取k =1时,则51(3)(2)(1)(2)0.5(1)66x x x r r -+=+,求得151 (3)36x =; 取k =2时,则51(4)(3)(2)(3)0.5(2)66x x x r r -+=+,求得941 (4)216 x =。 2. 古典法 (1) 零输入解 令输入为零,则得齐次方程 51 (2)(1)()066 x k x k x k +-++= (a) 根据差分方程定义的算子()()n d x k x k n =+,可得它的特征方程251 066 d d -+= 求得特征根为: 112d = ,21 3 d =

差分方程的解法

第三节 差分方程常用解法与性质分析 1、常系数线性差分方程的解 方程)(...110n b x a x a x a n k k n k n =+++-++ ( 8) 其中k a a a ,...,,10为常数,称方程(8)为常系数线性方程。 又称方程0...110=+++-++n k k n k n x a x a x a (9) 为方程(8)对应的齐次方程。 如果(9)有形如 n n x λ=的解,带入方程中可得: 0 ...1110=++++--k k k k a a a a λλλ (10) 称方程(10)为方程(8)、(9)的特征方程。 显然,如果能求出(10)的根,则可以得到(9)的解。 基本结果如下: (1) 若(10)有k 个不同的实根,则(9)有通解: n k k n n n c c c x λλλ+++=...2211, (2) 若(10)有m 重根λ,则通解中有构成项: n m m n c n c c λ )...(121----+++

(3)若(10)有一对单复根 βαλi ±=,令:?ρλi e ±=, αβ?βαρarctan ,22=+=,则(9)的通解中有构成项: n c n c n n ?ρ?ρsin cos 21--+ (4) 若有m 重复根:βαλi ±=,φρλi e ±=,则(9)的通项中有成 项: n n c n c c n n c n c c n m m m m n m m ?ρ?ρsin )...(cos )...(1221121---++---+++++++ 综上所述,由于方程(10)恰有k 个根,从而构成方程 (9)的通解中必有k 个独立的任意常数。通解可记为:-n x 如果能得到方程(8)的一个特解:*n x ,则(8)必有通解: =n x -n x +* n x (11) (1) 的特解可通过待定系数法来确定。 例如:如果)(),()(n p n p b n b m m n =为n 的多项式,则当b 不是特征 根时,可设成形如)(n q b m n 形式的特解,其中)(n q m 为m 次多项式;如 果b 是r 重根时,可设特解:r n n b )(n q m ,将其代入(8)中确定出系 数即可。

差分方程基本概念和方法

差分方程基本概念和方法 考察定义在整数集上的函数,(),,2,1,0,1,2, n x f n n ==-- 函数()n x f n =在n 时刻的一阶差分定义为: 1(1)()n n n x x x f n f n ?+=-=+- 函数()n x f n =在n 时刻的二阶差分定义为一阶差分的差分: 21212n n n n n n x x x x x x ???+++=-=-+ 同理可依次定义k 阶差分 k n x ? 定义1.含有自变量n ,未知函数n x 以及n x 的差分2,, n n x x ??的函数方程, 称为常 差分方程,简称为差分方程。出现在差分方程中的差分的最高阶数,称为差分方 程的阶。 k 阶差分方程的一般形式为 (,,, ,)0k n n n F n x x x ??= 其中(,,,,)k n n n F n x x x ??为,,, k n n n n x x x ??的已知函数,且至少k n x ?要在式中出 现。 定义2.含有自变量n 和两个或两个以上函数值1,, n n x x +的函数方程,称为(常) 差分方程,出现在差分方程中的未知函数下标的最大差,称为差分方程的阶。 k 阶差分方程的一般形式为 1(,,, ,)0n n n k F n x x x ++= 其中1(,,,,)n n n k F n x x x ++为1,,, n n n k n x x x ++的已知函数,且n x 和n k x +要在式中一定 要出现。 定义3.如果将已知函数()n x n ?=代入上述差分方程,使其对0,1,2, n =成为恒 等式,则称()n x n ?=为差分方程的解。如果差分方程的解中含有k 个独立的任意

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

第七章 常微分方程与差分方程 常微分方程是高等数学中理论性和应用性都较强的一部分,是描述客观规律的一种重要方法,是处理物理、力学、几何等应用问题的一个重要工具,微分和积分的知识是研究微分方程的基础。微分方程作为考试的重点内容,每年研究生考试均会考到。特别是微分方程的应用问题,既是重点,也是难点,在复习时必须有所突破。 【数学一大纲内容】常微分方程的基本概念;变量可分离的方程;齐次方程;一阶线性方程;伯努利(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 的通解为____________。

常微分方程数值解法

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 。

相关文档
最新文档