两点边值问题的差分方法及快速求解.pdfVIP

华南师范大学学报 (自然科学版) 2010年 11月 J0URNALOFSOUTH CHINANORMALUNIVERSn 2010年第4期 NOV．20lO (NATURALSCIENCEEDIT／ON) No．4，2010 文章编号：1000—5463(2010)04—0027—04 两点边值 问题的差分方法及快速求解 赵秉新 ，郑来运 (1．宁夏大学数学计算机学院，宁夏银川750021；2．宁夏大学机械工程学院，宁夏银川750021) 摘要：利用Richardson外推法，将求解两点边值问题的二阶中心格式和高精度差分格式推广到更高精度．通过数值算 例，提出了使用Excel电子表格软件进行快速求解的方法，并与编制的计算程序的计算结果进行了对比．结果表明， 使用Excel软件求解，能够快速得到满意的结果． 关键词：边值问题；差分方法；快速求解 中图分类号：O175．8；O241．8 文献标志码 ：A 在应用科学与工程技术领域中经常遇到常微分 解区间[a，6]进行网格剖分 将区间[a，6]作m等 方程的边值问题¨ J，由于边值问题解的存在性和 分，记h=(b—a)／m，=口=ih(0≤≤m)， = 1 唯一性远比初值问题复杂，因此求解边值问题不论 0≤≤m}．称h为网格步长，i为网格结点， 为网 在理论上还是数值方法上都要比初值问题困难．目 格，称定义在网格 上的函数为网格函数 设V= 前采用的数值解法有打靶法及其改进方法 和有 l0≤≤m}为 上的网格函数，记 限差分法 -s]，KANTH等 使用三次样条方法获得 1 1 二阶两点边值问题的近似解．求解该问题的方法虽 8xvi+l儿音( )，t，音(xUi+l一 -)， 然较多，但就其本质来说，最终都要编程实现，这些 1 1 并非所有工程技术人员都能做到，而利用已有的软 D (-一)／i)，D一 音(--gti_I)， 件或工具包可极大方便科技工作者 。。．本文利用 = = 6 =—— 一 =一h(龇6 …／2：—6 f_。I／2 ‘ Richardson外推法分别提高了二阶中心差分格式 (thesecond—orderCentralDifferenceSchcme--CDS) 应用Taylor公式能迅速得到二阶中心差分格式 和四阶精度差分格式 (Fourth—orderaccuracyDiffer— (CDS)： enceScheme--FDS)H 的精度；给出了一种利用 +p()gui+g( u= i)(1≤≤，n一1)，(3) Excel软件的自动迭代功能代替编程进行求解的新 Ⅱ0： ， = (4) 方法，同时编制FORTRAN程序，将比较和分析二者 式 (3)和式 (4)中共有m+1个方程，m+1个待 的计算结果． 定量 。， 一， ，局部截断误差具有二阶精度． 1．2 四阶紧致差分格式 1 差分方法 文献 [7]和文献[8]从中心差分公式出发，拆诸 差分余项反向修正，对两点边值问题得到了一种四 考虑如下定解问题： 阶精度的差分