三维地电场数值计算中的LANCZOS迭代过程及其改进算法
详细信息
- 责任者：宛新林^c安徽建筑工业学院^d1966^e博士 ; 席道瑛^c中国科学技术大学
- 刊名：物探化探计算技术
- 出版年：2009
- 卷期：31^b3总125
- 关键词：地电场 ; 三维模型20^a正演问题 ; 数值模拟 ; 电法勘探30^a线性方程组
- 页码：197～201
- 索书号：P630.6/363-14
- 图表：^c表1参9插图
- 分类号：2100
- 文摘期号：2009/11
摘要
针对三维地电场正演数值计算过程中形成的超大规模稀疏线性方程组,在分析此类线性方程组的一般解法基础上,着重阐述一种适宜求解此类方程组的Lanczos迭代过程与算法原理。同时,当地下介质的电性差异较大时,形成系数矩阵A的条件数就很大,可对算法进行适当改进。讨论采用不完全Cholesky分解方法进行预条件处理,经过条件数改善后,形成新的线性方程组系数矩阵,就会变为一个近似的单位矩阵。经改进后的Lanczos算法,将提高数值计算稳定性,从而加快迭代收敛速度,为提高反演质量提供基础。