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

一维导热方程 有限差分 北京化工大学 图形处理 matlab

10

3. 给出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图形仍为抛物线，但随着时间的增加，极值变小，图像变得平缓。

4. 用计算数据说明，空间、时间间隔对求解精度影响，并与有限差分法的计算结果做

比较。

调用前面做出来的真实值,跟pdepe做出来的值计算误差，再与有限差分法的误差比较。用pdepe函数求的误差远小于有限差分法，所以pdepe函数法更精确。

5. 用计算数据说明，有无稳定性要求，为什么？若有，如何对求解精度的影响。 不知道这个pdepe函数的稳定性要用什么检验。傅里叶级数检验不适用。