matlab验证中心切片定理

刚看到别人用matlab的程序验证中心切片定理，但是有两处不明，然后自己做了修改

clear all 
clc
I=phantom(256);%新建一个shapp_logan模型I
figure,imshow(I);%见于图1
H=imrotate(I,270);%将图像旋转270°，相当于θ=90°，
                  %使得X线以与X轴平行的方向进行投影
figure,imshow(H);%见于图2
p=sum(H,1);%进行投影
figure,plot(1:256,p);%见图1
f=fft(p);
fab=f.*conj(f);
%f3=fftshift(fab);%**修改处**
figure,plot(1:256,fab);%见图2
f2d=fft2(I);%此处是直接对原始图像进行二维傅立叶变换并进行验证
fab2d=f2d.*conj(f2d);
%f4=fftshift(fab2d);%**修改处**
figure,plot(fab2d(:,1));%见图3

傅里叶变换后不使用fftshift，直接输出傅里叶变换后的结果，进行比较