matlab有限差分法编程波导_一维导热方程-有限差分法-matlab实现

2

2

T

T

t

x













1.

用

Tylaor

展开法推导出

FTCS

格式的差分方程

2.

讨论该方程的相容性和稳定性，并说明稳定性要求对求解差分方程的影响。

3.

说明该方程的类型和定解条件，如何在程序中实现这些定解条件。

4.

编写

M

文件求解上述方程，并用适当的文字对程序做出说明。

(部分由网络搜索得

到，添加，修改后得到。

)

function rechuandaopde

%

以下所用数据，除了

t

的范围我根据题目要求取到了

20000

，其余均从

pdf

中得来

a=0.00001;%a

的取值

xspan=[0 1];%x

的取值范围

tspan=[0 20000];%t

的取值范围

ngrid=[100 10];%

分割的份数，前面的是

t

轴的，后面的是

x

轴的

f=@(x)0;%

初值

g1=@(t)100;%

边界条件一

g2=@(t)100;%

边界条件二

[T,x,t]=pdesolution(a,f,g1,g2,xspan,tspan,ngrid);%

计算所调用的函数

[x,t]=meshgrid(x,t);

mesh(x,t,T);%

画图，并且把坐标轴名称改为

x

，

t

，

T

xlabel('x')

ylabel('t')

zlabel('T')

T%

输出温度矩阵

dt=tspan(2)/ngrid(1);%t

步长

h3000=3000/dt;