星敏感器姿态确定仿真综合实验

航空航天大学专业综合实验报告学院宇航学院班级111514学号11151146姓名高荣荣指导老师王海涌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.星敏感器定姿基本原理nnfYs光轴ZsXsOsOuvpn(un,vn)w透镜第n颗星pun(un,0)OsXsYsZs—星敏感器坐标系Ouv—CCD成像面坐标系OsO之间距离f为光学透镜的焦距由图中的几何关系可得：第n颗星的单位矢量在星敏感器坐标系中的分量列阵：tantan/cosnnnnnufvfarctanarctan/cosnnnnnufvf或，根据星像点质心坐标直接计算得单位矢量：式中为星敏感器测量误差矢量。星敏视场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。n3时，可以采用最小二乘法得到计算式：Tis=(CTC)-1CTS姿态矩阵Tsi=TsbTblTli，则载体相对于发射点惯性系l系的姿态矩阵Tbl=(Tsb)-1(Tis)-1(Tli)-1=(TliTisTsb)-1，这里Tsb为星敏感器固联常量安装矩阵，发射点惯性坐标系l系相对于赤道惯性系i系的转移矩阵：其中：A为轴xl的方位角；sincoscoscossinnnnnnnnsnnXsYVZ2221nnnnsnnnnXusYfVuvfZv111111111213222222212223313233nnnnnnXYZUVWtttXYZUVWttttttXYZUVWcossincos()sinsin()coscos()sinsincos()cossin()cossin()sincos()coscoscossin()sinsinsinsin()coscos()sincoslileeiASASTTTSASASASASASASASAS为发射时刻的格林尼治恒星时，或发射时刻的春分点的格林时角GHAΥ；(λ,φ)为发射点经纬度。如果安装矩阵与火箭本体坐标系方向重合，则安装矩阵Tbs=单位阵I那么：Tbl=(TliTis)-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.3156cossin0sincos0001ieGHAGHACGHAGHAsinsincoscoscoscossinsincossin0cossinetC111213212223313233bltttTtttttt()()()yxzRRRlbSScoscoscossinsinsinsincoscossinsinsinsincoscossincoscossincossinsincossinsinsincoscoscosblT2312131133arctan();arcsin();arctan()ttttt00000sincossincoscos100etieiiiCCzyx由此得RA=300.2372Dec=40.5041指定γ=02开启ESS，进行参数设置，选择“星图模拟”“单帧重复”菜单命令，将（，，）计算值填入ESS的人机界面，观察“发送的星图”，此时ESS持续不断地将模拟星图数字视频传输给星敏感器PC仿真平台。设置：指定仿真星等阈值：6峰值饱和灰度对应星等：8得到模拟星图：其中，星模软件使用方法见附录；3.在PC仿真平台，存储接收的模拟星图，并选定六颗颗亮星。结果如图所示。得到亮星的坐标，星等，赤经和赤纬：310.358041845.280358371.3243.0664099873.2034322305.557124840.256699942.3514.3430786718.3319065311.552858233.970255722.6810.0597353993.0520476296.243895845.130739033.0272.523317369.3284183303.407916546.741358354.0190.4634936621.7853499299.076604335.083419324.0783.1009519466.73437074.运行虚拟天文馆软件stellarium.exe，选定（1）中地址和观星时间，目视对比模拟星图和stellarium界面，找出这六颗选定亮星的匹配星，记录其星号、赤经、赤纬和星等属性信息，计算在春分点赤道惯性系i系下的恒星方位矢量。运行虚拟天文馆结果如图所示：根据经纬度关系，并借助主观肉眼判断，匹配情况如下：星名/星等赤经/（J2000）赤纬天津四/1.25（脉动变星）310.25/20h41m25s45.16天津一/2.20（恒星）305.52/20h22m13.7s40.17天津九/2.45（恒星）311.50/20h16m41s34.01δCyg-18Cyg-HIP97165A/2.9（双星）296.25/19h45m20s45.07O1Cyg-31Cyg-V0695/3.8（食双星系统）303.404/20h13m37.9s46.46ηCyg-21Cyg-HIP98110/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)pixel。(2).阵列平面坐标系和星敏感器本体坐标系的关系如图所示：天津四天津一天津九δCyg-18ηCyg-21O1Cyg-31O(u0,v0)XYfZ成像阵列遮光罩Q(un,vn)(0,0)u,row(1024,1024)v,col光轴主点funn0u-tann0/cosv-tanfvnn其中星敏感器本体系s系下的星光矢量由下式计算得到sincoscoscossinnnnnnnnsnnXsYVZ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.35804180.100305.52/20h22m13.7s305.55712480.030311.50/20h16m41s311.55285820.050296.25/19h45m20s296.24389580.010303.404/20h13m37.9s303.40791650.003299.011/19h56m45s299.07660430.065分析：最大误差为0.1°，因此，在合理误差围之，说明星星确认正确；（2）赤