一维热传导方程数值解法及matlab实现分离变量法和有限差分法.pdfVIP

一维热传导方程数值解法及matlab实现

分离变量法和有限差分法

一维热传导方程的Matlab解法：分离变量法和有限差分

法。

问题描述：

本实验旨在利用分离变量法和有限差分法解决热传导方程

问题，并使用Matlab进行建模，构建图形，研究不同情况下

采用何种方法从更深层次上理解热量分布与时间、空间分布关

系。

实验原理：

分离变量法：

利用分离变量法，将热传导方程分解为两个方程，分别只

包含变量x和变量t，然后将它们相乘并求和，得到一个无穷

级数的解。通过截取该级数的前n项，可以得到近似解。

有限差分法：

利用有限差分法，将空间和时间分别离散化，将偏导数用

差分代替，得到一个差分方程组。通过迭代求解该方程组，可

以得到近似解。

分离变量法实验：

采用Matlab编写代码，利用分离变量法求解热传导方程。

首先设定x和t的范围，然后计算无穷级数的前n项，并将其

绘制成三维图形。

代码如下：

matlab

x=0:0.1*pi:pi;

y=0:0.04:1;

x。t]=meshgrid(x。y);

s=0;

m=length(j);

fori=1:m

s=s+(200*(1-(-1)^i))/(i*pi)*(sin(i*x).*exp(-i^2*t));

end

surf(x。t。s);

xlabel(x)。XXX(t)。zlabel(T);

title(分离变量法（无穷）);

axis([0pi010100]);

得到的三维热传导图形如下：

有限差分法实验：

采用Matlab编写代码，利用有限差分法求解热传导方程。

首先初始化一个矩阵，用于存储时间t和变量x。然后计算稳

定性系数S，并根据边界条件和初始条件，迭代求解差分方程

组，并将其绘制成三维图形。

代码如下：

matlab

u=zeros(10.25);

s=(1/25)/(pi/10)^2;

fprintf(稳定性系数S为:

);

disp(s);

fori=2:9

u(i。1)=100;

end

forj=1:25

u(1.j)=0;

u(10.j)=0;

end

forj=1:24

fori=2:9

u(i。j+1)=s*u(i+1.j)+(1-2*s)*u(i。j)+s*u(i-1.j);

end

end

x。t]=meshgrid(1:25.1:10);

surf(x。t。u);

xlabel(t)。XXX(x)。zlabel(T);

title(有限差分法解);

得到的热传导图形如下：

改写后的分离变量法实验：

采用Matlab编写代码，利用分离变量法求解热传导方程。

首先设定x和t的范围，然后计算无穷级数的前100项，并将

其绘制成三维图形。

代码如下：

matlab

x=0:0.1*pi:pi;

y=0:0.04:1;

x。t]=meshgrid(x。y);

s=0;

fori=1:100

s=s+(200*(1-(-1)^i))/(i*pi)*(sin(i*x).*exp(-i^2*t));

end

surf(x。t。s);

xlabel(x)。XXX(t)。zlabel(T);

title(分离变量法);

根据给定的坐标轴范围，使用分离变量法和有限差分法得

到的热传导图形基本相同，说明两种方法都适用于本题。

使用分离变量法时，可以考虑两种情况。第一种情况是取

无穷项，此时可以通过加上disp(s(:,6))来得到第六列（即

x=pi/2）处温度随时间的变化情况。根据计算结果，当温度低

于50度时