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

维波动方程的有限差 分法 -------------------------------------------------------------------------- 作者 : _____________ -------------------------------------------------------------------------- 日期 : _____________ 学 生 实 验 报 告 实验课程名称 偏微分方程数值解 开课实验室 数统学院 学 院 数 统 年级 2013 专业班 信计 02 班 学 生 姓 名 学 号 开 课 时 间 2015 至 2016 学年第 2 学期 总 成 绩 教师签名 2 数学与统计学院制 开课学院、实验室： 数统学院 实验时间 ： 2016年 6 月 20 日 实验项目 实验项目类型 一维波动方程的有限差分法 名 称 验证 演示 综合 设计 其他 指导教 成 绩 师 曾芳 是 一．实验目的 通过该实验，要求学生掌握求解一维波动方程的有限差分法，并能通过计算机语言编程实现。 二．实验内容 考虑如下的初值问题： 2u 2u , x 0,1 , t 0, 2 t 2 x 2 u x,0 sin x, u x,0 0 (1) t u 0, t u 1,t 0, t 0,2 1．在第三部分写出问题（ 1）三层显格式。 2．根据你写出的差分格式，编写有限差分法程序。将所写程序放到第四部分。 3．取 h 0.1, 0.1h ，分别将 t 0.5,1.0,1.5,2.0 时刻的数值解画图显示。 该问题的解析解为 u x,t cos t sin x ，将四个时刻的数值解的误差画图显示，对数值结果进行简单的讨论。 三．实验原理、方法（算法）、步骤 1、三层显格式建立 由于题中 h 0.1, 0.1h ， x 0,1 , t 0,2 ，取 N 10, M 200 ，故令网比 r 0.1 ， xj j h, j 0,1,2,L 10 ， tk k ,k 0,1,L , 200 ，在 x j , tk 内网个点处，利用二 h 阶中心差商得到如下格式： ukj 1 2ukj ukj 1 o 2 ukj 1 2ukj u kj 1o h2 （ 2） 2 h2 略去误差项得到： 3 ukj 1 r 2ukj 1 2 1 r 2 ukj r 2 ukj 1 ukj 1 （3） 其中 j 1,2,L 9, k 1,2,L ,199 ，局部截断误差为 o 2 h2 。 对于初始条件 u x,0 sin x ，建立差分格式为： u0j sin xj sin jh , j 0,1,2,L 10 （ 4） 对于初始条件 t u x,0 0 ，利用中心差商，建立差分格式为： u1j 1 uj 0, 即 u1j =u j 1， j 0,1,2,L 10 （ 5） 2 对于边界条件 u 0,t u 1, t 0, t 0,2 ，建立差分格式为： u0k uNk 0, k 0,1,L , 200 （ 6） 将差分格式延拓使 k 0 为内点，代入（ 3）得到的式子再与（ 5）联立消去 u j 1 后整理得到： u1j 1 r 2u0j 1 1 r 2 u 0j 1 r 2 u0j 1 （ 7） 2 2 综上（ 3）、（ 4）、（ 6）、（ 7）得到三层显格式如下：（局部截断误差为 o 2 h2 ） k 1 r 2 k 2 1 r 2 k 2 k k 1 , j L 9, k 1,2, L ,199 u j u j 1 u j r uj 1 u j 1,2, 0 sin x j sin jh , j 0,1,2, L u j 10 （8） u1j 1 r 2u0j 1 1 r 2 u 0j 1 r 2u0j 1, j 1,2,L 9 2 u0k uNk 2 0,k 0,1,L , 200 其中 r 0.1。 h 四．实验环境（所用软件、硬件等）及实验数据文件 Matlab 三层显格式程序如下： %一维波动方程，三层显格式求解法 h=0.1;tau=0.1*h; r=tau/h;N=1/h;M=2/tau; x=0:h:1;t=0:tau:2; u=sin(pi*x);% 计算 t=0 时刻的 u 值 u(1,11)=0; for j=2:N u(2,j)=0.5*r^2*u(1,j+1)+(1-r^2)*u(1,j)+0.5*r^2*u(1,j-1); end %定义 x=0 边界上的数值 for k=1:M+1 u(k,1)=0; end 4 %定义 x=1 边界上的数值 for k=1:M+1 u(k,N+1)=0; end %迭代计算