CN110780322A - 基于变步长梯度下降法的速镜面反射点计算方法及系统 - Google Patents

Abstract

本发明公开了基于变步长梯度下降法的镜面反射点计算方法及系统，所述方法包括：步骤1)计算第n次迭代的反射信号传播距离对镜面反射点的纬度导数和经度导数；步骤2)计算第n次迭代的曲率步长参数；步骤3)根据前后梯度关系计算第n次迭代的纬度步长参数和经度步长参数；步骤4)根据前三个步骤的计算结果，计算第n+1次迭代的镜面反射点的纬度和经度；步骤5)计算反射信号传播路径变化值，判断该值是否小于设定阈值，如果为是则进入步骤6)，否则，n加1后转入步骤1)；步骤6)将第n+1次迭代的镜面反射点的纬度和经度作为镜面反射点的位置。本发明的方法计算精度高，计算速度快，满足GNSS‑R接收机对精度和实时性的要求。

Description

基于变步长梯度下降法的速镜面反射点计算方法及系统

技术领域

本发明涉及全球导航卫星系统反射信号测量(Global Navigation SatelliteSystem Reflectometry,GNSS-R)领域，具体涉及基于变步长梯度下降法的镜面反射点计算方法及系统。

背景技术

GNSS-R接收机通过接收地球表面反射的导航卫星信号进行遥感探测。由于反射信号微弱，需要通过开环技术进行跟踪。镜面反射点计算是其核心算法。传统的计算方法忽略了地球长轴和短轴的差异计算，导致镜面反射点位置计算误差在公里量级。

现有的镜面反射点计算算法是TDS-1和CYGNSS卫星上使用的QSE(Quasi-Spherical Earth，近似圆地球)算法，该算法利用球型地球模型进行镜面反射点求解，然后通过坐标变化以及伸缩投影将结果转换到椭圆地球模型下。该方法的缺点在于模型误差大，计算得到的镜面反射点位置精度低，镜面反射点位置误差在公里级。

发明内容

本发明的目的在于克服上述技术缺陷，提出了基于变步长梯度下降法的镜面反射点计算方法，该方法基于采样梯度下降法求解镜面反射点位置，并结合变步长技术加快算法收敛速度，以实现米级误差的镜面反射点位置的快速计算。

为实现上述目的，本发明提出了一种基于变步长梯度下降法的镜面反射点计算方法，所述方法包括：

步骤1)计算第n次迭代的反射信号传播距离对镜面反射点的纬度导数和经度导数；

步骤2)计算第n次迭代的曲率步长参数；

步骤3)根据前后梯度关系计算第n次迭代的纬度步长参数和经度步长参数；

步骤4)根据前三个步骤的计算结果，计算第n+1次迭代的镜面反射点的纬度和经度；

步骤5)计算反射信号传播路径变化值，判断该值是否小于设定阈值，如果为是则进入步骤6)，否则，n加1后转入步骤1)；

步骤6)将第n+1次迭代的镜面反射点的纬度和经度作为镜面反射点的位置。

作为上述方法的一种改进，所述计算第n次迭代的反射信号传播距离对镜面反射点的纬度导数和经度导数；具体为：

第n次迭代的反射信号传播距离对镜面反射点的纬度导数

为：

第n次迭代的反射信号传播距离对镜面反射点的经度导数

为：

其中，Rⁿ为第n次迭代的反射信号传播距离，h是镜面反射点距离WGS84椭球面的高度，和

分别表示第n次迭代的镜面反射点的纬度和经度，(x_s ⁿ，y_s ⁿ，z_s ⁿ)是WGS84坐标系的镜面反射点三个坐标，通过

和

计算得到，(x_g，y_g，z_g)是WGS84坐标系的导航卫星三个坐标，(x_r，y_r，z_r)表示WGS84坐标系的接收机三个坐标，R_sr为镜面反射点与接收机的距离，R_sg为镜面反射点与导航卫星的距离；

E＝0.00669437990141；N＝a×k，a是地球长半轴6378137m。

作为上述方法的一种改进，所述计算第n次迭代的曲率步长参数，具体为：

第n次迭代的曲率步长参数σ_n为：

作为上述方法的一种改进，所述根据前后梯度关系计算第n次迭代的纬度步长参数和经度步长参数；具体包括：

为第n-1次迭代的纬度导数除以第n次迭代的纬度导数：

是第n-1次迭代的经度导数除以第n次迭代的经度导数：

n＝1时，镜面反射点的纬度

和经度

为接收机的纬度和经度；

判断

是否成立，

如果成立，则判断

是否成立，

如果成立，则

β_n＝α₃β_n-1，

否则，

β_n＝β_n-1；

否则，判断是否成立

如果成立，则

β_n＝α₄β_n-1，

否则，

β_n＝β_n-1；

其中，为第n次迭代的纬度步长参数；th1和th2是参数更新的两个阈值，α₁，α₂，α₃和α₄是经验参数，β_n为步长调整参数，β₀为已知的步长调整参数。

计算第n次迭代的纬度步长参数

的过程为：

判断

是否成立，

如果成立，则判断

是否成立，

如果成立，则

β_n＝α₃β_n-1，

否则，

β_n＝β_n-1；

否则，判断

是否成立

如果成立，则

β_n＝α₄β_n-1，

否则，

β_n＝β_n-1；

其中，

为第n次迭代的经度步长参数。

作为上述方法的一种改进，所述计算第n+1次迭代的镜面反射点的纬度和经度；具体为：

第n+1次迭代的镜面反射点的纬度和经度

其中，f_lat和f_lon为与初始误差相关的维度步长参数和经度步长参数；

作为上述方法的一种改进，所述纬度步长参数和经度步长参数的计算步骤为：

通过下式计算初始误差R_sa：

其中R_rg是接收机到导航卫星的水平距离，h_r是接收机的高度，h_g是导航卫星的高度，el是导航卫星相对于接收机的仰角；

然后通过初始误差R_sa估计纬度方向和经度方向的步长参数f_lat和f_lon：

其中K₁和K₂是经验参数，az表示导航卫星相对于接收机的方位角。

作为上述方法的一种改进，所述计算反射信号传播路径变化值，具体为：

根据第n+1次迭代的镜面反射点的纬度

和经度

计算其WGS84坐标系的三维坐标(x_s ⁿ⁺¹，y_s ⁿ⁺¹，z_s ⁿ⁺¹)：

根据下式计算第n+1次迭代的反射信号传播距离：

计算反射信号传播路径变化值ΔR：

ΔR＝|Rⁿ⁺¹-Rⁿ|。

本发明还提供了一种基于变步长梯度下降法的镜面反射点计算系统，所述系统包括：

导数计算模块，用于计算第n次迭代反射信号传播距离对镜面反射点的纬度导数和经度导数；

曲率步长参数计算模块，用于计算第n次迭代的曲率步长参数；

经纬度步长参数计算模块，用于根据前后梯度关系计算第n次迭代的纬度步长参数和经度步长参数；

镜面反射点经纬度计算模块，用于计算第n+1次迭代的镜面反射点的纬度和经度；

判断模块，用于计算反射信号传播路径变化值，判断该值是否小于设定阈值，如果为是则进入输出模块，否则，n加1后，进入导数计算模块；

输出模块，用于将第n+1次迭代的镜面反射点的纬度和经度为镜面反射点的位置并输出。

本发明的优势在于：

1、本发明的方法根据GNSS(全球导航卫星系统)接收机的定位信息预测GNSS信号在地球表面镜面反射点的位置，从而使接收机接收到微弱的反射信号；

2、本发明的方法计算精度高，计算速度快，满足GNSS-R接收机对精度和实时性的要求。

附图说明

图1为本发明的基于变步长梯度下降法的镜面反射点计算方法的流程图；

图2为本发明的第n次迭代的步长参数的计算流程图。

具体实施方式

下面结合附图对本发明的技术方案进行详细的说明。

镜面反射点的定义是使反射信号传播路径最短的反射点。梯度下降法是一种求解函数极小值的算法。通过在WGS84(世界大地坐标系1984)坐标系下建立反射信号传播路径的函数并在球坐标系中利用梯度下降法求解最小值从而得到镜面反射点的位置坐标。

梯度下降法稳定，计算精度高的特性。但是梯度下降法在接近最优值时收敛速度变慢，影响到算法的实时性。通过自适应步长调整技术，可以使算法快速收敛从而满足实时性的要求。步长的自适应调整依据以下三个参数：初始误差、经度圈的曲率、前后梯度的大小关系。

在迭代算法中，GNSS-R接收机的星下点被算作迭代计算的初始值，通过近似后的几何关系可以估计迭代算法的初始误差，从而估计所需要的调整量的大小。

由于算法在球坐标系中进行计算，在高位地区，经度圈的曲率变大，导致算法收敛变慢。通过经度圈曲率的平方来调整步长的大小可以有效地加快算法的收敛。

当算法迭代到最优值附近时，梯度接近0，收敛速度变慢，通过前后梯度的比值调整步长大小可以有效解决这个问题。

如图1所示，本发明的实施例1提出了一种基于变步长梯度下降法的镜面反射点计算方法，包括：

步骤1)将接收机的经纬度设定为镜面反射点的初始位置；

步骤2)利用几何关系估计与初始误差相关的步长参数f_lat和f_lon；

首先通过下式计算初始误差R_sa：

其中R_rg是接收机到导航卫星的水平距离，h_r是接收机的高度，h_g是导航卫星的高度，el是导航卫星相对于接收机的仰角；

然后通过初始误差R_sa估计经度和维度方向的步长参数f_lat和f_lon：

其中K₁和K₂是经验参数；az表示导航卫星相对于接收机的方位角。

步骤3)计算第n次迭代的反射信号传播距离对镜面反射点的纬度导数

和经度导数

第n次迭代的反射信号传播距离对镜面反射点的纬度导数

为：

第n次迭代的反射信号传播距离对镜面反射点的经度导数为：

其中，Rⁿ为第n次迭代的反射信号传播距离，h是镜面反射点距离WGS84椭球面的高度，和

分别表示第n次迭代的镜面反射点的纬度和经度，(x_s ⁿ，y_s ⁿ，z_s ⁿ)是WGS84坐标系的镜面反射点三个坐标，通过

和

计算得到，(x_g，y_g，z_g)是WGS84坐标系的导航卫星三个坐标，(x_r，y_r，z_r)表示WGS84坐标系的接收机三个坐标，R_sr为镜面反射点与接收机的距离，R_sg为镜面反射点与导航卫星的距离；

E＝0.00669437990141；N＝a×k，a是地球长半轴6378137m。

步骤4)计算第n次迭代的曲率步长参数σ_n：

步骤5)根据前后梯度关系计算第n次迭代的步长参数

和

为第n-1次迭代的纬度导数除以第n次迭代的纬度导数：

是第n-1次迭代的经度导数除以第n次迭代的经度导数：

n＝1时，镜面反射点的纬度和经度

为接收机的纬度和经度；

如图2所示，计算第n次迭代的纬度步长参数

的过程为：

判断

是否成立，

如果成立，则判断

是否成立，

如果成立，则

β_n＝α₃β_n-1，

否则，

β_n＝β_n-1；

否则，判断

是否成立

如果成立，则

β_n＝α₄β_n-1，

否则，

β_n＝β_n-1；

其中，

为第n次迭代的纬度步长参数；th1和th2是参数更新的两个阈值，α₁，α₂，α₃和α₄是经验参数，β_n为步长调整参数，β₀为已知的步长调整参数。

计算第n次迭代的纬度步长参数

的过程为：

判断

是否成立，

如果成立，则判断

是否成立，

如果成立，则

β_n＝α₃β_n-1，

否则，

β_n＝β_n-1；

否则，判断

是否成立

如果成立，则β_n＝α₄β_n-1，

否则，

β_n＝β_n-1；

其中，

为第n次迭代的经度步长参数。

步骤6)计算第n+1次迭代的镜面反射点的经度和纬度：

步骤7)计算反射信号传播路径变化值，如果小于设定阈值，则结束，否则，转入步骤3)。

计算反射信号传播路径变化值的步骤如下：

根据第n+1次迭代的镜面反射点的纬度和经度

计算其WGS84坐标系的三维坐标(x_s ⁿ⁺¹，y_s ⁿ⁺¹，z_s ⁿ⁺¹)：

根据下式计算第n+1次迭代的反射信号传播距离Rⁿ⁺¹：

计算反射信号传播路径变化值ΔR：

ΔR＝|Rⁿ⁺¹-Rⁿ|。

本发明的实施例2提出了一种基于变步长梯度下降法的镜面反射点计算系统，所述系统包括：

导数计算模块，用于计算第n次迭代反射信号传播距离对镜面反射点的纬度导数和经度导数；

曲率步长参数计算模块，用于计算第n次迭代的曲率步长参数；

经纬度步长参数计算模块，用于根据前后梯度关系计算第n次迭代的纬度步长参数和经度步长参数；

镜面反射点经纬度计算模块，用于计算第n+1次迭代的镜面反射点的纬度和经度；

判断模块，用于计算反射信号传播路径变化值，判断该值是否小于设定阈值，如果为是则进入输出模块，否则，n加1后，进入导数计算模块；

输出模块，用于将第n+1次迭代的镜面反射点的纬度和经度为镜面反射点的位置并输出

最后所应说明的是，以上实施例仅用以说明本发明的技术方案而非限制。尽管参照实施例对本发明进行了详细说明，本领域的普通技术人员应当理解，对本发明的技术方案进行修改或者等同替换，都不脱离本发明技术方案的精神和范围，其均应涵盖在本发明的权利要求范围当中。

Claims (8)

1.一种基于变步长梯度下降法的镜面反射点计算方法，所述方法包括：

步骤1)计算第n次迭代的反射信号传播距离对镜面反射点的纬度导数和经度导数；

步骤2)计算第n次迭代的曲率步长参数；

步骤3)根据前后梯度关系计算第n次迭代的纬度步长参数和经度步长参数；

步骤4)根据前三个步骤的计算结果，计算第n+1次迭代的镜面反射点的纬度和经度；

步骤5)计算反射信号传播路径变化值，判断该值是否小于设定阈值，如果为是则进入步骤6)，否则，n加1后转入步骤1)；

步骤6)将第n+1次迭代的镜面反射点的纬度和经度作为镜面反射点的位置。

2.根据权利要求1所述的方法，其特征在于，所述计算第n次迭代的反射信号传播距离对镜面反射点的纬度导数和经度导数；具体为：

第n次迭代的反射信号传播距离对镜面反射点的纬度导数为：

第n次迭代的反射信号传播距离对镜面反射点的经度导数

为：

其中，Rⁿ为第n次迭代的反射信号传播距离，h是镜面反射点距离WGS84椭球面的高度，

和

分别表示第n次迭代的镜面反射点的纬度和经度，(x_s ⁿ，y_s ⁿ，z_s ⁿ)是WGS84坐标系的镜面反射点三个坐标，通过

和

计算得到，(x_g，y_g，z_g)是WGS84坐标系的导航卫星三个坐标，(x_r，y_r，z_r)表示WGS84坐标系的接收机三个坐标，R_sr为镜面反射点与接收机的距离，R_sg为镜面反射点与导航卫星的距离；

E＝0.00669437990141；N＝a×k，a是地球长半轴6378137m。

3.根据权利要求2所述的方法，其特征在于，所述计算第n次迭代的曲率步长参数，具体为：

第n次迭代的曲率步长参数σ_n为：

4.根据权利要求3所述的方法，其特征在于，所述根据前后梯度关系计算第n次迭代的纬度步长参数和经度步长参数；具体包括：

为第n-1次迭代的纬度导数除以第n次迭代的纬度导数：

是第n-1次迭代的经度导数除以第n次迭代的经度导数：

n＝1时，镜面反射点的纬度

和经度

为接收机的纬度和经度；

判断

是否成立，

如果成立，则判断

是否成立，

如果成立，则

β_n＝α₃β_n-1，

否则，β_n＝β_n-1；

否则，判断

是否成立

如果成立，则

β_n＝α₄β_n-1，

否则，β_n＝β_n-1；

其中，

为第n次迭代的纬度步长参数；th1和th2是参数更新的两个阈值，α₁，α₂，α₃和α₄是经验参数，β_n为步长调整参数，β₀为已知的步长调整参数；

计算第n次迭代的纬度步长参数

的过程为：

判断是否成立，

如果成立，则判断是否成立，

如果成立，则

β_n＝α₃β_n-1，

否则，

β_n＝β_n-1；

否则，判断

是否成立

如果成立，则

β_n＝α₄β_n-1，

否则，

β_n＝β_n-1；

其中，

为第n次迭代的经度步长参数。

5.根据权利要求4所述的方法，其特征在于，所述计算第n+1次迭代的镜面反射点的纬度和经度；具体为：

第n+1次迭代的镜面反射点的纬度

和经度

其中，f_lat和f_lon为与初始误差相关的维度步长参数和经度步长参数。

6.根据权利要求5所述的方法，其特征在于，所述纬度步长参数和经度步长参数的计算步骤为：

通过下式计算初始误差R_sa：

其中R_rg是接收机到导航卫星的水平距离，h_r是接收机的高度，h_g是导航卫星的高度，el是导航卫星相对于接收机的仰角；

然后通过初始误差R_sa估计纬度方向和经度方向的步长参数f_lat和f_lon：

其中K₁和K₂是经验参数，az表示导航卫星相对于接收机的方位角。

7.根据权利要求6所述的方法，其特征在于，所述计算反射信号传播路径变化值，具体为：

根据第n+1次迭代的镜面反射点的纬度和经度计算其WGS84坐标系的三维坐标(x_s ⁿ⁺¹，y_s ⁿ⁺¹，z_s ⁿ⁺¹)：

根据下式计算第n+1次迭代的反射信号传播距离：

计算反射信号传播路径变化值ΔR：

ΔR＝|Rⁿ⁺¹-Rⁿ|。

8.一种基于变步长梯度下降法的镜面反射点计算系统，其特征在于，所述系统包括：

导数计算模块，用于计算第n次迭代反射信号传播距离对镜面反射点的纬度导数和经度导数；

曲率步长参数计算模块，用于计算第n次迭代的曲率步长参数；

经纬度步长参数计算模块，用于根据前后梯度关系计算第n次迭代的纬度步长参数和经度步长参数；

镜面反射点经纬度计算模块，用于计算第n+1次迭代的镜面反射点的纬度和经度；

判断模块，用于计算反射信号传播路径变化值，判断该值是否小于设定阈值，如果为是则进入输出模块，否则，n加1后，进入导数计算模块；

输出模块，用于将第n+1次迭代的镜面反射点的纬度和经度为镜面反射点的位置并输出。

Images