矩形波如何傅立叶展开_步进频连续波生命体回波处理——以MUSIC算法为例

此前，在MUSIC算法在微弱生命体回波处理中的应用？一文中，我们已经详细讲述了MUSIC算法的基本原理，并通过一个单频雷达波回波仿真实例证明了MUSIC算法的有效性，相比于传统的基于傅立叶变换的频谱分析方法，MUSIC算法可以克服谐波对信号的谱图的不良影响（谱图中不会出现谐波的尖峰）。MUSIC算法本质上是利用噪声子空间与信号有效成分子空间的正交性的实现对生命体探测的一种谱分析方法。然而，受限于单频雷达的工作模式，我们只能从回波信号中分析出信号中是否存在与生命体有关的特征频率成分，而无法实现生命体与雷达之间的测距。而步进频连续波雷达的出现，则为测距提供了可能。下面，我们将主要介绍步进频连续波雷达的工作模式，并给出仿真实例。

步进频连续波雷达

步进频，顾名思义，就是雷达的工作频率会随着时间而发生周期性的变化，对于有

个工作频率的步进频连续波雷达，其信号频率随时间发生的周期性变化的示意图如图1所示，由图可知，雷达的扫频周期为

。

图 1 步进频连续波雷达信号频率

假设该步进频雷达的起始频率为

，频率间隔为

，则每个频率的发射信号的时域表达式如式（1）所示：

其中，

是虚数单位，“rect”是矩形波函数。在上式中加入微弱生命体信号的时延

，且反射系数为

（R为反射系数的幅值），则步进频连续波的微弱生命体回波信号的表达式为：

由于

远大于

，且雷达接收机会对回波信号作

I/Q解调，上式变为：

至此，我们得到了步进频连续波雷达的回波信号的表达式，就可以实现步进频连续波雷达的生命体探测模拟了。

仿真实例

雷达工作的起始频率为1GHz，终止频率为3GHz，频率间隔为10 MHz，步进时间为100

，步进频率数为200。在距离雷达1 m处存在一个障碍物，5 m处存在一微弱生命体，微弱生命体的呼吸频率为0.3 Hz，振幅为4 mm，心跳频率为1.2 Hz，振幅为2 mm。雷达接收信号的采样点数为1000，采样率为50 Hz（

步进周期为0.02 s）。采取MUSIC算法分析回波信号（ 对于步进频连续波雷达，回波信号是一矩阵），最后得到谱分析结果，如图2所示。

图 2 步进频连续波雷达回波信号的谱分析结果

从图2可以看出，在5 m处存在两个尖峰，分别对应呼吸频率和心跳频率。需要注意的是回波矩阵频率轴的的步进频率索引

与距离

存在以下关系式：

其中，

是真空中电磁波的速度。

代码

SteppedFrequencyMusic.m

%% Simulation of weak life detection based on stepped frequency continuous wave radar

---

PS：完整的代码已经整理到本人的百度企业云盘中，感兴趣的同学可以从下方链接（有效期7天）下载哦~~

SFR_MUSIC​eyun.baidu.com

参考资料

[1]刘国菊. 微弱生命体回波处理算法研究[D].电子科技大学,2013.