一维抛物线偏微分方程数值解法(2)(3页).doc

《一维抛物线偏微分方程数值解法(2)(3页).doc》由会员分享，可在线阅读，更多相关《一维抛物线偏微分方程数值解法(2)(3页).doc（3页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

-一维抛物线偏微分方程数值解法(2)-第 3 页 一维抛物线偏微分方程数值解法（2）上一篇文章请参看 一维抛物线偏微分方程数值解法（1）解一维抛物线型方程(理论书籍可以参看孙志忠：偏微分方程数值解法）Ut-Uxx=0, 0x1,0t0) U(x,0)=ex, 0=x=1,U(0,t)=et,U(1,t)=e(1+t), 0t title(误差曲面); plot(t,e)误差较之前的欧拉向前差分格式 增长了两倍u p e x t=pwxywxh(0.1,0.05,10,20);plot(t,e)u p e x t=pwxywxh(0.01,0.05,100,20);plot(t,e)u p e x t=pwxywxh(0.01,0.01,100,100);plot(t,e)u p e x t=pwxywxh(0.01,0.005,100,200);plot(x,e)u p e x t=pwxywxh(0.01,0.005,100,200);plot(t,e)u p e x t=pwxywxh(0.005,0.005,200,200);plot(x,e)X=1时，出现了误差？? 不是边界条件吗？不能理解 这方法还是比前一种方法误差大呀不过可以随便改变时间、空间步长