我正在尝试通过MATLAB查找高斯拟合曲线的调制传递函数。高斯曲线如下:
x轴是距离(形式为-15mm to 15mm
),y轴是计数(幅度)。我使用以下代码从高斯曲线中找到傅立叶变换
FFT_y = fft(y); %take fourier transform
FF_mag = abs(FFT_y )/(length(FFT_y )); %find magnitude
FF_mag = (FF_mag-min(FF_mag))./(max(FF_mag)-min(FF_mag)); %normalize magnitude
我使用以下代码裁剪了FF_mag
FF_mag_nw = FF_mag(1:(length(FF_y)/32));
plot(FF_mag_nw);
我在上面的代码中使用32来获取图形的主要部分,并且得到了如下的MTF图:
我对X轴感到困惑。 lines per mm
中X轴的范围是多少?谁能帮助我给出一个计算MTF图X轴的想法?
谢谢!Manu
您对实际空间高斯使用了多少点。如果此数字为N,则要将fft x轴转换为(1 / mm),应除以N。