Matlab PDE工具箱有限元法求解偏微分方程

在科学技术各领域中，有很多问题都可以归结为偏微分方程问题。在物理专业的力学、热学、电学、光学、近代物理课程中都可遇见偏微分方程。

偏微分方程，再加上边界条件、初始条件构成的数学模型，只有在很特殊情况下才可求得解析解。随着计算机技术的发展，采用数值计算方法，可以得到其数值解。

偏微分方程基本形式

而以上的偏微分方程都能利用PDE工具箱求解。

PDE工具箱

PDE工具箱的使用步骤体现了有限元法求解问题的基本思路，包括如下基本步骤：

1) 建立几何模型

2) 定义边界条件

3) 定义PDE类型和PDE系数

4) 三角形网格划分

5) 有限元求解

6) 解的图形表达

以上步骤充分体现在PDE工具箱的菜单栏和工具栏顺序上，如下

具体实现如下。

打开工具箱

输入pdetool可以打开偏微分方程求解工具箱，如下

首先需要选择应用模式，工具箱根据实际问题的不同提供了很多应用模式，用户可以基于适

当的模式进行建模和分析。

在Options菜单的Application菜单项下可以做选择，如下

或者直接在工具栏上选择，如下

列表框中各应用模式的意义为：

① Generic Scalar：一般标量模式（为默认选项）。

② Generic System：一般系统模式。

③ Structural Mech.，Plane Stress：结构力学平面应力。

④ Structural Mech.，Plane Strain：结构力学平面应变。

⑤ Electrostatics：静电学。

⑥ Magnetostatics：电磁学。

⑦ Ac Power Electromagnetics：交流电电磁学。

⑧ Conductive Media DC：直流导电介质。

⑨ Heat Tranfer：热传导。

⑩ Diffusion：扩散。

可以根据自己的具体问题做相应的选择，这里要求解偏微分方程，故使用默认值。此外，对于其他具体的工程应用模式，此工具箱已经发展到了Comsol Multiphysics软件，它提供了更强大的建模、求解功能。

另外，可以在菜单Options下做一些全局的设置，如下

l Grid：显示网格

l Grid Spacing…：控制网格的显示位置

l Snap：建模时捕捉网格节点，建模时可以打开

l Axes Limits…：设置坐标系范围

l Axes Equal：同Matlab的命令axes equal命令

建立几何模型

使用菜单Draw的命令或使用工具箱命令可以实现简单几何模型的建立，如下

各项代表的意义分别为

l 绘制矩形或方形；

l 绘制同心矩形或方形；

l 绘制椭圆或圆；

l 绘制同心椭圆或圆；

l 绘制多义线。

这里只绘制一个圆如下

定义边界条件

选择Boundary菜单下的Specify Boundary Conditions…，如下

定义PDE类型和PDE系数

选择PDE菜单下的PDE Specifica tions…，如下

三角形网格划分

选择Mesh菜单下的Initialize Mesh初始化三角形网格，再选择Refine Mesh改进初始网格并细化网格，如下

初始化网格

细化网格

另外还可以进一步选择Jiggle Mesh微调网格。最后可以选择Display Triangle Quality 显示三角形网格的质量图，其中1表示质量最好，0表示最差，如下

有限元求解

选择Solve菜单下的Solve PDE选项进行PDE问题的求解，如下

解的图形表达

选择Plot菜单下的Parameters…可以设置显示的效果，如下

显示结果如下

比较数值解与精确解的误差：

可见数值解的精度是很高的。

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

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

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

Matlab求解微分方程(组)及偏微分方程(组)

第四讲 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

偏微分方程数值解法

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

Matlab求解微分方程(组)及偏微分方程(组)

第四讲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就是显示微分方程在积分区间tspan上从到用初始条件求解、 (3)如果要获得微分方程问题在其她指定时间点上得解,则令tspan(要求就是单调得)、 (4)因为没有一种算法可以有效得解决所有得ODE问题,为此,Matlab提供了多种求解器solver,对于不同得ODE问题,采用不同得solver、 表1 Matlab中文本文件读写函数

MatlabPDE工具箱有限元法求解偏微分方程

在科学技术各领域中，有很多问题都可以归结为偏微分方程问题。在物理专业的力学、热学、电学、光学、近代物理课程中都可遇见偏微分方程。 偏微分方程，再加上边界条件、初始条件构成的数学模型，只有在很特殊情况下才可求得解析解。随着计算机技术的发展，采用数值计算方法，可以得到其数值解。 偏微分方程基本形式 而以上的偏微分方程都能利用PDE工具箱求解。 PDE工具箱 PDE工具箱的使用步骤体现了有限元法求解问题的基本思路，包括如下基本步骤： 1) 建立几何模型 2) 定义边界条件 3) 定义PDE类型和PDE系数 4) 三角形网格划分

5) 有限元求解 6) 解的图形表达 以上步骤充分体现在PDE工具箱的菜单栏和工具栏顺序上，如下 具体实现如下。 打开工具箱 输入pdetool可以打开偏微分方程求解工具箱，如下 首先需要选择应用模式，工具箱根据实际问题的不同提供了很多应用模式，用户可以基于适

当的模式进行建模和分析。 在Options菜单的Application菜单项下可以做选择，如下 或者直接在工具栏上选择，如下 列表框中各应用模式的意义为： ① Generic Scalar：一般标量模式（为默认选项）。 ② Generic System：一般系统模式。 ③ Structural Mech.，Plane Stress：结构力学平面应力。

④ Structural Mech.，Plane Strain：结构力学平面应变。 ⑤ Electrostatics：静电学。 ⑥ Magnetostatics：电磁学。 ⑦ Ac Power Electromagnetics：交流电电磁学。 ⑧ Conductive Media DC：直流导电介质。 ⑨ Heat Tranfer：热传导。 ⑩ Diffusion：扩散。 可以根据自己的具体问题做相应的选择，这里要求解偏微分方程，故使用默认值。此外，对于其他具体的工程应用模式，此工具箱已经发展到了Comsol Multiphysics软件，它提供了更强大的建模、求解功能。 另外，可以在菜单Options下做一些全局的设置，如下 l Grid：显示网格 l Grid Spacing…：控制网格的显示位置 l Snap：建模时捕捉网格节点，建模时可以打开 l Axes Limits…：设置坐标系范围 l Axes Equal：同Matlab的命令axes equal命令

有限差分法求解偏微分方程MATLAB教学教材

有限差分法求解偏微分方程M A T L A B

南京理工大学 课程考核论文 课程名称：高等数值分析 论文题目：有限差分法求解偏微分方程姓名：罗晨 学号： 115104000545 成绩： 有限差分法求解偏微分方程

一、主要内容 1.有限差分法求解偏微分方程，偏微分方程如一般形式的一维抛物线型方程： 22(,)()u u f x t t x αα??-=??其中为常数 具体求解的偏微分方程如下： 22001 (,0)sin()(0,)(1,)00 u u x t x u x x u t u t t π???-=≤≤?????? =??? ==≥??? 2.推导五种差分格式、截断误差并分析其稳定性； 3.编写MATLAB 程序实现五种差分格式对偏微分方程的求解及误差分析； 4.结论及完成本次实验报告的感想。 二、推导几种差分格式的过程： 有限差分法（finite-difference methods ）是一种数值方法通过有限个微分方程近似求导从而寻求微分方程的近似解。有限差分法的基本思想是把连续的定解区域用有限个离散点构成的网格来代替；把连续定解区域上的连续变量的函数用在网格上定义的离散变量函数来近似；把原方程和定解条件中的微商用差商来近似，积分用积分和来近似，于是原微分方程和定解条件就近似地代之以代数方程组，即有限差分方程组，解此方程组就可以得到原问题在离散点上的近似解。 推导差分方程的过程中需要用到的泰勒展开公式如下： ()2 100000000()()()()()()()......()(()) 1!2!! n n n f x f x f x f x f x x x x x x x o x x n +'''=+-+-++-+- (2-1) 求解区域的网格划分步长参数如下：

有限差分法求解偏微分方程MATLAB

南京理工大学 课程考核论文 课程名称：高等数值分析 论文题目：有限差分法求解偏微分方程姓名：罗晨 学号： 成绩： 有限差分法求解偏微分方程

一、主要内容 1.有限差分法求解偏微分方程，偏微分方程如一般形式的一维抛物线型方程： 22(,)()u u f x t t x αα??-=??其中为常数 具体求解的偏微分方程如下： 22001 (,0)sin()(0,)(1,)00 u u x t x u x x u t u t t π???-=≤≤?????? =??? ==≥??? 2.推导五种差分格式、截断误差并分析其稳定性； 3.编写MATLAB 程序实现五种差分格式对偏微分方程的求解及误差分析； 4.结论及完成本次实验报告的感想。 二、推导几种差分格式的过程： 有限差分法（finite-difference methods ）是一种数值方法通过有限个微分方程近似求导从而寻求微分方程的近似解。有限差分法的基本思想是把连续的定解区域用有限个离散点构成的网格来代替；把连续定解区域上的连续变量的函数用在网格上定义的离散变量函数来近似；把原方程和定解条件中的微商用差商来近似，积分用积分和来近似，于是原微分方程和定解条件就近似地代之以代数方程组，即有限差分方程组，解此方程组就可以得到原问题在离散点上的近似解。 推导差分方程的过程中需要用到的泰勒展开公式如下： ()2100000000()()()()()()()......()(()) 1!2!! n n n f x f x f x f x f x x x x x x x o x x n +'''=+-+-++-+- (2-1) 求解区域的网格划分步长参数如下：

偏微分方程数值解法

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

内 容 简 介 本书试图用较少的篇幅描述偏微分方程的几种数值方法. 主要内容包括：Sobolev空间初步, 椭圆边值问题的变分问题, 椭圆问题的有限差分方法, 抛物型方程的有限差分方法, 双曲型方程的有限差分方法, 椭圆型方程的有限元方法, 抛物及双曲方程的有限元方法, 椭圆型方程的混合有限元方法, 谱方法等. 本书内容丰富, 深入浅出, 尽可能地用简单的方法来描述一些理论结果, 并根据作者对有限差分、有限元、混合有限元、谱方法的理解和研究生教学要求, 全面、客观地评价各种数值计算方法，并列举一些数值计算的例子, 阐述许多新的学术观点. 本书可作为高等学校数学系高年级本科生和研究生的教材或参考书, 也可作为计算数学工作者和从事科学与工程计算的科研人员的参考书. 图书在版编目(CIP)数据 偏微分方程数值解法/陈艳萍, 鲁祖亮, 刘利斌编著. —北京：科学出版社, 2015.1 (信息与计算科学丛书67) ISBN 978-7-03-000000-0 Ⅰ. ①偏… Ⅱ. ①陈… ②鲁… ③刘… Ⅲ. ① Ⅳ.① 中国版本图书馆CIP数据核字(2014) 第000000号 责任编辑: 王丽平／责任校对: 彭涛 责任印制: 肖钦／封面设计: 陈敬 出版 北京东黄城根北街16号 邮政编码: 100717 https://www.wendangku.net/doc/0319338267.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. 偏微分方程求解问题的描述 教材P653[12.1.1]椭圆型 教材P653[12.1.2] 教材P664[12.2.1]双曲型 教材P665[12.2.4]拉普拉斯泊松 对流 波动 教材P684[12.3.1]抛物型 教材P685[12.3.6]扩散 对流扩散 教材P686[12.3.8]二维扩散 教材P678[12.2.23]二维对流

??????????????????????≥≤≤==≥≤≤==≤≤=>≥≤≤≤≤???? ????+??=??0,0, ),(),,(),(),0,(0,0,),(),,(),(),,0(,0,),()0,,(0,0 , 0 , 0 21212222t L x t x v t L x u t x v t x u t L y t y t y L u t y t y u L y x y x y x u b t L y L x y u x u b t u μμ?Ω 求解域初值条件边值条件) ,,(t y x u 未知函数

????? ? ????????????????????≥<<-==≥<<==≥≤≤-==≥≤≤==≤≤==≤≤≤≤≤≤???? ????+??=??0 , 50 , sin 255sin ),(),5,(0 , 50 , 0),(),0,(0 , 50 , 5sin sin 25),(),,5(0 , 50 , 0),(),,0(5,0,0),()0,,( 10000 , 50 , 50 001.022********t x x x t x v t x u t x t x v t x u t y y y t y t y u t y t y t y u y x y x y x u t y x y u x u t u μμ?Ω 求解域初值条件边值条件以具体问题为例演示具体的求解过程) ,,(t y x u 未知函数

Matlab求解微分方程(组)及偏微分方程(组)

第四讲 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 上的解，则令 tspan 012[,,,]f t t t t =(要求是单调的). (4)因为没有一种算法可以有效的解决所有的ODE 问题，为此，Matlab 提供

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

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

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

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

偏微分方程数值解法答案

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),

有限差分法解偏微分方程

有限差分法解偏微分方程综述 绪论 有限元方法最早应用于结构力学，后来随着计算机的发展慢慢用于流体力学的数值模拟。在有限元方法中，把计算域离散剖分为有限个互不重叠且相互连接的单元，在每个单元内选择基函数，用单元基函数的线形组合来逼近单元中的真解，整个计算域上总体的基函数可以看为由每个单元基函数组成的，则整个计算域内的解可以看作是由所有单元上的近似解构成。 有限差分方法(FDM)是计算机数值模拟最早采用的方法，至今仍被广泛运用。该方法将求解域划分为差分网格，用有限个网格节点代替连续的求解域。有限差分法以Taylor 级数展开等方法，把控制方程中的导数用网格节点上的函数值的差商代替进行离散，从而建立以网格节点上的值为未知数的代数方程组。该方法是一种直接将微分问题变为代数问题的近似数值解法，数学概念直观，表达简单，是发展较早且比较成熟的数值方法。 对于有限差分格式，从格式的精度来划分，有一阶格式、二阶格式和高阶格式。 从差分的空间形式来考虑，可分为中心格式和逆风格式。 考虑时间因子的影响，差分格式还可以分为显格式、隐格式、显隐交替格式等。 目前常见的差分格式，主要是上述几种形式的组合，不同的组合构成不同的差分格式。差分方法主要适用于有结构网格，网格的步长一般根据实际地形的情况和柯朗稳定条件来决定。构造差分的方法有多种形式， 目前主要采用的是泰勒级数展开方法。其基本的差分表达式主要有三种形式：一阶向前差分、一阶向后差分、一阶中心差分和二阶中心差分等，其中前两种格式为一阶计算精度，后两种格式为二阶计算精度。通过对时间和空间这几种不同差分格式的组合，可以组合成不同的差分计算格式。 有限元方法的基础是变分原理和加权余量法，其基本求解思想是把计算域划分为有限个互不重叠的单元，在每个单元内，选择一些合适的节点作为求解函数的插值点，将微分方程中的变量改写成由各变量或其导数的节点值与所选用的插值函数组成的线性表达式，借助于变分原理或加权余量法，将微分方程离散求解。采用不同的权函数和插值函数形式，便构成不同的有限元方法。 有限差分法求解偏微分方程 在采用数值计算方法求解偏微分方程时，若将每一处导数由有限差分近似公式替代，从而把求解偏微分方程的问题转换成求解代数方程的问题，即所谓的有限差分法。有限差分法求解偏微分方程的步骤如下： 1、区域离散化，即把所给偏微分方程的求解区域细分成由有限个格点组成的网格； 2、近似替代，即采用有限差分公式替代每一个格点的导数； 3、逼近求解。换而言之，这一过程可以看作是用一个插值多项式及其微分来代替偏微分方程的解的过程 有限差分法的应用 抛物型方程的差分方法 1. 简单差分法

c4d教程CAD教案

CAD教案 第一课基本操作 1、AutoCAD： AutoCAD介绍：AutoCAD是由美国Autodesk公司于1982年推出的一款计算机辅助设计软件，AutoCAD简称CAD，版本从R1.0开始发展，到98年推出CAD R14，2000年起以年份命名，分别有2000版R15、02版R15.1、04/06版R16、07—09版R17， CAD2004版是目前行业中的主流版本；除AutoCAD外还有国产的中旺CAD（中旺CAD的操作习惯和命令与AutoCAD相同）、另外还有天正CAD（天正CAD是在AutoCAD基础上的二次开发命令，相当于一个外挂，主要用于建筑行业）、燕秀工具箱（用 于模具设计）； AutoCAD的功能主要有二维绘图(.dwg)和基本三维造型，它是一个用途非常广泛的软件，无论你是从事文职类、广告设计类、工业制造类、还是建筑装潢类的工作都需要使用AutoCAD；AutoCAD 是学习设计的基础课程，只有学好AutoCAD才能为后续学习广告设计/产品开发/模具设计/建筑装潢设计等提供坚实的基础；AutoCAD总共分为15节课，学习时间计划为20天；今天我们先学习第一课，首选我们在桌面上找一下AutoCAD的程序图标，再图标上双击左键即可打开CAD软件； 2、界面组成：标题栏、菜单栏、工具栏、绘图区、命令窗口、状态栏； 3、操作规定：左手控制键盘右手控制鼠标，空格键为确定； 4、TO工具栏：标准、绘图、修改； 5、启动命令：快捷命令（常用）、工具栏按钮、菜单命令； 6、L直线：U放弃上一段线，C终点与当次命令起点闭合； 7、空格作用：启动命令、结束命令（ESC）、延用命令； 8、F3捕捉开关：用于精确绘图，绘图时必须打开； 9、OS捕捉设置：一般全部选取；

偏微分方程数值解法

《偏微分方程数值解法》 课程设计 题目: 六点对称差分格式解热传导方程的初边 值问题 姓名: 王晓霜 学院: 理学院 专业: 信息与计算科学 班级: ０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 是自然数,

有限差分法求解偏微分方程

有限差分法求解偏微分方程 摘要：本文主要使用有限差分法求解计算力学中的系统数学模型，推导了有限差分法的理论基础，并在此基础上给出了部分有限差分法求解偏微分方程的算例验证了推导的正确性及操作可行性。 关键词：计算力学，偏微分方程，有限差分法 Abstract：This dissertation mainly focuses on solving the mathematic model of computation mechanics with finite-difference method. The theoretical basis of finite-difference is derived in the second part of the dissertation, and then I use MATLAB to program the algorithms to solve some partial differential equations to confirm the correctness of the derivation and the feasibility of the method. Key words：Computation Mechanics, Partial Differential Equations, Finite-Difference Method

1 引言 机械系统设计常常需要从力学观点进行结构设计以及结构分析，而这些分析的前提就是建立工程问题的数学模型。通过对机械系统应用自然的基本定律和原理得到带有相关边界条件和初始条件的微分积分方程，这些微分积分方程构成了系统的数学模型。 求解这些数学模型的方法大致分为解析法和数值法两种，而解析法的局限性众所周知，当系统的边界条件和受载情况复杂一点，往往求不出问题的解析解或近似解。另一方面，计算机技术的发展使得计算更精确、更迅速。因此，对于绝大多数工程问题，研究其数值解法更具有实用价值。对于微分方程而言，主要分为差分法和积分法两种，本论文主要讨论差分法。 2 有限差分法理论基础 2.1 有限差分法的基本思想 当系统的数学模型建立后，我们面对的主要问题就是微分积分方程的求解。基本思想是用离散的只含有限个未知量的差分方程组去近似地代替连续变量的微分方程和定解条件，并把差分方程组的解作为微分方程定解问题的近似解。将原方程及边界条件中的微分用差分来近似，对于方程中的积分用求和或及机械求积公式来近似代替，从而把原微分积分方程和边界条件转化成差分方程组。有限差分法求解偏微分方程的步骤主要有以下几步： ?区域离散，即把所给偏微分方程的求解区域细分成由有限个格点组成的网格， 这些离散点称作网格的节点； ?近似替代，即采用有限差分公式替代每一个格点的导数； ?逼近求解，换而言之，这一过程可以看作是用一个插值多项式及其微分来代 替偏微分方程的解的过程。 从原则上说，这种方法仍然可以达到任意满意的计算精度。因为方程的连续数值解可以通过减小独立变量离散取值的间格，或者通过离散点上的函数值进行插值计算来近似得到。理论上，当网格步长趋近于零时，差分方程组的解应该收敛于精确解，但由于机器字节的限制，网格步长不可能也没有必要取得无限小，

浩辰CAD教程模具之顶针绘制

浩辰CAD教程燕秀模具2011之八：顶针绘制 为了帮助设计师更快速、更准确的完成常见模具图素的设计，浩辰CAD燕秀模具2011软件提供了强大的模具标准件库，包含模胚、顶针等常见的标准件。在今天的CAD教程中，我们就着重来介绍下浩辰CAD燕秀模具2011中顶针的绘制、快速复制与复制对象还原等功能。 要绘制顶针并对顶针类型进行选择，我们需要在软件中依次点击菜单栏【燕秀工具箱】—【模具标准件】—【顶针】，调出顶针功能对话框（如图1） 图1 在对话框中，我们可以选择标准的顶针类型、规格、制式、剖面线、尺寸等参数，而且还可以为公司定制标准件的外挂，方便重复调用。 顶针在模具图的绘制过程中有很多方便快捷的绘制方法，比如快速绘制与复制对象还原。针对在快速复制之后顶针方向、直径标注等都会有一定角度的改变的问题，软件提供了复制对象还原的功能。 下面，我们通过具体案例来看下其功能的应用：首先，绘制顶针并快速复制（如图2）。

图2 左侧顶针为使用快速旋转复制功能来绘制的，可以看到顶针标注也旋转了。这里，我们就要用到复制对象还原命令了：依次点击【燕秀工具箱】—【块工具】—【块旋转还原】，然后窗选左侧顶针就可以快速将旋转对象还原（如图3）。 图3 以上命令的配合使用将会给提高设计师的设计效率，同时也增加了绘图的正确率，希望大家能在绘制实践中进行体会。 另外，浩辰CAD燕秀模具2011软件免费在线培训进行中，相关信息请点击浩辰CAD 官网并关注我们。https://www.wendangku.net/doc/0319338267.html,/onlinelist_1_124.html 若您对浩辰CAD燕秀模具2011软件有任何问题，可加入浩辰CAD技术QQ群：3382031。