一种基于RCS测量数据的自旋卫星姿态估计方法及系统技术方案

技术编号:18457915 阅读:43 留言:0更新日期:2018-07-18 12:21
本发明专利技术属于自旋卫星姿态估计技术,具体涉及一种基于RCS测量数据的自旋卫星姿态估计方法及系统,具体公开了一种基于RCS测量数据的自旋卫星姿态估计方法,包括如下步骤:S1:监测模块监测卫星在测站坐标系下卫星的位置矢量;S2:获取目标RCS序列;S3:根建立目标RCS序列库;S4:采用处理器内置的算法模块实现RCS序列的周期估计;S5:对自旋稳定卫星姿态估计进行参数设置;S6:设置粒子群库内的粒子群种群参数,初始化粒子群种群;S7:评估粒子群种群,获取最优粒子,若迭代次数达到设定的阈值,则结束反演;否则,转S8。本发明专利技术基于多频段RCS的自旋卫星姿态估计模型,采用改进粒子群算法进行优化求解,实现了自旋卫星在轨姿态估计。

A method and system for attitude estimation of spinning satellites based on RCS measurement data

The invention belongs to the spin satellite attitude estimation technique, which specifically involves a spin satellite attitude estimation method and system based on the RCS measurement data, and specifically discloses a spin satellite attitude estimation method based on the RCS measurement data, including the following steps: S1: the monitoring module monitors the position vector of the satellite in the station coordinate system. Quantity; S2: obtain the target RCS sequence; S3: the root establishes the target RCS sequence library; S4: uses the processor's built-in algorithm module to realize the periodic estimation of the RCS sequence; S5: the parameter setting for the spin stabilized satellite attitude estimation; S6: setting the particle swarm population parameters in the particle swarm library, initializing the particle swarm population; S7: evaluating grain The optimal particle is obtained by subgroup population. If the number of iterations reaches the set threshold, the inversion will be completed; otherwise, S8 will be transferred. Based on the multi band RCS's spin satellite attitude estimation model, the improved particle swarm optimization (PSO) algorithm is used to achieve the attitude estimation of the spin satellite in orbit.

【技术实现步骤摘要】
一种基于RCS测量数据的自旋卫星姿态估计方法及系统
本专利技术属于自旋卫星姿态估计技术,具体涉及一种基于RCS测量数据的自旋卫星姿态估计方法及系统。
技术介绍
随着各国航天事业的快速发展,在轨运行卫星的种类和数量迅速增加,空间环境越来越恶劣,对探测空间目标的在轨状态提出了更高的要求。空间目标的在轨状态主要有自旋稳定、自旋稳定、重力梯度稳定、翻滚等几类,而周期是衡量自旋稳定、翻滚空间目标在轨状态的一个重要参数。目前,利用非合作形式测量得到的特性数据估计空间翻滚目标周期的一种重要手段。实际测量中空间翻滚目标雷达散射截面(RadarCrossSection,RCS)受空间目标形状、姿态和雷达性能参数等因素影响,周期估计多根据RCS测量值的变化特征利用人工经验判别。文献[W-JZhong,J.-SWangandW-JJi.Theattitudeestimationofthree-axisstabilizedsatellitesusinghybridparticleswarmoptimizationcombinedwithradarcrosssectionpreciseprediction.ProcIMechEPartG:JournalofAerospaceEngineering2015;1-13.]提出变区间分组检验相乘积累进动周期估计方法对进动锥体目标RCS特性数据的统计分析,该方法存在倍频分频问题且运算量较大;文献[LuNL,JinJM.Applicationoffastmultipolemethodtofiniteelementboundary-integralsolutionofscatteringproblems.IEEETransactionsonantennasandpropagation1996;44:781-786.]利用方差分析法估计了空间目标RCS序列的数据周期,该方法同样存在倍频分频问题,且敏感于RCS序列整体的升降趋势;文献[RobertAM,JenniferLD,JohnLB.Aninformationinfrastructureforcoordinateearthscienceoberservations.Proceedingsofthe2ndIEEEInternationalConferenceonSpaceMissionChallengesforInformationTechnology2006;397-404.]利用循环自相关和循环平均幅度差函数相结合的方法估计导弹目标的进动周期,由于自相关与平均幅度差函数本身都存在倍频分频问题,两者结合仅仅改变了各分量的绝对大小,相对大小并没有改变,因此仍然存在错估的问题;文献[HarrubggtonRF.Fieldcomputationbymomentmethods.NewYork:Macmillan1968.]提出了基于非参数秩方差检验的经验模态分解的周期估计方法,该方法能有效克服虚假周期影响,且能改善翻滚目标周期估计精度。以上方法均是直接对测站采集的RCS序列进行分析获得空间目标的RCS,忽略了空间目标旋转、姿态等因素对利用RCS序列进行周期估计的影响。
技术实现思路
有鉴于此,本专利技术的主要目的是提供一种基于RCS测量数据的自旋卫星姿态估计方法及系统。本专利技术采用的技术方案是:一种基于RCS测量数据的自旋卫星姿态估计方法,包括如下步骤:S1:至少设置一个测站,依据测站设置的监测模块监测卫星在测站坐标系下卫星的位置矢量;S2:根据测站坐标系下卫星的位置矢量,获取卫星在测站坐标系下的俯仰角和方位角,依据俯仰角和方位角获得测站相对卫星的可见弧段,在俯仰角的弧段对测站进行引导跟踪卫星,获取目标RCS序列;S3:根据测站的监测模块跟踪的目标RCS序列,建立目标RCS序列库;S4:采用处理器内置的算法模块实现RCS序列的周期估计;S5:根据获取的先验信息,对自旋稳定卫星姿态估计进行参数设置;S6:设置粒子群库内的粒子群种群参数,初始化粒子群种群;S7:评估粒子群种群,获取最优粒子,若迭代次数达到设定的阈值,则结束反演;否则,转S8;S8:更新粒子群种群,否则转S7。进一步地,在S7中,最优粒子的获取方法如下:根据监测模块获取的轨道根数以及雷达视线角得到目标本体下入射角;根据入射角序列利用处理器内的计算仿真模块对RCS序列进行处理;根据计算模块获取单个粒子的适应度值;优选最优粒子。进一步地,所述测站坐标系下卫星在t时刻的位置矢量为r上述式中,ρ=(M)[(HG)r-Rb]HG=(B2)(B1)(N)(A)M为地球固定坐标系至J2000.0惯性坐标系的转换矩阵;Rb为站心在地球固定坐标下的位置矢量,B2为准地球固定坐标系至地球固定坐标下的转换矩阵;B1为瞬时真赤道坐标系至准地球固定坐标系下的转换矩阵;N为瞬时平赤道坐标系至瞬时真赤道坐标系下的转换矩阵;A为J2000.0惯性坐标系至瞬时平赤道坐标系下的转换矩阵;当俯仰角E>0时,即可满足测站对卫星的可见条件,获取卫星外推弹道每一时刻在测站坐标系下俯仰角和方位角,可获得测站相对卫星的可见弧段,即可在俯仰角E>0的弧段对测站进行引导跟踪卫星。进一步地,在S5中,自旋稳定卫星姿态估计为对自旋稳定卫星在卫星本体系下的主轴姿态进行反演,根据测站监测模块测量RCS序列同计算仿真模块RCS序列的残差的目标函数为式中:RCS_meas(m,n)是第m个测站监测模块测量的第n个弧段的RCS测量序列;RCS_simu(m,n)是第m个监测模块测量的第n个弧段的RCS仿真序列;Nstation是应用于卫星姿态反演的测站数;Narc是每个测站的采用的测量弧段数。因而,自旋稳定卫星的姿态确定的优化模型如下上述式中,卫星绕主轴旋转的角速度为ω,本专利技术还提供了一种基于RCS测量数据的自旋卫星姿态估计系统,包括设置至少一个跟踪卫星的测站,所述测站设置有监测卫星在测站坐标系下卫星位置矢量的监测模块,根据所述测站坐标系下卫星的位置矢量,获取目标RCS序列;设置一个RCS序列库,根据目标RCS序列库优选周期特征明显的RCS序列;依据处理器内置的算法模块实现优选RCS序列的周期估计。本专利技术基于多频段RCS的自旋卫星姿态估计模型。采用电磁场数值算法[5-6]快速计算卫星的雷达散射截面,推导了自旋模式下卫星本体坐标系下电磁波入射角的计算公式,利用RCS匹配获得卫星可跟踪弧段的理论RCS序列。采用改进粒子群算法进行优化求解,实现了自旋卫星在轨姿态估计。附图说明图1为本专利技术中流程图;图2为本专利技术中测站跟踪卫星的原理图;图3(a)、3(b)、3(c)、3(d)为本专利技术的仿真算法验证结果图;图4为本专利技术中卫星的单站RCS三维图;图5为本专利技术中PS-1A卫星RCS反射图;图6(a)和图6(b)基于RCS序列的卫星旋转周期估计;图7粒子群算法的收敛特性曲线;图8粒子群算法的收敛特性曲线。具体实施方式下面将结合附图以及具体实施例来详细说明本专利技术,在此本专利技术的示意性实施例以及说明用来解释本专利技术,但并不作为对本专利技术的限定。参照图1至图2,一种基于RCS测量数据的自旋卫星姿态估计方法,包括如下步骤:S1:至少设置一个测站,依据测站设置的监测模块监测卫星在测站本文档来自技高网
...

【技术保护点】
1.一种基于RCS测量数据的自旋卫星姿态估计方法,其特征在于,包括如下步骤:S1:至少设置一个测站,依据测站设置的监测模块监测卫星在测站坐标系下卫星的位置矢量;S2:根据测站坐标系下卫星的位置矢量,获取卫星在测站坐标系下的俯仰角和方位角,依据俯仰角和方位角获得测站相对卫星的可见弧段,在俯仰角的弧段对测站进行引导跟踪卫星,获取目标RCS序列;S3:根据测站的监测模块跟踪的目标RCS序列,建立目标RCS序列库;S4:采用处理器内置的算法模块实现RCS序列的周期估计;S5:根据获取的先验信息,对自旋稳定卫星姿态估计进行参数设置;S6:设置粒子群库内的粒子群种群参数,初始化粒子群种群;S7:评估粒子群种群,获取最优粒子,若迭代次数达到设定的阈值,则结束反演;否则,转S8;S8:更新粒子群种群,否则转S7。

【技术特征摘要】
1.一种基于RCS测量数据的自旋卫星姿态估计方法,其特征在于,包括如下步骤:S1:至少设置一个测站,依据测站设置的监测模块监测卫星在测站坐标系下卫星的位置矢量;S2:根据测站坐标系下卫星的位置矢量,获取卫星在测站坐标系下的俯仰角和方位角,依据俯仰角和方位角获得测站相对卫星的可见弧段,在俯仰角的弧段对测站进行引导跟踪卫星,获取目标RCS序列;S3:根据测站的监测模块跟踪的目标RCS序列,建立目标RCS序列库;S4:采用处理器内置的算法模块实现RCS序列的周期估计;S5:根据获取的先验信息,对自旋稳定卫星姿态估计进行参数设置;S6:设置粒子群库内的粒子群种群参数,初始化粒子群种群;S7:评估粒子群种群,获取最优粒子,若迭代次数达到设定的阈值,则结束反演;否则,转S8;S8:更新粒子群种群,否则转S7。2.根据权利要求1所述的基于RCS测量数据的自旋卫星姿态估计方法,其特征在于,在S7中,最优粒子的获取方法如下:根据监测模块获取的轨道根数以及雷达视线角得到目标本体下入射角;根据入射角序列利用处理器内的计算仿真模块对RCS序列进行处理;根据计算模块获取单个粒子的适应度值;优选最优粒子。3.根据权利要求1所述的基于RCS测量数据的自旋卫星姿态估计方法,其特征在于,所述测站坐标系下卫星在t时刻的位置矢量为r上述式中,M为地球固定坐标系至J2000.0惯性坐标系的转换矩阵;Rb为站心在地球固定坐标下的位置矢量,B2为准地球固定坐标系至地球固定坐标...

【专利技术属性】
技术研发人员:牛威
申请(专利权)人:西安四方星途测控技术有限公司北京四方继保自动化股份有限公司
类型:发明
国别省市:陕西,61

网友询问留言 已有0条评论
  • 还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。

1