基于matlab的fbp算法程序

function [imgfilter]=fbp2（rthetaxpimgrowsimgcolsfilter_type）;
[nk]=size（r）;
xp_offset=abs（min（xp））+1;
width=2^nextpow2（n）;
tmpimg=zeros（imgrowsimgcols）;
ramlak=2*[0:（width/2-1）width/2:-1:1]‘/width;
if（filter_type==0）
    filter=ones（width1）;
elseif（filter_type==1）
    filter=ramlak;
elseif（filter_type==2）
    Sinc = abs（sinc（2*（0:width-1）/2:-1:1））;
    temp = [ sinc（1:（width/2）） sinc（width/2:-1:1）];
    filter=ramlak.*temp‘;
elseif（filter_type==3）
    cosine=abs（cos（2*pi*（0:width-1）/（2*width）））;
    filter=ramlak.*cosine;
elseif（filter_type==4）
    hamming=0.54-0.46*cos（2*pi*（width-1）/width）;
    temp=[hamming（width/2:width） hamming（1:widyh/2-1）];
    filter=ramlak.*temp‘;
elseif（filter_type=