别怕数学！用Python的Scipy.fft给你的传感器数据做个‘降噪SPA’

当你的温湿度传感器传回的数据像被静电干扰的收音机信号，或是振动监测设备记录到不明来源的机械杂波时，频域处理就像给数据做了一次深度清洁护理。不同于时域中难以分离的信号与噪声，频域视角能让我们像调音师一样精确切除干扰频率——而这一切只需要几行Python代码。

1. 传感器数据降噪的本质：频域美容术

工业现场常见的温度传感器噪声往往集中在50Hz工频干扰，而加速度计采集的振动信号可能混入电机基频谐波。传统移动平均滤波会模糊真实波动，但傅立叶变换的独特价值在于：

• 选择性过滤 ：像用镊子拔除杂毛般精准清除特定频段
• 相位保留 ：逆变换后的信号时间戳与原始数据严格对齐
• 非破坏性处理 ：保留有用频段振幅不受影响

# 典型传感器噪声特征频率
noise_profiles = {
    '温度传感器': [50, 100],    # 工频及其二次谐波
    '振动传感器': [25, 75, 150], # 电机基频及轴承缺陷频率
    '麦克风': [60, 12000]       # 交流哼声与高频啸叫
}

实际案例：某光伏电站的逆变器温度采样数据显示，每15分钟出现一次规律性毛刺，经FFT分析发现是巡检无人机通信频段干扰

2. 实战四步曲：从原始信号到纯净波形

2.1 数据准备与可视化诊断

导入一段实际Arduino采集的温湿度数据：

import pandas as pd
import matplotlib.pyplot as plt

sensor_data = pd.read_csv('environment.csv', 
                         parse_dates=['timestamp'])
temp = sensor_data['temperature'].values
time = sensor_data['timestamp'].astype('int64')//1e9  # 转为Unix时间戳

plt.figure(figsize=(12,6))
plt.plot(time, temp, 'b-', alpha=0.5, label='Raw')
plt.title('Temperature Sensor with 50Hz Noise')
plt.xlabel('Timestamp (s)'); plt.ylabel('°C')

2.2 快速傅立叶变换实施

使用 scipy.fft 的rfft实现高效计算：

from scipy.fft import rfft, rfftfreq

n = len(temp)
sampling_rate = 1/(time[1]-time[0])  # 计算实际采样率
yf = rfft(temp - temp.mean())         # 去直流分量
xf = rfftfreq(n, 1/sampling_rate)

plt.plot(xf, np.abs(yf), 'r-')
plt.xlim(0, sampling_rate/2)          # 显示有效频段
plt.axvline(50, color='k', linestyle='--')

关键参数对照表：

参数	说明	典型取值
n	采样点数	建议1024的整数倍
sampling_rate	采样频率	根据传感器规格
xf	频率坐标	自动计算

2.3 动态阈值降噪算法

不同于固定阈值，动态算法根据信号特征自动调整：

def adaptive_threshold(spectrum):
    median = np.median(np.abs(spectrum))
    mad = 1.4826 * np.median(np.abs(spectrum - median))
    return median + 3*mad  # 基于统计学离群值检测

noise_threshold = adaptive_threshold(yf)
clean_mask = np.abs(yf) > noise_threshold
yf_clean = clean_mask * yf

2.4 逆变换与效果验证

from scipy.fft import irfft

temp_clean = irfft(yf_clean) + temp.mean()
plt.figure(figsize=(12,6))
plt.plot(time, temp, 'b-', alpha=0.3, label='Raw')
plt.plot(time, temp_clean, 'r-', lw=2, label='Denoised')
plt.legend()

3. 高级调参：避免过度滤波的陷阱

3.1 采样率与频率分辨率

两者关系满足： Δf = 采样率/采样点数 。举例说明：

采样率(Hz)	采样点数	分辨率(Hz)	适用场景
1000	1024	0.98	振动分析
10	512	0.02	温度监测

3.2 窗函数选择指南

常见窗函数对频谱泄露的影响：

窗类型	主瓣宽度	旁瓣衰减(dB)	适用场景
矩形窗	0.89	-13	瞬态信号
汉宁窗	1.44	-31	常规使用
平顶窗	3.77	-44	幅值精度要求高

from scipy.signal import windows

hanning = windows.hann(len(temp))
yf_windowed = rfft(temp * hanning)

4. 特殊场景处理方案

4.1 非平稳信号处理

对于随时间变化的噪声频率，建议采用短时傅立叶变换(STFT)：

from scipy.signal import spectrogram

f, t, Sxx = spectrogram(temp, fs=sampling_rate, 
                       nperseg=256, noverlap=128)
plt.pcolormesh(t, f, 10*np.log10(Sxx))
plt.ylabel('Frequency [Hz]')
plt.xlabel('Time [sec]')

4.2 实时处理架构

嵌入式设备上的实现方案：

# 伪代码示例
class RealTimeFFTFilter:
    def __init__(self, window_size=512):
        self.buffer = np.zeros(window_size)
        
    def update(self, new_sample):
        self.buffer = np.roll(self.buffer, -1)
        self.buffer[-1] = new_sample
        
        if time_to_process():
            spectrum = rfft(self.buffer)
            # ...处理逻辑...
            return irfft(clean_spectrum)

在树莓派上的实测性能：处理512点FFT仅需1.2ms，满足多数实时需求