CN110887509A - 一种仿复眼偏振传感器多方向标定方法 - Google Patents

一种仿复眼偏振传感器多方向标定方法 Download PDF

Info

Publication number
CN110887509A
CN110887509A CN201911251994.3A CN201911251994A CN110887509A CN 110887509 A CN110887509 A CN 110887509A CN 201911251994 A CN201911251994 A CN 201911251994A CN 110887509 A CN110887509 A CN 110887509A
Authority
CN
China
Prior art keywords
polarization
coordinate system
polarization sensor
compound eye
imitating
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.)
Granted
Application number
CN201911251994.3A
Other languages
English (en)
Other versions
CN110887509B (zh
Inventor
郭雷
刘鑫
杨健
王善澎
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beihang University
Original Assignee
Beihang University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Beihang University filed Critical Beihang University
Priority to CN201911251994.3A priority Critical patent/CN110887509B/zh
Publication of CN110887509A publication Critical patent/CN110887509A/zh
Application granted granted Critical
Publication of CN110887509B publication Critical patent/CN110887509B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C25/00Manufacturing, calibrating, cleaning, or repairing instruments or devices referred to in the other groups of this subclass

Landscapes

  • Engineering & Computer Science (AREA)
  • Manufacturing & Machinery (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Navigation (AREA)

Abstract

本发明提出一种仿复眼偏振传感器多方向标定方法,包括步骤:获取不同方向观测点的偏振角信息;利用载体所在的地理位置信息以及时间信息计算出地理坐标系下单位太阳矢量;结合高精度导航模块提供的姿态转换矩阵,天顶偏振传感器提供的偏振矢量,由观测点确定观测矢量,建立仿复眼偏振传感器坐标系与载体系的安装误差量测方程,并估计出两者之间的安装误差;利用加速度计获取重力信息,计算仿复眼偏振传感器坐标系下的太阳矢量信息;利用地理系下的太阳矢量与重力矢量的夹角作为量测,建立不同方向偏振传感器模块坐标系与仿复眼偏振传感器坐标系的安装误差量测方程,估计出不同方向偏振传感器的模块坐标系与仿复眼偏振传感器坐标系的安装误差。

Description

一种仿复眼偏振传感器多方向标定方法
技术领域
本发明涉及一种仿复眼偏振传感器多方向标定方法,用于估计仿复眼偏振传感器与高精度导航模块以及仿复眼偏振传感器中不同方向偏振传感器与天顶偏振传感器的位置安装误差,提高传感器的定姿、定位精度。
背景技术
自然界的生物具备高超的感知和获取空间运动信息的能力,经过研究发现生物可以通过感知天空偏振光分布,获取空间运动信息。因此,研究生物的导航机制,提出了偏振光导航的方法,由于偏振光导航具有自主性、抗干扰性以及环境适应性等优点,近几年受到国内外广泛的关注,而偏振传感器作为感知仪器,是信息获取的源头,导航的眼睛。因此,偏振传感器的研制成为仿生偏振导航系统的前端核心部分。
近几年,生物复眼的导航原理及结构成为研究热点,论文“多方向偏振光实时定位样机的设计与搭建”设计了5个方向的仿复眼偏振传感器,但是并未对仿复眼偏振传感器中不同方向偏振传感器间的位置关系进行标定和补偿。专利“一种基于扩展卡尔曼滤波的仿生偏振传感器多源误差标定方法”(申请号:201810129372.2),“一种基于无迹卡尔曼滤波的仿生偏振传感器多源误差标定方法”(申请号:201810129371.8),分析了单个传感器的误差标定及补偿方法,同样未涉及仿复眼偏振传感器中不同角度偏振传感器间的位置标定和补偿。随着对导航精度要求的提高,对仿复眼偏振传感器标定的准确性、全面性要求也提高,特别是对仿复眼偏振传感器中不同角度偏振传感器间位置关系的标定和补偿的研究还有待研究。
发明内容
为了解决上述问题,克服现有技术的不足,本发明提出一种仿复眼偏振传感器多方向标定方法,结合重力信息、天空偏振信息、太阳矢量信息、观测矢量信息、高精度导航模块,所述高精度是指姿态精度0.1°,航向精度0.1°,估计仿复眼偏振传感器与高精度导航模块以及仿复眼偏振传感器中不同方向偏振传感器与天顶偏振传感器的位置安装误差,提高传感器的定姿、定位精度。
本发明的技术解决方案为:一种仿复眼偏振传感器多方向标定方法,实现步骤如下:
步骤(1)、载体上安装有仿复眼偏振传感器,通过天文年历获取地理坐标系下的单位太阳矢量Sn,利用仿复眼偏振传感器的天顶偏振传感器测量天空相应点的单位偏振矢量
Figure BDA0002309302980000021
结合地理系下的太阳矢量信息,仿复眼偏振传感器坐标系下的观测矢量信息、偏振矢量信息,高精度导航模块提供的姿态转换矩阵信息,建立仿复眼偏振传感器坐标系与载体系,即高精度导航模块所在的坐标系的安装误差状态模型与量测模型,安装误差表示为φ1
Figure BDA0002309302980000022
为三维安装误差角,上标为数字1;
步骤(2)估计仿复眼偏振传感器坐标系与载体系的位置安装误差φ1,最后得到两者间的位置安装矩阵;
步骤(3)利用仿复眼偏振传感器中不同方向偏振传感器测量的偏振矢量
Figure BDA0002309302980000023
其中i=2~9,并结合天顶偏振传感器测量的偏振矢量,计算仿复眼偏振传感器坐标系下的太阳矢量;
步骤(4)根据载体系下的重力矢量,结合步骤(3)得到的仿复眼偏振传感器坐标系下的太阳矢量以及步骤(2)估计出的位置安装矩阵,并利用地理系与载体系两个不同坐标系下,太阳矢量与重力矢量夹角不变的原理,建立不同方向偏振传感器模块坐标系,即每个单独偏振传感器的坐标系与仿复眼偏振传感器坐标系的安装误差状态模型与量测模型,估计出不同方向偏振传感器的模块坐标系与仿复眼偏振传感器坐标系的位置安装误差φi
Figure BDA0002309302980000024
表示不同方向偏振传感器的模块坐标系与仿复眼偏振传感器坐标系的三维安装误差角,i=2~9;
步骤(5),经过步骤(1)~(4),完成仿复眼偏振传感器坐标系与载体系安装误差角以及不同方向偏振传感器的模块坐标系与仿复眼偏振传感器坐标系安装误差角的估计,即完成对仿复眼偏振传感器安装误差的标定。
进一步的,所述步骤(1)中,载体上安装有仿复眼偏振传感器,通过天文年历获取地理坐标系下的单位太阳矢量Sn,利用仿复眼偏振传感器的天顶偏振传感器测量天空相应点的单位偏振矢量
Figure BDA0002309302980000025
结合地理系下的太阳矢量信息,仿复眼偏振传感器坐标系下的观测矢量信息、偏振矢量信息,高精度导航模块提供的姿态转换矩阵信息,建立仿复眼偏振传感器坐标系与载体系的安装误差状态模型与量测模型,所述高精度导航模块的姿态精度0.1°,航向精度0.1°,安装误差表示为φ1
Figure BDA0002309302980000026
为三维安装误差角;仿复眼偏振传感器坐标系为m1系,天顶偏振传感器测得的偏振方位角为
Figure BDA0002309302980000027
得到仿复眼偏振传感器坐标系m1系下的单位偏振矢量
Figure BDA0002309302980000028
为:
Figure BDA0002309302980000031
Figure BDA0002309302980000032
中的下标为数字1,
Figure BDA0002309302980000033
与m1的下标均为数字1;
选取仿复眼偏振传感器坐标系与载体坐标系的安装误差φ1为状态量,建立状态模型:
Figure BDA0002309302980000034
量测模型:由于在同一坐标系下,偏振矢量垂直于太阳矢量与观测矢量,故:
Figure BDA0002309302980000035
其中,
Figure BDA0002309302980000036
Figure BDA0002309302980000037
为仿复眼偏振传感器坐标系与载体系之间已知的安装角度,
Figure BDA0002309302980000038
为已知的安装矩阵,
Figure BDA0002309302980000039
为待估计的安装误差角,M为待估计的安装误差矩阵,符号右上方的T表示转置;k1=±1,由偏振传感器测得的初始偏振矢量确定,
Figure BDA00023093029800000310
为观测矢量,
Figure BDA00023093029800000311
由高精度导航模块提供,sn表示地理系下的太阳矢量,V1表示天顶偏振传感器测量噪声;根据上述模型便估计出仿复眼偏振传感器坐标系与载体系的安装误差
Figure BDA00023093029800000312
进一步的,所述步骤(2),估计仿复眼偏振传感器坐标系与载体系的位置安装误差φ1,最后得到两者间的位置安装矩阵:
Figure BDA00023093029800000313
Figure BDA00023093029800000314
为已知的安装矩阵,M为待估计的安装误差矩阵。
进一步的,所述步骤(3),利用仿复眼偏振传感器中不同方向偏振传感器测量的偏振矢量
Figure BDA00023093029800000315
并结合天顶偏振传感器测量的偏振矢量,计算仿复眼偏振传感器坐标系下的太阳矢量;不同安装方向的偏振传感器中,每个偏振传感器固连于1个传感器模块坐标系即m系,天顶偏振传感器坐标系,即仿复眼偏振传感器坐标系为m1系,除天顶外,其他方向偏振传感器坐标系为mi系,每个传感器模块坐标系mi系(i=2~9)到仿复眼偏振传感器坐标系m1系的位置矩阵为
Figure BDA0002309302980000041
表示为:
Figure BDA0002309302980000042
其中,
Figure BDA0002309302980000043
为不同方向偏振传感器与天顶偏振传感器间已知的安装角度,
Figure BDA0002309302980000044
表示不同方向偏振传感器与天顶偏振传感器间等效旋转轴的方位角与高度角,
Figure BDA0002309302980000045
表示绕等效旋转轴旋转的角度;
Figure BDA0002309302980000046
为待估计的安装误差,
Figure BDA0002309302980000047
代表等效旋转轴,
Figure BDA0002309302980000048
表示等效旋转轴的第一个分量,
Figure BDA0002309302980000049
表示等效旋转轴的第二个分量,
Figure BDA00023093029800000410
表示等效旋转轴的第三个分量。
第i个偏振传感器测量得到偏振方位角
Figure BDA00023093029800000418
得到传感器模块坐标系mi下的单位偏振矢量
Figure BDA00023093029800000411
为:
Figure BDA00023093029800000412
进一步得到m1坐标系下第i个单位偏振矢量
Figure BDA00023093029800000413
为:
Figure BDA00023093029800000414
根据同一坐标系下太阳矢量与偏振矢量垂直的关系,得到m1系下的太阳矢量:
Figure BDA00023093029800000415
进一步的,所述步骤(4)中,根据载体系下的重力矢量,结合步骤(3)得到的仿复眼偏振传感器坐标系下的太阳矢量以及步骤(2)估计出的位置安装矩阵,并利用地理系与载体系两个不同坐标系下,太阳矢量与重力矢量夹角不变的原理,建立不同方向偏振传感器模块坐标系,即每个单独偏振传感器的坐标系与仿复眼偏振传感器坐标系的安装误差状态模型与量测模型,估计出不同方向偏振传感器的模块坐标系与仿复眼偏振传感器坐标系的位置安装误差φi
Figure BDA00023093029800000416
表示不同方向偏振传感器的模块坐标系与仿复眼偏振传感器坐标系的三维安装误差角,i=2~9,其中,状态方程,选取状态量为:
Figure BDA00023093029800000417
状态模型:
Figure BDA0002309302980000051
量测模型:由于太阳矢量与重力矢量在地理系与载体系下夹角相同,因此,选择地理系下太阳矢量与重力矢量的夹角作为量测量,建立量测模型:
Figure BDA0002309302980000052
其中,gn为地理系下的重力矢量,符号右上方的T表示转置,gb为载体系下的重力矢量,由加速度计测得,ki=±1,由太阳位置确定,Vi表示量测噪声;根据上述量测模型便估计出不同方向模块坐标系与仿复眼偏振传感器坐标系的安装误差
Figure BDA0002309302980000053
亦能够得到不同方向模块坐标系与仿复眼偏振传感器坐标系的位置安装矩阵
Figure BDA0002309302980000054
本发明与现有技术相比的优点在于:
首次提出了仿复眼偏振传感器与高精度导航模块(姿态精度0.1°,航向精度0.1°)以及仿复眼偏振传感器中不同方向偏振传感器与天顶偏振传感器的安装误差估计模型与方法,可进一步提高仿复眼偏振传感器的定向、定位精度。
附图说明
图1为本发明的方法流程图;
图2为仿复眼偏振传感器结构示意图;图中(1)为天顶偏振传感器,所在的坐标系为仿复眼偏振传感器坐标系;(2)为第二层不同方向的偏振传感器;(3)为第三层不同方向的偏振传感器;(4)为高精度导航模块(姿态精度0.1°,航向精度0.1°),所在的坐标系为载体系。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整的描述,显然,所描述的实施例仅为本发明的一部分实施例,而不是全部的实施例,基于本发明中的实施例,本领域的普通技术人员在不付出创造性劳动的前提下所获得的所有其他实施例,都属于本发明的保护范围。
本发明的载体可以是无人车、无人机、海洋舰船等,本发明的仿复眼偏振传感器是一种模仿生物复眼曲面多目结构的导航设备,该设备可用于解决无人车、无人机、深远海舰船等在卫星干扰、拒止以及地磁干扰环境下的自主导航、定位难题,可应用于电网巡线、环境监测、智能导航等多种领域。本发明提出的一种仿复眼偏振传感器多方向标定方法,可以补偿由于偏振传感器中偏振片制造工艺、安装以及不同方向传感器模块之间的安装问题引起的安装误差,这将大大提高仿复眼偏振传感器在工程应用中的定姿、定位精度。
根据本发明的一个实施例,仿复眼偏振传感器中包括高精度导航模块,安装于最下端;9个偏振传感器模块,分布于不同层的不同方向。其中每个偏振传感器为mi,其所在的坐标系为mi系,天顶偏振传感器为m1,分布于半球的第一层,m2,m4,m6,m8均匀分布于半球的第二层,与m1夹角为30°,m3,m5,m7,m9均匀分布于半球的第三层,与m1夹角为45°,高精度导航模块与m1平行安装,在不考虑高精度导航模块与m1安装误差的情况下,两者坐标系重合。
如图1所示,本发明的具体实现步骤如下:
步骤1、通过天文年历获取地理坐标系下的单位太阳矢量Sn,利用仿复眼偏振传感器的天顶偏振传感器测量天空相应点的单位偏振矢量
Figure BDA00023093029800000611
结合地理系下的太阳矢量信息,仿复眼偏振传感器坐标系下的观测矢量信息、偏振矢量信息,高精度导航模块(姿态精度0.1°,航向精度0.1°)提供的姿态转换矩阵信息,建立仿复眼偏振传感器坐标系与载体系的安装误差状态模型与量测模型,安装误差表示为φ1
Figure BDA0002309302980000061
为三维安装误差角。仿复眼偏振传感器坐标系为m1,天顶偏振传感器测得的偏振方位角为
Figure BDA0002309302980000062
得到仿复眼偏振传感器坐标系m1下的单位偏振矢量
Figure BDA0002309302980000063
为:
Figure BDA0002309302980000064
选取状态量为:
Figure BDA0002309302980000065
建立状态模型:
Figure BDA0002309302980000066
量测模型:由于在同一坐标系下,偏振矢量垂直于太阳矢量与观测矢量,可得:
Figure BDA0002309302980000067
其中,
Figure BDA0002309302980000068
Figure BDA0002309302980000069
为仿复眼偏振传感器坐标系与载体系之间已知的安装角度,
Figure BDA00023093029800000610
为已知的安装矩阵,
Figure BDA0002309302980000071
为待估计的安装误差角,M为待估计的安装误差矩阵,符号右上方的T表示转置。k1=±1,可以由偏振传感器测得的初始偏振矢量确定,
Figure BDA0002309302980000072
为观测矢量,
Figure BDA0002309302980000073
由高精度导航模块提供,sn表示地理系下的太阳矢量,V1表示天顶偏振传感器测量噪声。根据上述模型便可估计出仿复眼偏振传感器坐标系与载体系的安装误差
Figure BDA0002309302980000074
步骤2、估计仿复眼偏振传感器坐标系与载体系的位置安装误差φ1,最后得到两者间的位置安装矩阵:
Figure BDA0002309302980000075
步骤3、利用仿复眼偏振传感器中不同方向偏振传感器测量的偏振矢量
Figure BDA0002309302980000076
并结合天顶偏振传感器测量的偏振矢量,计算仿复眼偏振传感器坐标系下的太阳矢量。不同安装方向的偏振传感器中,每个偏振传感器固连于1个传感器模块坐标系即m系,天顶偏振传感器坐标系(即仿复眼偏振传感器坐标系)为m1系,除天顶外,其他方向偏振传感器其坐标系为mi系,每个传感器模块坐标系mi系(i=2~9)到仿复眼偏振传感器坐标系(m1系)的位置矩阵为
Figure BDA0002309302980000077
可表示为:
Figure BDA0002309302980000078
其中,
Figure BDA0002309302980000079
为不同方向偏振传感器与天顶偏振传感器间已知的安装角度,
Figure BDA00023093029800000710
表示不同方向偏振传感器与天顶偏振传感器间等效旋转轴的方位角与高度角,
Figure BDA00023093029800000711
表示绕等效旋转轴旋转的角度。
Figure BDA00023093029800000712
为待估计的安装误差,
Figure BDA00023093029800000713
代表等效旋转轴,
Figure BDA00023093029800000714
表示等效旋转轴的第一个分量,
Figure BDA00023093029800000715
表示等效旋转轴的第二个分量,
Figure BDA00023093029800000716
表示等效旋转轴的第三个分量。
第i个偏振传感器测量得到偏振方位角
Figure BDA00023093029800000717
得到传感器模块坐标系mi系下的单位偏振矢量
Figure BDA00023093029800000718
为:
Figure BDA00023093029800000719
进一步得到m1坐标系下第i个单位偏振矢量
Figure BDA00023093029800000720
为:
Figure BDA0002309302980000081
根据同一坐标系下太阳矢量与偏振矢量垂直的关系,得到m1系下的太阳矢量:
Figure BDA0002309302980000082
步骤4、根据载体系下的重力矢量,结合步骤3得到的仿复眼偏振传感器坐标系下的太阳矢量以及步骤2估计出的位置安装矩阵,并利用地理系与载体系两个不同坐标系下,太阳矢量与重力矢量夹角不变的原理,建立不同方向偏振传感器模块坐标系与仿复眼偏振传感器坐标系的安装误差状态模型与量测模型,估计出不同方向偏振传感器的模块坐标系,即每个单独偏振传感器的坐标系与仿复眼偏振传感器坐标系的位置安装误差φi
Figure BDA0002309302980000083
Figure BDA0002309302980000084
表示不同方向偏振传感器的模块坐标系与仿复眼偏振传感器坐标系的三维安装误差角,i=2~9,其中,状态方程,选取状态量为:
Figure BDA0002309302980000085
状态模型:
Figure BDA0002309302980000086
量测模型:由于太阳矢量与重力矢量在地理系与载体系下夹角相同。因此,选择地理系下太阳矢量与重力矢量的夹角作为量测量,建立量测模型:
Figure BDA0002309302980000087
其中,gn为地理系下的重力矢量,符号右上方的T表示转置,gb为载体系下的重力矢量,可由加速度计测得,ki=±1,可由太阳位置确定,Vi表示量测噪声。根据上述模型便可估计出不同方向模块坐标系与仿复眼偏振传感器坐标系的安装误差
Figure BDA0002309302980000088
亦可得到不同方向模块坐标系与仿复眼偏振传感器坐标系的位置安装矩阵
Figure BDA0002309302980000089
步骤(5)经过(1)~(4)步,完成仿复眼偏振传感器坐标系与载体系安装误差角以及不同方向偏振传感器的模块坐标系与仿复眼偏振传感器坐标系安装误差角的估计,即完成对仿复眼偏振传感器安装误差的标定。
尽管上面对本发明说明性的具体实施方式进行了描述,以便于本技术领域的技术人员理解本发明,且应该清楚,本发明不限于具体实施方式的范围,对本技术领域的普通技术人员来讲,只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内,这些变化是显而易见的,一切利用本发明构思的发明创造均在保护之列。

Claims (5)

1.一种仿复眼偏振传感器多方向标定方法,其特征在于,实现步骤如下:
步骤(1)、载体上安装有仿复眼偏振传感器,通过天文年历获取地理坐标系下的单位太阳矢量Sn,利用仿复眼偏振传感器的天顶偏振传感器测量天空相应点的单位偏振矢量
Figure FDA0002309302970000011
结合地理系下的太阳矢量信息,仿复眼偏振传感器坐标系下的观测矢量信息、偏振矢量信息,高精度导航模块提供的姿态转换矩阵信息,建立仿复眼偏振传感器坐标系与载体系,即高精度导航模块所在的坐标系的安装误差状态模型与量测模型,安装误差表示为φ1
Figure FDA0002309302970000012
Figure FDA0002309302970000017
为三维安装误差角,上标为数字1;
步骤(2)估计仿复眼偏振传感器坐标系与载体系的位置安装误差φ1,最后得到两者间的位置安装矩阵;
步骤(3)利用仿复眼偏振传感器中不同方向偏振传感器测量的偏振矢量
Figure FDA0002309302970000013
其中i=2~9,并结合天顶偏振传感器测量的偏振矢量,计算仿复眼偏振传感器坐标系下的太阳矢量;
步骤(4)根据载体系下的重力矢量,结合步骤(3)得到的仿复眼偏振传感器坐标系下的太阳矢量以及步骤(2)估计出的位置安装矩阵,并利用地理系与载体系两个不同坐标系下,太阳矢量与重力矢量夹角不变的原理,建立不同方向偏振传感器模块坐标系,即每个单独偏振传感器的坐标系与仿复眼偏振传感器坐标系的安装误差状态模型与量测模型,估计出不同方向偏振传感器的模块坐标系与仿复眼偏振传感器坐标系的位置安装误差φi
Figure FDA0002309302970000014
Figure FDA0002309302970000016
表示不同方向偏振传感器的模块坐标系与仿复眼偏振传感器坐标系的三维安装误差角,i=2~9;
步骤(5)、经过步骤(1)~(4),完成仿复眼偏振传感器坐标系与载体系安装误差角以及不同方向偏振传感器的模块坐标系与仿复眼偏振传感器坐标系安装误差角的估计,即完成对仿复眼偏振传感器安装误差的标定。
2.根据权利要求1所述的一种仿复眼偏振传感器多方向标定方法,其特征在于:
所述步骤(1)中,载体上安装有仿复眼偏振传感器,通过天文年历获取地理坐标系下的单位太阳矢量Sn,利用仿复眼偏振传感器的天顶偏振传感器测量天空相应点的单位偏振矢量
Figure FDA0002309302970000015
结合地理系下的太阳矢量信息,仿复眼偏振传感器坐标系下的观测矢量信息、偏振矢量信息,高精度导航模块提供的姿态转换矩阵信息,建立仿复眼偏振传感器坐标系与载体系的安装误差状态模型与量测模型,所述高精度导航模块的姿态精度0.1°,航向精度0.1°,安装误差表示为φ1
Figure FDA0002309302970000021
Figure FDA00023093029700000219
为三维安装误差角;仿复眼偏振传感器坐标系为m1系,天顶偏振传感器测得的偏振方位角为
Figure FDA00023093029700000218
得到仿复眼偏振传感器坐标系m1系下的单位偏振矢量
Figure FDA0002309302970000022
为:
Figure FDA0002309302970000023
Figure FDA0002309302970000024
中的下标为数字1,
Figure FDA0002309302970000025
与m1的下标均为数字1;
选取仿复眼偏振传感器坐标系与载体坐标系的安装误差φ1为状态量,建立状态模型:
Figure FDA0002309302970000026
量测模型:由于在同一坐标系下,偏振矢量垂直于太阳矢量与观测矢量,故:
Figure FDA0002309302970000027
其中,
Figure FDA0002309302970000028
Figure FDA0002309302970000029
Figure FDA00023093029700000210
为仿复眼偏振传感器坐标系与载体系之间已知的安装角度,
Figure FDA00023093029700000211
为已知的安装矩阵,
Figure FDA00023093029700000212
为待估计的安装误差角,M为待估计的安装误差矩阵,符号右上方的T表示转置;k1=±1,由偏振传感器测得的初始偏振矢量确定,
Figure FDA00023093029700000213
为观测矢量,
Figure FDA00023093029700000214
由高精度导航模块提供,sn表示地理系下的太阳矢量,V1表示天顶偏振传感器测量噪声;根据上述模型便估计出仿复眼偏振传感器坐标系与载体系的安装误差
Figure FDA00023093029700000215
3.根据权利要求1所述的一种仿复眼偏振传感器多方向标定方法,其特征在于:
所述步骤(2),估计仿复眼偏振传感器坐标系与载体系的位置安装误差φ1,最后得到两者间的位置安装矩阵:
Figure FDA00023093029700000216
Figure FDA00023093029700000217
为已知的安装矩阵,M为待估计的安装误差矩阵。
4.根据权利要求1所述的一种仿复眼偏振传感器多方向标定方法,其特征在于:
所述步骤(3),利用仿复眼偏振传感器中不同方向偏振传感器测量的偏振矢量
Figure FDA0002309302970000031
i=2~9,并结合天顶偏振传感器测量的偏振矢量,计算仿复眼偏振传感器坐标系下的太阳矢量;不同安装方向的偏振传感器中,每个偏振传感器固连于1个传感器模块坐标系即m系,天顶偏振传感器坐标系,即仿复眼偏振传感器坐标系为m1系,除天顶外,其他方向偏振传感器坐标系为mi系,每个传感器模块坐标系mi系(i=2~9)到仿复眼偏振传感器坐标系m1系的位置矩阵为
Figure FDA0002309302970000032
Figure FDA0002309302970000033
表示为:
Figure FDA0002309302970000034
Figure FDA0002309302970000035
Figure FDA0002309302970000036
Figure FDA0002309302970000037
其中,
Figure FDA0002309302970000038
为不同方向偏振传感器与天顶偏振传感器间已知的安装角度,
Figure FDA0002309302970000039
表示不同方向偏振传感器与天顶偏振传感器间等效旋转轴的方位角与高度角,
Figure FDA00023093029700000310
表示绕等效旋转轴旋转的角度;
Figure FDA00023093029700000311
为待估计的安装误差,
Figure FDA00023093029700000312
代表等效旋转轴,
Figure FDA00023093029700000313
表示等效旋转轴的第一个分量,
Figure FDA00023093029700000314
表示等效旋转轴的第二个分量,
Figure FDA00023093029700000315
表示等效旋转轴的第三个分量;
第i个偏振传感器测量得到偏振方位角
Figure FDA00023093029700000316
得到传感器模块坐标系mi下的单位偏振矢量
Figure FDA00023093029700000317
为:
Figure FDA00023093029700000318
进一步得到m1坐标系下第i个单位偏振矢量
Figure FDA00023093029700000319
为:
Figure FDA00023093029700000320
根据同一坐标系下太阳矢量与偏振矢量垂直的关系,得到m1系下的太阳矢量:
Figure FDA00023093029700000321
5.根据权利要求1所述的一种仿复眼偏振传感器多方向标定方法,其特征在于:
所述步骤(4)中,根据载体系下的重力矢量,结合步骤(3)得到的仿复眼偏振传感器坐标系下的太阳矢量以及步骤(2)估计出的位置安装矩阵,并利用地理系与载体系两个不同坐标系下,太阳矢量与重力矢量夹角不变的原理,建立不同方向偏振传感器模块坐标系,即每个单独偏振传感器的坐标系与仿复眼偏振传感器坐标系的安装误差状态模型与量测模型,估计出不同方向偏振传感器的模块坐标系与仿复眼偏振传感器坐标系的位置安装误差φi
Figure FDA0002309302970000041
Figure FDA0002309302970000047
表示不同方向偏振传感器的模块坐标系与仿复眼偏振传感器坐标系的三维安装误差角,i=2~9,其中,状态方程,选取状态量为:
Figure FDA0002309302970000042
状态模型:
Figure FDA0002309302970000043
量测模型:由于太阳矢量与重力矢量在地理系与载体系下夹角相同,因此,选择地理系下太阳矢量与重力矢量的夹角作为量测量,建立量测模型:
Figure FDA0002309302970000044
其中,gn为地理系下的重力矢量,符号右上方的T表示转置,gb为载体系下的重力矢量,由加速度计测得,ki=±1,由太阳位置确定,Vi表示量测噪声;根据上述量测模型便估计出不同方向模块坐标系与仿复眼偏振传感器坐标系的安装误差
Figure FDA0002309302970000045
亦能够得到不同方向模块坐标系与仿复眼偏振传感器坐标系的位置安装矩阵
Figure FDA0002309302970000046
CN201911251994.3A 2019-12-09 2019-12-09 一种仿复眼偏振传感器多方向标定方法 Active CN110887509B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911251994.3A CN110887509B (zh) 2019-12-09 2019-12-09 一种仿复眼偏振传感器多方向标定方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911251994.3A CN110887509B (zh) 2019-12-09 2019-12-09 一种仿复眼偏振传感器多方向标定方法

Publications (2)

Publication Number Publication Date
CN110887509A true CN110887509A (zh) 2020-03-17
CN110887509B CN110887509B (zh) 2021-09-07

Family

ID=69751154

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911251994.3A Active CN110887509B (zh) 2019-12-09 2019-12-09 一种仿复眼偏振传感器多方向标定方法

Country Status (1)

Country Link
CN (1) CN110887509B (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112782201A (zh) * 2021-01-26 2021-05-11 中国科学院上海技术物理研究所 一种基于x射线聚焦图像的龙虾眼光学器件标定方法
CN114018257A (zh) * 2021-11-05 2022-02-08 北京航空航天大学杭州创新研究院 一种偏振/惯性安装误差无依托自标定方法
CN117308926A (zh) * 2023-11-30 2023-12-29 北京航空航天大学 一种基于太阳敏感器和偏振传感器的太阳矢量寻优方法
CN117433564A (zh) * 2023-12-05 2024-01-23 北京航空航天大学 一种基于太阳敏感器辅助的偏振传感器误差在线标定方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102455522B (zh) * 2010-10-15 2014-11-12 香港科技大学 图案化偏振转换器
CN106441372A (zh) * 2016-11-24 2017-02-22 北京航空航天大学 一种基于偏振与重力信息的静基座粗对准方法
CN108225336A (zh) * 2018-01-24 2018-06-29 北京航空航天大学 一种基于置信度的偏振自主组合导航方法
CN108375381A (zh) * 2018-02-08 2018-08-07 北方工业大学 一种基于扩展卡尔曼滤波的仿生偏振传感器多源误差标定方法
CN108388720A (zh) * 2018-02-08 2018-08-10 北方工业大学 一种基于无迹卡尔曼滤波的仿生偏振传感器多源误差标定方法
CN108759820A (zh) * 2018-06-11 2018-11-06 北京航空航天大学 一种基于仿复眼多通道偏振传感器的太阳矢量解算方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102455522B (zh) * 2010-10-15 2014-11-12 香港科技大学 图案化偏振转换器
CN106441372A (zh) * 2016-11-24 2017-02-22 北京航空航天大学 一种基于偏振与重力信息的静基座粗对准方法
CN108225336A (zh) * 2018-01-24 2018-06-29 北京航空航天大学 一种基于置信度的偏振自主组合导航方法
CN108375381A (zh) * 2018-02-08 2018-08-07 北方工业大学 一种基于扩展卡尔曼滤波的仿生偏振传感器多源误差标定方法
CN108388720A (zh) * 2018-02-08 2018-08-10 北方工业大学 一种基于无迹卡尔曼滤波的仿生偏振传感器多源误差标定方法
CN108759820A (zh) * 2018-06-11 2018-11-06 北京航空航天大学 一种基于仿复眼多通道偏振传感器的太阳矢量解算方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
褚金奎等: "多方向偏振光实时定位样机的设计与搭建", 《光学精密工程》 *

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112782201A (zh) * 2021-01-26 2021-05-11 中国科学院上海技术物理研究所 一种基于x射线聚焦图像的龙虾眼光学器件标定方法
CN114018257A (zh) * 2021-11-05 2022-02-08 北京航空航天大学杭州创新研究院 一种偏振/惯性安装误差无依托自标定方法
CN114018257B (zh) * 2021-11-05 2024-05-24 北京航空航天大学杭州创新研究院 一种偏振/惯性安装误差无依托自标定方法
CN117308926A (zh) * 2023-11-30 2023-12-29 北京航空航天大学 一种基于太阳敏感器和偏振传感器的太阳矢量寻优方法
CN117308926B (zh) * 2023-11-30 2024-01-30 北京航空航天大学 一种基于太阳敏感器和偏振传感器的太阳矢量寻优方法
CN117433564A (zh) * 2023-12-05 2024-01-23 北京航空航天大学 一种基于太阳敏感器辅助的偏振传感器误差在线标定方法
CN117433564B (zh) * 2023-12-05 2024-03-01 北京航空航天大学 一种基于太阳敏感器辅助的偏振传感器误差在线标定方法

Also Published As

Publication number Publication date
CN110887509B (zh) 2021-09-07

Similar Documents

Publication Publication Date Title
CN110887509B (zh) 一种仿复眼偏振传感器多方向标定方法
CN103822629B (zh) 基于多方向偏振光导航传感器的定位系统及其定位方法
CN106767752B (zh) 一种基于偏振信息的组合导航方法
CN111323050B (zh) 一种捷联惯导和多普勒组合系统标定方法
CN109459027B (zh) 一种基于偏振-地磁矢量紧组合的导航方法
CN106643709B (zh) 一种海上运载体的组合导航方法及装置
CN106679645A (zh) 基于多方向偏振光的实时导航装置
CN110231025B (zh) 一种基于捷联式偏振光罗盘的动态定向方法及系统
CN108759819A (zh) 一种基于全天域偏振度信息的偏振导航实时定位方法
CN108759820B (zh) 一种基于仿复眼多通道偏振传感器的太阳矢量解算方法
CN110887473B (zh) 一种基于偏振度加权的仿生偏振自主组合导航方法
CN107727101B (zh) 基于双偏振光矢量的三维姿态信息快速解算方法
CN109556631A (zh) 一种基于最小二乘的ins/gnss/偏振/地磁组合导航系统对准方法
CN105973268B (zh) 一种基于共基座安装的传递对准精度定量评估方法
CN112066979A (zh) 一种基于偏振位姿信息耦合迭代自主导航定位方法
CN112923924B (zh) 一种锚泊船舶姿态与位置监测方法及系统
CN108225336A (zh) 一种基于置信度的偏振自主组合导航方法
CN103727937A (zh) 一种基于星敏感器的舰船姿态确定方法
CN110132308A (zh) 一种基于姿态确定的usbl安装误差角标定方法
CN112129288B (zh) 一种基于偏振光/地磁航向约束的位姿估计方法及系统
CN110887476B (zh) 基于偏振-天文夹角信息观测的自主航向与姿态确定方法
CN106767928A (zh) 一种自适应快速传递对准方法
CN116448145A (zh) 一种基于偏振矢量空间差分的航姿确定方法
CN104501809A (zh) 一种基于姿态耦合的捷联惯导/星敏感器组合导航方法
CN110398702B (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