使用MATLAB进行噪声信号的频谱分析源码步骤

你想了解的是在MATLAB中针对含噪声信号进行频谱分析的具体方法，核心需求是如何有效抑制噪声干扰、准确提取信号的真实频率特征，这也是工程实践中最常见的频谱分析场景。

一、噪声信号频谱分析的核心难点与解决思路

含噪声信号的频谱分析最大问题是：噪声会掩盖真实信号的频率峰值，且易出现频谱泄漏。MATLAB中解决这个问题的核心思路是：

1. 加窗处理：抑制频谱泄漏，让频率峰值更集中；
2. 功率谱分析（Welch方法）：通过分段、重叠、平均降低噪声方差，提升频谱信噪比；
3. 平滑处理：进一步滤除高频噪声毛刺。

二、完整实操：MATLAB噪声信号频谱分析（分步讲解）

场景说明

以“50Hz基波+120Hz谐波+高斯白噪声”的混合信号为例，对比“直接FFT”和“Welch功率谱”两种方法的分析效果，让你清晰看到噪声抑制的差异。

完整可运行代码

% ==================== 1. 基础参数设置 ====================Fs=500;% 采样频率(Hz)，需满足奈奎斯特采样定理（>2*最高信号频率）T=1/Fs;% 采样周期L_total=2000;% 总采样点数t=(0:L_total-1)*T;% 时间轴% ==================== 2. 生成含噪声信号 ====================% 纯净信号：50Hz(幅值1.5) + 120Hz(幅值0.8)x_clean=1.5*sin(2*pi*50*t)+0.8*sin(2*pi*120*t);% 加入高斯白噪声（噪声幅值0.6，模拟工程中的随机噪声）x_noisy=x_clean+0.6*randn(size(t));% ==================== 3. 方法1：直接FFT（无降噪） ====================Y_direct=fft(x_noisy);% 直接FFT变换L_half=floor(L_total/2);% 取采样点数的一半f=Fs*(0:L_half)/L_total;% 构建频率轴（0~Fs/2）% 计算单边幅值（修正：除总点数，单边频谱需乘以2）P1_direct=abs(Y_direct/L_total);P1_direct=P1_direct(1:L_half+1);P1_direct(2:end-1)=2*P1_direct(2:end-1);% ==================== 4. 方法2：加窗FFT（基础降噪） ====================win=hann(L_total);% 生成汉宁窗（抑制频谱泄漏）x_win=x_noisy.*win;% 信号加窗Y_win=fft(x_win);P1_win=abs(Y_win/L_total);P1_win=P1_win(1:L_half+1);P1_win(2:end-1)=2*P1_win(2:end-1);% ==================== 5. 方法3：Welch功率谱（最优降噪） ====================% 关键参数：分段长度=256，重叠长度=128，使用汉明窗[Pxx,f_welch]=pwelch(x_noisy,256,128,[],Fs);% ==================== 6. 结果可视化对比 ====================figure('Color','w','Position',[100,100,800,600]);% 子图1：含噪声信号时域波形subplot(2,2,1);plot(t,x_noisy);xlabel('时间 (s)');ylabel('幅值');title('含噪声信号时域波形');xlim([00.1]);% 仅显示前0.1秒，便于观察细节grid on;% 子图2：直接FFT频谱（噪声明显）subplot(2,2,2);plot(f,P1_direct);xlabel('频率 (Hz)');ylabel('幅值');title('直接FFT频谱（未降噪）');xlim([0150]);% 聚焦0~150Hz（覆盖信号频率范围）grid on;% 子图3：加窗FFT频谱（降噪效果一般）subplot(2,2,3);plot(f,P1_win);xlabel('频率 (Hz)');ylabel('幅值');title('加汉宁窗FFT频谱（基础降噪）');xlim([0150]);grid on;% 子图4：Welch功率谱（最优降噪）subplot(2,2,4);plot(f_welch,10*log10(Pxx));% 转换为分贝，更易观察xlabel('频率 (Hz)');ylabel('功率谱密度 (dB/Hz)');title('Welch功率谱（最优降噪）');xlim([0150]);grid on;

代码关键部分解释

1. 噪声生成：randn(size(t))生成高斯白噪声，乘以0.6控制噪声强度，贴近工程实际；
2. 加窗处理：hann(L_total)生成汉宁窗，通过.*逐点相乘给信号加窗，核心作用是减少频谱泄漏，让频率峰值更尖锐；
3. Welch功率谱：pwelch是MATLAB专为噪声信号设计的函数，原理是将长信号分段、加窗、做FFT后取平均，能大幅降低噪声方差。参数说明：
   • 第一个参数：待分析的噪声信号；
   • 第二个参数：每段的长度（256）；
   • 第三个参数：段与段之间的重叠长度（128，重叠50%效果最佳）；
   • 第四个参数：FFT点数（空则默认等于分段长度）；
   • 第五个参数：采样频率Fs。
4. 分贝转换：10*log10(Pxx)将功率谱转换为分贝刻度，能更清晰区分信号峰值和噪声基底。

运行结果说明

• 直接FFT频谱：50Hz、120Hz峰值被噪声掩盖，毛刺多，难以准确识别；
• 加窗FFT频谱：噪声略有降低，峰值更清晰，但仍有明显噪声毛刺；
• Welch功率谱：噪声基底大幅下降，50Hz、120Hz峰值突出，是噪声信号分析的最优选择。

三、工程实战技巧

1. 窗函数选择：
   • 汉宁窗/汉明窗：通用型，兼顾频谱分辨率和泄漏抑制，优先选择；
   • 布莱克曼窗：泄漏抑制更强，但分辨率稍低，适合噪声极多的场景；
   • 矩形窗：无加窗，仅适合无噪声、整周期采样的理想信号。
2. Welch参数调整：
   • 分段长度越长，频率分辨率越高，但降噪效果越弱；
   • 重叠率建议设为50%（如分段256，重叠128），平衡分辨率和降噪；
3. 噪声进一步抑制：可先对信号做低通/带通滤波（fir1/butter函数），再做频谱分析，效果更佳。