CN109459059A - 一种星敏感器外场转换基准测定系统及方法 - Google Patents
一种星敏感器外场转换基准测定系统及方法 Download PDFInfo
- Publication number
- CN109459059A CN109459059A CN201811392442.XA CN201811392442A CN109459059A CN 109459059 A CN109459059 A CN 109459059A CN 201811392442 A CN201811392442 A CN 201811392442A CN 109459059 A CN109459059 A CN 109459059A
- Authority
- CN
- China
- Prior art keywords
- theodolite
- coordinate system
- star sensor
- star
- rib body
- 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
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C25/00—Manufacturing, calibrating, cleaning, or repairing instruments or devices referred to in the other groups of this subclass
Abstract
本发明涉及一种星敏感器外场转换基准的测定系统及方法,所述的系统包括星敏感器1、天文引北经纬仪2、经纬仪A3、经纬仪B4、GPS天线5、GPS接收机6、气压温度传感器7和控制计算机8;所述的方法,使用本发明所述的系统;本方法应用恒星作为测量参考,测定星敏感器测量坐标系与基准棱体坐标系的转换关系。避免了传统实验室测定方法中,在星敏感器标定过程中,主点等参数受参数耦合影响而产生的估计误差。有效地提升了转换基准的测定精度。本方法广泛适用于各类型星敏感器外部引出基准的测定。
Description
技术领域
本发明属于天文导航领域,具体涉及一种星敏感器外场转换基准测定系统及方法。
背景技术
关于星敏感器转换基准测定方法,尚未见有公开报道的专利。
星敏感器是具有高度自主性的导航设备,通过视场中星点的观测矢量与天文星表中的参考矢量进行比对,测得载体的姿态信息。在实际测量过程中,星敏感器直接得到的是其测量坐标系相对参考坐标系的关系。测量坐标系建立在光电探测部件上,没有明确的外化特征。
一般使用六面棱体将星敏感器的姿态测量结果传递给其它导航设备。各台星敏感器出厂前,棱体坐标系相对星敏感器测量坐标系的转换基准需要事先测定并写入姿态解算程序中。实际工作时,星敏感器直接输出棱体坐标系相对参考坐标系的姿态信息。转换基准对于星敏感器的精度性能至关重要,其测定误差将导致姿态测量结果引入固定偏差。
目前星敏感器转换基准一般采用实验室方法测定:星敏感器安装在高精度转台上对长焦平行光管的模拟单星进行观测,转台发生不同角度模拟星点从不同方向入射。通过参数标定求解,可得出星敏感器测量坐标系与转台坐标系间的关系。转台个转轴上贴有基准棱镜,利用经纬仪可测得转台坐标系与星敏感器棱体坐标系间的关系。结合两次测量结果可测得星敏感器转换基准。
实验室标定过程中存在参数耦合问题:标定采集质心数据中的误差将使得各个标定参数相互耦合,估计偏离真值。由于星敏感器光学镜头的畸变较小,因此成像模型接近于针孔模型,表征测量坐标系的主点参数受耦合的影响最为严重。根据仿真结果分析,主点位置误差可达到角分级别。该误差将直接引入转换基准的最终测定结果中,造成巨大的转换误差。
发明内容
本发明的目的在于:克服传统星敏感器转换基准测定方法受实验标定过程参数耦合问题影响在测定精度方面的不足,提出一种星敏感器外场转换基准测定方法。该方法利用天文引北的方式确定棱体坐标系与当地惯性坐标系的关系。再结合星敏感器姿态信息、导航卫星信息和气象数据确定星敏感器测量坐标系与当地惯性坐标系的关系,最终测定转换基准。本方法测定精度高,广泛适用于各类星敏感器。
本发明的技术方案如下:一种星敏感器外场转换基准的测定系统,包括星敏感器、天文引北经纬仪、经纬仪A、经纬仪B、GPS天线、GPS接收机、气压温度传感器和控制计算机;
所述的星敏感器、GPS接收机和气压温度传感器分别与控制计算机相连,用于将相关数据信息传输至控制计算机;
所述的GPS接收机与GPS天线相连,用于接收导航卫星传入的地理位置和时间信息;
所述的天文引北经纬仪依据北极星建立当地惯性坐标系,为经纬仪A、经纬仪B提供基准信息,从而确定棱体坐标系与当地惯性坐标系的相互关系;
所述的经纬仪A、经纬仪B分别对星敏感器棱体上相互正交的两个面的法向进行测定;
所述的星敏感器对视场内的恒星进行观测,并将各帧星图数据送入控制计算机中;所述的控制计算机在收到星图数据后,结合导航卫星传入的地理位置和时间信息以及气压温度传感器传入的气压温度信息,确定星敏感器在当地惯性坐标系下的姿态;
结合棱体坐标系与当地惯性坐标系的关系和星敏感器在当地惯性坐标系下的姿态,进而求得星敏感器的转换基准。
一种星敏感器外场转换基准的测定方法,其特征在于本发明所述的系统,包括以下步骤:
步骤一、将经纬仪A、经纬仪B经过精调平后,对星敏感器棱体上相互正交的两个面的法向进行测定;得到星敏感器棱体坐标系的两个轴向与各自瞄准经纬仪坐标系的相互关系;
步骤二、天文引北经纬仪经过精调平后,以北极星为方位基准,确定当地惯性坐标系;
步骤三、通过天文引北经纬仪分别与经纬仪A和经纬仪B的互瞄,将天文引北经纬仪测得的基准传给经纬仪A和经纬仪B,并得到棱体坐标系与当地惯性坐标系的关系;
步骤四、通过对星敏感器视场中的恒星进行星图识别,得到观测恒星对应天文星表中的参考矢量,将这些矢量建立在天球坐标系下;
步骤五、根据GPS接收机传入的地理经纬度信息和时间信息,测算出天球坐标系与当地惯性坐标系间的转换关系;
步骤六、利用步骤五得到的天球坐标系与当地惯性坐标系间的转换关系,将天文星表中的参考矢量转换至当地惯性坐标系下;
步骤七、由于大气对星光矢量存在折射影响,当地惯性坐标系下的参考矢量需要进行蒙气差补偿:根据传感器传入的温度信息、气压信息建立蒙气差模型,对参考矢量的天顶距进行修正;
步骤八、将星敏感器测量得的观测矢量与经过步骤七的蒙气差补偿后的参考矢量联立,结算出星敏感器在地面惯性坐标系下的姿态;
步骤九、结合步骤三求得的棱体坐标系与当地惯性坐标系的关系和步骤八求得的星敏感器在地面惯性坐标系下的姿态,求得星敏感器转换基准;为进一步消除大气路径对星敏感器姿态测量结果的影响,进行多次测量,运用QUEST算法进行优化求解。
进一步的,所述的星敏感器的测量坐标系为OSS-XSSYSSZSS,星敏感器直接测得姿态的参考系;其中,XSS轴:在主视场像面上,正方向自左向右;YSS轴:在主视场像面上,正方向自下向上;ZSS轴:与主视场光轴重合,指向天空;
所述的棱体坐标系OLT-XLTYLTZLT,以星敏感器基准棱体各面法向建立的参考系;经纬仪A坐标系OTR1-XTR1YTR1ZTR1和经纬仪B坐标系OTR2-XTR2YTR2ZTR2,作为测量棱体两个正交面的经纬仪建立的参考系,XTR1ZTR1平面和XTR2ZTR2平面为水平面;
所述的当地惯性参考系ONTE-XNTEYNTEZNTE:XNTE轴,在水平面上,指向站心北向;YNTE轴,站心垂向上,指向天顶;ZNTE轴,在水平面上,指向站心东向。
进一步的,步骤一中,经纬仪A3放置在能够瞄向棱体-XLT轴向并且能够与天文引北经纬仪互瞄的位置,且精确调平;通过前者,经纬仪ATR1在经纬度方向与棱体坐标系的-XLT轴向重合;
在天文引北经纬仪和经纬仪A均精确调平的情况下,经纬仪A坐标系TR1与北天东坐标系仅相差一个偏航角的转动;记互瞄时,天文引北设备的经度是θref1,中间经纬仪的经度是θtr1;则棱体XLT轴向在北天东坐标系下的经纬度表示为:
棱体XLT轴向在北天东坐标系下的方向矢量表示为:
经纬仪B放置在能够瞄向棱体-YLT轴向并且能够与天文引北设备经纬仪互瞄的位置,且精确调平;通过前者,经纬仪BTR2在经纬度方向与棱体坐标系的-YLT轴向重合;
进一步的,步骤二中,通过精确调平和北极星引北,天文引北经纬仪坐标系与北天东坐标系重合。
进一步的,步骤三中,对互瞄进行过程必要的说明,在天文引北经纬仪和经纬仪B均精确调平的情况下,经纬仪B坐标系TR2与北天东坐标系仅相差一个偏航角的转动。记互瞄时,天文引北设备的经度是θref2,中间经纬仪的经度是θtr2;则棱体-YLT轴向在北天东坐标系下的经纬度表示为:
棱体-YLT轴向在北天东坐标系下的方向矢量表示为:
计算棱体坐标系至当地惯性坐标系关系根据矢量叉积定义,得到棱体ZLT轴向在北天东坐标系下的方向矢量表示为:
根据双矢量定姿,求得棱体坐标系至当地惯性坐标系的旋转矩阵为:
则当地惯性坐标系下的北向XNTE、天向YNTE和东向ZNTE在棱体坐标系的矢量表达为:
进一步的,步骤五中,星图识别后,需利用天球坐标系至当地惯性坐标系的转换关系对导航星库中的参考星矢量进行参考系转换;该转换关系的求解过程先后包括:岁差补偿、章动补偿、自转补偿、极移补偿;其中,自转补偿需利用导航卫星提供的地理经纬度信息和授时信息,且授时信息需与星敏感器1的采样时刻严格同步。
进一步的,步骤七中,经过参考系转换后的恒星参考矢量需要进行蒙气差补偿;蒙气差ρ采用中国天文年历中蒙气差修正模型:
ρ=(1+αtAt+B)ρ0 (8)
式中,ρ0为蒙气差常数,可直接利用下式计算:
其中Z为天顶距,Z=90°-δ。αt为气温变差乘数修订系数,当仰角为90°时,αt=1,当δ≤45°时可由下式计算:
At为气温变差乘数,与温度t(℃)有关:
B为气压变差乘数,与测站附近气压P有关:
进一步的,步骤八中,得到当地惯性坐标系至星敏感器测量坐标系关系后,求出当地惯性坐标系下的北向XNTE、天向YNTE和东向ZNTE在星敏感器测量坐标系的矢量表达为:
进一步的,步骤九中,为进一步消除大气路径对最终测定精度的影响,可更换星敏感器摆放位置,重复进行前序的操作,得出多组式(7)和式(13)的求解结果,利用QUEST姿态求解算法得到转换基准
本发明的显著效果在于:本转换基准关系测定方法应用天文引北的方式确定了星敏感器棱体坐标系在当地惯性坐标系下的准确指向。利用导航卫星接收机得到的地理位置信息和时间信息将星敏感器姿态测量的参考系从天球坐标系转换至当地惯性系下。此外,利用传感器传来的气压信息和温度信息对星光矢量进行补偿,避免了大气蒙气差产生的星光折射对测量结果的影响。
综上,本方法应用恒星作为测量参考,测定星敏感器测量坐标系与基准棱体坐标系的转换关系。避免了传统实验室测定方法中,在星敏感器标定过程中,主点等参数受参数耦合影响而产生的估计误差。有效地提升了转换基准的测定精度。本方法广泛适用于各类型星敏感器外部引出基准的测定。
附图说明
图1为本发明所述的一种星敏感器外场转换基准测定系统的结构示意图;
图2为星敏感器外场基准转换基准测定流程示意图;
图中:1.星敏感器、2.天文引北经纬仪、3.经纬仪A、4.经纬仪B、5.GPS天线、6.GPS接收机、7.气压温度传感器、8.控制计算机。
具体实施方式
下面结合附图及具体实施例对本发明所述的一种星敏感器外场转换基准测定系统及方法作进一步详细说明。
如图1所示,一种星敏感器外场转换基准的测定系统,包括星敏感器1、天文引北经纬仪2、经纬仪A3、经纬仪B4、GPS天线5、GPS接收机6、气压温度传感器7和控制计算机8;
所述的星敏感器1、GPS接收机6和气压温度传感器7分别与控制计算机8相连,用于将相关数据信息传输至控制计算机8;
所述的GPS接收机6与GPS天线5相连,用于接收导航卫星传入的地理位置和时间信息;
所述的天文引北经纬仪2依据北极星建立当地惯性坐标系,为经纬仪A3、经纬仪B4提供基准信息,从而确定棱体坐标系与当地惯性坐标系的相互关系;
所述的经纬仪A3、经纬仪B4分别对星敏感器1棱体上相互正交的两个面的法向进行测定;
所述的星敏感器1对视场内的恒星进行观测,并将各帧星图数据送入控制计算机8中;所述的控制计算机8在收到星图数据后,结合导航卫星传入的地理位置和时间信息以及气压温度传感器7传入的气压温度信息,确定星敏感器1在当地惯性坐标系下的姿态;
结合棱体坐标系与当地惯性坐标系的关系和星敏感器1在当地惯性坐标系下的姿态,进而求得星敏感器1的转换基准。
如图2所示,一种星敏感器外场转换基准的测定方法,其特征在于使用本发明所述的系统,包括以下步骤:
步骤一、将经纬仪A3、经纬仪B4经过精调平后,对星敏感器1棱体上相互正交的两个面的法向进行测定;得到星敏感器1棱体坐标系的两个轴向与各自瞄准经纬仪坐标系的相互关系;
步骤二、天文引北经纬仪2经过精调平后,以北极星为方位基准,确定当地惯性坐标系;
步骤三、通过天文引北经纬仪2分别与经纬仪A3和经纬仪B4的互瞄,将天文引北经纬仪2测得的基准传给经纬仪A3和经纬仪B4,并得到棱体坐标系与当地惯性坐标系的关系;
步骤四、通过对星敏感器1视场中的恒星进行星图识别,得到观测恒星对应天文星表中的参考矢量,将这些矢量建立在天球坐标系下;
步骤五、根据GPS接收机6传入的地理经纬度信息和时间信息,测算出天球坐标系与当地惯性坐标系间的转换关系;
步骤六、利用步骤五得到的天球坐标系与当地惯性坐标系间的转换关系,将天文星表中的参考矢量转换至当地惯性坐标系下;
步骤七、由于大气对星光矢量存在折射影响,当地惯性坐标系下的参考矢量需要进行蒙气差补偿:根据传感器传入的温度信息、气压信息建立蒙气差模型,对参考矢量的天顶距进行修正;
步骤八、将星敏感器1测量得的观测矢量与经过步骤七的蒙气差补偿后的参考矢量联立,结算出星敏感器1在地面惯性坐标系下的姿态;
步骤九、结合步骤三求得的棱体坐标系与当地惯性坐标系的关系和步骤八求得的星敏感器1在地面惯性坐标系下的姿态,求得星敏感器1转换基准;为进一步消除大气路径对星敏感器1姿态测量结果的影响,进行多次测量,运用QUEST算法进行优化求解。
进一步的,所述的星敏感器1的测量坐标系为OSS-XSSYSSZSS,星敏感器1直接测得姿态的参考系;其中,XSS轴:在主视场像面上,正方向自左向右;YSS轴:在主视场像面上,正方向自下向上;ZSS轴:与主视场光轴重合,指向天空;
所述的棱体坐标系OLT-XLTYLTZLT,以星敏感器1基准棱体各面法向建立的参考系;经纬仪A3坐标系OTR1-XTR1YTR1ZTR1和经纬仪B4坐标系OTR2-XTR2YTR2ZTR2,作为测量棱体两个正交面的经纬仪建立的参考系,XTR1ZTR1平面和XTR2ZTR2平面为水平面;
所述的当地惯性参考系ONTE-XNTEYNTEZNTE:XNTE轴,在水平面上,指向站心北向;YNTE轴,站心垂向上,指向天顶;ZNTE轴,在水平面上,指向站心东向。
进一步的,步骤一中,经纬仪A3放置在能够瞄向棱体-XLT轴向并且能够与天文引北经纬仪2互瞄的位置,且精确调平;通过前者,经纬仪A3TR1在经纬度方向与棱体坐标系的-XLT轴向重合;
在天文引北经纬仪2和经纬仪A3均精确调平的情况下,经纬仪A3坐标系TR1与北天东坐标系仅相差一个偏航角的转动;记互瞄时,天文引北设备的经度是θref1,中间经纬仪的经度是θtr1;则棱体XLT轴向在北天东坐标系下的经纬度表示为:
棱体XLT轴向在北天东坐标系下的方向矢量表示为:
经纬仪B4放置在能够瞄向棱体-YLT轴向并且能够与天文引北设备经纬仪互瞄的位置,且精确调平;通过前者,经纬仪B4TR2在经纬度方向与棱体坐标系的-YLT轴向重合;
进一步的,步骤二中,通过精确调平和北极星引北,天文引北经纬仪2坐标系与北天东坐标系重合。
进一步的,步骤三中,对互瞄进行过程必要的说明,在天文引北经纬仪2和经纬仪B4均精确调平的情况下,经纬仪B4坐标系TR2与北天东坐标系仅相差一个偏航角的转动。记互瞄时,天文引北设备的经度是θref2,中间经纬仪的经度是θtr2;则棱体-YLT轴向在北天东坐标系下的经纬度表示为:
棱体-YLT轴向在北天东坐标系下的方向矢量表示为:
计算棱体坐标系至当地惯性坐标系关系根据矢量叉积定义,得到棱体ZLT轴向在北天东坐标系下的方向矢量表示为:
根据双矢量定姿,求得棱体坐标系至当地惯性坐标系的旋转矩阵为:
则当地惯性坐标系下的北向XNTE、天向YNTE和东向ZNTE在棱体坐标系的矢量表达为:
进一步的,步骤五中,星图识别后,需利用天球坐标系至当地惯性坐标系的转换关系对导航星库中的参考星矢量进行参考系转换;该转换关系的求解过程先后包括:岁差补偿、章动补偿、自转补偿、极移补偿;其中,自转补偿需利用导航卫星提供的地理经纬度信息和授时信息,且授时信息需与星敏感器1的采样时刻严格同步。
进一步的,步骤七中,经过参考系转换后的恒星参考矢量需要进行蒙气差补偿;蒙气差ρ采用中国天文年历中蒙气差修正模型:
ρ=(1+αtAt+B)ρ0 (8)
式中,ρ0为蒙气差常数,可直接利用下式计算:
其中Z为天顶距,Z=90°-δ。αt为气温变差乘数修订系数,当仰角为90°时,αt=1,当δ≤45°时可由下式计算:
At为气温变差乘数,与温度t(℃)有关:
B为气压变差乘数,与测站附近气压P有关:
进一步的,步骤八中,得到当地惯性坐标系至星敏感器1测量坐标系关系后,求出当地惯性坐标系下的北向XNTE、天向YNTE和东向ZNTE在星敏感器1测量坐标系的矢量表达为:
进一步的,步骤九中,为进一步消除大气路径对最终测定精度的影响,可更换星敏感器1摆放位置,重复进行前序的操作,得出多组式(7)和式(13)的求解结果,利用QUEST姿态求解算法得到转换基准
Claims (10)
1.一种星敏感器外场转换基准的测定系统,其特征在于:包括星敏感器(1)、天文引北经纬仪(2)、经纬仪A(3)、经纬仪B(4)、GPS天线(5)、GPS接收机(6)、气压温度传感器(7)和控制计算机(8);
所述的星敏感器(1)、GPS接收机(6)和气压温度传感器(7)分别与控制计算机(8)相连,用于将相关数据信息传输至控制计算机(8);
所述的GPS接收机(6)与GPS天线(5)相连,用于接收导航卫星传入的地理位置和时间信息;
所述的天文引北经纬仪(2)依据北极星建立当地惯性坐标系,为经纬仪A(3)、经纬仪B(4)提供基准信息,从而确定棱体坐标系与当地惯性坐标系的相互关系;
所述的经纬仪A(3)、经纬仪B(4)分别对星敏感器(1)棱体上相互正交的两个面的法向进行测定;
所述的星敏感器(1)对视场内的恒星进行观测,并将各帧星图数据送入控制计算机(8)中;所述的控制计算机(8)在收到星图数据后,结合导航卫星传入的地理位置和时间信息以及气压温度传感器(7)传入的气压温度信息,确定星敏感器(1)在当地惯性坐标系下的姿态;
结合棱体坐标系与当地惯性坐标系的关系和星敏感器(1)在当地惯性坐标系下的姿态,进而求得星敏感器(1)的转换基准。
2.一种星敏感器外场转换基准的测定方法,其特征在于使用权利要求1所述的系统,包括以下步骤:
步骤一、将经纬仪A(3)、经纬仪B(4)经过精调平后,对星敏感器(1)棱体上相互正交的两个面的法向进行测定;得到星敏感器(1)棱体坐标系的两个轴向与各自瞄准经纬仪坐标系的相互关系;
步骤二、天文引北经纬仪(2)经过精调平后,以北极星为方位基准,确定当地惯性坐标系;
步骤三、通过天文引北经纬仪(2)分别与经纬仪A(3)和经纬仪B(4)的互瞄,将天文引北经纬仪(2)测得的基准传给经纬仪A(3)和经纬仪B(4),并得到棱体坐标系与当地惯性坐标系的关系;
步骤四、通过对星敏感器(1)视场中的恒星进行星图识别,得到观测恒星对应天文星表中的参考矢量,将这些矢量建立在天球坐标系下;
步骤五、根据GPS接收机(6)传入的地理经纬度信息和时间信息,测算出天球坐标系与当地惯性坐标系间的转换关系;
步骤六、利用步骤五得到的天球坐标系与当地惯性坐标系间的转换关系,将天文星表中的参考矢量转换至当地惯性坐标系下;
步骤七、由于大气对星光矢量存在折射影响,当地惯性坐标系下的参考矢量需要进行蒙气差补偿:根据传感器传入的温度信息、气压信息建立蒙气差模型,对参考矢量的天顶距进行修正;
步骤八、将星敏感器(1)测量得的观测矢量与经过步骤七的蒙气差补偿后的参考矢量联立,结算出星敏感器(1)在地面惯性坐标系下的姿态;
步骤九、结合步骤三求得的棱体坐标系与当地惯性坐标系的关系和步骤八求得的星敏感器(1)在地面惯性坐标系下的姿态,求得星敏感器(1)转换基准;为进一步消除大气路径对星敏感器(1)姿态测量结果的影响,进行多次测量,运用QUEST算法进行优化求解。
3.如权利要求2所述的一种星敏感器外场转换基准的测定方法,其特征在于:所述的星敏感器(1)的测量坐标系为OSS-XSSYSSZSS,星敏感器(1)直接测得姿态的参考系;其中,XSS轴:在主视场像面上,正方向自左向右;YSS轴:在主视场像面上,正方向自下向上;ZSS轴:与主视场光轴重合,指向天空;
所述的棱体坐标系OLT-XLTYLTZLT,以星敏感器(1)基准棱体各面法向建立的参考系;经纬仪A(3)坐标系OTR1-XTR1YTR1ZTR1和经纬仪B(4)坐标系OTR2-XTR2YTR2ZTR2,作为测量棱体两个正交面的经纬仪建立的参考系,XTR1ZTR1平面和XTR2ZTR2平面为水平面;
所述的当地惯性参考系ONTE-XNTEYNTEZNTE:XNTE轴,在水平面上,指向站心北向;YNTE轴,站心垂向上,指向天顶;ZNTE轴,在水平面上,指向站心东向。
4.如权利要求3所述的一种星敏感器外场转换基准的测定方法,其特征在于:步骤一中,经纬仪A(3)放置在能够瞄向棱体-XLT轴向并且能够与天文引北经纬仪(2)互瞄的位置,且精确调平;通过前者,经纬仪A(3)TR1在经纬度方向与棱体坐标系的-XLT轴向重合;
在天文引北经纬仪(2)和经纬仪A(3)均精确调平的情况下,经纬仪A(3)坐标系TR1与北天东坐标系仅相差一个偏航角的转动;记互瞄时,天文引北设备的经度是θref1,中间经纬仪的经度是θtr1;则棱体XLT轴向在北天东坐标系下的经纬度表示为:
棱体XLT轴向在北天东坐标系下的方向矢量表示为:
经纬仪B(4)放置在能够瞄向棱体-YLT轴向并且能够与天文引北设备经纬仪互瞄的位置,且精确调平;通过前者,经纬仪B(4)TR2在经纬度方向与棱体坐标系的-YLT轴向重合。
5.如权利要求3所述的一种星敏感器外场转换基准的测定方法,其特征在于:步骤二中,通过精确调平和北极星引北,天文引北经纬仪(2)坐标系与北天东坐标系重合。
6.如权利要求3所述的一种星敏感器外场转换基准的测定方法,其特征在于:步骤三中,对互瞄进行过程必要的说明,在天文引北经纬仪(2)和经纬仪B(4)均精确调平的情况下,经纬仪B(4)坐标系TR2与北天东坐标系仅相差一个偏航角的转动。记互瞄时,天文引北设备的经度是θref2,中间经纬仪的经度是θtr2;则棱体-YLT轴向在北天东坐标系下的经纬度表示为:
棱体-YLT轴向在北天东坐标系下的方向矢量表示为:
计算棱体坐标系至当地惯性坐标系关系根据矢量叉积定义,得到棱体ZLT轴向在北天东坐标系下的方向矢量表示为:
根据双矢量定姿,求得棱体坐标系至当地惯性坐标系的旋转矩阵为:
则当地惯性坐标系下的北向XNTE、天向YNTE和东向ZNTE在棱体坐标系的矢量表达为:
7.如权利要求3所述的一种星敏感器外场转换基准的测定方法,其特征在于:步骤五中,星图识别后,需利用天球坐标系至当地惯性坐标系的转换关系对导航星库中的参考星矢量进行参考系转换;该转换关系的求解过程先后包括:岁差补偿、章动补偿、自转补偿、极移补偿;其中,自转补偿需利用导航卫星提供的地理经纬度信息和授时信息,且授时信息需与星敏感器(1)的采样时刻严格同步。
8.如权利要求3所述的一种星敏感器外场转换基准的测定方法,其特征在于:步骤七中,经过参考系转换后的恒星参考矢量需要进行蒙气差补偿;蒙气差ρ采用中国天文年历中蒙气差修正模型:
ρ=(1+αtAt+B)ρ0 (8)
式中,ρ0为蒙气差常数,可直接利用下式计算:
其中Z为天顶距,Z=90°-δ。αt为气温变差乘数修订系数,当仰角为90°时,αt=1,当δ≤45°时可由下式计算:
At为气温变差乘数,与温度t(℃)有关:
B为气压变差乘数,与测站附近气压P有关:
9.如权利要求3所述的一种星敏感器外场转换基准的测定方法,其特征在于:步骤八中,得到当地惯性坐标系至星敏感器(1)测量坐标系关系后,求出当地惯性坐标系下的北向XNTE、天向YNTE和东向ZNTE在星敏感器(1)测量坐标系的矢量表达为:
10.如权利要求3所述的一种星敏感器外场转换基准的测定方法,其特征在于:步骤九中,为进一步消除大气路径对最终测定精度的影响,可更换星敏感器(1)摆放位置,重复进行前序的操作,得出多组式(7)和式(13)的求解结果,利用QUEST姿态求解算法得到转换基准
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811392442.XA CN109459059B (zh) | 2018-11-21 | 2018-11-21 | 一种星敏感器外场转换基准测定系统及方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811392442.XA CN109459059B (zh) | 2018-11-21 | 2018-11-21 | 一种星敏感器外场转换基准测定系统及方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109459059A true CN109459059A (zh) | 2019-03-12 |
CN109459059B CN109459059B (zh) | 2022-08-19 |
Family
ID=65611108
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811392442.XA Active CN109459059B (zh) | 2018-11-21 | 2018-11-21 | 一种星敏感器外场转换基准测定系统及方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109459059B (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112212856A (zh) * | 2020-09-21 | 2021-01-12 | 中国科学院长春光学精密机械与物理研究所 | 一种近地面白天星敏感器的高精度姿态测量方法 |
CN112833878A (zh) * | 2021-01-05 | 2021-05-25 | 上海航天控制技术研究所 | 近地面多源天文自主导航方法 |
CN113607188A (zh) * | 2021-08-02 | 2021-11-05 | 北京航空航天大学 | 基于经纬仪叉丝成像的多视场星敏感器的标定系统及方法 |
CN113916216A (zh) * | 2020-09-09 | 2022-01-11 | 蓝箭航天空间科技股份有限公司 | 惯性测量信息与星光测量信息数据同步方法及系统 |
CN114001756A (zh) * | 2021-11-03 | 2022-02-01 | 中国科学院光电技术研究所 | 一种小视场星敏感器外场地面寻星方法 |
CN114623833A (zh) * | 2022-01-28 | 2022-06-14 | 北京控制工程研究所 | 一种空间指向测量仪器视场标定系统和方法 |
CN114740801A (zh) * | 2022-03-21 | 2022-07-12 | 成都飞机工业(集团)有限责任公司 | 一种用于数控设备群协同生产线安装的基坐标系创建方法 |
Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS5973719A (ja) * | 1982-10-20 | 1984-04-26 | Natl Space Dev Agency Japan<Nasda> | 人工衛星搭載恒星センサ信号処理方式 |
CN102279002A (zh) * | 2011-06-27 | 2011-12-14 | 哈尔滨工业大学 | 星敏感器测量坐标系与载体坐标系转换矩阵的标定方法 |
CN102288199A (zh) * | 2011-06-22 | 2011-12-21 | 哈尔滨工业大学 | 一种星敏感器的地面测试方法 |
CN102426017A (zh) * | 2011-11-03 | 2012-04-25 | 北京航空航天大学 | 一种基于星敏感器确定载体相对于地理坐标系姿态的方法 |
CN102506894A (zh) * | 2011-10-11 | 2012-06-20 | 北京航空航天大学 | 基于精密计时的静基座平台绝对空间姿态基准的建立方法 |
CN102538825A (zh) * | 2011-12-29 | 2012-07-04 | 北京控制工程研究所 | 一种星敏感器探头组合体光轴指向标定方法 |
CN102901485A (zh) * | 2012-10-31 | 2013-01-30 | 中国科学院长春光学精密机械与物理研究所 | 一种光电经纬仪快速自主定向的方法 |
CN103913180A (zh) * | 2014-03-26 | 2014-07-09 | 中国科学院长春光学精密机械与物理研究所 | 一种船载大视场高精度星敏感器的安装角标定方法 |
CN105318891A (zh) * | 2014-07-25 | 2016-02-10 | 北京航天计量测试技术研究所 | 一种星敏感器基准立方镜安装误差的标定装置 |
CN105387874A (zh) * | 2015-12-14 | 2016-03-09 | 中国科学院长春光学精密机械与物理研究所 | 船载高精度星敏感器安装角标定方法 |
CN107883946A (zh) * | 2017-10-31 | 2018-04-06 | 北京航天计量测试技术研究所 | 一种三角形匹配式星敏感器星库的构建方法 |
-
2018
- 2018-11-21 CN CN201811392442.XA patent/CN109459059B/zh active Active
Patent Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS5973719A (ja) * | 1982-10-20 | 1984-04-26 | Natl Space Dev Agency Japan<Nasda> | 人工衛星搭載恒星センサ信号処理方式 |
CN102288199A (zh) * | 2011-06-22 | 2011-12-21 | 哈尔滨工业大学 | 一种星敏感器的地面测试方法 |
CN102279002A (zh) * | 2011-06-27 | 2011-12-14 | 哈尔滨工业大学 | 星敏感器测量坐标系与载体坐标系转换矩阵的标定方法 |
CN102506894A (zh) * | 2011-10-11 | 2012-06-20 | 北京航空航天大学 | 基于精密计时的静基座平台绝对空间姿态基准的建立方法 |
CN102426017A (zh) * | 2011-11-03 | 2012-04-25 | 北京航空航天大学 | 一种基于星敏感器确定载体相对于地理坐标系姿态的方法 |
CN102538825A (zh) * | 2011-12-29 | 2012-07-04 | 北京控制工程研究所 | 一种星敏感器探头组合体光轴指向标定方法 |
CN102901485A (zh) * | 2012-10-31 | 2013-01-30 | 中国科学院长春光学精密机械与物理研究所 | 一种光电经纬仪快速自主定向的方法 |
CN103913180A (zh) * | 2014-03-26 | 2014-07-09 | 中国科学院长春光学精密机械与物理研究所 | 一种船载大视场高精度星敏感器的安装角标定方法 |
CN105318891A (zh) * | 2014-07-25 | 2016-02-10 | 北京航天计量测试技术研究所 | 一种星敏感器基准立方镜安装误差的标定装置 |
CN105387874A (zh) * | 2015-12-14 | 2016-03-09 | 中国科学院长春光学精密机械与物理研究所 | 船载高精度星敏感器安装角标定方法 |
CN107883946A (zh) * | 2017-10-31 | 2018-04-06 | 北京航天计量测试技术研究所 | 一种三角形匹配式星敏感器星库的构建方法 |
Non-Patent Citations (2)
Title |
---|
郭敬明,等: "基于双星敏感器的船体姿态测量系统设计", 《光电子技术》 * |
郭敬明,等: "船载高精度星敏感器安装角的标定", 《光学精密工程》 * |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113916216A (zh) * | 2020-09-09 | 2022-01-11 | 蓝箭航天空间科技股份有限公司 | 惯性测量信息与星光测量信息数据同步方法及系统 |
CN112212856A (zh) * | 2020-09-21 | 2021-01-12 | 中国科学院长春光学精密机械与物理研究所 | 一种近地面白天星敏感器的高精度姿态测量方法 |
CN112833878A (zh) * | 2021-01-05 | 2021-05-25 | 上海航天控制技术研究所 | 近地面多源天文自主导航方法 |
CN113607188A (zh) * | 2021-08-02 | 2021-11-05 | 北京航空航天大学 | 基于经纬仪叉丝成像的多视场星敏感器的标定系统及方法 |
CN114001756A (zh) * | 2021-11-03 | 2022-02-01 | 中国科学院光电技术研究所 | 一种小视场星敏感器外场地面寻星方法 |
CN114001756B (zh) * | 2021-11-03 | 2023-09-19 | 中国科学院光电技术研究所 | 一种小视场星敏感器外场地面寻星方法 |
CN114623833A (zh) * | 2022-01-28 | 2022-06-14 | 北京控制工程研究所 | 一种空间指向测量仪器视场标定系统和方法 |
CN114740801A (zh) * | 2022-03-21 | 2022-07-12 | 成都飞机工业(集团)有限责任公司 | 一种用于数控设备群协同生产线安装的基坐标系创建方法 |
CN114740801B (zh) * | 2022-03-21 | 2023-09-29 | 成都飞机工业(集团)有限责任公司 | 一种用于数控设备群协同生产线安装的基坐标系创建方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109459059B (zh) | 2022-08-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109459059A (zh) | 一种星敏感器外场转换基准测定系统及方法 | |
CN106767714B (zh) | 提高卫星图像定位精度的等效失配模型多级标校方法 | |
CN105737858B (zh) | 一种机载惯导系统姿态参数校准方法与装置 | |
CN107607127B (zh) | 一种基于外场的星敏感器内部参数标定及精度快速验证系统 | |
CN104502926A (zh) | 一种基于精密机械臂的室外gnss接收机绝对天线相位中心校正方法 | |
CN105160125B (zh) | 一种星敏感器四元数的仿真分析方法 | |
CN102261921B (zh) | 一种修正大气折射对星敏感器精度影响的方法 | |
CN104655135B (zh) | 一种基于地标识别的飞行器视觉导航方法 | |
CN105444780B (zh) | 星载摆扫式光学相机图像定位验证系统及处理方法 | |
CN104848874A (zh) | 一种光电经纬仪的外场星校方法 | |
CN105758400B (zh) | 静止卫星成像导航与配准恒星敏感东西参数提取方法 | |
CN108362268A (zh) | 一种基于视频测量的自动天文测量方法和测量系统 | |
CN105444778B (zh) | 一种基于成像几何反演的星敏感器在轨定姿误差获取方法 | |
CN109269525B (zh) | 一种空间探测器起飞或着陆过程光学测量系统及方法 | |
CN102901485B (zh) | 一种光电经纬仪快速自主定向的方法 | |
CN111238540A (zh) | 一种基于恒星拍摄的珞珈一号相机-星敏安装标定方法 | |
CN104807477A (zh) | 一种基于靶标控制点的卫星ccd阵列影像几何检校方法 | |
CN113218577A (zh) | 一种星敏感器星点质心位置精度的外场测量方法 | |
CN106679695B (zh) | 一种基于星敏感器的焦平面热变形的测试方法 | |
CN106996779B (zh) | 基于gnss的紫外敏感器系统误差在轨标定方法 | |
CN105424060B (zh) | 一种飞行器星敏感器与捷联惯组安装误差的测量方法 | |
Liu et al. | Attitude calculation method based on full-sky atmospheric polarization mode | |
CN112097794B (zh) | 遥感卫星载荷平台标定方法及系统 | |
CN101324435A (zh) | 一种获得方向信息的方法与装置 | |
US20230243655A1 (en) | Method for joint estimation of stellar atmospheric refraction and star tracker atitude |
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 |