clc;
clear?all;
close?all;
N?=?180;?%圖像大小
N2?=?N^2;
I?=?phantom（N）;%?產生頭模型圖像
theta?=?linspace（0180181）;
theta?=?theta（1:180）;
%%?=?=?=?=?=?=?產生投影數據?=?=?=?=?=?=?%%
P_num?=?260;?%探測器通道個數
P?=?ParallelBeam（theta?N?P_num）;?%產生投影數據
%P?=?radon（Itheta）;
%%?=?=?=?=?=?=?獲取投影矩陣?=?=?=?=?=?=%%
delta?=?1;%?網格大小，角度增量
[W_indW_dat]?=?SystemMatrix（thetaNP_numdelta）;
%%?=?=?=?=?=?=?進行ART迭代?=?=?=?=?=?=?%%
F?=?zeros（N21）;??%初始圖像向量
lambda?=?0.25;??%松弛因子
c?=?0;??%迭代計數器
irt_num?=?5;??
while（c????for?j?=?1:length（theta）
????????for?i?=?1:1:P_num
????????????%?取得一條射線所穿過的網格編號和長度
????????????u?=?W_ind（（j-1）*P_num?+?i:）;??%?編號
????????????v?=?W