R?=?xlsread（‘A題附件.xls‘2）;?
W?=?512;???
theta?=?1:180;?
fourier?=?fft（R?W）;?
filter?=?2*[0:（W/2-1）?W/2:-1:1]‘/W;???
lvbo?=?zeros（W180）;??
for?i?=?1:180???
?????lvbo（:i）?=?fourier（:i）.*filter;???
end

inversse?=?real（ifft（lvbo））;???
N?=?360;
final?=?zeros（N）;?
for?i?=?1:180???
rad?=?theta（i）*pi/180;?
?for?x?=?（-N/2+1）:N/2???
?for?y?=?（-N/2+1）:N/2?
t?=?round（x*cos（rad+pi/2）+y*sin（rad+pi/2））;???
final（x+N/2y+N/2）=final（x+N/2y+N/2）+inversse（t+round（size（R1）/?2）i）;
?end
?end
end
final=final/180;
save?date?final;
imshow（final）title（‘final‘）
load?date?final

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件?????4231680??2017-08-09?10:09??A題附件.xls
?????文件?????????594??2018-08-28?21:14??wen1.m
?????文件????????1197??2018-08-28?21:15??wen2.m
?????文件????????1225??2018-08-28?21:16??wen3.m