ART算法代码程序

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