%SIRT?算法，可以用于CT重建
%做得比較簡單，只用了四個像素，六條投影路徑，具體情況見附圖
%像素的實際值是A=[97?50?;66?90]
W=[0?0.1?0.1?0;0.1?0.1?0?0;0?0?0.1?0.1;0.1?0?0?0.1;0?0.1?0?0.1;0.1?0?0.1?0];
%W是6*4的參數矩陣，元素代表某條路徑在某個像素中的長度，這里做的很簡略，把所有經過的長度都記為0.1，實際上應該好好算算的
P=[12?;15?;17?;20;?15?;17];
%TOF是投影值
count=100;
ng=4;%像素數
delta=1e-9;

%以下是SIRT迭代
X=ones（ng1）;%slowness
Xn=ones（ng1）;
Lie=sum（W）;%W的每列之和，是1*ng的行向量
k=count;
while?k>0
????P1=W*X;
????for?j=1:ng
??????A=W（:j）.*（P-P1）;
??????B=sum（A）;
??????C=0.0001+sum（Lie2）;
??????Xn（j）=B/（C*Lie（j））;
??????X（j）=X（j）+Xn（j）;
????end
????e=sqrt（sum（Xn.^2））;
????if?e????????k=0;
????else
????????k=k-1;
????end
end
X

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件????????778??2014-11-05?11:10??SIRT4grids.m

?????文件?????972054??2014-11-05?11:07??投影示意圖.bmp

-----------?---------??----------?-----??----

???????????????972832????????????????????2