二维波动方程的有限差分法.pdfVIP

学 生 实 验 报 告 实验课程名称 偏微分方程数值解 开课实验室 数统学院 学 院 数 统 年级 2013 专业班 信计02 班 学 生 姓 名 学 号 开 课 时 间 2015 至2016 学年第 2 学期 总 成 绩 教师签名 ,. 数学与统计学院制 开课学院、实验室： 数统学院 实验时间 ： 2016 年6 月20 日 实验项目 实验项目类型 二维波动方程的有限差分法 名 称 验证 演示 综合 设计 其他 指导教师 曾芳 成 绩 是 一．实验目的 通过该实验，要求学生掌握求解二维波动方程的有限差分法，并能通过计算机语言编 程实现。精品文档放心下载 二．实验内容 考虑如下的初值问题：       2u 2u 2u , x, y  0,1 , t  0,1.4   2 t 2 x2 y2 2 (1)         u x, y,0  sin x sin  y, t u x, y,0 0, x, y  0,1 u x, y, t  0,x,    精品文档放心下载 y  , t  0,1.4   1．在第三部分写出问题（1）三层显格式。 2 ．根据你写出的差分格式，编写有限差分法程序。将所写程序放到第四部分。精品文档放心下载 3 ．取h 0.1, 0.1h ，分别将t 0.5,1.0,1.4 时刻的数值解画图显示。谢谢阅读   4. 该问题的解析解为 u x, y, t cos 2t sin  x sin  y ，将四个时刻的数值解的误差画 图显示，对数值结果进行简单的讨论。谢谢阅读 三．实验原理、方法（算法）、步骤 ,. 网格划分 h  0.1h ，故 N 1 1.4 140 ，x ih, y jh, i, j 0,1,L ,10 ， 0.1, 10, M