CN114353831B - 一种船载usb测角误差修正模型及其实现方法 - Google Patents
一种船载usb测角误差修正模型及其实现方法 Download PDFInfo
- Publication number
- CN114353831B CN114353831B CN202111646453.8A CN202111646453A CN114353831B CN 114353831 B CN114353831 B CN 114353831B CN 202111646453 A CN202111646453 A CN 202111646453A CN 114353831 B CN114353831 B CN 114353831B
- Authority
- CN
- China
- Prior art keywords
- usb
- model
- angle
- error correction
- measurement error
- 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
- 238000005259 measurement Methods 0.000 title claims abstract description 44
- 238000012937 correction Methods 0.000 title claims abstract description 41
- 238000000034 method Methods 0.000 title claims abstract description 19
- 238000012886 linear function Methods 0.000 claims description 5
- 238000012545 processing Methods 0.000 claims description 5
- 238000012935 Averaging Methods 0.000 claims description 3
- 230000000903 blocking effect Effects 0.000 claims description 3
- 238000012360 testing method Methods 0.000 claims description 2
- 238000012795 verification Methods 0.000 claims 1
- 230000003068 static effect Effects 0.000 abstract description 2
- 230000003287 optical effect Effects 0.000 description 3
- 239000011159 matrix material Substances 0.000 description 2
- XOFYZVNMUHMLCC-ZPOLXVRWSA-N prednisone Chemical compound O=C1C=C[C@]2(C)[C@H]3C(=O)C[C@](C)([C@@](CC4)(O)C(=O)CO)[C@@H]4[C@@H]3CCC2=C1 XOFYZVNMUHMLCC-ZPOLXVRWSA-N 0.000 description 2
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 230000000452 restraining effect Effects 0.000 description 1
- 238000007665 sagging Methods 0.000 description 1
- 101150114545 secE gene Proteins 0.000 description 1
- 238000000638 solvent extraction Methods 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Landscapes
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明涉及一种船载USB测角误差修正模型及其实现方法,具体公开了一种新的船载S频段统一测控***(简称“USB”)测角误差修正模型,旨在提供一种低成本高精度的海上USB设备误差标定方法。与现有USB轴系参数模型相比,新模型考虑了在动态条件下USB天线本身的变化对测角的影响,提高了USB测角误差修正精度及可靠性;模型建立标定在动态或静态情况都可开展,降低了模型建立外界客观条件的限制;利用卫星星历和船的位置信息获取理论方位俯仰角度,无需其它设备辅助跟踪即可建立USB测角误差修正模型,实用价值高。
Description
技术领域
本发明属于航天测控技术领域,涉及一种船载USB测角误差修正模型及其实现方法。
背景技术
USB统一测控***是测量船重要的测控设备,其测角精度及稳定性将直接影响最终的卫星测控精度,是衡量测量船测控服务能力的重要指标。理想的情况下USB天线实际指向和理论的机械指向是一致的,但是由于外界环境、信号参数、天线变形等因素的影响,通常情况下USB天线实际指向和理论的机械指向有一定的差异。特别是测量船在动态条件下开展角度测量,受外界环境、天线变形等因素的影响更加显著。
为了实现USB***海上高精度的角度测量,船载测控***每次在进行海上测控之前需要通对设备进行标校。目前,USB测角误差修正模型主要采用7参数或9参数轴系参数模型,即:一种是通过跟踪过境卫星的方式直接得到7个电轴修正参数;另一种方法是先以经纬仪为基准,通过同步测恒星法得到USB天线标校电视光轴修正参数,然后再以光轴为基准,通过放球的方式进行光电偏差的标定,得到电轴的参数。
通过大量历史数据的分析发现,简单的轴系参数模型不能很好体现USB***在不同方位、俯仰上的测角误差变化情况,导致实际误差修正的效果并不理想。因此,针对当前轴系参数模型的不完整,有必要给出一种能够考虑USB***在不同方位、俯仰上的测角误差变化情况的误差修正模型。
发明内容
本发明所要解决的技术问题是针对上述现有技术提供了一种基于北斗***的船载USB测角误差修正模型及其修正方法,通过抑制海上复杂条件下USB测角误差变化的影响,提高USB测角精度与可靠性。
本发明解决上述问题所采用的技术方案为:一种船载USB测角误差修正模型的实现方法,所述方法包括以下步骤:
步骤1:根据计划提前安排跟踪时间段和船舶所在位置;
步骤2:获目标卫星星历数据;
步骤3:开展目标卫星可见性预报,并判断该时间段是否满足USB分段线性格网模型参数估计观测条件;如果满足最优观测条件,可在预定时间段和海域开展USB轴系参数标校;
步骤4:利用目标星历数据和测量船船位信息生成USB***引导数据;
步骤5:根据引导数据引导USB设备跟踪目标卫星;
步骤6:根据新的船载USB测角误差修正模型,结合步骤5中记录的相关数据,获得USB测角误差修正模型参数;
步骤7:跟踪其他卫星或在其他时间段跟踪该卫星对USB测角误差修正模型进行测试验证。
优选地,步骤5中在跟踪过程中实时记录USB跟踪角度测量数据、船舶姿态数据、船舶位置数据、船体变形数据和气象数据。
优选地,所述步骤6中的修正模型具体为一种基于分段线性函数的格网模型,假设模型F(A,E)由不同的网格点构成,每个格网均由4个待求误差点组成落入该格网内的F(A,E)值可通过以方位角和俯仰角为自变量的分段线性函数模型F(A,E)表示为
式中F(Ai,Ej)、F(Ai,Ej+1)、F(Ai+1,Ej)、F(Ai+1,Ej+1)分别为对应的四个格的格网模型值,F(A,E)为方位角A、俯仰角E处对应的模型改正值,Ej和Ej+1为俯仰角E最近两个格网处的俯仰角,Ai和Ai+1为方位角A最近两个格网点处的方位角。
优选地,所述步骤6具体包括以下流程:
a.残差分块:将所有天的俯仰角残差按照方位角和俯仰角等间隔地划分为残差小块;
b.统计求取F(A,E)模型格网点的值,将方位角A和俯仰角E处周围小块内的所有残差值取平均,即获得F(A,E)模型格网点的值:
式中,res(A,E)表示在方位角A和俯仰角E处理论俯仰角和实测俯仰角之间的差值,dA和dE分别表示方位角A和俯仰角E方向的格网尺度,N表示在区域内理论俯仰角和实测俯仰角之间的差值数据的总个数;
c.迭代:考虑到在处理过程中F(A,E)会被受偶然误差的影响,进行多次迭代来不断修正F(A,E)模型,直到F(A,E)各网点模型系数不再有明显的变化停止迭代。
与现有技术相比,本发明的优点在于:
本发明提供一种新的船载USB测角误差修正模型,新模型考虑了在动态条件下USB天线本身的变化对测角的影响,提高了USB测角误差修正精度及可靠性;模型建立标定在动态或静态情况都可开展,降低了模型建立外界客观条件的限制;利用卫星星历和船的位置信息获取理论方位俯仰角度,无需其它设备辅助跟踪即可建立USB测角误差修正模型,实用价值高。
附图说明
图1船载USB测角误差修正模型网格点线性插值示意图。
图2船载USB测角误差修正模型网格点系数获取示意图。
图3一种新的船载USB测角误差修正模型估计流程图。
具体实施方式
以下结合附图实施例对本发明作进一步详细描述。
如图1所示,本实施例提供一种新的船载USB测角误差修正模型,该模型为基于分段线性函数的格网模型,假设模型F(A,E)由不同的网格点构成,每个格网均由4个待求误差点组成,如附图1所示,落入该格网内的F(A,E)值可通过以方位角和俯仰角为自变量的分段线性函数模型F(A,E)可表示为
式中F(Ai,Ej)、F(Ai,Ej+1)、F(Ai+1,Ej)、F(Ai+1,Ej+1)分别为对应的四个格的格网模型值,F(A,E)为方位角A、俯仰角E处对应的模型改正值,Ej和Ej+1为俯仰角E最近两个格网处的俯仰角,Ai和Ai+1为方位角A最近两个格网点处的方位角。
基于上述修正模型的具体建立包括以下步骤:
1)利用船载USB设备跟踪指定的目标(以北斗卫星为例),获得USB观测数据,对观测数据进行处理可获得USB跟踪目标的方位俯仰角:
E=Ec-βm·cos(A-Am)-Ce-ΔEZ-ΔEg·cosE (2)
A=Ac-βm·tanE·sin(A-Am)-δm·tanE-Sb·secE-CS·secE-ΔAZ·secE (3)
式中,E轴系误差修正后的俯仰角,Ec地面零值修正后的俯仰角,A轴系误差修正后的方位角,Ac地面零值修正后的方位角,βm大盘不水平最大倾斜量,Am大盘不水平最大倾斜方位;CS天线光电轴横向不匹配,Ce天线光电轴纵向不匹配,ΔAZ方位动态滞后,ΔEZ俯仰动态滞后,ΔEg重力下垂量,δm俯仰轴和方位轴不正交,Sb俯仰轴和光轴不正交,Ac=Ausb-A0,Ec=Eusb-E0,Ausb和Eusb为USB方位角和俯仰角实测值,A0和E0为方位零位和俯仰零位。
2)利用测量船和目标的位置信息可计算出在目标当地空间直角坐标系下的位置矢量enu:
其中,x,y,z为卫星在空间直角坐标系下的位置,x0,y0,z0为测量船在空间直角坐标系下位置,λ0、为测量船在大地测量坐标系下的经度、纬度。enu表示一种当地空间直角坐标系,其坐标原点为测量船位置x0,y0,z0,e表示卫星在enu坐标系东向分量(指向东为正),n表示卫星在enu坐标系东向分量(指向北为正),u表示卫星在enu坐标系东向分量(指向天顶为正)。
3)利用测量船姿态等信息将目标在当地空间直角坐标下的坐标转换到USB测量坐标系:
式中,表示USB设备与惯导之间的变形欧拉角旋转矩阵,/>表示船体姿态欧拉角旋转矩阵,dx,dy,dz表示USB设备与惯导之间的位置偏差。
4)计算USB跟踪目标相应时刻的理论方位角At、俯仰角Et:
5)计算实测方位俯仰角和理论方位俯仰角差值:
6)对实测俯仰角和理论俯仰角差值进行分块建模,求解模型系数。在实际处理过程中,采用残差法估计模型系数F(Ai,Ej)。首先需要设定模型的格网大小,然后再根据格网的尺度大小进行分块,最后使用统计的方法统计获取各格网点的F(A,E)值。先以俯仰角误差模型系数求解为例,估计分段线性格网模型系数F(Ai,Ej),估计流程如下:
a.残差分块。将所有天的俯仰角残差按照方位角和俯仰角等间隔地划分为残差小块,如中的红线包络区域。
b.统计求取F(A,E)模型格网点的值。将方位角A和俯仰角E处周围小块(红线包络区域)内的所有残差值取平均,即可获得F(A,E)模型格网点的值
式中,res(A,E)表示在方位角A和俯仰角E处理论俯仰角和实测俯仰角之间的差值,dA和dE分别表示方位角A和俯仰角E方向的格网尺度,N表示在区域内理论俯仰角和实测俯仰角之间的差值数据的总个数。
c.迭代。考虑到在处理过程中F(A,E)会被受偶然误差的影响,因此需要进行多次迭代来不断修正F(A,E)模型,直到F(A,E)各网点模型系数不再有明显的变化停止迭代,通常需要迭代多次。
7)使用实测方位角和理论方位角数据重复步骤6)获取方位角误差修正模型系数。
8)跟踪S频段其他卫星,使用公式(1)验证模型的正确性。
除上述实施例外,本发明还包括有其他实施方式,凡采用等同变换或者等效替换方式形成的技术方案,均应落入本发明权利要求的保护范围之内。
Claims (2)
1.一种船载USB测角误差修正模型的实现方法,其特征在于:所述方法包括以下步骤:
步骤1:根据计划提前安排跟踪时间段和船舶所在位置;
步骤2:获目标卫星星历数据;
步骤3:开展目标卫星可见性预报,并判断该跟踪时间段是否满足USB分段线性格网模型参数估计观测条件;如果满足最优观测条件,可在预定时间段和海域开展USB轴系参数标校;
步骤4:利用目标卫星星历数据和测量船船位信息生成USB***引导数据;
步骤5:根据引导数据引导USB设备跟踪目标卫星,实时记录USB跟踪角度测量数据、船舶姿态数据、船舶位置数据、船体变形数据和气象数据;
步骤6:根据新的船载USB测角误差修正模型,结合步骤5中记录的相关数据,获得USB测角误差修正模型参数;
步骤7:跟踪其他卫星或在其他时间段的该目标卫星,对USB测角误差修正模型进行测试验证;
所述步骤6中的修正模型具体为一种基于分段线性函数的格网模型,假设模型F(A,E)由不同的格网点构成,每个格网均由4个待求误差点组成;落入该格网内的F(A,E)值可通过以方位角和俯仰角为自变量的分段线性函数模型表示为
式中F(Ai,Ej)、F(Ai,Ej+1)、F(Ai+1,Ej)、F(Ai+1,Ej+1)分别为对应的四个格的格网模型值,F(A,E)为方位角A、俯仰角E处对应的模型改正值,Ej和Ej+1为俯仰角E最近的两个格网处的俯仰角,Ai和Ai+1为方位角A最近的两个格网点处的方位角。
2.根据权利要求1所述的一种船载USB测角误差修正模型的实现方法,其特征在于:所述步骤6具体包括以下流程:
a.残差分块:将所有天的俯仰角残差按照方位角和俯仰角等间隔地划分为残差小块;
b.统计求取F(A,E)模型格网点的值,将方位角A和俯仰角E处周围小块内的所有残差值取平均,即获得F(A,E)模型格网点的值:
res(A,E)表示在方位角A和俯仰角E处理论俯仰角和实测俯仰角之间的差值,dA和dE分别表示方位角A和俯仰角E方向的格网尺度,N表示在区域内理论俯仰角和实测俯仰角之间的差值数据的总个数;
c.迭代:考虑到在处理过程中F(A,E)模型会受偶然误差的影响,需要进行多次迭代来不断修正F(A,E)模型,直到F(A,E)各网点模型系数不再有明显的变化停止迭代。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111646453.8A CN114353831B (zh) | 2021-12-29 | 2021-12-29 | 一种船载usb测角误差修正模型及其实现方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111646453.8A CN114353831B (zh) | 2021-12-29 | 2021-12-29 | 一种船载usb测角误差修正模型及其实现方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114353831A CN114353831A (zh) | 2022-04-15 |
CN114353831B true CN114353831B (zh) | 2024-05-28 |
Family
ID=81103504
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111646453.8A Active CN114353831B (zh) | 2021-12-29 | 2021-12-29 | 一种船载usb测角误差修正模型及其实现方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114353831B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116519011B (zh) * | 2023-03-11 | 2024-03-01 | 中国人民解放军国防科技大学 | 基于Psi角误差修正模型的长航时双惯导协同标定方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR906489A (fr) * | 1944-02-23 | 1946-01-08 | Soc Ind Des Procedes Loth | Procédé et dispositif automatique de correction des erreurs quadrantales des cadres radiogoniométriques à bord des navires ou des aéronefs |
CN101782370A (zh) * | 2010-03-09 | 2010-07-21 | 哈尔滨工业大学 | 基于usb摄像头的测量定位方法及测量运动物体的运动轨迹的方法 |
CN106052717A (zh) * | 2016-07-15 | 2016-10-26 | 中国西安卫星测控中心 | 一种利用精密弹道的电波折射修正效果标定方法 |
KR20190139032A (ko) * | 2018-06-07 | 2019-12-17 | 한국항공우주산업 주식회사 | 이동체용 스트랩다운 장치의 정렬 오차 보정장치, 보정시스템 및 보정방법 |
CN210664749U (zh) * | 2019-09-06 | 2020-06-02 | 中国石油大学(华东) | 一种基于角度修正的称重仪 |
CN111580053A (zh) * | 2020-05-30 | 2020-08-25 | 中国人民解放军海军航空大学 | 大型船舶目标时变雷达散射中心测量误差联合计算及修正方法 |
JP2020193954A (ja) * | 2019-05-30 | 2020-12-03 | 住友電気工業株式会社 | 位置補正サーバ、位置管理装置、移動体の位置管理システム及び方法、位置情報の補正方法、コンピュータプログラム、車載装置並びに車両 |
-
2021
- 2021-12-29 CN CN202111646453.8A patent/CN114353831B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR906489A (fr) * | 1944-02-23 | 1946-01-08 | Soc Ind Des Procedes Loth | Procédé et dispositif automatique de correction des erreurs quadrantales des cadres radiogoniométriques à bord des navires ou des aéronefs |
CN101782370A (zh) * | 2010-03-09 | 2010-07-21 | 哈尔滨工业大学 | 基于usb摄像头的测量定位方法及测量运动物体的运动轨迹的方法 |
CN106052717A (zh) * | 2016-07-15 | 2016-10-26 | 中国西安卫星测控中心 | 一种利用精密弹道的电波折射修正效果标定方法 |
KR20190139032A (ko) * | 2018-06-07 | 2019-12-17 | 한국항공우주산업 주식회사 | 이동체용 스트랩다운 장치의 정렬 오차 보정장치, 보정시스템 및 보정방법 |
JP2020193954A (ja) * | 2019-05-30 | 2020-12-03 | 住友電気工業株式会社 | 位置補正サーバ、位置管理装置、移動体の位置管理システム及び方法、位置情報の補正方法、コンピュータプログラム、車載装置並びに車両 |
CN210664749U (zh) * | 2019-09-06 | 2020-06-02 | 中国石油大学(华东) | 一种基于角度修正的称重仪 |
CN111580053A (zh) * | 2020-05-30 | 2020-08-25 | 中国人民解放军海军航空大学 | 大型船舶目标时变雷达散射中心测量误差联合计算及修正方法 |
Also Published As
Publication number | Publication date |
---|---|
CN114353831A (zh) | 2022-04-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
WO2020233591A1 (zh) | 一种面向三维地表形变估计的InSAR和GNSS定权方法 | |
CN107367716B (zh) | 一种高精度星载sar几何定标方法 | |
CN110220491B (zh) | 一种无人机的光学吊舱安装误差角估算方法 | |
CN110940354B (zh) | 一种光电跟踪***捷联惯导安装姿态的标定方法 | |
CN103675773B (zh) | 一种定标器与卫星指向对准的确定方法 | |
CN105548976A (zh) | 船载雷达海上精度鉴定方法 | |
CN108562279B (zh) | 一种无人机测绘方法 | |
CN103913180A (zh) | 一种船载大视场高精度星敏感器的安装角标定方法 | |
CN106052717B (zh) | 一种利用精密弹道的电波折射修正效果标定方法 | |
CN111060943B (zh) | 一种符合实际的gnss观测值定权方法 | |
CN102346033B (zh) | 基于卫星观测角误差估计的直接定位方法及*** | |
CN105387874A (zh) | 船载高精度星敏感器安装角标定方法 | |
CN107490364A (zh) | 一种大角度倾斜成像航空相机对地目标定位方法 | |
CN109085655B (zh) | 一种水下平台重力测量方案与验证方法 | |
CN113624197B (zh) | 一种基于无人机平台的测控天线大盘不水平度测量方法 | |
CN114353831B (zh) | 一种船载usb测角误差修正模型及其实现方法 | |
KR20190049086A (ko) | 무인기 고속비행 중 정밀위치 영상 획득 장치 및 그를 이용한 정밀위치 획득 방법 | |
CN103322969B (zh) | 一种船载无线电测量设备天线变形测量方法 | |
CN117192578A (zh) | 一种跟踪无人机的船载测控天线轴系参数标定方法 | |
CN109708667B (zh) | 一种基于激光陀螺的双动态目标跟踪引导方法 | |
CN113900069A (zh) | 一种基于干涉成像高度计的垂线偏差计算方法及其*** | |
CN110579744B (zh) | 一种基于dom卫片的场景建模仿真方法 | |
CN114355396A (zh) | 一种基于北斗***的usb轴系参数标定方法 | |
CN113483739B (zh) | 海上目标位置测量方法 | |
CN113625311B (zh) | Gnss-r镜面反射点的瞬时海洋反射面模型修正定位方法和*** |
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 |