fft频谱分析需四步校正：①建立物理频率轴f=(0:n-1)*fs/n；②幅值归一化p2=abs(y)/l并单边补偿×2；③加窗（如hanning）抑制泄漏，归一化用sum(win)；④零填充+插值提升频率估计精度。
一、理解FFT输出的物理频率映射
如果您在MATLAB中执行fft函数后得到一串复数值，却无法判断每个索引对应的实际频率，那么频谱图将失去物理意义。FFT结果的下标与真实赫兹频率之间存在严格换算关系，该关系由采样率Fs和FFT点数N共同决定。不建立正确的频率轴，所有幅值和相位分析均不可靠。
1、确定信号采样率Fs（单位：Hz），该值必须大于信号最高频率的两倍以满足奈奎斯特采样定理。
2、获取或设定FFT点数N，推荐使用2的整数次幂（如1024、2048）以启用快速算法。
3、计算物理频率向量：f = (0:N-1)*Fs/N，此向量长度为N，覆盖0至接近Fs的数字频率范围。
4、对实信号，取单边谱时仅保留前N/2+1个点，并对除直流和奈奎斯特点外的其余幅值乘以2，以补偿被舍弃的负频率能量：P1(2:end-1) = 2*P1(2:end-1)。
二、修正FFT幅值以匹配真实物理幅值
原始fft输出的复数模值未经归一化，直接绘图会导致幅值严重失真。MATLAB中fft(x)返回的是DFT定义下的未缩放结果，必须依据信号长度进行幅度校正，否则50Hz正弦波在频谱上可能显示为远高于其理论幅值2的结果。
1、计算双边幅值谱：P2 = abs(Y)/L，其中Y = fft(x)，L为原始信号长度。
2、提取单边谱：P1 = P2(1:L/2+1)，适用于实信号且仅关注0至Fs/2区间。
3、执行能量守恒修正：对P1中第2个至倒数第二个元素执行乘2操作，确保总能量与原始信号能量一致。
4、若使用补零后的N点FFT（N > L），应以N而非L进行归一化：P2 = abs(Y)/N，避免因补零引入虚假幅值衰减。
三、消除频谱泄漏：窗函数应用法
矩形截断（即默认无窗处理）会强制信号在边界处突变，导致能量扩散至邻近频率，形成虚假“毛刺”和主瓣展宽。这是新手频谱图出现多峰、分辨率下降的最常见原因。窗函数通过平滑截断边界来抑制泄漏，但会以牺牲频率分辨率为代价换取动态范围提升。
1、选择合适窗函数：对单频检测优先用矩形窗；对多频分离或弱信号检测推荐汉宁窗（hanning）或汉明窗（hamming）。
2、生成窗向量：win = hanning(L);，确保长度L与信号x相同。
3、加窗处理：x_win = x .* win;
4、执行FFT并归一化：Y_win = fft(x_win); P2_win = abs(Y_win)/sum(win); 使用窗系数和而非信号长度进行归一化，保证幅值准确性。
四、解决频率定位偏差：零填充与插值法
当信号实际频率介于两个相邻FFT谱线之间时，能量将分散到两侧谱线，造成峰值频率偏移和幅值低估。这种偏差并非误差，而是DFT固有栅栏效应所致。零填充可在不增加真实信息的前提下提高谱线密度，辅助更精确地估计峰值位置。
1、对原始信号x补零至更高点数Nz（如Nz = 4096），保持采样率Fs不变：x_zp = [x, zeros(1, Nz-length(x))];
2、执行Nz点FFT：Y_zp = fft(x_zp);
3、构建高密度频率轴：f_zp = (0:Nz-1)*Fs/Nz;
4、使用findpeaks函数定位主峰：[pks, locs] = findpeaks(abs(Y_zp(1:Nz/2+1))/Nz, 'MinPeakHeight', 0.1);
5、通过二次插值估算真实频率：f_est = f_zp(locs) + (abs(Y_zp(locs+1)) - abs(Y_zp(locs-1))) / (2*(2*abs(Y_zp(locs)) - abs(Y_zp(locs+1)) - abs(Y_zp(locs-1)))) * (Fs/Nz); 该公式可将频率估计精度提升至Fs/Nz量级以下。
五、应对相位失真：解缠绕与参考点对齐
FFT输出的angle()函数返回[-π, π]范围内的主值相位，当真实相位跨越±π边界时会出现跳变，导致相位谱呈现不连续锯齿状。此外，若信号起始时间未对齐周期整数倍，还会引入与时间偏移相关的线性相位斜率，掩盖真实系统相位响应。
1、对原始相位向量phase_raw = angle(Y(1:L/2+1))执行相位解缠绕：phase_unwrap = unwrap(phase_raw);
2、去除线性趋势以提取系统相位：phase_detrend = detrend(phase_unwrap, 'linear');
3、若需校正时间偏移影响，可先对信号做循环移位使其峰值对齐t=0：[~, idx_max] = max(abs(x)); x_shift = circshift(x, -idx_max+1);
4、对移位后信号重做FFT，并提取相位：此时相位谱更准确反映信号内在频率成分的相对延迟特性。
