偏微分实验报告四张伟.pdf

重 庆 大 学 学 生 实 验 报 告 实验课程名称 偏微分方程数值解 开课实验室 数统学院 学 院 数统 年级 2012专业班 信计1班 学 生 名 张伟 学 号 开 课 时 间 2014 至 2015 学年第 2 学期 总 成 绩 教师签名 数学与统计学院制 开课学院、实验室： 数统学院 实验时间 ： 2015年 5月20 日 实验项目 实验项目类型 热传导方程的有限差分法 验证 演示 综合 设计 其他 名 称 指导教师 李茂军 成 绩 是 一．实验目的 通过该实验，要求学生掌握求解常系数热传导方程的有限差分法，并能通过计算机语言 编程实现。 二．实验内容 考虑如下的初值问题：  2 u  u  , x 0,1 , t 0,0.2  2     t x u x,0 sinx, x 0,1      (1) u 0, t u 1, t 0, t 0,0.2         2    t   问题的解析解为ux,t e sin x 。 1．在第三部分写出问题 （1）最简显格式，最简隐格式和C-N格式，并给出这些格式的局部 截断误差。 2．根据你写出的差分格式，编写有限差分法程序。将所写程序放到第四部分。 3．取 h0.1 ，将t0.2 时刻的所有网点的数值解，精确解和误差列表显示或画图显示在 第五部分。 4. 取h0.1, 0.5h2 ，将t0.2 时刻的所有网点的数值解，精确解和误差列表显示或画图 显示在第五部分。 5. 在同一种算法下，比较两种网格的数值结果。 要求：将程序放到实验结 部分。将电子版的实验报告发送到邮箱 limj@cqu.edu.cn,邮件 标题为：偏微分方程实验报告2-姓名-学号。纸质版的实验报告在下次上机实验时提交。 三．实验原理、方法（算法）、步骤 1、矩形网格剖分区域 h  jk, x ,t  jhk,  , j0,1, , N,k 0,1, ,M 取步长为 ，时间为 ，节点 表示点 j k    1 1 其中N  , M  ，显然,我们需要求解这 N1  M 1 个点对应的函数值。事实上由已     h  知初边界条件附近的点可直接得到,所以只要确定微分方程的解在其它点上的取值即可。沿 2 k 用记号u u x ,t 。 j j k 2、建立差分格式