基于VMD算法的电子侦察信号预处理方法技术

一种基于变分模态分解VMD算法的电子侦察信号预处理方法，将实时的电子侦察频域信号，利用采样频率计算信号的中心频率，再通过变分模态分解VMD算法，得到信三个频域参数的更新公式，之后利用三个参数的更新公式更新每个信道的脉冲信号、信道的中心频率和拉格朗日乘子，最后得到预处理后的电子侦察脉冲信号。由于克服了现有技术将电子侦察信号信道中心频率均匀选取所造成的脉冲交叠与丢失的问题，使得本发明专利技术提高了电子侦察脉冲信号的分解精度。本发明专利技术简化了信号分解步骤，具有较高的信号分解效率，降低了所需处理电子侦察脉冲信号的时间复杂度。

【技术实现步骤摘要】
基于VMD算法的电子侦察信号预处理方法
本专利技术属于通信
，更进一步涉及雷达信号处理
中的一种基于变分模态分解VMD(VariationalModeDecomposition)算法的电子侦察信号预处理方法。本专利技术可用于电子侦察、电子支援和威胁告警系统中，对所接收到的电子侦察脉冲信号进行非均匀分解，便于进一步提取脉冲信号时频域特征，以供分析脉冲来波方向的目标信息。
技术介绍
电子侦察脉冲信号预处理是雷达电子对抗中一个重要组成部分，在电子支援和威胁告警系统中的发挥中重要的作用。信道化处理是电子侦察预处理重要组成之一，它对接收到的电子侦察脉冲信号在频域上进行分割，便于后续的识别与分选。非均匀的信道化处理是将电子侦察脉冲信号在频域上进行不均匀的分割，分解信号与噪声效果较好。吕卫祥等人在其发表的论文“一种改进的数字信道化结构设计及检测方法”(现代防御技术2018年第3期P73-79页)中提出了一种改进的信道化结构处理电子侦察脉冲信号。该方法利用半带滤波器插值后频谱重复复制的特点，将信号分解为频谱互补且无交叠的2路信号。然后采用奇偶划分的分析滤波器组对2路信号进行均匀信道化，在实现结构上将抽取因子与插值的半带滤波器互换位置。最后通过检测门限依据信道环境自适应的变检测长度的检测机制进行信道检测。该方法存在的不足之处是，使用均匀信道化处理电子侦察脉冲信号会造成频域信号的交叠和丢失，在处理复杂电子侦察脉冲信号时精度不高。北京无线电计量测试研究所在其申请的专利文献“一种机场跑道异物检测雷达信号预处理方法”(专利申请号：201611129998.0，申请公开号：106814359A)中提出了一种机场跑道异物检测电子侦察信号预处理方法。该方法首先建立理想中频信号的信号模型，并引入相位误差函数和相位误差的指数函数建立实测中频信号的信号模型，根据中频信号相位误差计算得到相位误差函数值，并根据相位误差函数值得到相位误差的指数函数值，之后根据实测中频信号的信号模型和相位误差函数的指数函数值计算得到电子侦察发射信号非线性校正后的中频信号，同时在距离频域对雷达发射信号非线性校正后的中频信号和相位误差的指数函数值进行匹配滤波，最后根据匹配滤波后的中频信号和相位误差的指数函数值计算得到频率偏移校正后的中频信号。该方法存在的不足之处是，对电子侦察脉冲信号进行预处理时将接收到的信号与构建的中频信号进行比对分解电子侦察脉冲信号，构建信号所用步骤繁杂，所需的时间复杂度较高。
技术实现思路
本专利技术的目的是针对上述现有技术的不足，提出了一种基于变分模态分解VMD算法的电子侦察预处理方法。实现本专利技术目的的思路是：将实时的电子侦察脉冲信号，通过模拟数字转换器转换模拟信号为数字信号，将数字信号经过离散傅里叶变换转换为频域信号，之后将信号分为若干个信道，同时利用采样频率计算信号的中心频率，再通过变分模态分解VMD算法，得到信号脉冲、中心频率和拉格朗日乘子三个参数的更新公式，之后利用三个参数的更新公式更新每个信道的脉冲信号、信道的中心频率和拉格朗日乘子，最后得到预处理后的电子侦察脉冲信号。为实现上述目的，本专利技术的具体步骤包括如下：(1)获取频域信号：将实时接收的电子侦察脉冲信号，通过模拟数字转换器转换模拟信号为数字信号，将数字信号经过离散傅里叶变换转换为频域信号；(2)利用中心频率计算公式，计算电子侦察脉冲信号中每个信道的中心频率值；(3)通过变分模态分解VMD算法，得到三个频域信号参数的更新公式：(3a)利用下述的带宽计算公式，得到由信道中心频率值和脉冲信号决定的信道带宽的最小值，其中每个信道的脉冲信号总和为接收到的频域信号：其中，Bk表示第k个电子侦察雷达信道的带宽，k＝1,...,K，表示带宽的最小值由uk和ωk决定，uk表示第k个电子侦察雷达信道的脉冲信号，ωk表示第k个雷达信道的中心频率，α表示电子侦察信号的噪声方差，||||2表示二范数操作，j表示虚数单位符号，ω表示待处理的电子侦察脉冲信号角频率，st.表示约束信道带宽最小值的条件，f表示待处理的频域信号，∑表示求和操作；(3b)利用下述的等价变换公式，得到带拉格朗日乘子的带宽计算公式：其中，表示积分操作，+∞表示正无穷大符号，||表示取模操作符号，λ表示拉格朗日乘子；(3c)分别对带拉格朗日乘子的带宽计算公式中的信号脉冲、信道中心频率和拉格朗日乘子求偏导，并且令三个偏导数等于零，得到信号脉冲、中心频率和拉格朗日乘子三个参数的更新公式；(4)利用电子侦察信号脉冲更新公式，更新当前每个电子侦察信道的脉冲信号；(5)利用电子侦察信道中心频率更新公式，更新当前每个电子侦察信道信道的中心频率；(6)判断是否更新完所有信道的电子侦察信道脉冲信号，若是，则执行步骤(7)；否则，执行步骤(4)；(7)利用拉格朗日乘子更新公式，更新当前的拉格朗日乘子；(8)判断更新后的每个信道的电子侦察脉冲信号是否满足收敛条件，若是，则执行步骤(9)；否则，用更新后的电子侦察脉冲信号替换当前的电子侦察脉冲信号后执行步骤(4)；(9)输出预处理后的电子侦察脉冲信号。本专利技术与现有技术相比较具有以下优点：第一，由于本专利技术通过变分模态分解VMD算法，得到的电子侦察信道中心频率更新公式，利用电子侦察信道中心频率更新公式更新中心频率，克服了现有技术将电子侦察脉冲信号信道中心频率均匀选取，出现脉冲交叠与丢失的问题，使得本专利技术解决了脉冲交叠和丢失的问题，提高了电子侦察脉冲信号的分解精度。第二，由于本专利技术通过变分模态分解VMD算法，得到的电子侦察信道脉冲信号更新公式，利用脉冲信号更新公式将更新后的电子侦察信道信号实时替换当前的电子侦察信道的频域信号，克服了现有技术构建信号步骤复杂，分解信号与噪声效果不佳的问题，使得本专利技术简化了分解步骤，具有较高的信号分解效率，降低了所需处理电子侦察脉冲信号的时间复杂度。附图说明图1是本专利技术的流程图；图2是本专利技术的仿真图。具体实施方式下面结合附图对本专利技术做进一步的详细描述。参照附图1，对本专利技术的具体步骤做进一步的详细描述。步骤1，获取频域信号。将实时接收的电子侦察脉冲信号，通过模拟数字转换器转换模拟信号为数字信号，将数字信号经过离散傅里叶变换转换为频域信号。步骤2，利用中心频率计算公式计算每个电子侦察信道的中心频率的初始值。所述的中心频率计算公式如下：其中，ωk表示第k个电子侦察雷达信道的中心频率，k表示信道的序号，k＝1,...,K，K表示根据电子侦察脉冲信号预处理所需的信道总数，ω0表示接收机的采样频率。步骤3，通过变分模态分解VMD算法，得到三个频域信号参数的更新公式。利用下述的带宽计算公式，得到由信道中心频率值和脉冲信号决定的信道带宽的最小值，其中每个信道的脉冲信号总和为接收到的频域信号：其中，Bk表示第k个电子侦察雷达信道的带宽，k＝1,...,K，表示带宽的最小值由uk和ωk决定，uk表示第k个电子侦察雷达信道的脉冲信号，ωk表示第k个雷达信道的中心频率，α表示电子侦察信号的噪声方差，||||2表示二范数操作，j表示虚数单位符号，ω表示待处理的电子侦察脉冲信号角频率，st.表示约束信道带宽最小值的条件，f表示待处理的频域信号，∑表示求和操作；利用下述的等价变换公式，得到带拉格朗日乘子的本文档来自技高网...

【技术保护点】
1.一种基于变分模态分解VMD算法的电子侦察信号预处理方法，其特征在于，通过变分模态分解VMD算法，得到三个频域信号参数的更新公式，利用三个参数更新公式更新电子侦察信号；该方法的步骤包括如下：(1)获取频域信号：将实时接收的电子侦察脉冲信号，通过模拟数字转换器转换模拟信号为数字信号，将数字信号经过离散傅里叶变换转换为频域信号；(2)利用中心频率计算公式，计算电子侦察脉冲信号中每个信道的中心频率值；(3)通过变分模态分解VMD算法，得到三个频域信号参数的更新公式：(3a)利用下述的带宽计算公式，得到由信道中心频率值和脉冲信号决定的信道带宽的最小值，其中每个信道的脉冲信号总和为接收到的频域信号：

【技术特征摘要】
1.一种基于变分模态分解VMD算法的电子侦察信号预处理方法，其特征在于，通过变分模态分解VMD算法，得到三个频域信号参数的更新公式，利用三个参数更新公式更新电子侦察信号；该方法的步骤包括如下：(1)获取频域信号：将实时接收的电子侦察脉冲信号，通过模拟数字转换器转换模拟信号为数字信号，将数字信号经过离散傅里叶变换转换为频域信号；(2)利用中心频率计算公式，计算电子侦察脉冲信号中每个信道的中心频率值；(3)通过变分模态分解VMD算法，得到三个频域信号参数的更新公式：(3a)利用下述的带宽计算公式，得到由信道中心频率值和脉冲信号决定的信道带宽的最小值，其中每个信道的脉冲信号总和为接收到的频域信号：其中，Bk表示第k个电子侦察雷达信道的带宽，k＝1,...,K，表示带宽的最小值由uk和ωk决定，uk表示第k个电子侦察雷达信道的脉冲信号，ωk表示第k个雷达信道的中心频率，α表示电子侦察信号的噪声方差，||||2表示二范数操作，j表示虚数单位符号，ω表示待处理的电子侦察脉冲信号角频率，st.表示约束信道带宽最小值的条件，f表示待处理的频域信号，∑表示求和操作；(3b)利用下述的等价变换公式，得到带拉格朗日乘子的带宽计算公式：其中，表示积分操作，+∞表示正无穷大符号，||表示取模操作符号，λ表示拉格朗日乘子；(3c)分别对带拉格朗日乘子的带宽计算公式中的信号脉冲、信道中心频率和拉格朗日乘子求偏导，并且令三个偏导数等于零，得到信号脉冲、中心频率和拉格朗日乘子三个参数的更新公式；(4)利用电子侦察信号脉冲更新公式，更新当前每个电子侦察信道的脉冲信号；(5)利用电子侦察信道中心频率更新公式，更新当前每个电子侦察信道信道的中心频率；(6)判断是否更新完所有信道的电子侦察信道脉冲信号，若是，则...

【专利技术属性】
技术研发人员：鲍丹，相巳琪，李鹏，蔡晶晶，
申请(专利权)人：西安电子科技大学，
类型：发明
国别省市：陕西,61