2、傅里叶分析入门:周期信号的傅里叶级数、非周期信号的傅里叶变换、频谱的概念及其在导弹制导中的应用

2.1 为什么导弹工程师要懂傅里叶?

说实话,我刚开始接触导弹制导时,也觉得傅里叶变换是数学家的玩具。

直到有一次,我在外场测试中遇到了一个棘手的抖动问题——导弹在末端制导阶段,导引头输出的角度信号总是不稳定。时域波形看起来乱七八糟,根本找不到规律。后来我用频谱仪一看,才发现是50Hz的电源工频干扰混进来了。

你看,这就是傅里叶分析的价值。它能把时间上乱成一团的信号,拆解成不同频率的成分。说白了,就是给你一副「频率眼镜」,让你看到时域里看不到的东西。

核心思想:任何信号都可以看作不同频率正弦波的叠加。傅里叶分析就是找出这些正弦波的频率、幅度和相位。

2.2 周期信号的傅里叶级数

先聊聊周期信号。什么是周期信号?就是每隔一段时间T,波形重复一次的信号。比如导弹发动机的振动信号、舵机的周期性偏转指令,都属于这一类。

傅里叶级数告诉我们:一个周期为T的信号f(t),可以写成这样:

f(t) = a₀/2 + Σ[aₙ·cos(nω₀t) + bₙ·sin(nω₀t)]   (n=1,2,3,...)

其中ω₀ = 2π/T,是基波角频率。a₀是直流分量,aₙ和bₙ是各次谐波的系数。

我个人习惯把傅里叶级数理解成「拆积木」:

  • 直流分量a₀:信号的平均值,相当于积木的底座
  • 基波(n=1):频率等于信号重复频率的成分,是积木的主体
  • 谐波(n≥2):频率为基波整数倍的成分,是积木的细节装饰

我在项目中遇到过一件事:某型导弹的舵机指令信号,时域看起来是标准的方波。但用傅里叶级数一展开,发现高次谐波的能量占比很大。这些高频成分会激发舵机的结构共振,导致控制精度下降。后来我们在指令前端加了个低通滤波器,问题就解决了。

实战技巧:工程中很少需要计算无穷多项。一般取到5~7次谐波就能还原90%以上的信号能量。取太多反而会引入噪声。

2.3 非周期信号的傅里叶变换

周期信号好处理,但导弹飞行过程中,更多是非周期信号——比如目标机动产生的加速度突变、导弹发射瞬间的冲击响应。这些信号没有固定的周期,傅里叶级数就派不上用场了。

这时候就要请出傅里叶变换:

F(ω) = ∫f(t)·e^(-jωt) dt   (积分从-∞到+∞)

嗯,这个公式看起来有点吓人。但你想想看,它其实是在做一件事:把信号f(t)和不同频率的复指数信号做「内积」,看看每个频率成分的贡献有多大。

我建议你这样理解:

  • 时域f(t):信号随时间的变化过程
  • 频域F(ω):信号在不同频率上的能量分布
  • 傅里叶变换:就是时域到频域的「翻译器」

举个例子。导弹发射时的冲击信号,时域上是一个短暂的尖峰。它的傅里叶变换是什么?是一个很宽的频谱——从低频到高频都有能量。这就是为什么冲击信号容易干扰电子设备,因为它「覆盖」了所有频率。

注意:傅里叶变换要求信号绝对可积,即∫|f(t)|dt < ∞。实际工程中,我们处理的信号都满足这个条件。但如果信号包含直流分量或趋势项,变换结果会在ω=0处出现冲击函数,需要先做预处理。

2.4 频谱的概念

频谱,就是信号在频率域上的「画像」。它包含两个部分:

  • 幅度谱:|F(ω)|,表示各频率成分的强度
  • 相位谱:∠F(ω),表示各频率成分的时间偏移

我刚开始做信号处理时,总盯着幅度谱看,觉得相位谱不重要。直到有一次,我在做导引头信号重构时,发现幅度谱完全正确,但重构出来的波形和原始信号对不上。查了半天,才发现是相位信息丢了。

你想想看,两个幅度相同、相位不同的正弦波叠加,结果可能天差地别。所以,相位谱和幅度谱同等重要,千万别忽略。

下面这张图展示了频谱分析的核心逻辑:

频谱分析核心逻辑 时域信号 f(t) 时间轴上的波形 (如导引头角度、舵机指令) 傅里叶变换 频域信号 F(ω) 频率轴上的分布 (幅度谱 + 相位谱) 逆变换 重构 f'(t) 导弹制导中的典型应用 🔹 导引头信号滤波:识别并滤除50Hz工频干扰 🔹 舵机频率响应分析:避免结构共振 🔹 目标机动特征提取:从回波中分离目标速度信息 🔹 制导指令平滑:去除高频抖动,提高跟踪稳定性 图:从时域到频域再到应用的完整链路

2.5 在导弹制导中的具体应用

讲完理论,咱们来点实际的。频谱分析在导弹制导中到底怎么用?我挑三个最常见的场景说说。

2.5.1 导引头信号去噪

导引头输出的角度信号,经常混有各种噪声。比如发动机振动、电源纹波、电磁干扰等。这些噪声在时域里和有用信号混在一起,很难分离。

但在频域里,情况就清晰多了:

  • 有用信号:频率较低,一般在几Hz到几十Hz(取决于目标运动速度)
  • 振动噪声:集中在发动机的基频及其谐波上(比如100Hz、200Hz)
  • 电源干扰:50Hz及其倍频

我曾经处理过一个案例:某型导弹在飞行试验中,导引头输出信号总是有周期性波动。时域波形看起来像正弦波,但频率不稳定。用FFT一分析,发现是舵机反馈信号通过电源线耦合进了导引头。我们在电源入口加了个EMI滤波器,问题就解决了。

避坑指南:我曾经在滤波器设计上吃过亏——为了滤除噪声,把截止频率设得太低,结果把有用信号的高频成分也滤掉了,导致导引头响应变慢。后来我养成了一个习惯:先做频谱分析,搞清楚信号和噪声的频率范围,再设计滤波器参数。

2.5.2 舵机频率响应分析

舵机是导弹控制系统的执行机构。它的频率响应特性直接影响控制精度。怎么测?给舵机输入不同频率的正弦信号,记录输出幅度和相位变化,然后画出伯德图。

这就是频谱分析在系统辨识中的应用。我记得有一次,某型舵机在80Hz附近出现了谐振峰,导致控制系统在该频率附近失稳。我们用频谱分析定位了问题,然后通过调整舵机结构阻尼,把谐振峰压了下去。

2.5.3 目标特征提取

在雷达制导中,目标回波信号包含多普勒频移信息。这个频移量正比于目标相对于导弹的径向速度。通过傅里叶变换提取回波信号的频谱,就能计算出目标速度。

说白了,就是利用频谱分析来「听」出目标在动。我参与过的一个项目中,我们通过分析回波信号的频谱宽度,还能判断目标是否在机动——机动目标会产生频谱展宽效应。

2.6 工程实现中的几个要点

理论讲完了,最后分享几个我在工程实践中总结的要点:

要点 说明 我的建议
采样率选择 根据奈奎斯特定理,采样率至少为信号最高频率的2倍 实际工程中取3~5倍,留出裕量
频谱泄漏 非整周期截断会导致频谱能量扩散到相邻频率 加窗函数(汉宁窗、海明窗)可以缓解
频率分辨率 Δf = fs / N,N为采样点数 想提高分辨率,要么降低采样率,要么增加采样点数
实时性要求 导弹制导系统对延迟敏感,FFT计算时间要可控 优先使用基2-FFT算法,点数选2的幂次

嗯,傅里叶分析的内容就讲到这里。它不是什么高深莫测的数学工具,而是我们导弹工程师手里的一把「手术刀」——用来切开信号,看清里面的结构。下次你在调试中遇到奇怪的信号抖动,不妨先做个频谱分析,说不定问题就一目了然了。

一句话总结:时域让你看到「发生了什么」,频域让你看到「为什么发生」。两者结合,才是信号处理的完整视角。

专注资料整理