热传导方程的差分格式.doc

该【热传导方程的差分格式 】是由【吴老师】上传分享，文档一共【21】页，该文档可以免费在线阅读，需要了解更多关于【热传导方程的差分格式 】的内容，可以使用淘豆网的站内搜索功能，选择自己适合的文档，以下文字是截取该文章内的部分文字，如需要获得完整电子版，请下载此文档到您的设备，方便您编辑和打印。热传导方程的差分格式第20页一维抛物方程的初边值问题分别用向前差分格式、向后差分格式、六点对称格式,求解以下问题: 在时刻的数值解,并与解析解。1差分格式形式设空间步长,时间步长,,网比.〔1〕向前差分格式 该问题是第二类初边值问题〔混合问题〕,我们要求出所需次数的偏微商的函数,满足方程和初始条件,及边值条件。 在相应区域光滑,并且在与边值相容,使问题有唯一充分光滑的解。 取空间步长,和时间步长,其中都是正整数。用两族平行直线和将矩形域分割成矩形网络,网络格节点为。以表示网格内点集合,即位于矩形的网点集合;表示闭矩形的网格集合;是网格界点的集合。向前差分格式,即〔1〕热传导方程的差分格式第2页,其中,以表示网比。〔1〕式可改写成如下:此格式为显格式。其矩阵表达式如下:〔2〕向后差分格式向后差分格式,即 〔2〕其中〔2〕式可改写成此种差分格式被称为隐格式。其矩阵表达式如下:〔3〕六点对称格式六点差分格式:〔3〕热传导方程的差分格式第20页将〔3〕式改写成其矩阵表达式如下:2利用MATLAB求解问题的过程对每种差分格式依次取,,,,用MATLAB求解并图形比拟数值解与精确解,用表格列出不同剖分时的误差。向前差分格式:::热传导方程的差分格式第4页::热传导方程的差分格式第20页热传导方程的差分格式第6页::热传导方程的差分格式第20页:向后差分格式:热传导方程的差分格式第8页热传导方程的差分格式第20页热传导方程的差分格式第10页