写一个适用于360扫描方式的CT重建系统矩阵计算matlab程序

以下是一个适用于360扫描方式的CT重建系统矩阵计算的Matlab程序的示例：

% CT系统参数
N = 512;            % 图像大小
D = 0.5;            % 探测器宽度
S = 1.0;            % 探测器数量
T = 0.5;            % X射线管到探测器距离
R = 0.5;            % 旋转半径
theta = 0:1:359;    % 角度范围

% 计算系统矩阵
A = zeros(N*S,length(theta)*N);
for i = 1:length(theta)
    for j = 1:N*S
        for k = 1:N
            x = (k-1)*D-N*D/2;
            y = (j-1)*D;
            z = T+R*cosd(theta(i));
            s = R*sind(theta(i));
            l = sqrt(x^2+(y-s)^2+z^2);
            cosphi = (z/l);
            sinphi = ((y-s)/l);
            A((i-1)*N+j,(i-1)*N+k) = 1/l*cosphi;
            A((i-1)*N+j,length(theta)*N+(j-1)*S+ceil(s/D)) = -1/l*cosphi*sinphi*D;
        end
    end
end

% 显示系统矩阵
figure;imagesc(A);colormap(gray);title('CT系统矩阵');

这个程序可以计算360度的CT系统矩阵，并且假设探测器宽度为0.5，探测器数量为1.0，X射线管到探测器距离为0.5，旋转半径为0.5。程序中的N是图像大小，可以根据需要修改。程序中的theta是角度范围，可以根据需要修改。