如何在MATLAB中计算高斯曲线的调制传递函数？

Question

我正在尝试使用 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轴感到困惑。 X 轴的lines per mm范围是多少？ 谁能帮我提供一个计算 MTF 图的 X 轴的想法？

提前致谢！ 摩奴

Answer 1

matlab手册对计算非常具体。 例如，看看这个

dx = 0.1
x = -15:dx:15;      
Fs = 1/dx;          % Sampling frequency
L = length(x);      % Signal length
sigma = 5;
amp   = 437500;
y     = amp/sqrt(2*pi*sigma^2) * exp(-x^2/(2*sigma^2));

n = 2^nextpow2(L);
Y = fft(y,n);
f = Fs*(0:(n/2))/n;
P = abs(Y/n);
plot(f,P(1:n/2+1)) 
title('Gaussian Pulse in Frequency Domain')
xlabel('Frequency f [in units of 1/dx]')
ylabel('|P(f)|')

Answer 2

你为真实空间高斯使用了多少点。 如果这个数字是 N，要在 (1/mm) 中转换 fft x 轴，你应该除以 N。