计算方法微分方程

个人收集整理-ZQ

1 / 1 习题

试用三种方法导出线性二步方法

122+++=n n n hf y y

用展开法求三步四阶方法类，并确定三步四阶显式方法. 形如

∑=++=k i k n k j n j f h y 0

βα

地阶方法称为方法，试确定一个三步方法，并给出其截断误差主项. 试用显式法及改进地法

)],(),([211n n n n n n n hf y t f y t f h

y y +++=++

给出线性多步法

])13()3[(4)1(212n n n n n f f h

y y y +++=--++++αααα

为零稳定地条件，并证明该方法为零稳定时是二阶收敛地. 给出题()题中1=α时地公式地绝对稳定域.

指出方法

地相容阶，并给出由该方法以步长计算初值问题()地步骤. 试述刚性问题地基本特征，并给出级方法为稳定地条件. 设有???=='0

0)()

,(y x y y x f y ，试构造形如

)()(11011--++++=n n n n n f f h y y y ββα

地二阶方法，并推导其局部截断误差首项.

设有常微分方程初值问题???=='0

0)()

,(y x y y x f y 地

单步法)],(2),([3111+++++=n n n n n n y x f y x f h

y y ，证明该方法是无条件稳定地.

081数值计算方法—常微分方程(组)

科学计算—理论、方法及其基于MATLAB 的程序实现与分析微分方程（组）数值解法 §1 常微分方程初值问题的数值解法微分方程（组）是科学研究和工程应用中最常用的数学模型之一。如揭示质点运动规律的Newton 第二定律： ()()()?????'='==0 00022x t x x t x t F dt x d m （１）和刻画回路电流或电压变化规律的基尔霍夫回路定律等，但是，只有一些简单的和特殊的常微分方程及常微分方程组，可以求得用公式给出的所谓“解析解”或“公式解”，如一阶线性微分方程的初值问题： () ()0 0y y t f ay dt dy =+= （２）的解为： ()()()τττd f e y e t y t t a at ?-+=00 （３）但是，绝大多数在实际中遇到的常微分方程和常微分方程组得不到“解析解”，因此，基于如下的事实：

１、绝大多数的常微分方程和常微分方程组得不到（有限形式的）解析解；２、实际应用中往往只需要知道常微分方程（组）的解在（人们所关心的）某些点处的函数值（可以是满足一定精度要求的近似值）；如果只需要常微分方程（组）的解在某些点处的函数值，则没有必要非得通过求得公式解，然后再计算出函数值不可，事实上，我们可以采用下面将介绍的常微分方程（组）的初值问题的数值解法，就可以达到这一目的。一般的一阶常微分方程（组）的初值问题是指如下的一阶常微分方程（组）的定解问题： ()()0 00,y t y t t t y t F dt dy f =≤≤= （７）其中 ()()()()???? ?? ? ??=t y t y t y t y n 21 （８） ()()()()???? ?? ? ??=y t f y t f y t f y t F n ,,,,21 （９）常微分方程（组）的初值问题通常是对一动态过程（动态系统、动力系统）演化规律的描述，求解常微分方程（组）的初值问题就是要了解和掌握动态过程演化规律。 §１.1 常微分方程（组）的Cauch 问题数值解法概论

常微分方程在数学建模中的应用(免费版)

常微分方程在数学建模中的应用这里介绍几个典型的用微分方程建立数学模型的例子. 一、人口预测模型由于资源的有限性,当今世界各国都注意有计划地控制人口的增长,为了得到人口预测模型,必须首先搞清影响人口增长的因素,而影响人口增长的因素很多,如人口的自然出生率、人口的自然死亡率、人口的迁移、自然灾害、战争等诸多因素,如果一开始就把所有因素都考虑进去，则无从下手.因此,先把问题简化,建立比较粗糙的模型,再逐步修改,得到较完善的模型. 例1（马尔萨斯（Malthus ）模型）英国人口统计学家马尔萨斯（1766—1834）在担任牧师期间,查看了教堂100多年人口出生统计资料,发现人口出生率是一个常数,于1789年在《人口原理》一书中提出了闻名于世的马尔萨斯人口模型,他的基本假设是：在人口自然增长过程中,净相对增长（出生率与死亡率之差）是常数,即单位时间内人口的增长量与人口成正比,比例系数设为r ,在此假设下,推导并求解人口随时间变化的数学模型. 解设时刻t 的人口为)(t N ,把)(t N 当作连续、可微函数处理（因人口总数很大,可近似地这样处理,此乃离散变量连续化处理）,据马尔萨斯的假设,在t 到t t ?+时间段内,人口的增长量为 t t rN t N t t N ?=-?+)()()(, 并设0t t =时刻的人口为0N ,于是 ?????==．， 00)(d d N t N rN t N 这就是马尔萨斯人口模型,用分离变量法易求出其解为 )(00e )(t t r N t N -=, 此式表明人口以指数规律随时间无限增长. 模型检验：据估计1961年地球上的人口总数为9 1006.3?,而在以后7年中,人口总数以每年2%的速度增长,这样19610=t ,901006.3?=N ,02.0=r ,于是 ) 1961(02.09 e 1006.3)(-?=t t N . 这个公式非常准确地反映了在1700—1961年间世界人口总数.因为,这期间地球上的人口大约每35年翻一番,而上式断定34.6年增加一倍（请读者证明这一点）．但是,后来人们以美国人口为例,用马尔萨斯模型计算结果与人口资料比较,却发现有很大的差异,尤其是在用此模型预测较遥远的未来地球人口总数时,发现更令人不可思议的问题,如按此模型计算,到2670年,地球上将有36 000亿人口.如果地球表面全是陆地（事实上,地球表面还有80%被水覆盖）,我们也只得互相踩着肩膀站成两层了,这是非常荒谬的,因此,这一模型应该修改. 例2（逻辑Logistic 模型）马尔萨斯模型为什么不能预测未来的人口呢？这主要是地

常微分方程的实际应用

常微分方程的实际应用于萍摘要：常微分方程在当代数学中是极为重要的一个分支，它的实用价值很高，应用也很广泛，本文主要介绍常微分方程在几何、机械运动、电磁振荡方面的应用，并举例说明，体会常微分方程对解决实际问题的作用，在解决实际问题过程中通常是建立起实际问题的数学模型，也就是建立反映这个实际问题的微分方程，求解这个微分方程，用所得的数学结果解释实际问题，从而预测到某些物理过程的特定性质，以便达到能动地改造世界，解决实际问题的目的。关键字：常微分方程，几何，机械运动，电磁振荡，应用

Abstract: Nomal differential equation is an important part of math at it has a high practical value. This thesis shows the use in geometry, mechaics and electrothermal and makes some examples. Also, it summarizes the normal move of dealing with practical problems by the normal differential equation. Normal, we set up the maths matic model of the problem, solute the normal differentical equation make the use of the result to explain practical problems and make a forecast of some special character of physical process. Key: Normal differetial equation geometry mechanics electrothermal use

第十章-偏微分方程数值解法

第十章偏微分方程数值解法偏微分方程问题，其求解十分困难。除少数特殊情况外，绝大多数情况均难以求出精确解。因此，近似解法就显得更为重要。本章仅介绍求解各类典型偏微分方程定解问题的差分方法。 §1 差分方法的基本概念 1.1 几类偏微分方程的定解问题椭圆型方程：其最典型、最简单的形式是泊松（Poisson ）方程 ),(22 2 2y x f y u x u u =??+??=? 特别地，当0),(≡y x f 时，即为拉普拉斯（Laplace ）方程，又称为调和方程 22 22=??+??=?y u x u u Poisson 方程的第一边值问题为 ?? ???Ω ?=Γ=Ω∈=??+??Γ∈),(),(),() ,(),(22 22y x y x u y x y x f y u x u y x ?

其中 Ω为以Γ为边界的有界区域，Γ为分段光滑曲线，ΓΩY 称为定解区域，),(y x f ，),(y x ?分别为Ω，Γ上的已知连续函数。第二类和第三类边界条件可统一表示为 ),(),(y x u u y x ?α=??? ? ??+??Γ∈n 其中n 为边界Γ的外法线方向。当0=α时为第二类边界条件， 0≠α时为第三类边界条件。抛物型方程：其最简单的形式为一维热传导方程 2 20(0)u u a a t x ??-=>?? 方程可以有两种不同类型的定解问题：初值问题 ?? ???+∞ <<∞-=+∞<<-∞>=??-??x x x u x t x u a t u )()0,(,00 22? 初边值问题

2 212 00,0(,0)()0(0,)(),(,)()0u u a t T x l t x u x x x l u t g t u l t g t t T ????-=<<<

偏微分方程求解方法及其比较

偏微分方程求解方法及其比较发表时间：2008-12-11T09:32:01.530Z 来源：《科海故事博览科教创新》2008年第10期供稿作者：曹海洋吕淑娟王淑芬 [导读] 近些年来，无限维动力系统得到了很大的发展.随着对它研究的深入和计算能力的迅速提高，使得与之相关的数值研究越来越被人们关注.谱方法作为一种数值求解偏微分方程的方法，它具有无穷阶收敛性.因此，谱方法也就引起人们更多的关注. 摘要:近些年来，无限维动力系统得到了很大的发展.随着对它研究的深入和计算能力的迅速提高，使得与之相关的数值研究越来越被人们关注.谱方法作为一种数值求解偏微分方程的方法，它具有无穷阶收敛性.因此，谱方法也就引起人们更多的关注. 关键词：谱方法；偏微分；收敛；逼近； 1偏微分方程及其谱方法的介绍偏微分方程主要借助于未知函数及其导数来刻画客观世界的物理量的一般变化规律。理论上，对偏微分方程解法的研究已经有很长的历史了。最初的研究工作主要集中在物理，力学，几何学等方面的具体问题，其经典代表是波动方程，热传导方程和位势方程（调和方程）。通过对这些问题的研究，形成了至今仍然使用的有效方法，例如，分离变量法，fourier变换法等。早期的偏微分方程研究主要集中在理论上，而在实际操作中其研究方法和研究结果都难以得到广泛的应用。求解的主要方法为：有限差分法，有限元法，谱方法。谱方法起源于Ritz-Galerkin方法，它是以正交多项式（三角多项式，切比雪夫多项式，勒让得多项式等）作为基函数的Galerkin方法、Tau 方法或配置法，它们分别称为谱方法、Tau方法或拟谱方法（配点法），通称为谱方法。谱方法是以正交函数或固有函数为近似函数的计算方法。从函数近似角度看．谱方法可分为Fourier方法．Chebyshev或Legendre方法。前者适用于周期性问题，后两者适用于非周期性问题。而这些方法的基础就是建立空间基函数。下面介绍几种正交多项式各种节点的取值方法及权重。 1) Chebyshev-Gauss: 2) Chebyshev-Gauss-Radau: x0 =1, 3) Chebyshev-Gauss-Lobatto: x0 =1, xN =1, 4)Legendre-Gauss: xj 是的零点且 5） Legendre-Gauss-Radau: xj 是的N+1个零点且 6） Legendre-Gauss-Lobatto: x0=-1,xN=1其它N-1个点是的零点且下面介绍谱方法中最重要的Jacobi正交多项式其迭代公式为：其中： Jacobi正交多项式满足正交性：而Chebyshev多项式是令时Jacobi多项式的特殊形式，另外Legendre多项式是令时Jacobi多项式的特殊形式。 2 几种典型的谱方法谱方法是以正交函数或固有函数为近似函数的计算方法。谱近似可以分为函数近似和方程近似两种近似方式。从函数近似角度看．谱方法可分为Fourier方法．Chebyshev或Legendre方法。前者适用于周期性问题，后两者适用于非周期性问题。从方程近似角度看，谱方法可分为在物理空间离散求解的Collocation法、在谱空间进行离散求解的Galerkin法，以及先在物理空间离散求积，再变换到谱空间求解的Pseudo-spectral法。Collocation法适用于非线性问题．Galerkin法适用于线性问题，而Pseudo-spectral法适用于展开方程时的非线性项的处理。谱方法的特点是对光滑函数指数性逼近的谱精度；以较少的网格点得到较高的精度；无相位误差；适合多尺度的波动性问题；计算精度高于其他方法。快速傅立叶变化的提出大大促进了谱方法的发展，迄今已有各种的谱方法计算格式被提出．并被应用于天文学、电磁学、地理学等各种问题的计算。下面介绍一下应用于各个区域的几种谱方法： 1）以Fourier谱方法为例介绍谱方法解方程的主要过程以一阶波动方程为例：其中u(x,t)为方程的解，L是包含u和u关于空间变量的导数的算子，除了方程以有初始条件和适当的边界条件。故可设其中为试探空间的基函数，ak(t)为展开系数，对于傅立叶谱方法中的共轭有：其中从而利用其正交性和周期性可以减少工作量，另外再结合边界条件就可以求出来。 2) Galerkin方法是谱方法中十分经典的解偏微分方程的方法，但还有其局限性，而利用Hermite谱方法中依赖时间的权函数对经典的Galerkin方法进行拓展后的新的方法能适用范围扩大了很多。它能很好的应用在微分方程最优控制问题有限元方法的分析中，并且如果能够灵活运用利用Chebyshev方法、Galerkin方法和配置方法，则会形成更强的计算方法。如将Tau方法的思想成功地应用于奇数阶微分方程Petrov－Galerkin谱方法。 3)在无界区域上谱方法和拟谱方法发展了以Hermite函数和Laguerre函数为基函数的正交逼近和插值理论，在这些结果的基础上发展了全空间和半空间上数理方程的谱方法和拟谱方法，从而形成一种新的能更好解决误解区域问题的方法，此种方法被很好的应用于统计物理、量子力学和流体力学中。 4) 我们利用非一致带权Sobolev空间中的Jacobi多项式正交逼近和Jacobi-Gauss型插值理论，提出以Jacobi多项式为基函数的Jacobi谱方法和拟谱方法用来解决一些奇异问题和计算某些特定的无界区域问题。 5）有限谱方法是基于有限点、有限项的局域谱方法。这种方法要求近似函数应具有等同隔网格和非周期性的性质。有限谱方法分为基于非

偏微分方程数值解法

一、问题用有限元方法求下面方程的数值解 2 u u u f t ?-?+=? in (]0,T Ω? 0u = on []0,T ?Ω? ()00,u x u = in Ω 二、问题分析第一步利用Green 公式，求出方程的变分形式变分形式为：求()()21 00,;u L T H ∈Ω，使得 ()())(2 ,,,,u v u v u v f v t ???+??+= ???? ()10v H ?∈Ω （*）以及 ()00,u x u =. 第二步对空间进行离散，得出半离散格式对区域Ω进行剖分，构造节点基函数，得出有限元子空间：()12,,,h NG V span ???=???，则（*）的Galerkin 逼近为： []0,t T ?∈，求()()1 0,h h u t x V H ∈?Ω，使得 ()()()()() () )(2 ,,,,h h h h h h h d u t v u t v u t v f v dt +??+= h h v V ?∈ （**）以及()0,0h h u u =，0,h u 为初始条件0u 在h V 中的逼近，设0,h u 为0u 在h V 中的插值. 则0t ?≥，有()()1 N G h i i i u t t ξ? == ∑，0,h u =01 N G i i i ξ?=∑，代人（**）即可得到一常微分方程组. 第三步进一步对时间进行离散，得到全离散的逼近格式对 du dt 用差分格式.为此把[]0,T 等分为n 个小区间[]1,i i t t -，其长度1i i T t t t n -?=-= ,n t T =. 这样把求i t 时刻的近似记为i h u ，0 h u 是0u 的近似.这里对（**）采用向后的欧拉格式，即 ()()() () )(2 11 11 1 ,,,,i i i i h h h h h h h i h u u v u v u v f v t ++++-+??+ = ? h h v V ?∈ (***) i=0,1,2…,n-1. 0 h u =0,h u 由于向后欧拉格式为隐式格式且含有非线性项，故相邻两时间步之间采用牛顿迭代，即：

常微分方程在数学建模中的应用.

微分方程应用 1 引言常微分方程的形成与发展和很多学科有着密切的联系，例如力学、天文学、物理学等.数学的其他分支的快速发展，产生出很多新兴学科，这些新兴学科的产生都对常微分方程的发展有着深刻的影响，而且当前计算机的快速发展更是为常微分方程的应用及理论研究提供了非常有力的工具. 数学解决实际问题就必须建立模型，而数学建模就是把数学语言描述实际现象的过程.利用数学去解决各类实际问题时，建立数学模型是十分重要的一步，但是也是最困难的一步.建立数学模型的过程，是把错综复杂的实际问题简化、抽象为合理的数学结构的过程.要通过大量调查、收集相关数据资料，观察和研究实际对象的固有特征和内在规律，抓住问题的主要矛盾，建立起反映实际问题的数量关系，然后利用数学的理论和方法去分析和解决问题. 因此本文先简要介绍了如何建立微分方程模型，并通过具体的实例来简单地介绍了微分方程在数学建模中的应用. 2 数学模型简介通常我们把现实问题的一个模拟称为模型．如交通图、地质图、航空模型和建筑模型等．利用字母、数学及其它数学符号建立起来的等式或不等式以及图表、图象、框图等来模拟现实的模型称为数学模型．数学模型在实际生活中经常碰到，如求不规则图形的面积，可建立定积分的数学模型，求变化率的问题可建立导数模型，统计学中抽样调查，买彩票中奖的概率问题等等．学会建立数学模型对解决实际生活问题会有很大的帮助．建立数学模型是沟通摆在面前的实际问题与数学工具之间联系的一座必不可少的桥梁．随着科学技术的进步，特别是电子计算机技术的迅速发展，数学已经渗透到从自然科学技术到工农业生产建设，从经济生活到社会生活的各个领域．一般地说，当实际问题需要我们对所研究的现实对象提供分析、预报、决策、控制等方面的定量结果时，往往都离不开数学的应用，而建立数学模型则是这个过程的关键环节． 3 常微分方程模型 3．1 常微分方程的简介

偏微分方程的应用

偏微分方程在生物学上的应用刘富冲pb06007143 1偏微分方程的发展偏微分方程是反映有关的未知变量关于时间的导数和关于空间变量的导数之间制约关系的等式。许多领域中的数学模型都可以用偏微分方程来描述，物理学中的许多基本方程本身就是偏微分方程。早在微积分理论刚形成后不久，人们就开始用偏微分方程来描述、解释或预见各种自然现象，并将所得到的研究方法和研究成果运用于各门科学和工程技术中，不断地取得了显著的成效，显示了偏微分方程对于人类认识自然界基本规律的重要性。逐渐地，以物理、力学等各门科学中的实际问题为背景的偏微分方程的研究成为传统应用数学中的一个最主要的内容，它直接联系着众多自然现象和实际问题，不断地提出和产生出需要解决的新课题和新方法，不断地促进着许多相关数学分支(如泛函分析、微分几何、计算数学等)的发展，并从它们之中引进许多有力的解决问题的工具。偏微分方程已经成为当代数学中的一个重要的组成部分，是纯粹数学的许多分支和自然科学及工程技术等领域之间的一座重要的桥梁。在国外，对偏微分方程的应用发展是相当重视的。很多大学和研究单位都有应用偏微分方程的研究集体，并得到国家工业、科学部门及军方、航空航天等方面的大力资助。比如在国际上有重大影响的美国的Courant研究所、法国的信息与自动化国立研究所等都集中了相当多的偏微分方程的研究人员，并把数学模型、数学方法、应用软件及实际应用融为一体，在解决实际课题、推动学科发展及加速培养人才等方面都起了很大的作用。 2偏微分方程的应用在科技和经济发展中，很多重要的实际课题都需要求解偏微分方程，为相应的工程设计提供必要的数据，保证工程安全可靠且高效地完成任务。在很多的实际课题中，有不少课题(特别是国防课题)是不能或很难用工程试验的方法来进行研究的(一方面是危险系数大，另一方面是耗费大)，因此就需要尽可能地减少试验的次数或在试验前给出比较准确的预计。随着电子计算机的出现及计算技术的发展，电子计算机成为解决这些实际课题的重要工具。但是有效地利用电子计算机，必须具备如下先决条件：针对所考虑的实际问题建立合理的数学模型，而这些能精确描述问题的模型大都是通过偏微分方程给出的。对相应的偏微分方程模型进行定性的研究。根据所进行的定性研究，寻求或选择有效的求解方法。编制高效率的程序或建立相应的应用软件，利用电子计算机对实际问题进行模拟。因此，总体上来说，上述这些先决条件都属于偏微分方程应用的研究范围，这些问题解决的好坏直接影响到使用电子计算机所得结果的精确性及耗费的大小。如果解决得好，就会对整个问题的解决起到事半功倍的效果。到目前为止，偏微分方程已经在解决有关人口问题、传染病动力学、高速飞行、石油开发及城市交通等方面的实际课题中做出了重大的贡献。下面主要讲一下大家比较熟悉的人口问题及传染病动力学问题，详细阐述偏微分方程在解决实际问题中的应用。

常微分方程在高中物理中的应用

微分方程在高中物理中的应用高中阶段，我们经常会遇到一些需要定性分析的物理问题，其实如果我们应用高等数学的知识，可以把其中一些问题进行定量的分析。例如，质量为m 的物体从高度H 自由下落，所受阻力f 与速度v 成正比，g 为重力加速度这是我们平时常见的一类问题。但我们只知道速度V 最终会趋近于某一数值v0。下面我进行一下定量分析。根据题目所给信息，可列出动力学方程 mg-kv=ma ① a=dv/dt ② 结合①式可得mg-kv=mdv/dt 这里移项可得dt=mdv/(mg-kv)③ 两边同时积分便可的到 V=mg(ce*(-kt/m)+1)/k 又∵自由下落，可得t=0时v=.0 ∴v=mg(1-e*(-kt/m))/k ④ 由④式知，当t 趋近于正无穷时，e*(-kt/m)=0，此时v=mg/k ⑤ 若按照正常思路，当物体受力平衡时，mg=kv,此时也能得到⑤式的结论。而在高考中，更为常见的是在电磁场中的同类问题，我们不妨看一下下面这一道例题（2012·山东理综）如图所示，相距为L 的两条足够长的光滑平行金属导轨与水平面的夹角为θ，上端接有定值电阻，匀强磁场垂直于导轨平面，磁感应强度为B 。将质量为m 的导体棒由静止释放，当速度达到v 时开始匀速运动，此时对导体棒施加一平行于导轨向下的拉力，并保持拉力的功率为P ，导体棒最终以2v 的速度匀速运动。导体棒始终与导轨垂直且接触良好，不计导轨和导体棒的电阻，重力加速度为g ，下列选项正确的是 A ．P =2mg sin θ B ．P =3mg sin θ C ．当导体棒速度达到v /2时加速度为12 g sin θ D ．在速度达到2v 以后匀速运动的过程中，R 上产生的焦耳热等于拉力所做的功我们根据题目也可以列出动力学方程 Mgsin θ-B*2L*2V/R=ma ① a=dv/dt ② 同样可以解得v=(mgR sin θ/B*2L*2)(1-e*(-B*2L*2t/mR))③ 从③式可以看出当t 趋近于正无穷时，v=mgR sin θ/B*2L*2即B*2L*2v/R=mg sin θ转化而来。所以题目中所说当速度到达V 时开始匀速运动存在明显错误。应改为近似于做匀速直线运动。

第九章偏微分方程差分方法

170 第9章偏微分方程的差分方法含有偏导数的微分方程称为偏微分方程。由于变量的增多和区域的复杂性，求偏微分方程的精确解一般是不可能的，经常采用数值方法求方程的近似解。偏微分方程的数值方法种类较多，最常用的方法是差分方法。差分方法具有格式简单，程序易于实现，计算量小等优点，特别适合于规则区域上偏微分方程的近似求解。本章将以一些典型的偏微分方程为例，介绍差分方法的基本原理和具体实现方法。 9.1椭圆型方程边值问题的差分方法 9.1.1 差分方程的建立最典型的椭圆型方程是Poisson （泊松）方程 G y x y x f y u x u u ∈=??+??-≡?-),(),,()(2222 （9.1） G 是x ,y 平面上的有界区域，其边界Γ为分段光滑的闭曲线。当f (x ,y )≡0时，方程（9.1）称为Laplace(拉普拉斯）方程。椭圆型方程的定解条件主要有如下三种边界条件第一边值条件 ),(y x u α=Γ （9.2）第二边值条件 ),(y x n u β=??Γ （9.3）第三边值条件 ),()( y x ku n u γ=+??Γ （9.4）这里，n 表示Γ上单位外法向，α(x,y ),β(x,y ),γ(x,y )和k (x,y )都是已知的函数，k (x,y )≥0。满足方程（9.1）和上述三种边值条件之一的光滑函数u (x ,y )称为椭圆型方程边值问题的解。用差分方法求解偏微分方程，就是要求出精确解u (x ,y )在区域G 的一些离散节点（x i ，y i ）上的近似值u i ,j ≈(x i ，y i ）。差分方法的基本思想是，对求解区域G 做网格剖分，将偏微分方程在网格节点上离散化，导出精确解在网格节点上近似值所满足的差分方程，最终通过求解差分方程，通常为一个线性方程组，得到精确解在离散节点上的近似值。设G ={0

偏微分方程数值解法答案

1. 课本2p 有证明 2. 课本812,p p 有说明 3. 课本1520,p p 有说明 4. Rit2法，设n u 是u 的n 维子空间，12,...n ???是n u 的一组基底，n u 中的任一元素n u 可表为1n n i i i u c ?==∑ ，则,11 11()(,)(,)(,)(,)22j n n n n n n i j i j j i j j J u a u u f u a c c c f ???=== -=-∑∑是12,...n c c c 的二次函数，(,)(,)i j j i a a ????=，令 () 0n j J u c ?=?，从而得到12,...n c c c 满足1 (,)(,),1,2...n i j i j i a c f j n ???===∑，通过解线性方程组，求的i c ，代入1 n n i i i u c ?==∑，从而得到近似解n u 的过程称为Rit2法简而言之，Rit2法：为得到偏微分方程的有穷维解，构造了一个近似解，1 n n i i i u c ?== ∑，利用,11 11()(,)(,)(,)(,)22j n n n n n n i j i j j i j j J u a u u f u a c c c f ???===-=-∑∑确定i c ，求得近似解n u 的过程 Galerkin 法：为求得1 n n i i i u c ? == ∑形式的近似解，在系数i c 使n u 关于n V u ∈，满足(,)(,) n a u V f V =，对任意 n V u ∈或（取 ,1j V j n ?=≤≤） 1 (,)(,),1,2...n i j i j i a c f j n ???===∑的情况下确定i c ，从而得到近似解1 n n i i i u c ?==∑的过程称 Galerkin 法为 Rit2-Galerkin 法方程： 1 (,)(,)n i j i j i a c f ???==∑ 5. 有限元法：将偏微分方程转化为变分形式，选定单元的形状，对求解域作剖分，进而构造基函数或单元形状函数，形成有限元空间，将偏微分方程转化成了有限元方程，利用有效的有限元方程的解法，给出偏微分方程近似解的过程称为有限元法。 6. 解：对求解区间进行网格剖分，节点01......i n a x x x x b =<<<<=得到相邻节点1,i i x x -

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

郑州大学研究生课程（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 =′

数值分析_第五章_常微分方程数值解法

图５畅２令珔h ＝h λ，则y n ＋１＝１＋珔 h ＋１２珔h ２＋１６珔h ３＋１２４珔 h ４y n ．由此可知，绝对稳定性区域在珔h ＝h λ复平面上满足｜１＋珔 h ＋１２珔h ２＋１６珔h ３＋１２４珔h ４｜≤１的区域，也就是由曲线１＋珔h ＋１２珔h ２＋１６珔h ３＋１２４珔h ４＝e i θ 所围成的区域．如图５畅２所示．例22　用Euler 法求解 y ′＝－５y ＋x ，y （x ０）＝y ０，　x ０≤x ≤X ．从绝对稳定性考虑，对步长h 有何限制？解　对于模型方程y ′＝λy （λ＜０为实数）这里λ＝抄f 抄y ＝－５．由｜１＋h λ｜＝｜１－５h ｜＜１得到对h 的限制为：０＜h ＜０畅４．四、习题１畅取步长h ＝０畅２，用Euler 法解初值问题 y ′＝－y －x y ２， y （０）＝１．　（０≤x ≤０畅６），２畅用梯形公式解初值问题 y ′＝８－３y ，　（１≤x ≤２），

取步长h＝０畅２，小数点后至少保留５位．３畅用改进的Euler公式计算初值问题 y′＝１x y－１x y２， y（１）＝０畅５，　１＜x＜１畅５，取步长h＝０畅１，并与精确解y（x）＝ x １＋x比较．４畅写出用梯形格式的迭代算法求解初值问题 y′＋y＝０， y（０）＝１的计算公式，取步长h＝０畅１，并求y（０畅２）的近似值，要求迭代误差不超过１０－５．５畅写出用四阶经典Runge唱Kutta法求解初值问题 y′＝８－３y， y（０）＝２的计算公式，取步长h＝０畅２，并计算y（０畅４）的近似值，小数点后至少保留４位．６畅证明公式 y n＋１＝y n＋h９（２K１＋３K２＋４K３）． K１＝f（x n，y n）， K２＝f x n＋h２，y n＋h２K１， K３＝f x n＋３４h，y n＋３４h K２，至少是三阶方法．７畅试构造形如 y n＋１＝α（y n＋y n－１）＋h（β０f n＋β１f n－１）

二阶常微分方程的解法及其应用.

目录 1 引言 (1) 2 二阶常系数常微分方程的几种解法 (1) 2.1 特征方程法 (1) 2.1.1 特征根是两个实根的情形 (2) 2.1.2 特征根有重根的情形 (2) 2.2 常数变异法 (4) 2.3 拉普拉斯变化法 (5) 3 常微分方程的简单应用 (6) 3.1 特征方程法 (7) 3.2 常数变异法 (9) 3.3 拉普拉斯变化法 (10) 4 总结及意义 (11) 参考文献 (12)

二阶常微分方程的解法及其应用摘要：本文通过对特征方程法、常数变易法、拉普拉斯变换法这三种二阶常系数常微分方程解法进行介绍，特别是其中的特征方程法分为特征根是两个实根的情形和特征根有重根的情形这两种情况，分别使用特征值法、常数变异法以及拉普拉斯变换法来求动力学方程，现今对于二阶常微分方程解法的研究已经取得了不少成就，尤其在二阶常系数线性微分方程的求解问题方面卓有成效。应用常微分方程理论已经取得了很大的成就，但是，它的现有理论也还远远不能满足需要，还有待于进一步的发展，使这门学科的理论更加完善。关键词：二阶常微分方程;特征分析法；常数变异法；拉普拉斯变换

METHODS FOR TWO ORDER ORDINARY DIFFERENTIAL EQUATION AND ITS APPLICATION Abstract:This paper introduces the solution of the characteristic equation method, the method of variation of parameters, the Laplasse transform method the three kind of two order ordinary differential equations with constant coefficients, especially the characteristic equation method which is characteristic of the root is the two of two real roots and characteristics of root root, branch and don't use eigenvalue method, method of variation of constants and Laplasse transform method to obtain the dynamic equation, the current studies on solution of ordinary differential equations of order two has made many achievements, especially in the aspect of solving the problem of two order linear differential equation with constant coefficients very fruitful. Application of the theory of ordinary differential equations has made great achievements, however, the existing theory it is still far from meeting the need, needs further development, to make the discipline theory more perfect. Keywords:second ord er ordinary differential equation; Characteristic analysis; constant variation method; Laplasse transform 1 引言数学发展的历史告诉我们，300年来数学分析是数学的首要分支，而微分方程

微分方程在经济学中的应用

第四节微分方程在经济学中的应用微分方程在经济学中有着广泛的应用，有关经济量的变化、变化率问题常转化为微分方程的定解问题．一般应先根据某个经济法则或某种经济假说建立一个数学模型，即以所研究的经济量为未知函数，时间t 为自变量的微分方程模型，然后求解微分方程，通过求得的解来解释相应的经济量的意义或规律，最后作出预测或决策，下面介绍微分方程在经济学中的几个简单应用．一、供需均衡的价格调整模型在完全竞争的市场条件下，商品的价格由市场的供求关系决定，或者说，某商品的供给量S 及需求量D 与该商品的价格有关，为简单起见，假设供给函数与需求函数分别为 S =a 1+b 1P , D =a -bP , 其中a 1,b 1,a ,b 均为常数，且b 1＞0,b ＞0；P 为实际价格．供需均衡的静态模型为 ?? ???=+=-=).()(,,11P S P D P b a S bP a D 显然，静态模型的均衡价格为 P e =1 1b b a a +-．对产量不能轻易扩大，其生产周期相对较长的情况下的商品，瓦尔拉（Ｗalras ）假设：超额需求［D (P )-S (P )］为正时，未被满足的买方愿出高价，供不应求的卖方将提价，因而价格上涨；反之，价格下跌，因此，t 时刻价格的变化率与超额需求D -S 成正比，即 t P d d =k (D -S )，于是瓦尔拉假设下的动态模型为 ??? ????-=+=-=)].()([), (),(11P S P D k t P t P b a S t bP a D d d 整理上述模型得 t P d d =λ(P e -P ), 其中λ=k (b +b 1)＞0,这个方程的通解为 P (t )=P e +C e -λt ．假设初始价格为P (0)=P 0,代入上式得，C =P 0-P e ,于是动态价格调整模型的解为 P (t )=P e +(P 0-P e )·e -λt ，由于λ＞0，故 lim ()t P t →+∞=P e ．这表明，随着时间的不断延续，实际价格P (t )将逐渐趋于均衡价格P e ．二、索洛(Solow)新古典经济增长模型

求解偏微分方程三种数值方法

数值模拟偏微分方程的三种方法介绍（有限差分方法、有限元方法、有限体积方法） I．三者简介有限差分方法（Finite Difference Methods）是数值模拟偏微分方程最早采用的方法，至今仍被广泛使用。该方法包括区域剖分和差商代替导数两个步骤。首先将求解区域划分为差分网格，用有限个网格节点代替连续的求解区域。其次，利用Taylor级数展开等方法将偏微分方程中的导数项在网格节点上用函数值的差商代替进行离散，从而建立以网格节点上的值为未知量的代数方程组。该方法是一种直接将微分问题变为代数问题的近似数值解法，数学概念直观，表达简单，是发展较早且十分成熟的数值方法。差商代替导数后的格式称为有限差分格式，从格式的精度来考虑，有一阶格式、二阶格式和高阶格式。从差分的空间离散形式来考虑，有中心格式和迎风格式。对于瞬态方程，考虑时间方向的离散，有显格式、隐格式、交替显隐格式等。目前常见的差分格式，主要是以上几种格式的组合，不同的组合构成不同的差分格式。差分方法主要适用于结构网格，网格的大小一般根据问题模型和Courant 稳定条件来决定。有限元方法（Finite Element Methods）的基础是虚位移原理和分片多项式插值。该方法的构造过程包括以下三个步骤。首先，利用虚位移原理得到偏微分方程的弱形式，将计算区域划分为有限个互不重叠的单元（三角形、四边形、四面体、六面体等），在每个单元上选择合适的节点作为求解函数的插值点，将偏微分方程中的变量改写成由各变量或其导数的节点值与所选用的分片插值基函数组成的线性表达式，得到微分方程的离散形式。利用插值函数的局部支集性质及数值积分可以得到未知量的代数方程组。有限元方法有较完善的理论基础，具有求解区域灵活（复杂区域）、单元类型灵活（适于结构网格和非结构网格）、程序代码通用（数值模拟软件多数基于有限元方法）等特点。有限元方法最早应用于结构力学，随着计算机的发展已经渗透到计算物理、流体力学与电磁学等各个数值模拟领域。