对比抛物Radon正变换几种矩阵求解方法

• 作者: 姜岩 王维红
• 来源期刊: 物探化探计算技术
• 学科: 地球科学
• 关键词: 抛物Radon正变换 共轭梯度法 Levinson递推法 Cholesky分解法 矩阵求逆
• 年，卷（期）: 2007,（6）
• 所属期刊栏目: 计算技术与计算机应用
• 页码范围: 555-559
• 页数: 5页
• 分类号: P631.4
• 字数: 4202字
• 语种: 中文
• DOI: 10.3969/j.issn.1001-1749.2007.06.020

摘要

抛物Radon变换法(Parabolic Radon Transform)在地震资料处理中有广泛的应用.PRT可对不同频率的地震数据解耦处理,这一特点使得抛物Radon变换的计算效率比双曲Radon变换有数量级上的提高.在频率域求解时,需要对每一个频率成份求解同样大小的线性方程组.求解抛物Radon正变换的计算方法主要有Levinson递推法、共轭梯度法、Cholesky分解法和直接矩阵求逆法.最小平方抛物Radon正变换所形成的矩阵具有Toeplitz结构,可采用Levinson递推法进行计算.高分辨率抛物Radon正变换所形成矩阵的Toeplitz结构被破坏,一般采用共轭梯度法或Cholesky分解法进行求解.这里详细推导了复Toeplitz矩阵的Levinson递推算法,并分别对求解方程的四种方法进行了讨论,最后给出抛物Radon正变换求解的数值算例,并对所给出的四种方程求解方法的计算效率及计算精度进行了对比.

作者信息

序号	姓名	单位	发文数	被引次数	H指数	G指数
1	姜岩		22	259	8.0	15.0
3	王维红	大庆油田有限责任公司勘探开发研究院	17	244	9.0	15.0