SIRT算法其思想在于利用通过该像素的全部射线，其迭代过程对图像每个像素的更新量是对所有投影线的修正按照贡献因子取加权平均，然后反投影得到。与ART每条投影线都对图像更新一次不同，SIRT算法综合了所有投影线的贡献，可以避免一条投影线上的误差对重建结果带来过大影响，因而可以有效抑制重建图像中的噪声。 2、算法实现步骤 （1）对第 条射线，计算估计投影值 （2）计算实际投影与估计投影的误差 （3）反投影值 其中， 是所有投影角度下光线的集合. （4）对第 个像素点的值进行修正 （5）将上一轮的结果作为初值，重复（1）~（4）的过程，直到达到收敛要求或指定的迭代次数。 因此，SIRT算法的迭代公式为： 其中， 是松弛因子， 是迭代次数。 3、优缺点分析 由于SIRT算法对所有投影线的修正量进行了加权平均，显著地降低了迭代的收敛速度。另一方面，对每个像素更新时，需要计算好所有投影线的贡献，因此在实际计算中需要对各个投影线的贡献量进行存储，存储量至少比ART算法多一倍。因此，SIRT算法具有更好的稳定性，但是收敛速度慢、存储容量大，此两点成为影响其应用的主要问题。

地震走时层析成像实验 单元数：9  12 单元边长：3.0  5.0m 参数：速度（km/s) 实验一：单边激发单边接收 (左发右收) 数据：12  12

SIRT方法是一种基本的地球物理反演方法，是与SVD,LSQR反演方法一类的基本反演方法，一般用于地震正演数据的反演数据，从而得出地下地层的速度参数，最终得出地下地质体情况

基于GPU图形加速卡，编程实现多线程同步加速，适用于各类图像重建算法，尤其是迭代类计算量比较大的算法。ART,SIRT,SART 等。

本代码是基于matlab的层析反演SIRT算法过程，能够通过matlab来帮助大家实现整个过程

解线性方程组art and sirt 算法

基本的SIRT重建算法，已调试，可以运行。 做的比较简单，有兴趣的话可以在这个的基础上加点东西。欢迎交流讨论~

CT算法MATLAB，全部，含源代码，及测试数据，及使用说明书； 对CT算法学习极有帮助。。。SIRT, ART ,MART

本程序主要对SIRT算法进行Matlab实现，共包含三个m文件。运行时，首先执行GenerateSIRTSystemMatrix.m生成系统矩阵A。然后，每需要重建一副图像时运行一次SIRTOnce.m或SIRTOnceSTD.m。其中，SIRTOnce.m对应了固定步长的SIRT算法，SIRTOnceSTD.m为采用最速下降原理对迭代步长进行优化以后的SIRT算法，并采用FBP算法进行对比重建。