一阶格式和二阶格式误差对比.pdfVIP

导数边界值问题一阶格式和二阶格式的误差比较 陈永鑫 （一）问题 应用差分格式计算如下两点边值问题对比一阶格式和二阶格式误差  x u (x)u(x) e (sinx 2cosx)   u(0) 1,u() e （二）求解方法描述 x 首先该方程有解析解u(x) e sinx ，有了这个以后可以与我们 接下来的数值解进行比较，以此来判断求解方法是否精确。 A 一阶格式  u 2u  u  i1 i i1  qu f  h2 i i i  (u u ) h,u u h ,1 i m  1 0 m m1 h2 将 乘到该方程的两端可得到：  2 2 u [2h q(x )]u u h f (x )  i1 i i i1 i (u u ) ,u u h,1 i m  1 0 m m1 π 因为在这个题目中q(x)=1， 1, e 1 0,2 0所以 ， 将这个线性方程写成矩阵形式：  1 1  u0   h       2 2 1 2h 1  u1   h f (x )  1                2 2  1 2h 1um1  h f (x )  m1     π  1 1 u he   m    即可以化成方程Au=b，若A 非奇异则u=inv(A)*b，则可以得到解。 所得到u ,u ,u u 即为以h 为网格步长得到的网格结点的数值 0 1 2 m 解。 （三）程序的思路 首先在matlab 中建立两个函数文件