微分方程数值解法抛物型方程问题的差分格式.docVIP

微分方程数值解法抛物型方程问题的差分格式.doc

  1. 1、本文档共8页,可阅读全部内容。
  2. 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
  3. 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  4. 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
  5. 5、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
  6. 6、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们
  7. 7、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
  8. 8、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
微分方程数值解法抛物型方程问题的差分格式

目 录 一、问题的描述 5 二、算法设计及流程图 5 2.1 算法设计 5 2.2 流程图 6 三、算法的理论依据及其推导 6 3.1 截断误差分析 6 3.2 稳定性分析 7 四、数值结果及分析 7 五、总结 9 六、附件(源代码) 10 抛物型方程问题的差分格式 一、问题的描述 有限差分方法就是一种数值解法,它的基本思想是先把问题的定义域进行网格剖分,然后在网格点上,按适当的数值微分公式把定解问题中的微商换成差商,从而把原问题离散化为差分格式,进而求出数值解。此外,还要研究差分格式的解的存在性和唯一性、解的求法、解法的数值稳定性、差分格式的解与原定解问题的真解的误差估计、差分格式的解当网格大小趋于零时是否趋于真解(即收敛性),等等。 偏微分方程边值问题的差分法物理上的定常问题其定解问题为各种边值问题, 即要求解在某个区域内满足微分方程,在边界上满足给定的边界条件。方程的差分解法可归结为选取合理的差分网格,建立差分格式 (1.1) 二、算法设计及流程图 2.1 算法设计 运用加权隐式格式求解常系数扩散问题(1.1),(1.6) 步骤1 输入初始值,确定加权隐式格式的参数; 步骤2 定义向量A,把初边值条件离散,得到,j=0,1,…,J的值存入向量A 步骤3 利用加权隐式差分格式由第n层计算第n+1层,建立相应线性方程组,求解并且存入向量A; 步骤4 计算到t=1,输出u 2.2 流程图 三、算法的理论依据及其推导 3.1 截断误差分析 常系数扩散问题(1.1) ,(1.6) 其中, 把(1.6)改写为便于计算的形式 = (1.7) 其中。 下面来求出差分格式(1.6)的截断误差,设u(x,t)是方程(1.1)的充分光滑的解,在处进行Taylor级数展开并经化简有 . 由此可以看出,当时,截断误差为。特别引起注意的是的情况,此时差分格式的截断误差是,即差分格式是二阶精度的。 3.2稳定性分析 我们用Fourier方法分析差分格式(1.6)式的稳定性。容易求出(1.6)式的增长因子为 根据von Neumann 条件是稳定性的充要条件,因此只验证,右端一定成立,只需要考虑左端的不等式,即 因,因此只需要,这就是差分格式(1.6)的稳定性限制。此条件也可以写得更明确些,即加权隐式格式稳定的条件是 ,当; 无限制, 当 四、数值结果及分析 采用加权隐式差分格式求解 其精确解为。取轴方向的网格步长分别为轴方向的网格步长为。计算在的近似解。 (1)取0.5 在 Matlab 命令窗口键入:u = peParabWegImp(1,0.5,0.01,11,0,1,0,0,10) 输出结果为: u = Columns 1 through 9 0 0.1203 0.2248 0.3076 0.3607 0.3790 0.3607 0.3076 0.2248 Columns 10 through 11 0 (2) 取0 在 Matlab 命令窗口键入:u = peParabWegImp(1,0,0.01,11,0,1,0,0,10) 输出结果为: u = Columns 1 through 10 0 -0.8819 5.5939 -7.9892 9.9736 -9.5804 9.9736 -7.9892 5.5939 -0.8819 Column 11 0 (3) 取1 在 Matlab 命令窗口键入:u = peParabWegImp(1,1,0.01,11,0,1,0,0,10) 输出结果为: u = Columns 1 through 10 0 0.1215 0.2310 0.3180 0.3738 0.3930 0.3738 0.3180 0.2310 0.1215 Column 11 0 (4)取0.2 在 Matlab 命令窗口键入:u = peParabWegImp(1,0.2,0.01,11,0,1,0,0,10) 输出结果为: u = Columns 1 through 9 0 0.1208 0.2228 0.3017 0.3549 0.3708 0.3549 0.3017 0.2228 Columns 10 through 11 0.1208 0 实验结果

文档评论(0)

gz2018gz + 关注
实名认证
文档贡献者

该用户很懒,什么也没介绍

1亿VIP精品文档

相关文档