偏微分方程分类

偏微分方程数值解法试题与答案

一．填空（1553=?分） 1．若步长趋于零时，差分方程的截断误差0→lm R ，则差分方程的解lm U 趋近于微分方 程的解lm u . 此结论_______（错或对）； 2．一阶Sobolev 空间{} )(,,),()(21 Ω∈''=ΩL f f f y x f H y x 关于内积=1),( g f _____________________是Hilbert 空间； 3．对非线性（变系数）差分格式，常用 _______系数法讨论差分格式的_______稳定性； 4．写出3 x y =在区间]2,1[上的两个一阶广义导数：_________________________________， ________________________________________； 5．隐式差分格式关于初值是无条件稳定的. 此结论_______（错或对）。 二．（13分）设有椭圆型方程边值问题 用1.0=h 作正方形网格剖分 。 （1）用五点菱形差分格式将微分方程在内点离散化； （2）用截断误差为)(2 h O 的差分法将第三边界条件离散化； （3）整理后的差分方程组为 三．（12）给定初值问题 x u t u ??=?? ， ()10,+=x x u 取时间步长1.0=τ，空间步长2.0=h 。试合理选用一阶偏心差分格式（最简显格式）, 并以此格式求出解函数),(t x u 在2.0,2.0=-=t x 处的近似值。 1．所选用的差分格式是： 2．计算所求近似值： 四．（12分）试讨论差分方程 ()h a h a r u u r u u k l k l k l k l ττ + - = -+=++++11,111 1 逼近微分方程 0=??+??x u a t u 的截断误差阶R 。 思路一：将r 带入到原式，展开后可得格式是在点（l+1/2,k+1/2）展开的。 思路二：差分格式的用到的四个点刚好是矩形区域的四个顶点，可由此构造中心点的差分格 式。

(完整版)偏微分方程的MATLAB解法

引言 偏微分方程定解问题有着广泛的应用背景。人们用偏微分方程来描述、解释或者预见各种自然现象，并用于科学和工程技术的各个领域fll。然而，对于广大应用工作者来说，从偏微分方程模型出发，使用有限元法或有限差分法求解都要耗费很大的工作量，才能得到数值解。现在，MATLAB PDEToolbox已实现对于空间二维问题高速、准确的求解过程。 偏微分方程 如果一个微分方程中出现的未知函数只含一个自变量，这个方程叫做常微分方程，也简称微分方程；如果一个微分方程中出现多元函数的偏导数，或者说如果未知函数和几个变量有关，而且方程中出现未知函数对几个变量的导数，那么这种微分方程就是偏微分方程。 常用的方法有变分法和有限差分法。变分法是把定解问题转化成变分问题，再求变分问题的近似解；有限差分法是把定解问题转化成代数方程，然后用计算机进行计算；还有一种更有意义的模拟法，它用另一个物理的问题实验研究来代替所研究某个物理问题的定解。虽然物理现象本质不同，但是抽象地表示在数学上是同一个定解问题，如研究某个不规则形状的物体里的稳定温度分布问题，由于求解比较困难，可作相应的静电场或稳恒电流场实验研究，测定场中各处的电势，从而也解决了所研究的稳定温度场中的温度分布问题。 随着物理科学所研究的现象在广度和深度两方面的扩展，偏微分方程的应用范围更广泛。从数学自身的角度看，偏微分方程的求解促使数学在函数论、变分法、级数展开、常微分方程、代数、微分几何等各方面进行发展。从这个角度说，偏微分方程变成了数学的中心。

一、MATLAB方法简介及应用 1.1 MATLAB简介 MATLAB是美国MathWorks公司出品的商业数学软件，用于算法开发、数据可视化、数据分析以及数值计算的高级技术计算语言和交互式环境，主要包括MATLAB和Simulink两大部分。 1.2 Matlab主要功能 数值分析 数值和符号计算 工程与科学绘图 控制系统的设计与仿真 数字图像处理 数字信号处理 通讯系统设计与仿真 财务与金融工程 1.3 优势特点 1) 高效的数值计算及符号计算功能，能使用户从繁杂的数学运算分析中解脱出来； 2) 具有完备的图形处理功能，实现计算结果和编程的可视化； 3) 友好的用户界面及接近数学表达式的自然化语言，使学者易于学习和掌握； 4) 功能丰富的应用工具箱(如信号处理工具箱、通信工具箱等) ，

Maab求解微分方程组及偏微分方程组

第四讲 Matlab 求解微分方程（组） 理论介绍：Matlab 求解微分方程（组）命令 求解实例：Matlab 求解微分方程（组）实例 实际应用问题通过数学建模所归纳得到的方程，绝大多数都是微分方程，真正能得到代数方程的机会很少.另一方面，能够求解的微分方程也是十分有限的，特别是高阶方程和偏微分方程（组）.这就要求我们必须研究微分方程（组）的解法：解析解法和数值解法. 一．相关函数、命令及简介 1.在Matlab 中，用大写字母D 表示导数，Dy 表示y 关于自变量的一阶导数，D2y 表示y 关于自变量的二阶导数，依此类推.函数dsolve 用来解决常微分方程（组）的求解问题，调用格式为： X=dsolve(‘eqn1’,’eqn2’,…) 函数dsolve 用来解符号常微分方程、方程组，如果没有初始条件，则求出通解，如果有初始条件，则求出特解. 注意，系统缺省的自变量为t 2.函数dsolve 求解的是常微分方程的精确解法，也称为常微分方程的符号解.但是，有大量的常微分方程虽然从理论上讲，其解是存在的，但我们却无法求出其解析解，此时，我们需要寻求方程的数值解，在求常微分方程数值解方面，MATLAB 具有丰富的函数，我们将其统称为solver ，其一般格式为： [T,Y]=solver(odefun,tspan,y0) 说明：(1)solver 为命令ode45、ode23、ode113、ode15s 、ode23s 、ode23t 、ode23tb 、ode15i 之一. (2)odefun 是显示微分方程'(,)y f t y =在积分区间tspan 0[,]f t t =上从0t 到f t 用初始条件0y 求解. (3)如果要获得微分方程问题在其他指定时间点012,,,,f t t t t L 上的解，则令tspan 012[,,,]f t t t t =L (要求是单调的). (4)因为没有一种算法可以有效的解决所有的ODE 问题，为此，Matlab 提供了多种求解器solver ，对于不同的ODE 问题，采用不同的solver.

偏微分方程数值解法

一、 问题 用有限元方法求下面方程的数值解 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 由于向后欧拉格式为隐式格式且含有非线性项，故相邻两时间步之间采用牛顿迭代，即：

各种类型的微分方程及其相应解法教程文件

各种类型的微分方程及其相应解法 专业班级：交土01班 姓名：高云 学号：1201110102 微分方程的类型有很多种，解题时先判断微分方程是哪种类型，可以帮助我们更快解题，所以我们有必要归纳整理一下各类型（主要是一阶和二阶）的微分方程及其相应解法。 一、一阶微分方程的解法 1.可分离变量的方程 dx x f dy y g )()(=,或)()(y g x f dx dy = 其特点是可以把变量x 和y 只分别在等式的两边，解法关键是把变量分离后两边积分。 例1.求微分方程ydy dx y xydy dx +=+2的通解. 解 先合并dx 及dy 的各项,得dx y dy x y )1()1(2-=- 设,01,012≠-≠-x y 分离变量得 dx x dy y y 1112-=- 两端积分??-=-dx x dy y y 1112得 ||ln |1|ln |1|ln 2 112C x y +-=- 于是 2212)1(1-±=-x C y 记,21C C ±=则得到题设方程的通解 .)1(122-=-x C y 2.齐次方程 （1）)(x y f dx dy = （2） )(c by ax f dx dy ++=（a ，b 均不等于0） 例2求解微分方程.2222xy y dy y xy x dx -=+- 解 原方程变形为=+--=2222y xy x xy y dx dy ,1222?? ? ??+--??? ??x y x y x y x y 令,x y u =则,dx du x u dx dy +=方程化为,1222u u u u dx du x u +--=+ 分离变量得?? ????-+--??? ??--112212121u u u u ,x dx du = 两边积分得 ,ln ln ln 2 1)2ln(23)1ln(C x u u u +=----

偏微分方程数值解法

“十二五”国家重点图书出版规划项目 信息与计算科学丛书 67 偏微分方程数值解法 陈艳萍鲁祖亮刘利斌编著

内 容 简 介 本书试图用较少的篇幅描述偏微分方程的几种数值方法. 主要内容包括：Sobolev空间初步, 椭圆边值问题的变分问题, 椭圆问题的有限差分方法, 抛物型方程的有限差分方法, 双曲型方程的有限差分方法, 椭圆型方程的有限元方法, 抛物及双曲方程的有限元方法, 椭圆型方程的混合有限元方法, 谱方法等. 本书内容丰富, 深入浅出, 尽可能地用简单的方法来描述一些理论结果, 并根据作者对有限差分、有限元、混合有限元、谱方法的理解和研究生教学要求, 全面、客观地评价各种数值计算方法，并列举一些数值计算的例子, 阐述许多新的学术观点. 本书可作为高等学校数学系高年级本科生和研究生的教材或参考书, 也可作为计算数学工作者和从事科学与工程计算的科研人员的参考书. 图书在版编目(CIP)数据 偏微分方程数值解法/陈艳萍, 鲁祖亮, 刘利斌编著. —北京：科学出版社, 2015.1 (信息与计算科学丛书67) ISBN 978-7-03-000000-0 Ⅰ. ①偏… Ⅱ. ①陈… ②鲁… ③刘… Ⅲ. ① Ⅳ.① 中国版本图书馆CIP数据核字(2014) 第000000号 责任编辑: 王丽平／责任校对: 彭涛 责任印制: 肖钦／封面设计: 陈敬 出版 北京东黄城根北街16号 邮政编码: 100717 https://www.360docs.net/doc/6c18451071.html, 印刷 科学出版社发行 各地新华书店经销 * 2015年1月第一版开本: 720×1000 1/16 2015年1月第一次印刷印张: 14 字数: 280 000 定价: 88.00元 (如有印装质量问题, 我社负责调换)

偏微分方程组解法

偏微分方程组解法 某厚度为10cm 平壁原温度为20C ?，现其两侧面分别维持在20C ?和120C ?，试求经过8秒后平壁温度分布，并分析温度分布随时间的变化直至温度分布稳定为止。 22x t a t ??=??τ 式中a 为导温系数，/s m c 2；2=a 。 解： 模型转化为标准形式： 2 21x t t a ??=??τ 初始条件为： ()200,=x t 边界条件为： ()120,0=τt ，()20,1.0=τt 函数： pdefun.m %偏微分方程(一维动态传热) function [c,f,s]=pdefun(x,t,u,dudx) c=1/2e-4;f=dudx;s=0; icbun.m %偏微分方程初始条件(一维动态传热) function u0=icbun(x) u0=20; bcfun.m %偏微分方程边界条件(一维动态传热) function [pl,ql,pr,qr]=bcfun(xl,ul,xr,ur,t) pl=ul-120;ql=0;pr=ur-20;qr=0; 命令： x=linspace(0,10,20)*1e-2; t=linspace(0,15,16); sol=pdepe(0,pdefun,icfun,bcfun,x,t); mesh(x,t,sol(:,:,1)) %温度与时间和空间位置的关系图 %画1、2、4、6、8、15s 时刻温度分布图

plot(x,sol(2,:,1)) 1s时刻，（因为本题sol第一行为0时刻） hold on plot(x,sol(3,:,1)) plot(x,sol(5,:,1)) plot(x,sol(7,:,1)) plot(x,sol(9,:,1)) plot(x,sol(16,:,1)) 计算结果： %第8秒时温度分布 x sol(9,:,1) 经过8秒时的温度分布为： x/cm 0 0.5263 1.0526 1.5789 2.1053 2.6316 3.1579 t/C ?120.0000 112.5520 105.1653 97.8994 90.8100 83.9477 77.3562 x/cm 3.6842 4.2105 4.7368 5.2632 5.7895 6.3158 6.8421 t/C ?71.0714 65.1202 59.5200 54.2784 49.3930 44.8518 40.6338 x/cm 7.3684 7.8947 8.4211 8.9474 9.4737 10.0000 t/C ?36.7095 33.0419 29.5877 26.2982 23.1207 20.0000 或者求第8秒时，x=0,2,4,,6,8,10cm处的温度 [uout,duoutdx]=pdeval(0,x,sol(9,:,:),[0,2,4,6,8,10]*1e-2) 120.0000 92.2279 67.5007 47.5765 32.3511 20.0000

偏微分方程数值解法答案

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 -

偏微分方程数值解法试题与答案

x 1 ?若步长趋于零时，差分方程的截断误差 R m 0,则差分方程的解 U i m 趋近于微分方 程的解U m ?此结论 ________ (错或对)； 1 2.一 阶 Sobolev 空间 H ( ) f (x,y) f , f x , f y L ?() 关于内积(f,g )1 _____________________________________ 是Hilbert 空间； 3 ?对非线性(变系数)差分格式，常用 ____________ 系数法讨论差分格式的 ________ 稳定性； 4?写出y x 3在区间［1,2］上的两个一阶广义导数： ______________________________________ _____ ____ ______________ _ ____ ________ ; 5 ?隐式差分格式关于初值是无条件稳定的 ?此结论 ________ (错或对)。 (13分)设有椭圆型方程边值问题 0.1作正方形网格剖分 。 (1) 用五点菱形差分格式将微分方程在内点离散化； (2) 用截断误差为 O (h 2)的差分法将第三边界条件离散化; (3) 整理后的差分方程组为 U C 三.(12)给定初值问题 u x,0 x 1 取时间步长 0.1，空间步长h 0.2。试合理选用一阶偏心差分格式(最简显格式) 2 u ~2 x 2 u ~2 y 0 x 0.3 0.2 x 0.3 2y 1, — u n 2x y 0.2

并以此格式求出解函数u(x,t)在x 0.2,t 0.2处的近似值。 x

1.所选用的差分格式是： 2 .计算所求近似值： 1 a k 1 四.（12分）试讨论差分方程 u l 1 k k k 1 u | r u | 1 u | , r h a 1 h 逼近微分方程 u a u 0 t x 的截断误差阶R 。 思路一：将r 带入到原式，展开后可得格式是在点（ l+1/2,k+1/2 ）展开的。 思路二：差分格式的用到的四个点刚好是矩形区域的四个顶点，可由此构造中心点的差分格 式。 2 —2 ,考虑 Du Fort-Frankel 格式 X 试论证该格式是否总满足稳定性的 Von-Neumann 条件？ 六. （12分）（1 ）由Green 第一公式推导 Green 第二公式: （2） 对双调和方程边值问题 n 2 选择函数集合（空间）为： 推导相应的双线性泛函和线性泛函: A （u,v ） F （v ） 相应的虚功问题为： 极小位能问题为 七. （ 12分）设有常微分方程边值问题 y y f （x ） , a x b y a 1, y b 1 五.（12分） 对抛物型方程 U |k1 U |k 2 |k 1 (U |k1 U |k1) U |k 1 ) 2 (u)vdxdy G (u) u vdxdy :[v v u ]ds n f (x,y) (x,y) g 1(x , y), g 2(x, y) (x,y),

二阶线性偏微分方程的分类与小结

第六章 二阶线性偏微分方程的分类与小结 一 两个自变量的二阶线性方程 1 方程变换与特征方程 两个自变量的二阶线性偏微分方程总表示成 f cu u b u b u a u a u a y x yy xy xx =+++++212212112 ① 它关于未知函数u 及其一、二阶偏导数都是线性的，其中f u c b b a a a ，，，，，，，21221211都是自变量y x ,的已知函数，假设它们的一阶偏 导数在某平面区域D 内都连续，而且221211a a a ，，不全为0 。 设),(000y x M 是D 内给定的一点，考虑在0M 的领域内对方程进行简化。取自变量变换 ),(y x ξξ=,),(y x ηη= 其中它们具有二连续偏导数，而且在0M 处的雅可比行列式。 = ??),(),(y x ηξy x y x ηηξξ =x y y x ηξηξ- 根据隐函数存在定理，在0M 领域内存在逆变换, ),(ηξx x =,),(ηξy y = 因为 x x x u u u ηξξξ+=，y y y u u u ηξξξ+=

xx xx x x x x xx u u u u u u ηξηηξξηξηηξηξξ++++=222 yy yy y y y y yy u u u u u u ηξηηξξηξηηξηξξ++++=222 xy xy y x x y y x x x xy u u u u u u ηξηηηξηξξξηξηηξηξξ+++++=)( 将代入①使其变为 F Cu u B u B u A u A u A =+++++ηξηηξηξξ212212112 经过变换后，方程的阶数不会升高，由变换的可逆性，方程的阶数也不会降低，所以221211,,A A A 不全为0。并可验证 222112122211212))((x y y x a a a A A A ηξηξ--=- 这表明，在可逆变换下2 22112 12A A A -与22112 12 a a a -保持相同的正负号。 定理 在0M 的领域内，不为常数的函数),(y x ?是偏微分方程022*******=++y y x x a a a ????之解的充分必要条件是： C y x ≡),(?是常微分方程的 0)(2)(22212211=++dx a dxdy a dy a 通解。 2 方程的类型及其标准形式 根据以上结论简化方程的问题归结为寻求其特征曲线。为此将特征方程分解成两个方程： 11 22 11 2 12 12 a a a a a dx dy -+=，11 22 11 2 12 12 a a a a a dz dy --= (1) 若在0M 的邻域内022112 12>-a a a 时，方程可以化为

一维偏微分方程的pdepe(matlab)函数 解法

本文根据matlab帮助进行加工，根据matlab帮助上的例子，帮助更好的理解一维偏微分方程的pdepe函数解法，主要加工在于程序的注释上。 Examples Example 1.This example illustrates the straightforward formulation, computation, and plotting of the solution of a single PDE. This equation holds on an interval for times . The PDE satisfies the initial condition and boundary conditions It is convenient to use subfunctions to place all the functions required by pdepe in a single function. function pdex1 m = 0; x = linspace(0,1,20); %linspace(x1,x2,N)linspace是Matlab中的一个指令，用于产生x1,x2之间的N点行矢量。 %其中x1、x2、N分别为起始值、终止值、元素个数。若缺省N，默认点数为100 t = linspace(0,2,5); sol = pdepe(m,@pdex1pde,@pdex1ic,@pdex1bc,x,t);

% Extract the first solution component as u. u = sol(:,:,1); % A surface plot is often a good way to study a solution. surf(x,t,u) title('Numerical solution computed with 20 mesh points.') xlabel('Distance x') ylabel('Time t') % A solution profile can also be illuminating. figure plot(x,u(end,:)) title('Solution at t = 2') xlabel('Distance x') ylabel('u(x,2)') % -------------------------------------------------------------- function [c,f,s] = pdex1pde(x,t,u,DuDx) c = pi^2; f = DuDx; s = 0; % -------------------------------------------------------------- function u0 = pdex1ic(x) u0 = sin(pi*x); % -------------------------------------------------------------- function [pl,ql,pr,qr] = pdex1bc(xl,ul,xr,ur,t) pl = ul; ql = 0; pr = pi * exp(-t); qr = 1;

偏微分方程的数值解法

《偏微分方程数值解法》试题 (专业:凝聚态物理学号:2013201260 姓名:鄢建军) 1.考虑定解问题 （1）用迎风格式(P、45)求解 1,0 (,0) 0,0 t x u u x u x x += ? ? ≤ ? ? =? ?> ? ? 。 利用迎风格式编写Fortran程序语言,运行结果如下: Fig 1、迎风格式求解结果 （2）用Beam-Warming格式(P、51)求解。 利用Beam—Warming格式编写Fortran程序语言,运行结果如下 :

Fig 2、 Beam —Warming 格式求解结果 （3） 比较两种方法结果的异同。 将两种格式运行的结果绘制在一起,要求时间步长与空间步长在两种格式中都相同,运行结果如下图所示: Fig 3、 迎风格式与Beam-Warming 格式求解结果比较 从两种格式的运行结果来瞧,都存在边缘的误差现象,相比而言,Beam-Warming 格式的运行结果差一些。但就是理论上分析,迎风格式的截断误差为()h οτ+,而Beam-Warming 格式的截断误差为22()h h οττ++。稳定性上来分析,迎风格式的稳定性较好,要求1(/)a h λλτ≤=,Beam-Warming 格式的稳定性条件为2(/)a h λλτ≤=。 2. 考虑定解问题212 1110,04(,0)sin ,0(0,)(,)0u u a x l t t u x x x l l u t u l t π???-=<

偏微分方程数值解法

《偏微分方程数值解法》 课程设计 题目: 六点对称差分格式解热传导方程的初边 值问题 姓名: 王晓霜 学院: 理学院 专业: 信息与计算科学 班级: ０91１01２ 学号: ０９１10１218 指导老师：翟方曼 ２０12年12月14

日 一、题目 用六点对称差分格式计算如下热传导方程的初边值问题 222122,01,01(,0),01 (0,),(1,),01x t t u u x t t x u x e x u t e u t e t +???=<<<≤?????=≤≤??==≤≤??? 已知其精确解为 2(,)x t u x t e += 二、理论 1.考虑的问题 考虑一维模型热传导方程 （1.１) )(22x f x u a t u +??=??,T t ≤<0 其中a 为常数。)(x f 是给定的连续函数。（1.１)的定解问题分两类: 第一,初值问题（Cauch ｙ 问题)：求足够光滑的函数()t x u ,,满足方程（１.1）和初始条件： (1.２) ()()x x u ?=0,， ∞<<∞-x 第二,初边值问题（也称混合问题）:求足够光滑的函数()t x u ,，满足方程(１.１)和初始条件： ()13.1 ()()x x u ?=0,， l x l <<- 及边值条件 ()23.1 ()()0,,0==t l u t u ， T t ≤≤0 假定()x f 和()x ?在相应的区域光滑,并且于()0,0，()0,l 两点满足相容条件,则上述问题有唯一的充分光滑的解。 现在考虑边值问题（１.1),（1.３)的差分逼近 取 N l h = 为空间步长,M T =τ为时间步长,其中N ,M 是自然数,

matlab偏微分方程组求解

MATLAB学习（序列1）偏微分方程组的求解 ode23 解非刚性微分方程，低精度，使用Runge-Kutta法的二三阶算法。 ode45 解非刚性微分方程，中等精度，使用Runge-Kutta法的四五阶算法。 ode113 解非刚性微分方程，变精度变阶次Adams-Bashforth-Moulton PECE算法。 ode23t 解中等刚性微分方程，使用自由内插法的梯形法则。 ode15s 解刚性微分方程，使用可变阶次的数值微分（NDFs）算法。 ode23s 解刚性微分方程，低阶方法，使用修正的Rosenbrock公式。 ode23tb 解刚性微分方程，低阶方法，使用TR-BDF2方法，即Runger-Kutta公式的第一级采用梯形法则，第二级采用Gear法。 [t，YY]=solver('F',tspan,Yo 解算ODE初值问题的最简调用格式。 solver指上面的指令。 tspan=[0,30]; ％时域t的范围 y0=[1;0]; ％y（1）y(2的初始值 [tt,yy]=ode45(@DyDt,tspan,y0; plot(tt,yy(:,1,title('x(t' function ydot=DyDt(t,y ydot=[y(2; 2*(1-y(1^2*y(2-y(1] 刚性方程：刚性是指其Jacobian矩阵的特征值相差十分悬殊。在解的性态上表现为，其中一些解变化缓慢，另一些变化快，且相差较悬殊，这类方程常常称为刚性方程，又称为Stiff方程。 刚性方程和非刚性方程对解法中步长选择的要求不同。 刚性方程一般不适合由ode45这类函数求解，而应该采用ode15s等。 如果不能分辨是否是刚性方程，先试用ode45，再用ode15s。 [t，YY，Te，Ye,Ie] = solver('F',tspan,Yo,options,p1,p2,… 解算ODE初值问题的最完整调用格式。 为了能够解出方程，要用指令odeset确定求解的条件和要求。在MATLAB中，求解方程组的指令都有默认的求解的条件和要求（由结构数组options表示），但可以用odeset修改或重新建立，也可以用odeget去获取已有的“优化选项”的信息。指令odeset和odeget用法介绍如下： 语句格式如下： options=odeset(‘name1’,value1,’name2’,value2,…

《偏微分方程概述及运用matlab求解偏微分方程常见问题》.

北京航空航天大学 偏微分方程概述及运用matlab求解微分方 程求解常见问题 姓名徐敏 学号57000211 班级380911班 2011年6月

偏微分方程概述及运用matlab求解偏微分 方程常见问题 徐敏 摘要偏微分方程简介，matlab偏微分方程工具箱应用简介，用这个工具箱解方程的过程是：确定待解的偏微分方程；确定边界条件；确定方程所在域的几何形状；划分有限元；解方程 关键词MATLAB 偏微分方程程序 如果一个微分方程中出现的未知函数只含有一个自变量，这个方程叫做常微分方程，也简称微分方程：如果一个微分方程中出现多元函数的偏导数，或者说如果未知函数和几个变量有关，而且方程中出现未知函数对几个变量的导数，那么这种微分方程就是偏微分方程。 一，偏微分方程概述 偏微分方程是反映有关的未知变量关于时间的导数和关于空间变量的导数之间制约关系的等式。许多领域中的数学模型都可以用偏微分方程来描述，很多重要的物理、力学等学科的基本方程本身就是偏微分方程。早在微积分理论刚形成后不久，人们就开始用偏微分方程来描述、解释或预见各种自然现象，并将所得到的研究方法和研究成果运用于各门科学和工程技术中，不断地取得了显著的成效，显示了偏微分方程对于人类认识自然界基本规律的重要性。逐渐地，以物

理、力学等各门科学中的实际问题为背景的偏微分方程的研究成为传统应用数学中的一个最主要的内容，它直接联系着众多自然现象和实际问题，不断地提出和产生出需要解决的新课题和新方法，不断地促进着许多相关数学分支(如泛函分析、微分几何、计算数学等)的发展，并从它们之中引进许多有力的解决问题的工具。偏微分方程已经成为当代数学中的一个重要的组成部分，是纯粹数学的许多分支和自然科学及工程技术等领域之间的一座重要的桥梁。 在国外，对偏微分方程的应用发展是相当重视的。很多大学和研究单位都有应用偏微分方程的研究集体，并得到国家工业、科学部门及军方、航空航天等方面的大力资助。比如在国际上有重大影响的美国的Courant研究所、法国的信息与自动化国立研究所等都集中了相当多的偏微分方程的研究人员，并把数学模型、数学方法、应用软件及实际应用融为一体，在解决实际课题、推动学科发展及加速培养人才等方面都起了很大的作用。 在我国，偏微分方程的研究起步较晚。但解放后，在党和国家的大力号召和积极支持下，我国偏微分方程的研究工作发展比较迅速，涌现出一批在这一领域中做出杰出工作的数学家，如谷超豪院士、李大潜院士等，并在一些研究方向上达到了国际先进水平。但总体来说，偏微分方程的研究队伍的组织和水平、研究工作的广度和深度与世界先进水平相比还有很大的差距。因此，我们必须继续努力，大力加强应用偏微分方程的研究，逐步缩小与世界先进水平的差距 二，偏微分方程的内容

第二章 二阶线性偏微分方程的分类

第二章 二阶线性偏微分方程的分类 1．把下列方程化为标准形式： （1）02=+++++u cu bu au au au y x yy xy xx 解：因为 02 22112 12=?-=-a a a a a a 所以该方程是抛物型方程，其特征方程为 12 2 =-± =a a a a dx dy 。 它只有一族实的特征线 c x y =- 在这种情况下，我们设x y -=ξ，x =η（或令y =η，总之，此处η是与ξ无关的任一函数，当然宜取最简单的函数形式x =η或y =η）。 方法一：用抛物型方程的标准形式 ][12122 F Cu u B u B A +++- =ηξηηη 先算出： ? ??? ? ? ?? ? ? ?-====?+?+?+?+?=++++=?+-+?+?+?=++++==?+?+=++=b c C b c b a a a b b a a a B c b a a a b b a a a B a a a a a a a A y x yy xy xx y x yy xy xx y y x x 0F ,1010020 2 1)1(0020 2 002 2212212112 2122121112 221221122ηηηηηξξξξξηηηη ∴])[(1 u bu u c b a u +++--=ηξηη 即 01=+ + -+ u a u a b u a b c u ηξηη 方法二：应用特征方程，作自变量变换，求出 ??? ??=+-=+-=+--==+-= ,2 ,ξξηξξξηηξηξξηηηξξηξξξηξu u u u u u u u u u u u u u u u u u yy xy xx y x 代入原方程得，0)(=++-+u bu u b c au ηξξη

偏微分方程定解问题

第一章 偏微分方程定解问题 引言：在研究、探索自然科学和工程技术中，经常遇到各种微分方程。 如 牛顿定律 22d x dt m g = ------(1) 波动方程 222222222(,,,)f t x y z u u u u a t x y z ?? ? ???+????=++????------(2) 热传导方程 2222222(,,,)f t x y z u u u u a t x y z ?? ? ??? +????=++???? ------(3) 静电场位方程 2222 222(,,)f x y z u u u a x y z ?? ?=- ??? ???++??? ------(4) 激波方程 0u u u t x ??+=?? ------(5) 等等。 其中(1)为一维常微分方程；(2)----(4)为三维偏微分方程；(5)为一维偏微分方程。 这些数学中的微分方程均来自物理问题，有着各自的物理背景，从数量关系上反映着相应的物理规律，称为数学物理方程，简称数理方程。 数学物理方程是数学与物理学的交叉分支学科。从物理上讲它是理论物理的基本工具；在数学上属于应用数学的(偏)微分方程分支。 本课程主要研究和讨论三类数理方程(2)，(3)，(4)的建立(导出)以及几种常用的典型的求解方法。 为了下面研究和讨论的方便，先引入有关微分方程的几个基本概念

(术语)。 1. 常，偏微分方程 只含一个自变量，关于该变量的未知函数，以及未知函数对该变量的导数的微分方程为常微分方程，如(1)。 含有多个自变量，关于这些变量的未知函数，以及未知函数对这些变量的偏导数的微分方程为偏微分方程，如(2)----(5)。 2. 阶 上述(1)----(5)均可改写成如下形式 220d x m g dt -= ------(1’) 222 30u t a u f -???-= -------(2’) 230u t a u f -???-= ------(3’) 230a u f ?+= ------(4’) 0u u t x u +????= ------(5’) 其中 222 3222x y z ????=++???，x=x(t)，u=u(t,x,y,z)或u(x,y,z)，f=f(t,x,y,z) 或f(x,y,z)。 这些方程可归纳为如下形式 12 121212,,,,,,,,,,n m n m m m n n u u u u F x x x u x x x x x x ?? ? ?? ? ?????????????????????????=0， 其中12n m m m m =++???+为导数的最高阶数，成为方程的阶。 3. 线性、非线性偏微分方程

偏微分方程数值解法

偏微分方程数值解法 Prepared on 22 November 2020

第十章偏微分方程数值解法 偏微分方程问题，其求解十分困难。除少数特殊情况外，绝 大多数情况均难以求出精确解。因此，近似解法就显得更为重要。本章仅介绍求解各类典型偏微分方程定解问题的差分方法。 §1差分方法的基本概念 几类偏微分方程的定解问题 椭圆型方程：其最典型、最简单的形式是泊松（Poisson）方程 特别地，当 ) , (≡ y x f 时，即为拉普拉斯（Laplace）方程，又 称 为调和方程 Poisson方程的第一边值问题为 其中Ω为以Γ为边界的有界区域，Γ为分段光滑曲线， ΓΩ 称为定解区域， ) ,(y x f ， ) , (y x ? 分别为Ω，Γ上的已知 连续函数。 第二类和第三类边界条件可统一表示为 其中n 为边界Γ的外法线方向。当0 = α时为第二类边界条件，0 ≠ α时为第三类边界条件。 抛物型方程：其最简单的形式为一维热传导方程方程可以有两种不同类型的定解问题： 初值问题

初边值问题 其中 ) (x ? ， )( 1 t g ， )( 2 t g 为已知函数，且满足连接条件 边界条件 )( ),( ), ( ),0( 2 1 t g t l u t g t u= = 称为第一类边界条 件。 第二类和第三类边界条件为 其中 )( 1 ≥ t λ ， )( 2 ≥ t λ 。当 )( )( 2 1 ≡ =t tλ λ 时，为第二 类边界条件， 否则称为第三类边界条件。 双曲型方程: 最简单形式为一阶双曲型方程 物理中常见的一维振动与波动问题可用二阶波动方程 描述，它是双曲型方程的典型形式。方程的初值问题为 边界条件一般也有三类，最简单的初边值问题为 差分方法的基本概念 差分方法又称为有限差分方法或网格法，是求偏微分方程定解问题的数值解中应用最广泛的方法之一。 它的基本思想是：先对求解区域作网格剖分，将自变量的连续变化区域用有限离散点（网格点）集代替；将问题中出现的连续变量的函数用定义在网格点上离散变量的函数代替；通过用网格点上函数的差商代替导数，将含连续变量的偏微分方程定解问题化成只含有限个未知数的代数方程组（称为差分格式）。如果差分格式有解，且当网格无限变小时其解收敛于原微分方程定解

常见的几个偏微分方程的解法参考文献汇总

Boussinesq:0)(220=---xx xxxx xx tt u u u c u βα 1. Jacobi 椭圆函数展开法 [2]刘式适,傅遵涛,刘式达,赵强.Jacobi 椭圆函数展开法及其在求解非线性波动方程中的应用[J].物理学报,2001(11):2068-2073. ；[29]闻小永.Boussinesq 方程的Jacobi 椭圆函数精确解[J].北京机械工业学院学报,2007(01):23-26. 2. 三角函数法和吴文俊消元 [26]贺锋,郭启波,刘辽.用三角函数法获得非线性Boussinesq 方程的广义孤子解[J].物理学报,2007(08):4326-4330. 3. 双函数法、吴文俊消元 [27]黄文华,张解放,盛正卯.Boussinesq 方程的新显式精确行波解[J].浙江大学学报(理学版),2003(02):145-149. 4. 齐次平衡法、backlund 变换 [28]夏铁成,张鸿庆,李佩春.Boussinesq 方程精确解析解研究[J].大连理工大学学报,2003(04):393-396. 5.推广的Tanh 法、Jacobi 椭圆函数、双曲函数 [30]高亮,徐伟,申建伟,唐亚宁.Boussinesq 方程新的显式行波解[J].西南民族大学学报(自然科学版),2006(01):54-59. 6.试探方程法、齐次平衡法 [39]杨玉婷,崔泽建.用试探方程法求解Boussinesq 方程[J].重庆文理学院学报(自然科学版),2012,31(03):5- 7. 7.拓展的Jacobi 椭圆函数法 [44]钟太勇,钟远涛.用形变映射法求KdV 方程的显式精确行波解 [J].江汉大学学报(自然科学版),2009,37(03):10-12. 8.改进的试探函数法 [45]谢元喜,唐驾时.用改进的试探函数法求解Boussinesq 方程[J].安阳工学院学报,2005(06):73-76. 9.形变映射法 [46]方建平.形变映射法构造非线性Boussinesq 方程的行波解[J].丽水师范专科学校学报,2003(02):12-15. Sine -Gordon 方程：0sin =+-u u u tt xx 1. 直接积分法 [1]李志斌. 非线性数学物理方程的行波解[M]. 科学出版社, 2007. 2. 混合指数方法. [1]李志斌. 非线性数学物理方程的行波解[M]. 科学出版社, 2007. 3. F -展开法 齐次平衡法 [13]王明亮,聂惠,李向正.用F 展开法解Sine -Gordon 方程[J].河南科技大学学报(自然科学版),2005(01):79-82. ；[32]范建华,闫杰生.Sine -Gordon 方程的精确解[J].商丘职业技术学院学报,2004(06):11-13+21. 4.扩展的sinh -Gordon 方程展开法 [95]杨先林,唐驾时.非线性演化方程的新Jacobi 椭圆函数解[J].动力学与控制学报,2011,9(02):147-151. 5.双线性算子、齐次平衡 [99]杨琼芬,唐再良,罗守双.用双线性形式求得sine -Gordon 方程新的精确解[J].绵阳师范学院学报,2015,34(11):12-14+29. 6.Jacobi 椭圆函数展开法 [100]沈水金.利用Jacobi 椭圆函数展开法求解特殊类型的方程[J].上海大学学报(自然科学版),2010,16(04):383-386. Fisher 方程：0)1(=---u u u u xx t βα 1. 观察试凑法 [1]李志斌. 非线性数学物理方程的行波解[M]. 科学出版社, 2007. 2. 指数函数法 [51]王军帽,张睿,张文亮,张苗,韩家骅.Exp 函数法与Fisher 方程新的精确解[J].安徽大学学报(自然科学版),2009,33(01):53-56. ； [52]张桂戌,李志斌,段一士.非线性波方程的精确孤立波解 [J].中国科学(A 辑),2000(12):1103-1108. 3.正切函数变换 [53]张宏.Fisher 方程的新孤波解[J].青海师范大学学报(自然科学版),2006(03):37-38+67. 4.推广的tanh 函数法、复tanh 函数法、广义幂指函数法 [54]庄红波. 函数变换法求经典Fisher 方程的显示解[D].四川师范大学,2006.