中心切片定理的matlab实现

clear
sl=phantom;
the = 0:179;
[Rxp] = radon（slthe）;
[widthlength]=size（R）;      
width1=width;
length1=width;
rec=zeros（width1length1）;  
flou=zeros（widthlength）;
for i=1:length
flou（:i） =fftshift（R（:i））; 
flou（:i）=fft（flou（:i））;  
end
o=round（width1/2）;
 for i=1:width1
    for j=1:length1
        x=j-o;
        y=o-i;        
            if（x==0）       
                  if（y>=0）
                     theta=90;
                 else
                     theta=-90;
                  end
            else
                theta=（atan（y/x））*180/pi;
            end            
            if（x<=0）
                theta=180+theta;
            elseif（y<0）
                theta=2*180+theta;
            end            
           r=sqrt（x^2+y^2）;
            if（theta<180）