基于经验模态分解的飞机微多普勒特性分析

0 引 言

飞机上普遍都存在旋转部件，如直升机的主旋翼和尾旋翼、螺旋桨飞机的螺旋桨叶片和喷气式飞机发动机的压缩叶片。它们在飞机飞行的同时自身也存在周期性的旋转运动。美国著名学者VICTOR C CHEN将旋转部件的这种旋转运动称为微运动。在一定的目标姿态角范围内，这种微运动会引起雷达回波中的多普勒调制。他将这种微运动产生的多普勒调制称为微多普勒效应。[1-3]微多普勒效应的提出为低分辨雷达对飞机的分类提供了新的途径。

1998年，Huang等人首次提出一种用来分析非线性非平稳信号的自适应分解算法，称为经验模态分解(简称EMD)方法。[4-5]经验模态分解方法与传统的傅里叶变换和小波变换的不同之处在于该方法不需要预先设定一组基信号，而是根据信号自身的局部特征自适应地将其分解成一系列“基信号”。Flandrin等人[6]指出，EMD可以根据输入信号自适应地提取出不同频率的窄带信号。经验模态分解方法已广泛应用于各个领域，包括雷达领域。文献[7]中利用EMD方法将雷达信号中的车身分量和微多普勒分量进行了分离，较好地提取了雷达回波中的微多普勒信息。文献[8]中利用EMD方法将ISAR图像中的螺旋桨机身分量和螺旋桨产生的微多普勒分量进行了分量，消除了微多普勒分量对螺旋桨飞机ISAR成像的影响。本文研究经验模态分解方法在飞机微多普勒特征提取中的应用。

1 经验模态分解(EMD)方法

经验模态分解(EMD)算法在对信号进行分解时并不预先设定一组基，而是自适应地分解出“基信号”，并把该“基信号”称为本征模态函数(IMF)。一个IMF必须满足以下两个性质：

(1) 在整个时间范围内，函数的局部极值点和过零点的数目要么相等要么最多相差一个；

(2) 在任意时刻点，局部最大值的包络(上包络线)和局部最小值的包络(下包络线) 的平均必须为零。

实测直升机的回波信号如图3(a)所示。 对该直升机的复信号及该复信号经验模态分解后各分量进行FFT处理及平移处理。从图3(a)中可以看出，飞机目标频谱中含有一个700 Hz附近幅度很强的机身频谱、幅度较强的零频率处分量。对比原始信号的频谱，从图3(b)中可以看出，第一本征模态函数分量频谱在[-2 500 Hz,-500 Hz]和[500 Hz,2 500 Hz]的频率部分与原始信号频谱基本相同，零频率处的低频部分与原始信号频谱相比其幅度明显较弱。从图3(c)中可以看出，除第一本征模态函数分量外其余本征模态函数和的频谱中只含有[-1 000 Hz,1 000 Hz]低频分量。从图3(d)可以看出，分解得到的残差频谱中只有零频率处幅度较强的分量。

(1) 令中间变量信号

“云石？挖掘云石，那可是我的使命呢。”丁达插话道，“不过，我好像要先处理眼前这堆麻烦。”他皱了皱眉头，冲壶天晓大喊道，“壶天晓，把我的腿还给我！”

(2) 在信号 中搜索极大值点，对极大值点序列进行插值构造出信号 的上包络

(3) 在信号 中搜索极小值点，对极小值点序列进行插值构造出信号 的下包络

(4) 计算上包络和下包络的均值，得到均值信号：

(3) 点序列插值，得到φk方向的包络γφk；

(6) 令并重复步骤(2)～(5)直到m成为一个IMF。

其中，mi表示第i个IMF，qL为经过L次分解后的余项。

 

当得到一个IMF后，将该IMF从原信号中减去，得到余项q=s-m。将余项q视为新的待分解信号，对其重复上述“筛”过程，便能得到一系列IMF。最终，得到信号s的分解形式为

由经验模态分解(EMD)的分解过程可以看出，分解得到的每个IMF分量都来自于上一次分解后的剩余信号中，所以各个IMF的频率没有重叠，并且呈递减趋势。因此，借助经验模态分解方法可以将信号分解成一系列不同频率的子信号，并且信号的频率在逐渐减小。

Rilling等人将实信号EMD算法的思想扩展到复数域，提出了适用于复信号的EMD算法——CEMD。对于复信号，计算旋转的包络均值均需要搜索M个方向上的极值点。这些方向可以由φk=2πk/M表示，其中0≤k≤M-1。给定一个待分解信号s，复信号单一IMF计算的“筛”过程可以表示如下：

(1) s投影至方向φk上：Pφk=Re(e-jφks)；

(2) 影信号Pφk的极值点位置；

(5) 将均值信号从 中减去，得到细节信号：

(4) 重复步骤(1)～(3)，直到得到K个方向的包络；

(5) 计算包络均值：

在取得现场检测的采集视频资料后，通过管道检测软件对检测视频进行判读并作出量化分析，得出山前水库涵洞存在缺陷统计，见表1。

贝伐单抗多囊脂质体的处方优化及体外释放特性研究…………………………………………………… 王毅云等（7）：922

(6) 从原信号中减去均值得到细节项m=s-smean；

(7) 将细节项m作为新的信号，重复步骤(1)～(6)，直到m成为一个IMF。

3.3.1 社区建设规划 社区建设规划分为两部分：一是房屋建设；二是基础设施建设。房屋应与当地传统民居建筑风格相协调，结合现代美学原理进行改造，符合整体特色，提升居住的舒适度和建筑美感。社区基础设施应配备齐全，保证居民交通方便和生活舒适。

喜欢一个人=感觉他身上有wifi啊；暗恋一个人=不知道wifi密码是啥；失恋=本来能自动连接的wifi突然有一天连不上了；暧昧=wifi信号不稳定，有时候强，有时候弱；热恋=知道密码是啥，信号稳定；异地恋=知道wifi密码是啥，但是距离太远，连不上。

2 详细流程设计

实测喷气式飞机的回波信号如图1(a)所示。对该喷气式飞机的复信号及该复信号经验模态分解后各分量进行FFT处理及平移处理。从图1(b)中可以看出，飞机目标原始频谱中含有一个-2 000 Hz附近幅度很强的机身频谱，-1 500 Hz、300 Hz和800 Hz处上分布着幅度较强的微动调制分量。对比原始信号的频谱，从图1(b)中可以看出，第一本征模态函数分量频谱在[-2 500 Hz,-1 500 Hz]和[1 500 Hz,2 500 Hz]的较高频率部分与原始信号频谱基本相同，(0 Hz,1 000 Hz)的低频部分与原始信号频谱相比其幅度明显较弱，以及分解产生的1 200 Hz附近的较强分量。从图1(c)中可以看出，除第一本征模态函数分量外其余本征模态函数和的频谱中含有300 Hz和800 Hz附近幅度较强的微动调制分量，以及分解产生的1 200 Hz附近的较强分量。从图2(d)可以看出，分解得到的残差频谱中只有零频率处的微弱分量。

(1) 设定N个相位方向，令相位方向k=1；

在同类研究中,优质护理管理模式的护理效果为97.73%,护理满意度为100%,从上述结果中可以看出,本次研究结果与同类研究中的效果基本一致[3]。

(2) 计算复信号s在第k个方向上投影实信号p的极大值点；

分析师建议，种植户、储户要密切关注货源质量，做好储备工作；如果储存较大，则要有序出货。对于华北地区种植户来说，目前多数货源都已经入窖，建议密切关注货源质量，做好库内通风、保温等工作，防止货源质量出现问题。（信息来源：安徽农网）

(6) 计算细节信号m的极值点数与零点数之差、能量比值，并重复步骤(2)～(5)，循环迭代直到满足终止条件， m成为一个IMF分量；

(4) 重复步骤(1)～(3)，计算所有相位方向上包络线的均值；

(5) 从复信号s中减去均值得到细节信号m；

(3) 对极大值点进行插值计算，获取第k个相位方向上的包络，

(7) 从复信号s中减去m得到新的复信号s，重复步骤(2)～(6)，直到分解出所有预期的IMF分量。

3 测试与分析

3.1 喷气式飞机

根据实信号EMD算法和复信号CEMD算法的算法思想，设计算法的详细流程如下所示：

  

图1 喷气式飞机的经验模态分解效果

3.2 螺旋桨飞机

实测螺旋桨飞机的回波信号如图2(a)所示。 对该螺旋桨飞机的复信号及该复信号经验模态分解后各分量进行FFT处理及平移处理。从图2(b)中可以看出，飞机目标原始频谱中含有一个-1 000 Hz附近幅度很强的机身频谱和[-500 Hz,1 000 Hz]上分布的幅度较强的微动调制分量。对比原始信号的频谱，从图2(b)中可以看出，第一本征模态函数分量频谱在[-2 500 Hz,-500 Hz]和[500 Hz,2 500 Hz]的频率部分与原始信号频谱基本相同，(-500 Hz,500 Hz)的低频部分与原始信号频谱相比其幅度明显较弱。从图2(c)中可以看出，除第一本征模态函数分量外其余本征模态函数和的频谱中只含有[-700 Hz,700 Hz]低频分量，包含目标调制谱分量。从图2(d)可以看出，分解得到的残差频谱中只有零频率处的微弱分量。

分析四边自由的钢筋混凝土板：长(X方向)4 m，宽(Y方向)3 m，厚0.4 m，密度为2 450 kg/m3。钢筋与混凝土的弹性模量分别为210 GPa和34.3 GPa；钢筋的泊松比为0.3，混凝土的泊松比为0.167；基础板按Φ12@200配置钢筋，水平向配筋率为0.565%；考虑到箍筋为Φ8钢筋，竖向配筋率为0.126%，由此算出钢筋混凝土板的正交各向异性参数如表1所示。

3.3 直升机

经验模态分解(EMD)算法的分解过程类似于“筛”过程。给定一个实信号s，“筛”过程的步骤如下：

[3]周大风、刘春曙、王耀华、丁献芝：《中国大百科全书·音乐舞蹈卷·畲族音乐》，北京：中国大百科全书出版，1989年，第198页。

  

图2 螺旋桨飞机的经验模态分解效果

  

图3 直升机的经验模态分解效果

4 结束语

根据经验模态分解算法的核心思想和算法步骤，本文设计了详细的复信号的经验模态分解算法流程，并将经验模态分解方法应用于飞机雷达回波复信号的分解。通过对实测飞机数据的分析，可以看出经验模态分解方法得到的第一本征模态函数频谱很好地保留了飞机目标频谱的高频部分。除第一本征模态函数外的其余本征模态函数和的频谱较好地保留了飞机目标频谱的低频部分，分解得到的残差频谱近似于目标频谱的零频率杂波。进一步研究经验模态分解方法在飞机雷达回波复信号的分解可以为设计更加稳定的微多普勒特征提取方法打下坚实的基础。

参考文献:

区块链技术是使用块式和链式的存储结构来认证和保存数据，使用共识算法实现生成新区块，使用非对称加密算法保证数据在信道中的安全传输，使用智能合约来处理数据的新型分布式技术。区块链分为私有链、联盟链和公有链。从本质上讲，区块链就是一个去中心化的分布式数据库，任何用户都可以参与到区块链中。用户周围的路由器设备就是一个节点，每个节点都拥有一整套数据的备份，并且各个节点间使用相同的共识机制，通过竞争计算来生成或更新区块链。基于区块链结果的特点，如果任何一个节点失败，其他节点仍能进行正常的工作，且能分辨出是哪一个节点失败。因此，区块链技术解决了传统平台易受攻击或篡改的缺陷。

[1] Chen V C，Li F, Ho S, et al. Micro-Doppler Effect in Radar Phenomenon, Model and Simulation Study[J]. IEEE Trans. on AES, 2006, 42(1): 2-21.

要推进经济转型、产业升级，关键在于推进技术进步。以组织编写《路线图》的方式来推动技术进步，是中国印刷及设备器材工业协会一次全新的尝试。

[2] 张群, 罗迎, 何劲. 雷达目标微多普勒效应研究概述[J]. 空军工程大学学报(自然科学版)，2011，12(2): 22-26.

[3] 张琳, 盛卫星, 马晓峰. 基于微多普勒效应的雷达目标识别算法[J]. 现代雷达, 2007, 29(12): 35-39.

[4] Huang N E, Shen Z, Long S R, et al. Theempiri-cal mode decomposition and the Hilbert spectrum for non-linear and non-stationary time series analysis[J]. Proc. Royal

[5] Yanbing Li, Lan Du, Hongwei Liu. Hierarchical Cla-ssification of Moving Vehicles Based on Empirical Mode Decomposition of Micro-Doppler Signatures[J]. IEEE Trans. on Geoscience and remote sensing, 2012, 51(5): 3001-3013.

[6] Flandrin P, Rilling G, Goncalves P. Empirical mode decomposition as a filter bank[J]. IEEE Signal Processing Letters, 2014, 11(2):112-114.

[7] Cai C J, Liu W X, Fu J S, et al. Empirical mode decomposition of micro-Doppler signature[J]. In Proc. of IEEE Int. Radar Conf, 2005: 895-899.

[8] Bai X, Xing M, Zhou F, et al. Imaging of micro-motion targets with rotating parts based on empirical mode decomposition[J]. IEEE Trans. on Geosci. Remote Sens. 2008, 46(11): 3514-3523.

谈及恩师王选在心中的形象，肖建国做了如下的总结：第一，追求创新。在提出研究课题时，恩师经常会尝试许多新鲜的想法，提倡首创精神，在开辟新思路的前提下再确定课题的研究方向。第二，忘我工作。无论是研究新技术、开发新课题，还是在其他的学习工作中，恩师都非常刻苦认真，一旦进入工作状态，就会潜心钻研非常忘我。第三，务实精神。在科研创新方面，说到做到、稳扎稳打；在用人方面，恩师非常重视实干能力，“尤其在关键岗位，一定要选有真本事的人”。