地球物理反演概论 射线过 单元 射线不过 单元4.3 SIRT 联合代数 重建 技术为提高ART算法 精度,对其进 一步改 进。将过第j 个单元格的所 有射线 的修正 量都计 算出来,然后 取平均值作为 模型参 数的修 正量1.SIRT 联合代数 重建技术 基本原理 1 2 3 4 5 6 7 8 98t7t6t5t4t3t2t1t 第j 个单元格有几条射 线通过?答:等于 的 第j列非零 元素个数(1)行元素几个非零,射线过几个单元(2)列元素几个非零,单元格过几射线2.SIRT 迭代 重建算法流程(1)令。(2)for j=1,M,计算 第j 个单元 格经过 的射线总条数(中非零元素 个数)。射线过 单元 射线不过 单元(8)判断解是否 收敛,不收敛,令(4)for i=0,1,N,计算i 射线的 实际长 度(的非零元素之和乘以)。(5)令。SIRT算法优点:精度 高于ART 算法。缺点:计算量大于ART 算法。(3)for i=0,1,N,计算 第i 条射 线经过 的单元格总数(中非零 元素 个数)。(6)for i=0,1,N-1;j=1,2,M,计 算(7)forj=1,2,M,令