CN113267794B - 一种基线长度约束的天线相位中心校正方法及装置 - Google Patents
一种基线长度约束的天线相位中心校正方法及装置 Download PDFInfo
- Publication number
- CN113267794B CN113267794B CN202110819933.3A CN202110819933A CN113267794B CN 113267794 B CN113267794 B CN 113267794B CN 202110819933 A CN202110819933 A CN 202110819933A CN 113267794 B CN113267794 B CN 113267794B
- Authority
- CN
- China
- Prior art keywords
- antenna
- phase center
- gnss
- baseline
- base line
- 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
- 238000000034 method Methods 0.000 title claims abstract description 32
- 238000005457 optimization Methods 0.000 claims abstract description 15
- 230000008859 change Effects 0.000 claims abstract description 9
- 230000004927 fusion Effects 0.000 claims abstract description 8
- 230000003068 static effect Effects 0.000 claims abstract description 3
- 238000012360 testing method Methods 0.000 claims description 18
- 239000013598 vector Substances 0.000 claims description 16
- 230000006870 function Effects 0.000 claims description 10
- 238000001914 filtration Methods 0.000 claims description 4
- 238000010801 machine learning Methods 0.000 claims description 4
- 230000008569 process Effects 0.000 description 7
- 239000011159 matrix material Substances 0.000 description 5
- 230000001934 delay Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 238000009499 grossing Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/13—Receivers
- G01S19/23—Testing, monitoring, correcting or calibrating of receiver elements
- G01S19/235—Calibration of receiver components
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Variable-Direction Aerials And Aerial Arrays (AREA)
Abstract
Description
技术领域
本发明涉及天线测量技术领域,尤其是一种基线长度约束的天线相位中心校正方法及装置。
背景技术
全球卫星导航接收机通过天线接收卫星信号,天线的接收空间是个半球形或球形,卫星信号通过不同的入射角度进入时,会有不同的时间延迟。使用一个虚拟球面拟合所有方向的延迟,则该球面称为等相位面,该球心称为等相位中心或平均相位中心。这样,天线参考点和平均相位中心的偏差称为天线相位中心偏差;某个具体角度的入射信号的瞬时相位中心和平均相位中心的偏差称为天线相位中心变化。对于高精度GNSS定位,接收机天线的PCO、PCV是重要的误差来源,因此需要进行建模和校正。
对于某些固定基线长度的应用,例如GNSS定向、定姿应用,两个或多个天线安装在载体上面随着载体一起运动,天线间的基线长度存在约束,约束为固定不变,基线方向随载体运动而改变。对于这类应用无需分别估计两个天线的相位中心参数,而将两个天线相位中心作为一个整体进行建模和校正。
参考中国专利公开号为中国专利CN104502926A以及中国专利公开号为CN107290762B公开的一种室外绝对校正方法,这些方法的共同点是:参考天线固定不动,被测天线转动,从而实现对一个被测天线相位中心的校正。需要高精度的机械臂,机械臂每次运动后输出其末端的精确位置或机械臂坐标系,这样才能得到精确的基线估计值。由于需要高精度的机械臂,使得校正成本高,操作难度也增大。
发明内容
本发明解决了现有技术校正成本高和操作难度大的问题,提出一种基线长度约束的天线相位中心校正方法及装置,本发明针对基线长度固定的应用,将A、B两个天线一起固定在转动装置上运动,并且将两个天线相位中心作为一个整体建模和校正,另外,使用传感器和信息融合优化算法实现精确的基线最优估计值,因此无需高精度机械臂。
为实现上述目的,提出以下技术方案:
一种基线长度约束的天线相位中心校正方法,包括以下步骤:
S1,针对基线约束的GNSS应用,将GNSS的天线A和天线B的综合天线相位中心变化作为一个整体进行建模和校正;
S2,在静止状态下计算天线A和天线B的初始的状态估计值;
S3,保持天线A和天线B的基线长度不变,并以一定转动方向转动天线A和天线B,所述转动方向包括水平方向和倾斜方向;
本发明在基线约束的GNSS应用情况下,将两个GNSS天线都安装在转动装置上运动,天线A和天线B的转动依赖于转动装置,本发明增加传感器模块采集状态数据,然后由转动装置控制天线A和天线B转动并采集数据,转动过程中保持天线A和天线B的基线长度不变。使用传感器数据构建动力学方程计算被测设备的状态估计值。本发明将两个天线的相位变化作为一个整体进行建模和校正,整个过程中无需高精度机械臂,使用一般的机械臂、转台、可活动支架等转动装置,甚至也可以由人手工操作,具备成本低、效率高和精度高的优点。
作为优选,所述S1具体包括将天线A和天线B的综合天线相位中心误差函数作为一个整体进行建模和校正。
作为优选,所述水平方向转动范围为0~360°,所述倾斜方向转动范围为0~90°。
作为优选,所述转动为连续转动或以固定间隔进行转动。
作为优选,所述S4具体包括以下步骤:
所述信息融合和最优化算法包括最小二乘法、卡尔曼滤波算法、图优化算法和机
器学习算法,直接估计最优的基线最优估计值,或分步进行估计,先估计最优化的长度,然后估计最优化的单位方向,最后将两者合并得到基线最优估计值:
作为优选,所述GNSS差分模型包括单差模型、双差模型和三差模型。
作为优选,所述天线A和天线B综合的天线相位中心变化进行建模采用分段线性模型或二维线性插值模型或球谐函数模型。
一种基线长度约束的天线相位中心校正装置,适用于上述的一种基线长度约束的GNSS接收机天线相位中心校正方法,包括GNSS定向应用设备、测试基台和传感器模块,所述GNSS定向应用设备包括天线A和天线B,所述天线A和天线B的基线向量102的长度始终保持不变,所述传感器模块用于获取天线A和天线B的状态数据,所述测试基台上设有转动装置,所述转动装置与GNSS定向应用设备固定连接。
本发明的有益效果是:本发明将两个天线的相位变化作为一个整体进行建模和校正,整个过程中无需高精度机械臂,使用一般的机械臂、转台、可活动支架等转动装置,甚至也可以由人手工操作,具备成本低、效率高和精度高的优点。
附图说明
图1是实施例GNSS定向应用设备的结构示意图;
图2是实施例校正装置测试状态示意图;
图3是实施例校正方法的流程图;
101、GNSS定向应用设备 102、基线向量 103、天线B 104、天线A 105、传感器模块 106、载体 201、转动装置 202、测试基台 102a、t-1时刻的基线向量 102b、t时刻的基线向量。
具体实施方式
实施例1:
本实施例提出一种基线长度约束的GNSS应用形式,参考图1,在载体106上安装GNSS天线A 104、天线B 103。载体106为汽车、无人机或机械等,由于天线是固定在载体上,载体是刚体,所以两个天线的基线向量102的长度始终保持不变。但是随着载体的运动,例如汽车行驶中转弯,基线向量102的方向会随之改变,从而能够指示出载体的运动状态。由两个天线共同构成了一个GNSS定向应用设备101。另外还可以增加传感器模块105,包括加速度计、陀螺仪、磁力计等。
本发明只考虑天线相位中心误差情况下,单个GNSS天线A 104、天线B 103对卫星i的测量值:
其中:表示测量值,表示天线平均相位中心到卫星的真实距离,表示天
线A的天线相位中心误差函数,表示卫星的入射角(仰角和方位角),因为
天线A、B的基线长度远小于到卫星的距离,可认为同一颗卫星对两个天线的入射角度是相
同的。在进行相对定位时需要使用差分模型,包括单差、双差、三差等。考虑最简单的单差模
型:
实施例2:
本实施例在实施例1的基础上,提出一种基线长度约束的天线相位中心校正装置,包括GNSS定向应用设备101、测试基台202和传感器模块105,所述GNSS定向应用设备101包括天线A 104和天线B 103,所述天线A 104和天线B 103的基线向量102的长度始终保持不变,所述传感器模块105用于获取天线A和天线B的状态数据,所述测试基台202上设有转动装置201,所述转动装置201与GNSS定向应用设备101固定连接。在开阔环境中,在测试基台202上安装转动装置201,该转动装置201可以是精密机械臂、普通机械臂、转台、可转动的支架甚至是由人工操作,将基线长度固定的GNSS定向应用设备101安装在转动装置末端。
测试过程中控制转动装置201不断的转动, 图中显示了从 t-1 时刻到 t 时刻的
转动。 假设已知t-1时刻被测设备101的状态估计值,其中向量102a表
示基线估计值。对包含在设备101中传感器模块105中的加速度计、陀螺仪的输出在t-1
至t时刻进行积分得到两个时刻间的位置和姿态变化量,那么,被测设备在t时刻的状
态估计值,其中向量102b表示基线估计值,t时刻的状态估计值表示为:
其中加号表示广义加法,具体对于基线矢量而言表示矢量加法,对于姿态而言
是旋转矩阵乘法或四元数旋转变换,以此类推可以得到任意第i个时刻的被测设备的状态
估计值。对于被测设备的初始状态估计值可以在测试转动开始前计算:基线估计值
用GNSS卫星定位得到;姿态估计值用加速度计、磁力计计算得到。
注意在上述所有时刻中,真实的基线长度是相同的,即有如下约束:
因此利用“真实基线长度都相同”这个约束对估计基线进行优化。
实施例3:
本实施例在实施例1和实施例2的基础上提出一种基线长度约束的天线相位中心校正方法,参考图3,包括以下步骤:
S301,做好测试前的准备工作。在开阔环境中,将转动装置安装在测试基台上,该
转动装置可以是精密机械臂、普通机械臂、转台、可转动的支架甚至是由人工操作。将基线
长度固定的被测设备安装在转动装置末端。被测设备包括两个天线和传感器模块,计算被
测设备的初始状态估计值。
S302,转动测试并采集数据过程。转动包括两个方向:水平方向0~360°,倾斜方向0~90°,可以连续转动,也可以按固定间隔,如5°转动。采集的数据包括天线A和天线B对应的GNSS接收机卫星观测值、传感器模块的输出值。传感器模块包括加速度计、陀螺仪、磁力计,另外,根据所使用转动装置的不同,采集的数据还包括转动装置末端坐标、齿轮转数、分度盘角度、视频和图像等。
S303,使用信息融合最优化算法估计基线的最优估计值。图1中基线102是一个三
维向量,包括长度、单位方向,在最优化算法中直接估计最优的,或分
开估计,分别得到最优化的长度和最优化的单位方向,然后将两者合并得到基线最
优估计值。本实施例以分步最优估计为例介绍。
计算最优化基线长度。
根据公式:
当然也可以使用其他最优估计方法,比如最小二乘法、卡尔曼滤波(平滑)、图优化、机器学习等。
计算最优化基线方向。
计算整个过程中被测设备所有时刻的姿态估计值序列:
其中,表示传感器输出值,包括加速度计、陀螺仪、磁力
计。另外,还得到一些观测量,例如:由卫星信号定位得到两个天线的基线、转动装置末端
坐标、机械臂的齿轮转数、转台分度盘的转动角度、视频图像、检测到的特殊位置,如和之前
某个运动状态相同的位置、相垂直的位置等,将这些观测量统称为。由状态估计量和观测
量之间的差异表示为:
其中:
一个卫星可以构造一个单差载波观测值方程,天线A和B可以同时接收多个卫星,那么可以构造多个单差载波观测值方程,将它们合在一起写成矩阵形式如下:
因此,定义观测值残差估计值矩阵为:
S305,用观测值残差估计值对天线A和天线B的综合相位中心变化进行建模和校正,建模过程使用分段线性模型、二维线性插值模型或球谐函数模型。
本实施例使用一个简单的模型为例:只和仰角有关的一维多项式函数:
综合可得:
以最小二乘法为例求解:
Claims (7)
1.一种基线长度约束的天线相位中心校正方法,其特征是,包括以下步骤:
S1,针对基线约束的GNSS应用,将GNSS的天线A和天线B的综合天线相位中心变化作为一个整体进行建模和校正;具体包括将天线A和天线B的综合天线相位中心误差函数作为一个整体进行建模和校正;
S2,在静止状态下计算天线A和天线B的初始的状态估计值;
S3,保持天线A和天线B的基线长度不变,并以一定转动方向转动天线A和天线B,所述转动方向包括水平方向和倾斜方向;
2.根据权利要求1所述的一种基线长度约束的天线相位中心校正方法,其特征是,所述水平方向转动范围为0~360°,所述倾斜方向转动范围为0~90°。
3.根据权利要求1或2所述的一种基线长度约束的天线相位中心校正方法,其特征是,所述转动为连续转动或以固定间隔进行转动。
5.根据权利要求1所述的一种基线长度约束的天线相位中心校正方法,其特征是,所述GNSS差分模型包括单差模型、双差模型和三差模型。
6.根据权利要求1所述的一种基线长度约束的天线相位中心校正方法,其特征是,所述天线A和天线B综合的天线相位中心变化进行建模采用分段线性模型或二维线性插值模型或球谐函数模型。
7.一种基线长度约束的天线相位中心校正装置,适用于权利要求1所述的一种基线长度约束的GNSS接收机天线相位中心校正方法,其特征是,包括GNSS定向应用设备(101)、测试基台(202)和传感器模块(105),所述GNSS定向应用设备(101)包括天线A(104)和天线B(103),所述天线A(104)和天线B(103)的基线向量(102)的长度始终保持不变,所述传感器模块(105)用于获取天线A和天线B的状态数据,所述测试基台(202)上设有转动装置(201),所述转动装置(201)与GNSS定向应用设备(101)固定连接。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110819933.3A CN113267794B (zh) | 2021-07-20 | 2021-07-20 | 一种基线长度约束的天线相位中心校正方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110819933.3A CN113267794B (zh) | 2021-07-20 | 2021-07-20 | 一种基线长度约束的天线相位中心校正方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113267794A CN113267794A (zh) | 2021-08-17 |
CN113267794B true CN113267794B (zh) | 2021-10-29 |
Family
ID=77236851
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110819933.3A Active CN113267794B (zh) | 2021-07-20 | 2021-07-20 | 一种基线长度约束的天线相位中心校正方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113267794B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113899957B (zh) * | 2021-09-30 | 2024-03-29 | 武汉大学 | 基于双差相位观测值的室外天线绝对标定方法及系统 |
CN117075166B (zh) * | 2023-10-17 | 2023-12-15 | 智慧司南(天津)科技发展有限公司 | 一种船用卫星罗经艏向平滑方法 |
CN117724141B (zh) * | 2024-02-18 | 2024-04-16 | 中国民航大学 | 一种面向飞行器的北斗正交基线测向的角度优化方法 |
Family Cites Families (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7855681B2 (en) * | 2008-11-19 | 2010-12-21 | Harris Corporation | Systems and methods for determining element phase center locations for an array of antenna elements |
CN102323489B (zh) * | 2011-08-18 | 2013-09-11 | 北京华力创通科技股份有限公司 | 对接收机天线的相位中心进行校准的方法和装置 |
CN104502926B (zh) * | 2014-12-17 | 2017-10-24 | 武汉大学 | 一种基于精密机械臂的室外gnss接收机绝对天线相位中心校正方法 |
CN106443722B (zh) * | 2016-09-09 | 2019-05-17 | 中国人民解放军国防科学技术大学 | 一种天线相位中心偏差检测方法 |
CN107290762B (zh) * | 2016-09-22 | 2018-09-07 | 武汉大学 | 基于载波三次差分的室外gnss绝对天线相位中心标校方法 |
CN107315185B (zh) * | 2017-08-29 | 2023-07-04 | 中国地震局第一监测中心 | 连续旋转法gps天线相位中心检测系统及方法 |
CN109031359B (zh) * | 2018-06-11 | 2022-05-13 | 武汉大学 | 一种基于微波暗室的gnss接收天线绝对相位中心校正方法 |
CN111381256B (zh) * | 2020-03-10 | 2022-07-26 | 上海卫星工程研究所 | 主动遥感卫星天线相位中心偏移误差计算的方法和系统 |
CN112285747B (zh) * | 2020-10-16 | 2022-04-22 | 中国地震局第一监测中心 | Gnss接收机天线相位中心在线测量装置和测量方法 |
-
2021
- 2021-07-20 CN CN202110819933.3A patent/CN113267794B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN113267794A (zh) | 2021-08-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113267794B (zh) | 一种基线长度约束的天线相位中心校正方法及装置 | |
CN108375382B (zh) | 基于单目视觉的位置姿态测量系统精度检校方法和装置 | |
CN106643792B (zh) | 惯性测量单元和地磁传感器整体标定装置及标定方法 | |
CN107270893B (zh) | 面向不动产测量的杆臂、时间不同步误差估计与补偿方法 | |
Hol | Sensor fusion and calibration of inertial sensors, vision, ultra-wideband and GPS | |
CN106405670B (zh) | 一种适用于捷联式海洋重力仪的重力异常数据处理方法 | |
CN113311436B (zh) | 一种移动平台上激光测风雷达运动姿态测风订正方法 | |
JP5084303B2 (ja) | 移動体姿勢計測装置 | |
CN108363078B (zh) | 用于导航定位系统的动态定位误差测试装置、系统和方法 | |
CN110133692B (zh) | 惯导技术辅助的高精度gnss动态倾斜测量系统及方法 | |
JP2020528553A (ja) | 磁力計を較正する方法 | |
CN110926468A (zh) | 基于传递对准的动中通天线多平台航姿确定方法 | |
CN108106637A (zh) | 一种分布式pos的精度检校方法和装置 | |
CN110044377B (zh) | 一种基于Vicon的IMU离线标定方法 | |
CN110793518A (zh) | 一种海上平台的定位定姿方法及系统 | |
JP2019120587A (ja) | 測位システム及び測位方法 | |
CN105892498A (zh) | 基于三轴云台的目标凝视与扫视控制系统 | |
CN107807375B (zh) | 一种基于多gps接收机的无人机姿态追踪方法及系统 | |
CN113552599B (zh) | 基于信息融合的gnss接收机天线相位中心校正方法及装置 | |
CN117031513A (zh) | 一种道路及附属物的实时监测定位方法、系统及装置 | |
CN110954081A (zh) | 一种磁罗盘快速校准装置及方法 | |
CN110954080A (zh) | 一种消除载体磁干扰的磁罗盘校准方法 | |
CN112902950B (zh) | 一种面向低速运动载体中mems级imu的初始对准方法 | |
CN105928513B (zh) | 一种基于位置姿态测量系统的机载合成孔径雷达运动参数测量方法 | |
CN113064208B (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 |