- 1、本文档共25页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 5、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 6、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 7、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 8、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
抛物型偏微方程.ppt
* * 第七章 解抛物型偏微分方程的差分解法 稳定性讨论是差分法的重要问题。 常遇到的抛物型偏微分方程有: 与上一章的边值问题不同, 这些问题一般都是初值问题. 基本作法:按等间隔剖分,用各离散点上的差商来近似代替该点的偏微商,把偏微分方程转化为线性代数方程组,进而解出各离散点上的待求函数随时间的变化。 主要内容 §7.1简单差分法及其不稳定性 §7.2隐格式差分法与三对角方程 §7.3二维扩散方程的差分格式 §7.4含时间薛定格方程的差分格式 §7.1简单差分法及其不稳定性 一、差分格式 二、例 设 ,初条件和边界条件为: 容易检验,满足微分方程(7.4)和条件(7.8)的理论解为: 现在,我们利用递推式 作数值计算,并比较它与理论解(7.9)的误差。 取空间间隔,时间间隔,时间步数为20(t=0.015)和60(t=0.045)。 当取 时,数值解与理论解的误差很小。 但当取 时,结果就变得很糟。 DO 20 ITER=1,NITER POLD=0. DO 30 IX=1,NSTEP-1 PNEW=PHI(IX)+DTH*(POLD+PHI(IX+1)- * 2*PHI(IX)) EXACT(X,T)=GAUSS(X,T)-GAUSS(X-1.,T)- * GAUSS(X+1.,T) POLD=PHI(IX) PHI(IX)=PNEW 30 CONTINUE IF (MOD(ITER,10) .EQ. 0) THEN PRINT *, iteration = , ITER, * time = ,ITER*DT T=ITER*DT DO 40 IX=1,NSTEP-1 DIFF=PHI(IX)-EXACT(IX*H,T) PRINT *, phi = , PHI(IX), error = , DIFF 40 CONTINUE END IF 20 CONTINUE GOTO 50 END 71.FOR PARAMETER (NSTEP=25) DIMENSION PHI(0:NSTEP) GAUSS(X,T)=EXP(-20.*(X-.5)**2/ * (1.+80*T))/ SQRT(1+80*T) H=1./NSTEP 50 PRINT *, Enter time step and * total time (0 to stop) READ *,DT,TIME IF (DT .EQ. 0.) STOP NITER=TIME/DT DTH=DT/H**2 T=0. PHI(0)=0. PHI(NSTEP)=0. DO 10 IX=1,NSTEP-1 PHI(IX)=EXACT(IX*H,T) 10 CONTINUE 这个差分递推公式在超过一定条件时是不稳定的。 三、 稳定性分析。 利用该算符,差分格式(7.7)可改写成矩阵形式: (7.12) 令 为算符 的本征函数.相应的本征本值为 由于 是厄密算符, 所以本征值是实数. 由于 所以有精确解 而从差分格式 我们有 因此,只要步长选得足够小, 它近似于精确解. 但当步长过大, 以至于对于最大本征值的分量 这种方法就变得不稳定. 可以从数学上严格证明,算符 H 的本征值都是正值,最大值为 . 为了使计算稳定,则要求 当时 ,差分格式的稳定性条件要求 §7.2隐格式差分法与三对角方程 把方程中的空间导数换成 n+1 时刻有 用矩阵 H 并写成矩阵形式 形式解可写成 显然,对于任意的的本征值,都会使
您可能关注的文档
- 强电讲义.ppt
- 徐徐道来永定河--20170517.pptx
- 微分方程 8.ppt
- 微分方程differential-equations.ppt
- 微分方程与差分方程建模.pptx
- 微分方程建模的若干问题讲稿下.ppt
- 微分方程数值解法第二章.ppt
- 微分方程模型03--人口模型等.ppt
- 微分方程模型全.ppt
- 微分方程的数值解法 1.ppt
- 2025年山东省济南市单招(语文)测试模拟题库精编.docx
- 2025年安徽省蚌埠市选调生考试(行政职业能力测验)综合能力测试题及答案1套.docx
- 2025年安徽黄梅戏艺术职业学院单招语文测试题库最新.docx
- 2025年安徽省宣城市事业单位招聘考试(职业能力倾向测验)题库参考答案.docx
- 2025年安徽省阜阳市选调生考试(公共基础知识)综合能力题库完整版.docx
- 2025年宣化科技职业学院单招语文测试模拟题库及答案参考.docx
- 2025年安徽省芜湖市单招语文测试题库精编.docx
- 2025年初级经济师之初级经济师人力资源管理题库500道及1套参考答案.docx
- 2025年国家电网招聘之电网计算机题库500道附完整答案(典优).docx
- 2025年安徽省芜湖市单招(语文)测试题库学生专用.docx
文档评论(0)