一种基于可见光谱的空间碎片姿态分析方法技术

一种基于可见光谱的空间碎片姿态分析方法，包括如下步骤：步骤1，空间碎片光谱数据采集，记为观测谱；步骤2，光谱数据预处理；步骤3，建立空间碎片的三维模型，确定各面元的法线方向及材料谱；步骤4，利用在轨测量，确定空间碎片质心到测站的方向矢量和空间碎片质心到太阳的方向矢量；步骤5，记空间碎片的偏航角、俯仰角、滚动角，计算旋转矩阵R；步骤6，计算在本体参考系下的空间碎片指向测站的方向矢量和空间碎片指向太阳的方向矢量；步骤7，由朗伯散射模型计算预测谱；步骤8，计算预测谱与观测谱之差的欧氏范数(残余)作为姿态指向判别度量；步骤9，估计目标指向，以最小残余所对应的偏航角、俯仰角、滚动角作为空间碎片姿态的空间指向。本发明专利技术有效适用于深空环境、准确性更高。

【技术实现步骤摘要】

本专利技术涉及载人航天实施空间监测
，尤其是一种基于可见光谱的空间碎片姿态分析方法。
技术介绍
自1957年10月前苏联发射首颗人造地球卫星以来，美国、法国、日本也相继发射了人造地球卫星，中国也于1970年发射了第一颗自己的人造地球卫星“东方红一号”；据2013年统计，截止2013年5月31号，在轨运行的人造地球卫星总数为1071颗；而目前具体有多少颗卫星，并不确定，因为该数字几乎每天都会在变，有新的卫星发射，也有老卫星失效后被废弃；并且随着小(纳、皮)卫星的发射，其数量变化更快。人类几十年的空间探索活动产生了大量的空间碎片，对人类航天活动的安全造成严重威胁，成为空间环境的主要污染源，并在一定程度上对航天活动的正常开展产生了影响。传统空间碎片测量以位置信息测量为主，包括空间碎片的三维位置坐标、速度、加速度等参数，可衍生出各类地球轨道参数。为了提高空间碎片的监测预警能力，对测量系统除了要求获得其位置信息之外，更需要获得空间碎片的特征信息，如空间碎片姿态等特征信息，它对空间监测、预警尤为重要。光谱测量技术是空间碎片状态分析的一种重要方法，空间碎片一般本身并不发光，其亮度来自太阳光的反射，仅由亮度变化不足以分辨空间碎片的材料组成及空间指向，分析其光谱特征成为我们辨认空间碎片姿态的主要手段。在地基(地面观测站)条件下，光谱仪所获取光谱由空间碎片材料成分、太阳光谱、地球大气吸收谱，以及测量过程中产生的噪声等因素综合决定，因而由光谱推断空间碎片指向并没有一个解析解。本专利技术采用穷举搜索方法，对空间碎片姿态指向作判别。空间碎片的测量有雷达测量和光学测量，雷达测量可以克服天气、太阳和及地影的影响，能全天候全天时工作，但由于雷达测量时，其反射回波信号的强度与距离的四次方成反比，因而雷达测量比较适合于低地球轨道的空间碎片；而对于无源光学测量，信号反射强度与物体的距离的平方成反比，因而能探测高轨道空间碎片，另外它仅接收空间碎片对太阳光的反射，属于被动探测形式，可较好的适合于某些特殊场合。
技术实现思路
为了克服已有空间碎片姿态指向分析方法无法适用于深空环境、准确性较差的不足，本专利技术提供一种有效适用于深空环境、准确性更高的基于可见光谱的空间碎片姿态分析方法。本专利技术解决其技术问题所采用的技术方案是：一种基于可见光谱的空间碎片姿态分析方法，所述方法包括如下步骤：步骤1，空间碎片光谱数据采集：利用光谱仪，采集空间碎片光谱数据，记为观测谱Ytest；步骤2，光谱数据预处理：进行去噪处理；步骤3，建立空间碎片的3维模型，确定各面元的法线方向及材料谱；步骤4，利用在轨测量，确定空间碎片质心到测站的方向矢量o*和空间碎片质心到太阳的方向矢量s*；步骤5，记空间碎片的偏航角Yaw、俯仰角Pitch、滚动角Roll，计算旋转矩阵R，有：R＝Rx(Roll)*Ry(Pitch)Rz(Yaw)其中： R x ( R o l l ) = 1 0 0 0 c o s ( R o l l ) - s i n ( R o l l ) 0 s i n ( R o l l ) cos ( R o l l ) , R y ( P i t c h ) = c o s ( P i t c h ) 0 本文档来自技高网...

【技术保护点】
一种基于可见光谱的空间碎片姿态分析方法，其特征在于：所述方法包括如下步骤：步骤1，空间碎片光谱数据采集：利用光谱仪，采集空间碎片光谱数据，记为观测谱Ytest；步骤2，光谱数据预处理：进行去噪处理；步骤3，建立空间碎片的3维模型，确定各面元的法线方向及材料谱；步骤4，利用在轨测量，确定空间碎片质心到测站的方向矢量o*和空间碎片质心到太阳的方向矢量s*；步骤5，记空间碎片的偏航角Yaw、俯仰角Pitch、滚动角Roll，计算旋转矩阵R，有：R＝Rx(Roll)*Ry(Pitch)Rz(Yaw)其中：Rx(Roll)=1000cos(Roll)-sin(Roll)0sin(Roll)cos(Roll),Ry(Pitch)=cos(Pitch)0sin(Pitch)010-sin(Pitch)0cos(Pitch)]]>Rz(Yaw)=cos(Yaw)-sin(Yaw)0sin(Yaw)cos(Yaw)0001]]>步骤6，计算在本体参考系下的空间碎片指向测站的方向矢量o和空间碎片指向太阳的方向矢量ss=R*s*o=R*o*]]>步骤7，由朗伯散射模型计算预测谱：Ypred=A*a1<n1,o><n1,s>...ak<nk,o><nk,s>]]>其中，A为k个材料谱所组成矩阵，ai表示第i面元的面积，ni表示第i面元的法线矢量，其中，点积运算符<·,·>还要求满足，对任意实数y∈R，y表示由点积运算结果值：<y>=yy>00y≤0]]>步骤8，计算预测谱Ypred与观测谱Ytest之差的欧氏范数(残余)作为姿态指向判别度量，记残余：残余＝||Ypred‑Ytest||2步骤9，估计空间碎片指向，以最小残余所对应的偏航角Yaw、俯仰角Pitch、滚动角Roll作为空间碎片姿态的空间指向。...

【技术特征摘要】
1.一种基于可见光谱的空间碎片姿态分析方法，其特征在于：所述方
法包括如下步骤：
步骤1，空间碎片光谱数据采集：利用光谱仪，采集空间碎片光
谱数据，记为观测谱Ytest；
步骤2，光谱数据预处理：进行去噪处理；
步骤3，建立空间碎片的3维模型，确定各面元的法线方向及材
料谱；
步骤4，利用在轨测量，确定空间碎片质心到测站的方向矢量o*和
空间碎片质心到太阳的方向矢量s*；
步骤5，记空间碎片的偏航角Yaw、俯仰角Pitch、滚动角Roll，
计算旋转矩阵R，有：
R＝Rx(Roll)*Ry(Pitch)Rz(Yaw)
其中：
R x ( R o l l ) = 1 0 0 0 cos ( R o l l ) - sin ( R o l l ) 0 sin ( R o l l ) cos ( R o l l ) , R y ( P i t c h ) = cos ( P i t c h ) 0 sin ( P i t c h ) 0 1 0 - sin ( P i t c h ) 0 cos ( P i t c h ) ]]> R z ( Y a w ) = c o s ( Y a w ) - s i n ( Y a w ...

【专利技术属性】
技术研发人员：庄德文，唐轶峻，康湘华，夏晨阳，
申请(专利权)人：浙江工业大学，
类型：发明
国别省市：浙江;33