申请/专利权人：合肥工业大学

申请日：2022-02-17

公开（公告）日：2024-04-05

公开（公告）号：CN114488208B

主分类号：G01S19/21

分类号：G01S19/21;G06F18/10;G06F17/14

优先权：

专利状态码：有效-授权

法律状态：2024.04.05#授权;2022.05.31#实质审查的生效;2022.05.13#公开

摘要：本发明公开了一种联合经验小波和SPWVD变换的北斗信号抗干扰方法，包括：1、将含复杂干扰信号通过解析表示；2、将解析信号分解为按频率顺序排列的经验小波函数，不仅解决了常规小波基难以构造的缺陷，也有效避免复杂干扰信号分解时的模态混叠现象；3、对于经验小波变换后的信号，将有用信息和干扰信号进行分离和筛选；4、对于存在的相同频率信号分量之间的交叉项，采用平滑伪Wigner‑Ville分布SPWVD变换来进一步消除。本发明能有效地抑制交叉项的干扰和噪声的影响，并且保留了良好的时频定位特性，可广泛应用于GPS、GLONASS、Galileo、北斗导航卫星信号的处理，具有广阔的应用前景。

主权项：1.一种联合经验小波和SPWVD变换的北斗信号抗干扰方法，其特征是按如下步骤进行：步骤1：利用式1将北斗接收机接收到的原始连续时间信号进行解析变换，得到解析信号yat： 式1中，j是虚根单位，yt表示原始连续时间信号；表示yt经过希尔伯特变换后的信号；步骤2：将所述解析信号yat进行经验小波变换，得到频率不同、相互独立的各模态函数分量；步骤2.1：对所述解析信号yat进行傅里叶变换得到其对应的频谱yaω；步骤2.2：在所述频谱yaω上找到M个局部极大值点，并对其进行降序排列，得到排序后的局部极大值点；根据排序后的局部极大值点，在[0,π]上对所述频谱yaω划分为N个连续的部分，得到划分后的N个频谱段并形成N+1个端点；若MN，则保留排序后的前N-1个局部极大值点；若MN，则重置参数N，使得N＝M；令ωn表示相邻第n-1个和第n个频谱段的边界点，令初始边界点ω0＝0，末端边界点ωN＝π；n＝1,…,N；则第n个频谱段Λn＝[ωn-1,ωn]，从而得到其中，ωn-1表示相邻第n-2个和第n-1个频谱段的边界点；步骤2.3：利用式2和式3分别构造经验尺度函数的表达式和经验小波的小波函数表达式： 式2和式3中，φnω表示第n个经验尺度函数，ψnω表示第n个经验小波函数，γ表示0到1之间的常量，β·是一个值在[0,1]上的函数；ωn+1表示相邻第n个和第n+1个频谱段的边界点；ω是角频率；步骤2.4：利用式4和5分别定义细节相关系数和近似相关系数Wfε0,t： 式4中，“—”表示共轭运算，当n＝1时，φ1t由φnω傅里叶反变换得到函数值；ψnt由ψnω傅里叶反变换得到的函数，yat由yaω傅里叶反变换得到的函数；F-1·表示傅里叶反变换；τ是积分变量；步骤2.5：利用式6对解析信号yat进行重构，从而得到解析信号yat的各经验模态分量如式7和式8所示： 式6中，*表示卷积运算，是的傅里叶变换形式，表示n＝0时的值；式7中，y0t表示n＝0时的调幅-调频单分量，表示的傅里叶反变换形式；式8中，ynt表示不同频率的调幅-调频分量，n＝1,2,3,...,N；步骤3：利用式9对解析信号yat的各模态分量ynt进行平滑伪Wigner-Ville变换，得到平滑伪Wigner-Ville时频分布表达式 式9中，和表示ynt的瞬时相关函数，表示的共轭；gs-t是时域窗函数，hτ是频域窗函数，且g0＝h0＝1，s和τ均为积分变量；ω是角频率；步骤4：利用式10将各个模态分量ynt的平滑伪Wigner-Ville分析结果线性叠加，得到解析信号yat的平滑伪Wigner-Ville分布 步骤5：对所述解析信号yat的滑伪Wigner-Ville分布进行干扰检测，得到的干扰信号再采用陷波IIR滤波器进行滤波，从而得到滤波后的滑伪Wigner-Ville分布。

全文数据：

权利要求：

百度查询： 合肥工业大学 联合经验小波和SPWVD变换的北斗信号抗干扰方法

免责声明
1、本报告根据公开、合法渠道获得相关数据和信息，力求客观、公正，但并不保证数据的最终完整性和准确性。
2、报告中的分析和结论仅反映本公司于发布本报告当日的职业理解，仅供参考使用，不能作为本公司承担任何法律责任的依据或者凭证。

阅读全文 双屏查看 官方信息 专利公告 收藏专利 下载PDF 下载WORD