光束偏转技术的FBP-matlab滤波器设计

光束偏转技术的FBP-matlab滤波器设计,matlab,filter,tomography-reconstruction,Matlab,Filter,Tomography Reconstruction,我需要你的帮助来回答一个关于iradon过滤器设计的问题。实际上,在视线投影(光束偏转)中,与吸收等其他技术相反,ram lak滤波器被另一个离散响应为: 在频域中,其表示为: 受iradon的启发,我实现了以下过滤器: n = 0:(order/2); filtImpResp = zeros(1,(order/2)+1); filtImpResp(1) = 0; % Set the DC term filtImpResp(2:2:end) = 1./(pi^2.*n(2:2:en

我需要你的帮助来回答一个关于iradon过滤器设计的问题。实际上,在视线投影(光束偏转)中,与吸收等其他技术相反,ram lak滤波器被另一个离散响应为:

在频域中,其表示为:

受iradon的启发,我实现了以下过滤器:

n = 0:(order/2);

filtImpResp = zeros(1,(order/2)+1); 

filtImpResp(1) = 0; % Set the DC term 

filtImpResp(2:2:end) = 1./(pi^2.*n(2:2:end));

filtImpResp = [filtImpResp -filtImpResp(end-1:-1:2)]; 

filt =imag(fft(filtImpResp))

% but it doesn't give a good results.
你能看一下吗?我的实现准确吗?特别是我对原始iradon函数中的这一行有疑问:

filt = 2*real(fft(filtImpResp)); 
为什么它要乘以2?我该怎么办

非常感谢


致以最良好的祝愿首先,我设置
顺序=100

你提到的那个,我猜是过滤器的问题。您可以通过以下方式在MATLAB上检索它:

绘图(角度(fft移位(fft(filtempresp)))

以相同的方式,幅值响应:

绘图(abs(fftshift(fft(filtempresp)))

标量2可能用于规范化滤波器的大小,例如从0到1,即使它在代码上写为实数