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，直接輸出傅里葉變換後的結果，進行比較