技术特征：
1.一种基于加权窗求直射与反射信号snr干涉峰值的方法，其特征在于，所述方法包括：获取接收机接收到的直射信号与反射信号的干涉信号；对所述干涉信号进行单周期提取；对单周期提取后的干涉信号进行归一化处理，得到归一化序列；基于所述归一化序列计算窗的长度；基于所述窗的长度计算每次迭代的峰值，得到峰值序列；检测所述峰值序列的个数是否为1，若为1，则表示全部搜索完成，输出此时的峰值序列，若不为1，则重复上述步骤直到峰值序列的个数为1。2.根据权利要求1所述的基于加权窗求直射与反射信号snr干涉峰值的方法，其特征在于，所述干涉信号的表达式为：其中，a
d
表示直射信号辐射值，a
m
表示反射信号辐射值，表示直射信号与反射信号的相位差。3.根据权利要求1所述的基于加权窗求直射与反射信号snr干涉峰值的方法，其特征在于，对所述干涉信号进行单周期提取具体包括：计算所述干涉信号的均值；判断当前序列值是否大于均值；若否，序列左移，返回上一步骤；若是，开始信号截取序列左移，序列计数器加1；判断序列计数器是否大于周期t，若是，则输出截取序列；若否，则返回步骤“开始信号截取序列左移，序列计数器加1”。4.根据权利要求3所述的基于加权窗求直射与反射信号snr干涉峰值的方法，其特征在于，计算所述干涉信号的均值具体采用以下公式：其中，θ表示gnss卫星高度角，e[]表示求均值运算，snr为接收机接收到的干涉信号的信噪比。5.根据权利要求3所述的基于加权窗求直射与反射信号snr干涉峰值的方法，其特征在于，对单周期提取后的干涉信号进行归一化处理，得到归一化序列具体采用以下公式：其中，x(θ)表示干涉信号的半周期信号，x
p
(θ)表示归一化序列，表示干涉信号的均值。6.根据权利要求1所述的基于加权窗求直射与反射信号snr干涉峰值的方法，其特征在于，基于所述归一化序列计算窗的长度具体包括：给定窗的大小其中，[]为向上取整符号，n为归一化序列的长度；计算归一化的功率和γ，γ＝sum(x
p
(θ))；
当前序列和l置零；l＝l+当前序列x；判断[l]＝[γ/k]是否成立；若不成立，则序列左移，序列个数λ+1，返回步骤“l＝l+当前序列x”；若成立，则判断窗个数n是否等于窗的大小k；若是，则输出序列λ；若否，则窗的个数n+1，返回步骤“当前序列和l置零”。7.根据权利要求1所述的基于加权窗求直射与反射信号snr干涉峰值的方法，其特征在于，基于所述窗的长度计算每次迭代的峰值，得到峰值序列具体包括：建立矩阵m；所述矩阵m为λ
j
×
k矩阵；对矩阵m的每一行求和，得到一个l维列向量；计算所述l维列向量最小值对应的下标；取所述矩阵m的前r列，构成新矩阵m
js
；计算所述新矩阵m
js
的每一列的标准差σ
i
；所有σ
i
＝0的元素对应的下标为原始序列的峰值，记为λ＝[λ1,λ2,
…
λ
k
]。8.根据权利要求7所述的基于加权窗求直射与反射信号snr干涉峰值的方法，其特征在于，所述矩阵m的表达式如下：9.根据权利要求7所述的基于加权窗求直射与反射信号snr干涉峰值的方法，其特征在于，计算所述新矩阵m
js
的每一列的标准差σ
i
具体采用以下公式：其中，r表示前r列，表示矩阵元素。
技术总结
本发明涉及一种基于加权窗求直射与反射信号SNR干涉峰值的方法，所述方法包括：获取接收机接收到的直射信号与反射信号的干涉信号；对所述干涉信号进行单周期提取；对单周期提取后的干涉信号进行归一化处理，得到归一化序列；基于所述归一化序列计算窗的长度；基于所述窗的长度计算每次迭代的峰值，得到峰值序列；检测所述峰值序列的个数是否为1，若为1，则表示全部搜索完成，输出此时的峰值序列，若不为1，则重复上述步骤直到峰值序列的个数为1。本发明中的上述方法能够降低噪声对干涉信号所带来的的毛刺影响。所带来的的毛刺影响。所带来的的毛刺影响。
技术研发人员：杨东凯 李杰 汉牟田 王峰
受保护的技术使用者：北京航空航天大学
技术研发日：2021.08.20
技术公布日：2021/11/21