学 生 实 验 报 告

实验课程名称 偏微分方程 数统学院

学 院 数 统 年级 2013 专业班 信计02班

学 生 姓 名 学 号

开 课 时 间 2015 至 2016 学年第 2 学期

总 成 绩教师签名

数学与统计学院制

开课学院、实验室： 数统学院 2016年 6月20日

实验项目

名 称一维波动方程的有限差分法实验项目类型验证演示综合设计其他指导教师曾芳成 绩是一．实验目的

通过该实验，要求学生掌握法并能通过计算机语言编程实现。如下的初值问题：

(1)

1．在第三部分写出问题(1)三层显格式。

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

3．取，分别将时刻的数值解画图显示。

4. 该问题的解析解为将时刻的数值解的误差画图显示，的讨论，，取，故令网比，，，在内网个点处，利用二阶中心差商得到如下格式：

(2)

略去误差项得到：

(3)

其中，局部截断误差为。

对于初始条件，建立差分格式为：

(4)

对于初始条件，利用中心差商，建立差分格式为：

(5)

对于边界条件，建立差分格式为：

(6)

将差分格式延拓使为内点，代入(3)得到的式子再与(5)联立消去后整理得到：

(7)

综上(3)、(4)、(6)、(7)得到三层显格式如下：(局部截断误差为)

(8)

其中。

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

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

%定义x=1边界上的数值

for k=1:M+1

u(k,N+1)=0;

end

%迭代计算开始,差分格式

for k=2:M

for j=2:N

u(k+1,j)=r^2*u(k,j+1)+2*(1-r^2)*u(k,j)+r^2*u(k,j-1)-u(k-1,j);

end

end

u(201,:)=zeros(1,11);

%计算k=201行的数值解

u2(201,11)=0;

for j=2:N

u2(201,j)=r^2*u(200,j+1)+2*(1-r^2)*u(200,j)+r^2*u(200,j-1)-u(199,j);

end

u=u+u2;

u=rot90(u,2);%将矩阵u旋转180度赋值于u

%作出图像

[x,t]=meshgrid(0:0.1:1,0:0.01:2);%划分网格

%作出数值解的函数图像

subplot(2,2,1);

mesh(x,t,u);

title('u(x,t)数值解的函数图像');

xlabel('x变量');

ylabel('t变量');

zlabel('u值');

%作出精确解的函数图像

subplot(2,2,2);

u1=cos(pi*t).*sin(pi*x);

mesh(x,t,u1);

title('u(x,t)精确解的函数图像');

xlabel('x变量');

ylabel('t变量');

zlabel('u值');

%作出t=0.5，1.0，1.5, 2.0时刻的绝对误差图像

subplot(2,2,3);

wucha=abs(u-u1);

x=0:h:1;

plot(x,wucha(51,:),'g*-');

hold on

grid on

plot(x,wucha(101,:),'ro-');

hold on

plot(x,wucha(151,:),'ks-');

hold on

plot(x,wucha(201,:),'mp-');

title('t=0.5，1.0，1.5, 2.0时刻的绝对误差函数图像');

xlabel('x变量'); ylabel('绝对误差