CN113670953B - 一种基于加权窗求直射与反射信号snr干涉峰值的方法 - Google Patents
一种基于加权窗求直射与反射信号snr干涉峰值的方法 Download PDFInfo
- Publication number
- CN113670953B CN113670953B CN202110961464.9A CN202110961464A CN113670953B CN 113670953 B CN113670953 B CN 113670953B CN 202110961464 A CN202110961464 A CN 202110961464A CN 113670953 B CN113670953 B CN 113670953B
- Authority
- CN
- China
- Prior art keywords
- sequence
- window
- interference
- direct
- signal
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N22/00—Investigating or analysing materials by the use of microwaves or radio waves, i.e. electromagnetic waves with a wavelength of one millimetre or more
- G01N22/04—Investigating moisture content
Landscapes
- Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Radar Systems Or Details Thereof (AREA)
- Noise Elimination (AREA)
Abstract
本发明涉及一种基于加权窗求直射与反射信号SNR干涉峰值的方法,所述方法包括:获取接收机接收到的直射信号与反射信号的干涉信号;对所述干涉信号进行单周期提取;对单周期提取后的干涉信号进行归一化处理,得到归一化序列;基于所述归一化序列计算窗的长度;基于所述窗的长度计算每次迭代的峰值,得到峰值序列;检测所述峰值序列的个数是否为1,若为1,则表示全部搜索完成,输出此时的峰值序列,若不为1,则重复上述步骤直到峰值序列的个数为1。本发明中的上述方法能够降低噪声对干涉信号所带来的的毛刺影响。
Description
技术领域
本发明涉及遥感领域,特别是涉及一种基于加权窗求直射与反射信号SNR干涉峰值的方法。
背景技术
利用GNSS(Global Navigation Satellites System)反射信号进行地球陆地物理参数的反演是遥感领域的新型技术之一,具有信号源广、设备便携,受大气影响小等优势。利用GNSS反射信号对土壤湿度进行测量能够有效降低成本,提高测量覆盖面,将传统的单点测量转换到面的测量,减少了所需要的测量设备,使得测量数据更有效利用在农林、地下水循环等相关领域。
土壤湿度,又称为土壤含水量,是用来表示土壤干湿程度大小的物理量。土壤湿度在研究陆地能量循环和水循环过程中是一项非常重要的变量,除了可以影响净辐射能量向潜热和感热转换的比率,还可以影响大气降水如雨、雪等转变为蒸发、下渗以及地表径流的比例。尽管土壤含水量在全球水资源中只占到0.0005%,但是对农作物生长的重要性是不可忽略的。
在传统的卫星信号接收中,接收机收到的除了视距传播的直射信号外,还有经过多个路径反射后的信号,又称为多径反射信号。在导航定位应用中通常把多径反射信号当作干扰进行抑制,因为其会对定位准确性产生影响,并且难以完全消除。但实际上反射信号与反射面的物理特征如土壤湿度具有相关性。单一天线接收到的GNSS信号包括了直射分量和反射的多径分量,由于地基情况下直射跟反射信号频率近似相等,只是传播路径长度不同,这两个分量会在接收机天线处发生较为稳定的干涉现象,形成干涉信号。在低高度角时,干涉振荡现象十分明显。测绘级接收机会将这种干涉信号以信噪比(SNR,signal-to-noise ratio)的形式记录存储在设备记录中。干涉现象的变化决定于土壤表面的反射特性,因此,可以利用这种干涉现象反演反射面土壤湿度等一些重要物理属性。SNR的数学模型可用下式表示:
只有反射信号的特性与反射面的物理特性有关,所以需要去掉干涉信号中的直射信号,得到反射信号的SNRm。传统的方法没有考虑周期信号中噪声所产生的毛刺影响,往往会输出多个伪峰值。本发明根据归一化干涉幅度进行加权,设置一个比较窗,该窗长度单位为卫星仰角,对窗内的数据进行比较,得到该窗内的局部峰值,作为该窗的特征值,然后将所有同样单位窗的特征值作为新的序列,改变窗的大小重复上述步骤进行比较。
发明内容
本发明的目的是提供一种基于加权窗求直射与反射信号SNR干涉峰值的方法,降低噪声对干涉信号所带来的的毛刺影响。
为实现上述目的,本发明提供了如下方案:
一种基于加权窗求直射与反射信号SNR干涉峰值的方法,所述方法包括:
获取接收机接收到的直射信号与反射信号的干涉信号;
对所述干涉信号进行单周期提取;
对单周期提取后的干涉信号进行归一化处理,得到归一化序列;
基于所述归一化序列计算窗的长度;
基于所述窗的长度计算每次迭代的峰值,得到峰值序列;
检测所述峰值序列的个数是否为1,若为1,则表示全部搜索完成,输出此时的峰值序列,若不为1,则重复上述步骤直到峰值序列的个数为1。
可选的,所述干涉信号的表达式为:
可选的,对所述干涉信号进行单周期提取具体包括:
计算所述干涉信号的均值;
判断当前序列值是否大于均值;
若否,序列左移,返回上一步骤;
若是,开始信号截取序列左移,序列计数器加1;
判断序列计数器是否大于周期T,若是,则输出截取序列;
若否,则返回步骤“开始信号截取序列左移,序列计数器加1”。
可选的,计算所述干涉信号的均值具体采用以下公式:
其中,θ表示GNSS卫星高度角,E[]表示求均值运算,SNR为接收机接收到的干涉信号的信噪比。
可选的,对单周期提取后的干涉信号进行归一化处理,得到归一化序列具体采用以下公式:
可选的,基于所述归一化序列计算窗的长度具体包括:
计算归一化的功率和Γ,Γ=sum(xp(θ));
当前序列和L置零;
L=L+当前序列x;
判断[L]=[Γ/K]是否成立;
若不成立,则序列左移,序列个数Λ+1,返回步骤“L=L+当前序列x”;
若成立,则判断窗个数n是否等于窗的大小K;
若是,则输出序列Λ;
若否,则窗的个数n+1,返回步骤“当前序列和L置零”。
可选的,基于所述窗的长度计算每次迭代的峰值,得到峰值序列具体包括:
建立矩阵M;所述矩阵M为Λj×K矩阵;
对矩阵M的每一行求和,得到一个L维列向量;
计算所述L维列向量最小值对应的下标;
取所述矩阵M的前r列,构成新矩阵Mjs;
计算所述新矩阵Mjs的每一列的标准差σi;
所有σi=0的元素对应的下标为原始序列的峰值,记为λ=[λ1,λ2,…λK]。
可选的,所述矩阵M的表达式如下:
可选的,计算所述新矩阵Mjs的每一列的标准差σi具体采用以下公式:
根据本发明提供的具体实施例,本发明公开了以下技术效果:
本发明中的上述方法可以有效去除噪声所带来的毛刺影响;
可以提升每次迭代的计算效率,减少所消耗的时间;
可以通过改变设计窗的个数来决定计算精度和计算效率,更加灵活方便。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例一种基于加权窗求直射与反射信号SNR干涉峰值的方法流程图;
图2为本发明实施例序列截取流程图;
图3为本发明实施例窗大小计算流程图;
图4为本发明实施例局部峰值计算流程图;
图5为本发明实施例一种基于加权窗求直射与反射信号SNR干涉峰值的系统结构示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明的目的是提供一种基于加权窗求直射与反射信号SNR干涉峰值的方法,降低噪声对干涉信号所带来的的毛刺影响。
为使本发明的上述目的、特征和优点能够更加明显易懂,下面结合附图和具体实施方式对本发明作进一步详细的说明。
图1为本发明实施例一种基于加权窗求直射与反射信号SNR干涉峰值的方法流程图,如图1所示,所述方法包括:
步骤101:获取接收机接收到的直射信号与反射信号的干涉信号。
所述干涉信号的表达式如下:
步骤102:对所述干涉信号进行单周期提取。
具体流程图如图2所示,步骤如下:
计算所述干涉信号的均值;
判断当前序列值是否大于均值;其中,当前序列值指的t1时刻的值。
若否,序列左移,返回上一步骤;其中,序列左移相当于时间右移,即处理下一个时刻的值。
若是,开始信号截取序列左移,序列计数器加1;从第一个均值点开始,截取一个周期T的长度,即为一个单周期。
判断序列计数器是否大于周期T,若是,则输出截取序列;
若否,则返回步骤“开始信号截取序列左移,序列计数器加1”。
步骤103:对单周期提取后的干涉信号进行归一化处理,得到归一化序列。
步骤104:基于所述归一化序列计算窗的长度,具体流程如图3所示。
记上一步骤中的归一化的功率和为Γ=sum(xp(θ)),则第一个窗的大小表示为:L1=x1+x2+…+xi,其包含的序列个数为:Λ1=i-1+1,使得[L1]=[Γ/K],同理可求得第二个窗的大小:L2=xi+xi+1+…+xs,其包含的序列个数为:Λ2=s-i+1,使得[L2]=[Γ/K],记所求的所有窗包含的序列个数为:Λ=[Λ1,Λ2,…,Λk],具体步骤如下:
计算归一化的功率和Γ,Γ=sum(xp(θ));
当前序列和L置零;
L=L+当前序列x;
判断[L]=[Γ/K]是否成立;
若不成立,则序列左移,序列个数Λ+1,返回步骤“L=L+当前序列x”;
若成立,则判断窗个数n是否等于窗的大小K;
若是,则输出序列Λ;
若否,则窗的个数n+1,返回步骤“当前序列和L置零”;
步骤105:基于所述窗的长度计算每次迭代的峰值,得到峰值序列,具体流程如图4所示。
得到所有窗包含的序列个数后,接着建立一个Λj×K矩阵M;
其中,矩阵中每一个元素的定义如下:
计算所述L维列向量最小值对应的下标:μj=arg min(γk j);
取所述矩阵M的前r列,构成新矩阵Mjs;
所有σi=0的元素对应的下标为原始序列的峰值,记为λ=[λ1,λ2,…λK]。
步骤106:检测所述峰值序列的个数是否为1,若为1,则表示全部搜索完成,输出此时的峰值序列,若不为1,则重复上述步骤直到峰值序列的个数为1。
图5为本发明实施例一种基于加权窗求直射与反射信号SNR干涉峰值的系统示意图,如图5所示,该系统包括:
单周期提取模块、归一化功率处理模块、窗设计模块、局部峰值求解模块、迭代控制模块。首先需要给定一个窗的个数K,其个数决定了计算精度和计算效率。干涉信号首先通过单周期提取模块,产生一个周期内的干涉信号;然后通过归一化功率处理模块,得到基于单周期均值的归一化功率;接着窗设计模块根据前置模块所得到的的归一化功率来设计窗的大小,归一化功率越小的位置,所设计的窗大小越大,窗大小与归一化功率成反比;根据前置模块所设计的窗,局部峰值求解模块对窗内的序列利用二分法进行峰值的寻找,将寻找到的峰值存储起来,记为该窗的1阶特征值;迭代控制模块在收到所有窗内序列一次迭代全部计算完成的信号后,将该次迭代的全部特征值合并成一个新的序列,重复进行以上的工作,直到窗的加权值为归一化的1,即窗已经覆盖了全部的序列,此时该窗的特征值即为峰值。
所述的归一化处理模块首先计算单个周期信号的均值,然后对单个周期内的全部周期信号进行归一化,并存储。
所述的窗设计模块根据窗的个数K与归一化的功率值设计窗的大小,窗的个数为K个,窗的大小与归一化的功率值成反比。
所述的局部峰值求解模块在单个窗内进行峰值寻找,在单个窗内利用二分法不断逼近峰值,由于所涉及的信号为正弦函数,所以设计二分法的加权系数比值为正切函数的比值,即为1:tan(1/2),对单个窗函数所求解的峰值进行存储,记为n阶特征值,n为迭代的次数。
所述的迭代控制模块控制前述三个模块,当局部峰值求解模块的所有窗全部计算完成后,向迭代控制模块发送迭代完成的信号,收到此信号后,迭代控制模块将局部分值求解模块的全部n阶特征值记为新的序列,然后控制局部峰值求解模块进行n+1次迭代。若第n次的局部峰值求解模块窗的加权系数为1,即已经覆盖到了全部序列信号后,迭代结束,迭代控制模块输出此时的n阶特征值。
本说明书中各个实施例采用递进的方式描述,每个实施例重点说明的都是与其他实施例的不同之处,各个实施例之间相同相似部分互相参见即可。对于实施例公开的系统而言,由于其与实施例公开的方法相对应,所以描述的比较简单,相关之处参见方法部分说明即可。
本文中应用了具体个例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想;同时,对于本领域的一般技术人员,依据本发明的思想,在具体实施方式及应用范围上均会有改变之处。综上所述,本说明书内容不应理解为对本发明的限制。
Claims (5)
1.一种基于加权窗求直射与反射信号SNR干涉峰值的方法,其特征在于,所述方法包括:
获取接收机接收到的直射信号与反射信号的干涉信号;
对所述干涉信号进行单周期提取;
对所述干涉信号进行单周期提取具体包括:
计算所述干涉信号的均值;
判断当前序列值是否大于均值;
若否,序列左移,返回上一步骤;
若是,开始信号截取序列左移,序列计数器加1;
判断序列计数器是否大于周期T,若是,则输出截取序列;
若否,则返回步骤“开始信号截取序列左移,序列计数器加1”;
对单周期提取后的干涉信号进行归一化处理,得到归一化序列;
对单周期提取后的干涉信号进行归一化处理,得到归一化序列具体采用以下公式:
基于所述归一化序列计算窗的长度;
基于所述归一化序列计算窗的长度具体包括:
计算归一化的功率和Γ,Γ=sum(xp(θ));
当前序列和L置零;
L=L+当前序列x;
判断[L]=[Γ/K]是否成立;
若不成立,则序列左移,序列个数Λ+1,返回步骤“L=L+当前序列x”;
若成立,则判断窗个数n是否等于窗的大小K;
若是,则输出序列Λ;
若否,则窗的个数n+1,返回步骤“当前序列和L置零”;
基于所述窗的长度计算每次迭代的峰值,得到峰值序列;
基于所述窗的长度计算每次迭代的峰值,得到峰值序列具体包括:
建立矩阵M;所述矩阵M为Λj×K矩阵;
对矩阵M的每一行求和,得到一个L维列向量;
计算所述L维列向量最小值对应的下标;
取所述矩阵M的前r列,构成新矩阵Mjs;
计算所述新矩阵Mjs的每一列的标准差σi;
所有σi=0的元素对应的下标为原始序列的峰值,记为λ=[λ1,λ2,L λK];
检测所述峰值序列的个数是否为1,若为1,则表示全部搜索完成,输出此时的峰值序列,若不为1,则重复上述步骤直到峰值序列的个数为1。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110961464.9A CN113670953B (zh) | 2021-08-20 | 2021-08-20 | 一种基于加权窗求直射与反射信号snr干涉峰值的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110961464.9A CN113670953B (zh) | 2021-08-20 | 2021-08-20 | 一种基于加权窗求直射与反射信号snr干涉峰值的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113670953A CN113670953A (zh) | 2021-11-19 |
CN113670953B true CN113670953B (zh) | 2022-07-22 |
Family
ID=78544612
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110961464.9A Active CN113670953B (zh) | 2021-08-20 | 2021-08-20 | 一种基于加权窗求直射与反射信号snr干涉峰值的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113670953B (zh) |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2015135295A1 (zh) * | 2014-03-12 | 2015-09-17 | 中兴通讯股份有限公司 | 一种干扰环境下随机接入信号的检测方法、装置和系统 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7688878B2 (en) * | 2006-03-16 | 2010-03-30 | The Boeing Company | Method and device of peak detection in preamble synchronization for direct sequence spread spectrum communication |
US20090088116A1 (en) * | 2007-09-28 | 2009-04-02 | Electronics And Telecommunications Research Institute | Method and apparatus of successive interference cancellation for wireless communication system |
CN106093339B (zh) * | 2016-05-23 | 2019-02-19 | 北京航空航天大学 | 一种基于gnss干涉功率峰谷值的土壤湿度测量方法 |
CN106371110B (zh) * | 2016-08-16 | 2018-12-25 | 上海航天测控通信研究所 | 一种gnss-r双补时延干涉处理系统及方法 |
CN110058281B (zh) * | 2019-04-29 | 2021-09-17 | 湖南国科微电子股份有限公司 | 动态定位方法及装置 |
CN112880633A (zh) * | 2021-01-12 | 2021-06-01 | 上海海洋大学 | 基于伯格算法的海面高度测量方法 |
-
2021
- 2021-08-20 CN CN202110961464.9A patent/CN113670953B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2015135295A1 (zh) * | 2014-03-12 | 2015-09-17 | 中兴通讯股份有限公司 | 一种干扰环境下随机接入信号的检测方法、装置和系统 |
Also Published As
Publication number | Publication date |
---|---|
CN113670953A (zh) | 2021-11-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US10571544B2 (en) | Direction finding using signal power | |
CN106226761B (zh) | 一种高性能相干高频雷达多频探测方法 | |
CN108549059B (zh) | 一种复杂地形条件下的低空目标仰角估计方法 | |
CN108459312B (zh) | 基于复合多径因子估计的加权多频最大似然仰角估计方法 | |
CN105223560B (zh) | 基于杂波俯仰方位谱稀疏恢复的机载雷达目标检测方法 | |
CN104515909B (zh) | 一种基于相关法的大天线方向图测量方法 | |
CN112414554B (zh) | 海面盐度获取方法、装置、设备及介质 | |
CN104569625B (zh) | 一种基于可转动辅助天线的大型天线方向图测量方法 | |
CN111487478B (zh) | 基于深度神经网络的角度依赖型复杂阵列误差校准方法 | |
CN106199600B (zh) | 基于多普勒估计的方位多通道合成孔径雷达成像方法 | |
CN110221241A (zh) | 一种基于rbf神经网络的低仰角doa估计方法 | |
CN110058222B (zh) | 一种基于传感器选择的双层粒子滤波检测前跟踪方法 | |
CN110389327A (zh) | 接收站位置误差下多站多外辐射源雷达双基距定位方法 | |
CN109507635A (zh) | 利用两个未知方位辅助源的阵列幅相误差估算方法 | |
CN109188019A (zh) | 基于多重信号分类算法的三维风速风向测量方法 | |
CN101644760A (zh) | 一种适用于高分辨阵列的快速鲁棒的信源个数检测方法 | |
CN104199020A (zh) | 基于多帧信息融合的米波阵列雷达目标仰角测量方法 | |
CN112383370A (zh) | 一种卫星时变信道的建模仿真方法及系统 | |
Zhao et al. | Altitude measurement of low elevation target based on iterative subspace projection | |
CN111198387A (zh) | 一种抗欺骗干扰的空时采样导航定位方法 | |
CN111352083A (zh) | 一种高频地波雷达多接收通道增益自动校准方法及装置 | |
CN113670953B (zh) | 一种基于加权窗求直射与反射信号snr干涉峰值的方法 | |
CN108761384A (zh) | 一种抗差的传感器网络目标定位方法 | |
CN117471397A (zh) | 一种基于图信号处理的圆形阵列二维doa估计方法 | |
CN115826004B (zh) | 一种基于二维角度及时差联合的三星协同直接定位方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |