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