第九章偏微分方程差分方法汇总.pdf
《第九章偏微分方程差分方法汇总.pdf》由会员分享,可在线阅读,更多相关《第九章偏微分方程差分方法汇总.pdf(15页珍藏版)》请在咨信网上搜索。
1、170第第 9 9 章章 偏微分方程的差分方法偏微分方程的差分方法含有偏导数的微分方程称为偏微分方程。由于变量的增多和区域的复杂性,求偏微分方程的精确解一般是不可能的,经常采用数值方法求方程的近似解。偏微分方程的数值方法种类较多,最常用的方法是差分方法差分方法。差分方法具有格式简单,程序易于实现,计算量小等优点,特别适合于规则区域上偏微分方程的近似求解。本章将以一些典型的偏微分方程为例,介绍差分方法的基本原理和具体实现方法。9.19.1 椭圆型方程边值问题的差分方法椭圆型方程边值问题的差分方法9.1.19.1.1 差分方程的建立差分方程的建立最典型的椭圆型方程是 PoissonPoisson(
2、泊松)方程(泊松)方程 (9.1)Gyxyxfyuxuu),(),()(2222G 是 x,y 平面上的有界区域,其边界为分段光滑的闭曲线。当 f(x,y)0 时,方程(9.1)称为 Laplace(Laplace(拉普拉斯)方程拉普拉斯)方程。椭圆型方程的定解条件主要有如下三种边界条件 第一边值条件 (9.2)),(yxu 第二边值条件 (9.3)),(yxnu 第三边值条件 (9.4)),()(yxkunu这里,n 表示 上单位外法向,(x,y),(x,y),(x,y)和 k(x,y)都是已知的函数,k(x,y)0。满足方程(9.1)和上述三种边值条件之一的光滑函数 u(x,y)称为椭圆型
3、方程边值问题的解。用差分方法求解偏微分方程,就是要求出精确解 u(x,y)在区域 G 的一些离散节点(xi,yi)上的近似值 ui,j(xi,yi)。差分方法的基本思想是,对求解区域 G 做网格剖分,将偏微分方程在网格节点上离散化,导出精确解在网格节点上近似值所满足的差分方程,最终通过求解差分方程,通常为一个线性方程组,得到精确解在离散节点上的近似值。设 G=0 xa,0y0,B(x,y)Bmin 0,E(x,y)0。引进半节点,12121hxxii利用一阶中心差商公式,在节点(i,j)处可有,22121hyyii173)(2),1(),1(),()(),1(),(),(),1(1)(),21
4、)(),21)(1),)(211211,211,211211hOhjiujiujixuhOhjiujiuAhjiujiuAhhOjixuAjixuAhjixuAxjiji对类似处理,就可推得求解方程(9.9)的差分方程yuyuBy),((9.10)hjijijijijijijijijijiGjijifuauauauaua),(),(,1,1,1,1,1,1,1,1其中 (9.11)jijijijijijijijijijijijijijijijijijiEBBhAAhaDhBhaDhBhaChAhaChAha,21,21,22,21,2121,221,221,221,221,1,2121,1,1
5、,2121,1)()()2()2()2()2(显然,当系数函数 A(x,y)=B(x,y)=1,C(x,y)=D(x,y)=E(x,y)=0 时,椭圆型方程(9.9)就成为 Poisson 方程(9.1),而差分方程(9.10)就成为差分方程(9.6)。容易看出,差分方程(9.10)的截断误差为阶。)(2221hhO9.1.29.1.2 一般区域的边界条件处理一般区域的边界条件处理前面已假设 G 为矩形区域,现在考虑 G 为一般区域情形,这里主要涉及边界条件的处理。考虑 Poisson 方程第一边值问题 (9.12)),(),(),(),(yxyxuGyxyxfu其中 G 可为平面上一般区域,
6、例如为曲边区域。仍然用两组平行直线:x=x0+ih1,y=y0+jh2,i,j=0,1,对区域 G 进行矩形网格剖分,见图 9-3。174 如果一个内节点(i,j)的四个相邻节点(i+1,j),(i-1,j),(i,j+1)和(i,j-1)属于,则称其为正则内点正则内点,见图 9-3 中打“。”号者;如果一个节点 GG(i,j)属于且不为正则内点,则称其为非正则内点非正则内点,见图 9-3 中打“.”号者。G记正则内点集合为,非正则内点集合为。显然,当 G 为矩形区域时,hGh成立。hhhhGG,在正则内点(i,j)处,完全同矩形区域情形,可建立五点差分格式 hjijijijijijijiGj
7、ifuuuhuuuh),(,2121,1,1,22,1,121(9.13)在方程(9.13)中,当(i,j)点临近边界时,将出现非正则内点上的未知量,因此必须补充非正则内点处的方程。若非正则内点恰好是边界点,如图 9-4 中D 点,则利用边界条件可取 uD=(D)对于不是边界点的非正则内点,如图 9-4 中B 点,一般可采用如下两种处理方法。a.直接转移法直接转移法.取与点 B 距离最近的边界点(如图 9-4 中 E 点)上的 u 的值作为 u(B)的近似值 uB,即 uB=u(E)=(E)直接转移法的优点是简单易行,但精度较低,只为一阶近似。b.线性插值法线性插值法.取 B 点的两个相邻点(
8、如图 9-4 中边界点 A 和正则内点 C 作为插值节点对 u(B)进行线性插值175)()()()(21hOCuxxxxAuxxxxBuACABACBC则得到点 B 处的方程 ABCBxxuhAhhu,)(111线性插值法精度较高,为二阶近似。对每一个非正则内点进行上述处理,将所得到的方程与(9.13)式联立,就组成了方程个数与未知量个数相一致的线性代数方程组。求解此方程组就可得到一般区域上边值问题(9.12)的差分近似解。对于一般区域上二阶椭圆型方程(9.9)的第一边值问题,可完全类似处理。第二、三边值条件的处理较为复杂,这里不再讨论。9.29.2 抛物型方程的差分方法抛物型方程的差分方法
9、本节介绍抛物型方程的差分方法,重点讨论差分格式的构造和稳定性分析。9.2.19.2.1 一维问题一维问题 作为模型,考虑一维热传导的初边值问题 (9.14)Ttlxtxfxuatu0,0),(22 (9.15)lxxxu0),()0,((9.16)Tttgtlutgtu0),(),(),(),0(21其中 a 是正常数,都是已知的连续的函数。)()(),(),(21tgtgxtxf和现在讨论求解问题(9.14)-(9.18)的差分方法。首先对求解区域 G=0 xl,0tT进行网格剖分。取空间步长 h=l/N,时间步长=T/M,其中 N,M 是正整数,作两族平行直线 MkkttNjjhxxkj,
10、1,0,1,0,将区域 G 剖分成矩形网格,见图 9-5,网格交点(xj,tk)称为节点。176用差分方法求解初边值问题(9.14)-(9.16)就是要求出精确解 u(x,t)在每个节点(xj,tk)处的近似值。为简化记号,简记节点(xj,tk)),(kjkjtxuu=u(j,k)。利用一元函数的 Taylor 展开公式,可推出下列差商表达式 (9.17))(),()1,(),(Okjukjukjtu (9.18))()1,(),(),(Okjukjukjtu (9.19))(2)1,()1,(),(2Okjukjukjtu (9.20))(),1(),(2),1(),(2222hOhkjuk
11、jukjukjxu1.1.古典显格式古典显格式 在区域 G 的内节点(j,k)处,利用公式(9.17)和(9.20),可将偏微分方程(9.14)离散为)(),1(),(2),1(),()1,(22hOfhkjukjukjuakjukjukj其中。舍去高阶小项,就得到节点近似值(差分解)),(kikjtxff)(2hO所满足的差分方程kju (9.21)kjkjkjkjkjkjfhuuuauu21112显然,在节点(j,k)处,差分方程(9.21)逼近偏微分方程(9.14)的误差为,这个误差称为截断误差截断误差,它反映了差分方程逼近偏微分方程的精度。)(2hO177现将(9.21)式改写为便于计
12、算的形式,并利用初边值条件(9.15)与(9.16)补充上初始值和边界点方程,则得到 (9.22)MktgutguNjxuMkNjfruurruukkNkkjjkjkjkjkjkj,1,0),(),(1,2,1),(1,1,0,1,2,1)21(2100111其中称为网比网比。2har与时间相关问题差分方程的求解通常是按时间方向逐层进行的。对于差分方程(9.22),当第 k 层节点值已知时,可直接计算出第 k+1 层节点值。这kju1kju样,从第 0 层已知值开始,就可逐层求出各时间层的节点值。差分方程)(0ijxu(9.22)的求解计算是显式的,无须求解方程组,故称为古典显格式古典显格式。
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 第九 微分方程 方法 汇总
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【精****】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【精****】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。