星敏感器姿态确定仿真综合实验共14页.doc

精选优质文档-----倾情为你奉上航空航天大学专业综合实验报告学 院 宇航学院 班 级 学 号 姓 名 高荣荣指导老师 王海涌 2015年1月3日星敏感器姿态确定仿真综合实验摘要：通过对电子星图模拟器和星敏感器PC仿真平台的操作，实现星敏感器姿态确定，以及借助星象天文馆，来实现电子星图模拟器对星的标定完成星敏感器系统仿真关键词：星敏感器定姿星图模拟星图姿态矩阵一实验目的通过电子星图模拟器（ESS）和星敏感器PC仿真平台的操作，熟悉星图模拟和星敏感器姿态基本流程及各模块功能，完成星敏感器系统仿真二实验原理1.星图模拟原理（1）星图模拟系统是一种近似模拟星空的仿真系统为星敏感器算法调试、星敏感器产品测试及天文导航半物理仿真系统运行提供标准的星图输入，并提供已知参考星光矢量及星像中心的理想映射坐标分光学物理星模、电子星模和计算机软件星模三种2）星图模拟系统的实现旋转关系：OZ是光轴指向：星图模拟是根据光轴指向及旋角(,,γ)所确定的既定视场，将其围的星空目标映射到星敏感器CCD面阵上并模拟出目标图像的过程。

其中，需要确定第二赤道坐标系、航天器本体坐标系、星敏感器坐标系的转换矩阵星敏感器固联（安装矩阵为常数阵），那么只考虑第二赤道坐标系和星敏感器坐标系之间的转换关系令O-UVW为第二赤道坐标系，令O’-XYZ为星敏感器坐标系那么星光矢量在两个坐标系下的分量列阵的关系可以表示为：[X,Y,Z]T = Tsi[U,V,W]T，其中Tsi为转换矩阵2.星敏感器定姿基本原理OsXsYsZs—星敏感器坐标系Ouv— CCD 成像面坐标系OsO之间距离 f 为光学透镜的焦距由图中的几何关系可得：第 n 颗星的单位矢量在星敏感器坐标系中的分量列阵：或，根据星像点质心坐标直接计算得单位矢量：式中为星敏感器测量误差矢量星敏视场 n 颗恒星在星敏感器坐标系Sb中的单位矢量坐标分别为 (X1, Y1, Z1), (X2, Y2, Z2), …, (Xn, Yn, Zn)经过星图匹配，获知这 n 颗恒星在惯性空间i系中单位矢量坐标：[U1, V1, W1], [U2, V2, W2], … , [Un, Vn, Wn]，是由星表中赤经赤纬经过直角坐标转换求得则存在以下矩阵转换关系式：上式各矩阵赋予命名，对应关系式：S = CTis当观测星数 n 等于3颗时，星历矩阵C3×3为非奇异矩阵，那么转移矩阵Tis ＝ C-1S。

n > 3时，可以采用最小二乘法得到计算式：Tis = (CTC)-1CTS 姿态矩阵Tsi = TsbTblTli，则载体相对于发射点惯性系l系的姿态矩阵Tbl = (Tsb)-1 (Tis)-1(Tli)-1 =(Tli Tis Tsb)-1，这里Tsb为星敏感器固联常量安装矩阵，发射点惯性坐标系 l系相对于赤道惯性系i系的转移矩阵：其中：A 为轴xl的方位角；S 为发射时刻的格林尼治恒星时，或发射时刻的春分点的格林时角GHAΥ；(λ, φ) 为发射点经纬度如果安装矩阵与火箭本体坐标系方向重合，则安装矩阵Tbs = 单位阵 I那么：Tbl = (Tli Tis)-1欧拉角法表示火箭姿态发惯系l，原点在发射点L，轴yl在当地铅垂线向上，轴 xl 在当地水平，且在名义射击平面发射点惯性坐标系Lxlylzl转动 3 次转到本体坐标系oxbybzb转动顺序：那么，从Sl到Sb的坐标变换矩阵：矩阵元素对照得主值：姿态角θ和ψ的取值围都在[-90°,90°]；φ的取值围都在[-180°,180°]然后基于主值再获得真值三实验步骤1.设定观星地址和观星时间，采用静基座天顶观测方式，解算出相机光轴的赤经、赤纬和旋角。

1）设定观星地点为市延庆县马匹营村，北纬40°30′14.571″，东经116°04′20.815″（2）设定时间为2011年8月20日22h22min222s（3）2011年8月20日当天世界时0时对应的恒星时为12h16m40.6407s;=26.3156由此得RA=300.2372Dec=40.5041指定γ=02开启ESS，进行参数设置，选择“星图模拟”“单帧重复”菜单命令，将（，，）计算值填入ESS的人机界面，观察“发送的星图”，此时ESS持续不断地将模拟星图数字视频传输给星敏感器PC仿真平台设置：指定仿真星等阈值：6峰值饱和灰度对应星等：8得到模拟星图：其中，星模软件使用方法见附录；3.在PC仿真平台，存储接收的模拟星图，并选定六颗颗亮星结果如图所示得到亮星的坐标，星等，赤经和赤纬：310.45.1.3243.873.305.40.2.3514.718.311.33.2.6810.993.296.45.3.0272.369.303.46.4.0190.621.299.35.4.0783.466.4.运行虚拟天文馆软件stellarium.exe，选定（1）中地址和观星时间，目视对比模拟星图和stellarium界面，找出这六颗选定亮星的匹配星，记录其星号、赤经、赤纬和星等属性信息，计算在春分点赤道惯性系i系下的恒星方位矢量。

运行虚拟天文馆结果如图所示：δCyg-18O1 Cyg-31天津四ηCyg-21天津九天津一根据经纬度关系，并借助主观肉眼判断，匹配情况如下：星名/星等赤经/（J2000）赤纬天津四/1.25（脉动变星）310.25/20h41m25s45.16天津一/2.20（恒星）305.52/20h22m13.7s40.17天津九/2.45（恒星）311.50/20h16m41s34.01δCyg-18Cyg-HIP 97165 A/2.9（双星）296.25/19h45m20s45.07O1 Cyg-31 Cyg-V0695/3.8（食双星系统）303.404/20h13m37.9s46.46ηCyg-21 Cyg-HIP 98110/3.85(恒星)299.011/19h56m45s35.06注意事项：（1）在运用虚拟天文馆软件时，要注意调节观察视角，一般为20度左右事宜本次实验，调节观察视角读数为22.8°2）根据纬度关系显示，将虚拟天文馆的视场调节到指定的围，以便容易找到星星本次实验，大约根据纬度关系调节到45°的围3）在记录赤经、赤纬的信息时，有两种表示的方法：一种是J2000，另一种是J2011.2。

两者之间有非常细微的差距因为地轴的进动，使得春分点不断西退，而赤经赤纬又是依靠春分点的位置来确定，所以随着春分点位置的不断变化，赤经赤纬也会不断变化J2000版指的是以2000年春分点作为基准天体的坐标，日期版就是以现在的春分点作为基准的坐标差距是非常细微的，对一般的爱好者来说，忽略这个微小的差距也无妨，但是对于专业的观测来说，这点误差是要命的在本实验中，记录的是J2000下的赤经赤纬信息记录星等信息时，有目视星等和绝对星等两种星等信息绝对星等是假定把恒星放在距地球10秒差距（32.6光年）的地方测得的恒星的亮度，用以区别于视星等它反映天体的真实发光本领由于本实验中记录的是目视星等5.给定该相机焦距f的标定值，查找星敏感器PC平台下的日志文件获取6颗选定亮星的星象质心，是在星图预处理环节由灰度重心法计算获得的，计算星敏感器本体系s系下的星光矢量1).焦距f=2889.2643，单位：pixel，主点O的值(508.3, 513.3) pixel2).阵列平面坐标系和星敏感器本体坐标系的关系如图所示：其中星敏感器本体系s系下的星光矢量由下式计算得到6.基于i系下的恒星方位矢量和s系下的星光矢量，计算出观星时刻星敏感器相对于i系的姿态矩阵，并换算出星敏感器的光轴赤经赤纬和旋角，对比星敏感器PC仿真平台的计算结果。

Matlab程序见附录；得到赤经赤纬及旋角为：RA=300.2372Dec=40.5041四实验结论及误差分析：1.以观测者所在位置，运用查表法得到的赤经赤纬为Ra=300.2372，Dec=40.5041，Gamma=0；再运用模拟星图ESS和虚拟天文馆得到的赤经赤纬为RaDeter=300.2385，DecDeter=40.5180，GammaDeter=-0.3402;所以，经运算得到绝对误差和相对误差如下：绝对误差为：|Ra-RaDeter|=0.0013 |Dec-DecDeter|=0.0139 |Gamma-GammaDeter|=0.3402相对误差为：η(Ra)=0.004%η(Dec)=0.031%η(Gamma)=0.340%2.在星图模拟和虚拟天文馆的对比找星过程中，也存在一定的误差，主要是，星等误差，赤经和赤纬对比误差等等误差定量计算如下：（1）赤经绝对误差：绝对误差310.25/20h41m25s310.0.100305.52/20h22m13.7s305.0.030311.50/20h16m41s311.0.050296.25/19h45m20s296.0.010303.404/20h13m37.9s303.0.003299.011/19h56m45s299.0.065分析：最大误差为0.1°，因此，在合理误差围之，说明星星确认正确；（2）赤纬绝对误差：绝对误差45.1645.0.1240.1740.0.0834.0133.0.0445.0745.0.0646.4646.0.2835.0635.0.02分析：最大误差为0.28，计算相对误差为0.6%，误差总体较小，因此也在合理围之。

说明确认星星正确；（3）星等绝对误差绝对误差1.251.30.052.202.30.102.452.60.052.93.00.103.84.00.203.854.00.15分析：最大误差为0.20，计算相对误差一定比较大，但是不影响星星判定，星等的误差，大多还是由于虚拟天文馆的制作造成的，因此，从整体绝对误差上看，星星的误差还是在合理围之的因此，认为确定星星正确综上所述：经过赤经，赤纬，星等的误差对比，确认模拟星图所找的星星在虚拟天文馆中得到一一对应误差可视为合理误差五实验总结为期一周的综合实验到这里结束了，首先，要感王海涌老师的悉心教导，在王老师的帮助指导下，将天文导航这门课程的精髓在实验中得到了应用，收获颇丰。