两点边值问题的差分求解--课程设计_第1页
两点边值问题的差分求解--课程设计_第2页
两点边值问题的差分求解--课程设计_第3页
两点边值问题的差分求解--课程设计_第4页
全文预览已结束

下载本文档

版权说明：本文档由用户提供并上传，收益归属内容提供方，若内容存在侵权，请进行举报或认领

文档简介

两点边值问题的差分求解摘 要：本文给出了一种二阶微分方程的差分格式，通过 Taylor 级数展开给出其截断误差， 数值例子验证了其误差。 关键词：有限差分方法，两点边值问题，泰勒级数展开1. 引言考虑二阶常微分方程边值问题：最简型：变系型守恒型其中上的连续函数，和为给定的常数。2. 格式推导 网格剖分（区域剖分）取空间步长为，其中N为整数。用直线，将区间分割成网格，网格节点为。 截断误差（Taylor级数展开）在处有：将，在点Taylor展开令，由，得：(6)式+(7)式，得：由此可得：其中表示方括号内的函数在点取值。于是方程（1）可以写成：其中 即为截断误差。 差分格式（差分方程+定解条件）当h是足够小时，是h的二阶无穷小量，舍去可得差分方程：用数值解代替精确解，并给定方程的边值条件，可得如下差分格式： 格式求解可将差分格式改写为矩阵形式：其系数矩阵为：3. 数值实验以如下方程为例：该问题的精确解为：将作N等分，Matlab 程序求出 4 个结点处的精确解和取不同步长所得的数值解如表1所示；表2 给出了这些结点处取不同步长时所得数值解和精确解差的绝对值；表3给出了取不同步长时所得数值解的大误差；图1给出取不同步长时所得数值解的误差曲线图。hx/8/43/8/2/80.44581.21052.13884.7574/160.09981.48552.61384.8060/320.0948 1.57092.74964.9448/640.46571.76452.94514.3813表1 部分结点处的精确解和取不同步长时所得的数值解0.20.40.60.8/80.75140.59200.91270.7778/160.09740.31710.43760.7291/320.10240.23160.30180.5904/645.8229e-040.03800.10630.1539精确解0.19721.80263.05154.5351表2 不同步长下部分节点处数值解的误差的绝对值h/80.6914/160.4454/320.6214/640.8541表3 不同步长时数值解的最大误差图1 取不同步长时所得数值解的误差曲线图4. 结论由表1，表2对比可知，数值解与精确解高度越接近，由图1可知步长对误差的影响，因此我们可以认为本文接近二阶常微分方程的查分格式。参考文献：1 孙志忠. 偏微分方程数值解法M. 北京:科学出版社, 2005.2 李荣华. 偏微分方程数值解法M. 北京:
人人文库> 全部分类> 应用文书 > 事务文书

温馨提示

  • 1. 本站所有资源如无特殊说明，都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等，如果需要附件，请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸，网页内容里面会有图纸预览，若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间，仅对用户上传内容的表现方式做保护处理，对用户上传分享的文档内容本身不做任何修改或编辑，并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容，请与我们联系，我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

最新文档

评论

0/150

提交评论