偏微分实验报告五.docx

重庆大学

学生实验报告

实验课程名称

开课实验室 数统学院

学院数统年级2012专业班信计1班

学生姓名张伟学号开课时间2014至2015学年第2学期

总成绩

教师签名

数学与统计学院制

开课学院、实验室：数统学院

实验时间：2015年6月3日

实验项目

名称

一维双曲型方程的有限差分法

实验项目类型

验证

演示

综合

设计

其他

指导教师

李茂军

成绩

是

一.实验目的

通过该实验，要求学生掌握求解一维双曲型方程的有限差分法，并能通过计算机语言编程实现。

二.实验内容

考虑如下的初值问题:

du du 八

——+——=0,dt dx

u(x,0)=1,x0

u(x,0)=

0.5,x=0

0,x0

1.在第三部分写出问题（1）左偏心格式，Lax-Wendroff格式，并给出这些格

式的局部截断误差。

2.根据你写出的差分格式，编写有限差分法程序。将所写程序放到第四部分。

3.取h=0.1具=0.5h，分别将t=0.05,0.1,0.15,0.2时刻的数值解uk(j=0,±1,±2,...,±10)画

图显示。

4.对数值结果进行简单的讨论。

要求：将电子版的实验报告发送到邮箱limj@cqu.edu.cn,邮件标题为：偏微分方程实验报告5-姓名-学号。纸质版的实验报告在下次上机实验时提交。

实验原理、方法（算法）、步骤

1、左偏心格式

对方程⑴，利用差商代替微商的思想，采取先后差商替代告令土（网比）得:

(2)uk+1—uk uk—uk 。

(2)

t h

局部截断误差为o(r+h)，只能取正值，作以下平移：uk+i=(1-r)uk+ru

局部截断误差为o(r+h)，只能取正值，作以下平移：

uk+i=(1-r)uk+ruk (3)

,5：这里j?的前四个点是为了求解uk+1设定的点，因此最后j

其中j=1,2,,25,k=1,2,

得到的应该为一个5x21的矩阵A，从下到上分别代表k=0,1,2,3,4，从左到右分别代表

??????

j=5,6,,25，即分别表示t=0,0.05,0.1,0.15,0.2的时刻,xe(—1,1)的数值。(当然也可以利用

赋值思、想，将左边需用到的点赋值于1，得到的就是一个矩阵形式，不过这样需人为干涉)综上得到左偏心差分格式如下：

(1—r)uk+ruk,j=1,2,,25,k=1,2,,5

(4)1,j=14j—

(4)

0.5,j=15

0,j=16,17,,25

2、Lax-Wendroff格式

该格式由Taylor展开求得，即

(5)uk+1=—(r2—r)uk+(1—r2)uk+—(r2+r)uk

j2 j+1 j2 j-1

(5)

局部截断误差为o(T2+h2)，结合问题3，其中j=0,±1,,±10,k=0,1,2,3,4，考虑到下

标只能取正值，作以下平移： …

其中j=1,2,,29,k=1,2,,5，这里/的前四个点以及后四个点都是为了求解uk+1设定j

的点，因此最后得到的应该为一个5x21的矩阵A，从下到上分别代表k=0,1,2,3,4，从左到

??????

右分别代表j=5,6, ,25，即分别表示t=0,0.05,0.1,0.15,0.2的时刻,xe(—1,1)的数值。

综上，得到Lax-Wendroff差分格式如下：

???

uk+1=

j2 j+1】,j=1,2,140.5,j=150,j=16,17,,29—(r2—r)uk+(—r2)uk+—

uk+1=

j

2 j+1

】,j=1,2,14

0.5,j=15

0,j=16,17,,29

u1

u1=

j

实验环境(所用软件、硬件等)及实验数据文件

Matlab

1、左偏心格式程序代码

clc;

h=0.1;%定义步长

tau=0.5*h;a=1;r=tau/h;%定义网比

%分析得到欲求t=0.2时刻的解，需向左平移4个点，且下标采用正值，故取25个点%本程序未采取矩阵形式求解，主要是格式本身带来的麻烦，以及便于迭代分别画图%k=0,即t=0时刻的数值解u_0forj=1:25

ifj=1j=14

u(j,1)=1;

end

ifj==15

u(j,1)=0.5;

end

ifj=16j=25

u(j,1)=0;

end

end

u1=u;

u_0=u(5:25,1);

%k=1,即t=0