偏微分课程课件9_椭圆型方程的有限差分方法(I)_第1页
偏微分课程课件9_椭圆型方程的有限差分方法(I)_第2页
偏微分课程课件9_椭圆型方程的有限差分方法(I)_第3页
偏微分课程课件9_椭圆型方程的有限差分方法(I)_第4页
偏微分课程课件9_椭圆型方程的有限差分方法(I)_第5页
已阅读5页,还剩28页未读 继续免费阅读

下载本文档

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

文档简介

1、1第五章 椭圆型方程的差分方法 (一) Poisson方程(二)差分格式的性质(三)边界条件的处理(四)变系数方程(五)双调和方程(六)特征值问题2(一) Poisson方程区域以及边界离散4内点边界点 返回5 返回1.五点差分格式利用Taylor级数展开有则8 由于上差分方程中只出现u在(i,j)及其四个临点上的值, 故称为五点差分格式Possion方程第一边值问题的差分逼近为改写为代数方程, 仅考虑解:11 先仅下标 j 小的未知量放前面, 相同的 j 的再按 i 从小到大顺利排列11按自然顺序排列网点(i,j)按自然顺序排列网点(i,j)定义向量于是差分方程为:12分析系数矩阵H对于第一

2、个结点(1,1), 未知节点分析系数矩阵H对于第二个结点(2,1), 分析系数矩阵H对于前I个结点, 系数矩阵H为 五点差分格式等价于B,I为I维方阵, H为I*J维方阵对于第一个结点, 即例:17同理, 对于另外三个结点, 分别有联立, 有线性代数方程组18计算机实现五点差分格式的系数矩阵算法程序返回19五点差分格式的系数矩阵对possion方程第一边值问题建立五点差分格式将双下标数列ui,j按自然顺序排成单下标uk,那么20这样有代数方程组Hu=g, 其中21这样有代数方程组Hu=g, 其中H是稀疏矩阵, 每一行最多5个数不为0H是对角占优矩阵H是对称正定课堂练习23 微分方程差分格式真解u=u(x,y)真解u=un实用性分析 相容性收敛性 稳定性椭圆型方程收敛性分析区别于双曲与抛物方程对于双曲与抛物方程:24 相容性所以五点差分格式是相容的.2.九点差分格式27四阶精度3 极坐标形式的差分格式若求解域是圆环、环形域或扇形域,则采用极坐标是方便的,此时Poisson方程形如任意点用中心差商公式分别取等步长下面用有限体积法获得在点 的差分方程用r乘方程上积分并对得到的方程在用中心差商代替上面的微商,并用 除得利用中

温馨提示

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

评论

0/150

提交评论