1. FFT频谱分析的核心原理第一次接触FFT时我被这个能把时域信号变成频域信号的魔法工具惊艳到了。简单来说FFT快速傅里叶变换就像是一个信号翻译官它能把杂乱无章的波形信号翻译成我们看得懂的频率成分。想象你在听交响乐FFT就是那个能告诉你现在是小提琴在演奏还是大提琴在发声的神奇工具。采样定理是FFT分析的地基。我刚开始做实验时就踩过坑用100Hz采样率去分析一个80Hz的信号结果频谱完全失真。后来才明白采样频率必须大于信号最高频率的2倍奈奎斯特频率这个铁律。比如要分析1kHz的信号采样率至少得2kHz实际应用中我通常会选2.56倍以上。MATLAB中的fft函数用起来很简单Fs 1000; % 采样频率1kHz t 0:1/Fs:1-1/Fs; % 1秒时间向量 x cos(2*pi*100*t); % 100Hz余弦信号 X fft(x); % 执行FFT但这里有个关键细节FFT结果的前半部分才是有效的频率分量。因为FFT结果是对称的后半部分其实是前半部分的镜像。所以我在分析时通常只取前N/2个点N是采样点数。2. MATLAB实验的关键参数设置2.1 采样频率的选择策略记得有次做电机振动分析设置的采样率不够高结果高频噪声把真实信号完全淹没了。后来我总结出一个实用公式采样频率 目标最高频率 × 2.56 × 安全系数(1.2~1.5)比如要分析500Hz以内的振动我会设置Fs 500 * 2.56 * 1.2; % 约1.5kHz采样率2.2 变换区间N的黄金法则N的选择直接影响频率分辨率。我常用这个对比表来说明N值频率分辨率计算量适用场景256较高小实时监测1024中等中一般分析4096很高大精密测量在MATLAB中可以通过补零来提高频谱显示效果N 1024; % 实际采样点数 N_fft 4096; % 使用4096点FFT X fft(x, N_fft); % 补零后的FFT2.3 窗函数的选择技巧不加窗就做FFT就像不用三脚架拍长曝光照片。我最常用的是汉宁窗(Hanning)它对抑制频谱泄露效果很好win hanning(N); % 生成汉宁窗 x_windowed x .* win; % 加窗 X fft(x_windowed); % 再做FFT不同窗函数的对比矩形窗频率分辨率最高但泄露严重汉宁窗平衡分辨率和泄露通用首选平顶窗幅值精度最高但分辨率较低3. 频谱分析中的误差控制实战3.1 频谱泄露的识别与解决有一次分析50Hz工频信号时发现频谱上出现拖尾现象这就是典型的频谱泄露。解决方法有三步确保采样时长是信号周期的整数倍选择合适的窗函数增加采样点数MATLAB示例% 错误示范非整数倍周期 t 0:1/Fs:0.019; % 19ms采样不是20ms(50Hz周期) x sin(2*pi*50*t); X fft(x); % 正确做法 t 0:1/Fs:0.02-1/Fs; % 精确20ms x sin(2*pi*50*t);3.2 栅栏效应的应对方案栅栏效应就像透过栅栏看风景会错过某些细节。要减少这种效应增加FFT点数补零使用更高采样率采用Zoom FFT技术这里有个实用技巧在MATLAB中观察频谱时用interp1函数进行插值f (0:N/2-1)*Fs/N; % 原始频率轴 f_interp linspace(0,Fs/2,10*N); % 10倍插值 X_interp interp1(f, abs(X(1:N/2)), f_interp, spline);3.3 幅值校正的注意事项FFT得到的幅值需要校正才能反映真实值。对于不同窗函数校正系数也不同% 汉宁窗幅值校正 X_corrected 2 * abs(X(1:N/2)) / sum(hanning(N));实测数据对比正弦信号理论幅值1.0不加窗测量值0.99汉宁窗未校正0.5校正后1.0024. 进阶实战多信号分析与案例解析4.1 含噪声信号的频谱分析模拟一个含噪声的信号x_clean 0.5*sin(2*pi*50*t) sin(2*pi*120*t); noise 0.5*randn(size(t)); x_noisy x_clean noise;用平均功率谱密度(PSD)提高信噪比[Pxx,f] pwelch(x_noisy, hanning(256), 128, 1024, Fs); plot(f, 10*log10(Pxx))4.2 实际工程案例轴承故障诊断通过特征频率识别故障% 轴承外圈故障特征频率计算 BPFO (n/2) * (1 - d/D * cos(φ)) * rpm/60; % 在频谱中查找BPFO附近的峰值4.3 自动峰值检测技巧用findpeaks函数自动识别频谱峰值[peaks, locs] findpeaks(abs(X(1:N/2)), MinPeakHeight, 0.1); f_peaks (locs-1) * Fs/N; % 转换为实际频率最后分享一个调试心得频谱分析时一定要先目测时域信号是否正常。有次我花了半天调频谱后来发现是采集卡接触不良导致信号断续。记住好的频谱一定来自好的时域信号。