DFT变换 频率分辨率问题 子载波间隔

关于DFT后的频谱图横纵坐标单位，或DFT的频率分辨能力。
正变换：

逆变换：

MATLAB仿真程序如下：

clear all;close all;clc
Fs = 1000; % 采样频率
Ts = 1/Fs; % 采样周期
N = 1000; % 信号长度 由此知，频率分辨率为Fs/N = 1hz
t = (0:N-1)*Ts; % 时间向量
S = 0.7*exp(1i*2*pi*50*t) + exp(1i*2*pi*150*t); % 原始函数 按照 Fs = 1000进行采样。
figure;
subplot(2,1,1);plot(t,real(S));
subplot(2,1,2);plot(t,imag(S));

Y = fft(S);
P = abs(Y/N); % 每个频域分量的幅度值
f = (Fs/N).*(0:N-1); %频谱图横轴，单位Hz

figure;
stem(f,P)
title('Amplitude Spectrum of S(t)')
xlabel('f(Hz)')
ylabel('|P(f)|')

结合网上资料，个人理解，欢迎交流指正！