3、40(0.5,0.2)2.01342.01383.4361e-0460(0.5,0.3)2.22512.22554.1249e-0480(0.5,0.4)2.45912.45964.6788e-04100(0.5,0.5)2.71782.71835.2148e-04120(0.5,0.6)3.00363.00425.7794e-04140(0.5,0.7)3.31953.32016.3932e-04160(0.5,0.8)3.66863.66937.0677e-04180(0.5,0.9)4.05444.05527.8118e-04200(0.5,1.0)
4、4.48084.48178.6337e-04表3.2部分结点处数值解、精确解和误差的绝对值(h=l/10,T二1/100)k(x,t)数值解精确解1精确解-数值解11(0.5,0.01)1.66521.66536.8968e-052(0.5,0.02)1.68191.68201.3932e-043(0.5,0.03)1.69871・69892.1107e-044(0.5,0.04)1.71571.71602.8425e-04L□(0.5,0.05)1.73291.73333.5888e-046(0.5,0.06)1.75041.75072.7944e-04
5、7(0.5,0.07)1.76731.76839.7606e-048(0.5,0.08)1.78751.78601.4242e-039(0.5,0.09)1.79701.80407.0151e-0310(0.5,0.10)1.84231.82212.0150e-0211(0.5,0.11)1.77541.84046.5060e-0212(0.5,0.12)2.05461.85891.9565e~0113(0.5,0.13)1.28701.87765.9063e-0114(0.5,0.14)3.65191.89651.7554e+0015(0.5,0.15)
6、-3.27761.91555.1931e+0016(0.5,0.16)17.21141.93481.5277e+0117(0.5,0.17)-42.83711・95424.4791e+0118(0.5,0.18)132.94771・97391.3097e+02表3.3取不同步长时数值解的最大误差(“1/2)hTEgfji,t)Ee(2h,4T)/E^(htt)1/101/2008.6337e-04*1/201/8002.1748e-043.9698e+001/401/32005.4366e-054.0004e+001/801/128001.3591e-05
7、4.0001e+00图3.3I时的数值解曲线(h=1/10J=1/200)与粽确解曲线图3.453.x10'31：•0.9•0.8・t■时不同步长数值解的谋差曲线(r=1/2)h=1/10,T=1/200e—h=1/20,T=1/800——h=1/40,T=1/3200o3a7o654•••ooo_(L・x=xxLxn_2ao不同步长数值解的误差曲面(r=1/2)-(rx)xxrx)三Ih=1/10tT=1/200I]h^17201=1/8000.2四.小结从表3.3可以看出，当空间步长缩小到原来的1/2,时间步长缩小到原来的1/4时，最大误差约缩小到原来
8、的1/4.五.程序表3.1clc,clearformatshort
