一维导热方程有限差分法matlab实现

长度曲线画图,其图像就是三维图的截面Рfunction [c,f,s]=pdefun(x,t,T,DuDx)%PDE方程函数Рc=100000;Рf=DuDx;Рs=0;Рfunction u0=icfun(x)%初始条件函数Рu0=0;Рfunction [pl,ql,pr,qr]=bcfun(xl,Tl,xr,Tr,t)%边界条件函数Рpl = Tl-100;Рql = 0;Рpr = Tr-100;Рqr = 0;Р将温度随时间变化情况用曲线表示。Р给出3000、9000、15000三个时刻的温度分布情况。РT3000 =100.0000 67.1058 39.8611 21.1973 10.9885 7.8279 10.9885 21.1973 39.8611 67.1058 100.0000РT9000=100.0000 83.4839 68.6032 56.8191 49.2705 46.6732 49.2705 56.8191 68.6032 83.4839 100.0000РT15000=100.0000 90.8310 82.5601 75.9972 71.7845 70.3330 71.7845 75.9972 82.5601 90.8310 100.0000Р根据数据分析,在同一个x点上温度随时间的增加而增加,但增幅变小。x-T图形仍为抛物线,但随着时间的增加,极值变小,图像变得平缓。Р用计算数据说明,空间、时间间隔对求解精度影响,并与有限差分法的计算结果做比较。Р调用前面做出来的真实值,跟pdepe做出来的值计算误差,再与有限差分法的误差比较。用pdepe函数求的误差远小于有限差分法,所以pdepe函数法更精确。Р用计算数据说明,有无稳定性要求,为什么?若有,如何对求解精度的影响。Р不知道这个pdepe函数的稳定性要用什么检验。傅里叶级数检验不适用。