计算物理 ›› 1993, Vol. 10 ›› Issue (4): 413-421.

• 论文 • 上一篇    下一篇

常微分方程边值问题的高阶三对角OCI差分法

彭点云   

  1. 西南物理研究院, 成都 610041
  • 收稿日期:1992-04-07 修回日期:1993-01-26 出版日期:1993-12-25 发布日期:1993-12-25

HIGH ORDER TRIDIAGONAL OCI DIFFERENCE SCHEME FOR ORDINARY DIFFERENTIAL EQUATIONS WITH BOUNDARY VALUES

Peng Dianyun   

  1. Southwestern institute of Physics, Chengdu 610041
  • Received:1992-04-07 Revised:1993-01-26 Online:1993-12-25 Published:1993-12-25

摘要: 本文给出了二阶线性常微分方程两点边值问题(ODETPBVP)的高阶差分格式构造的基本思想,推导出六阶三对角OCI差分格式,并对端点有奇异性的方程进行了极限值处理,消去了奇异性,对边界层问题采用了非均匀网格上的六阶三对角OCI差分格式。通过大量的数值比较实验表明,这种高阶三对角OCI差分格式能很好地求解奇异性问题,固有不稳定性问题,奇异摄动问题,对生不稳定性问题和振荡性问题。

关键词: 三对角OCI差分格式, 极限值方法, 非均匀网格, 常微分方程两点边值问题, 追赶法

Abstract: This paper presents basic idea to construct high order difference schemes for linear two order ordinary differential equations with boundary values. A six order tridiagonal OCI difference scheme is deduced and the singularity at ends of interval are treated by means of limited values. A six order tridiagonal OCI difference scheme on non-uniform mesh is adopted for layer problems. A large number of numerical experimental results show that this high order tridiagonal OCI difference scheme can solve singular problems, inherent instability problems, dichotomous instability problems, singular perturbation problems and oscillating problems very well.

Key words: tridiagonal OCT difference scheme, limited value method, non-uniform mesh, ODETPBVP, double back substitution