matlab有限差分一维导热,一维导热方程-有限差分法-matlab实现.docx

此文档收集于网络，如有侵权请联系网站删除

第五次作业(前三题写在作业纸上)

一、用有限差分方法求解一维非定常热传导方程，初始条件和边界条件见说明.pdf文件，热扩散系数α=const，

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%输出温度矩阵