A novel satellite image encryption algorithm based on hyperchaotic systems and Josephus problem

With the rapid growth of the number of Earth observation satellite (EOS) supporting critical applications, it is required to improve the security techniques to protect the sensitive data and images during the transmission between the satellites and the ground stations. This paper introduces a new satellite image encryption algorithm based on the Linear Feedback Shift Register (LFSR) generator, SHA 512 hash function, hyperchaotic systems, and Josephus problem. LFSR generates a matrix that is used to construct the 512-bits value of the hash function. These bits are used to set the initial values and parameters of the proposed encryption algorithm. Firstly, the six dimensions (6-D) hyperchaotic system is divided into three parts, where every two equations are considered as one part. Secondly, the 1-D hyperchaotic logistic-tent system is considered as the controller to select one part. The selected part is used to generate a matrix that is XORed with the original image. Thirdly, the scrambling operation by Josephus sequences is applied to the output of the previous step by scrambling the rows and the columns according to the selected part to produce the pre-encrypted image. Finally, if the number of iterations is less than the required number which is considered as a parameter of the secret key, the previous operations will be repeated in the pre-encrypted image; otherwise, the pre-encrypted image is considered as the final cipher image. Experimental and analyses results show that the proposed algorithm has good performance in terms of high level of security, large enough key-space, tolerance to Single Event Upsets (SEU) as well as low time complexity.     

轨控大干扰下的系统角动量管理

摘要： 针对利用角动量交换装置系统实现轨控期间姿态维持的情况，研究当轨控推力器存在大力矩扰动时系统角动量管理方法.结合对地运行航天器轨道运动特性，分析轨控干扰力矩的积累角动量变化规律；考虑角动量管理装置系统的角动量存贮容量，基于角动量积累规律提出了一种结合偏置角动量建立与对称分布轨控位置选择的系统角动量管理方法，实现了对称轨控位置积累角动量互相抵消的自平衡效果，极大程度地提高轨控效率且避免系统角动量饱和现象.所提出的方法适用于喷气欠驱动控制系统，方法的有效性通过了数学仿真验证及在轨型号应用.     

Identification of seismic activity sources on the subsatellite track by ionospheric plasma disturbances detected with the Sich-2 onboard probes

Using a cylindrical Langmuir probe and the authors’ proprietary two-channel pressure transducer, ionospheric plasma parameter distributions along the orbit of the Sich-2 satellite (Ukraine, 2011–2012) were measured. This paper is concerned with identifying the space–time location of ionospheric plasma disturbance sources, including the epicenters of actual earthquakes (before or during the satellite flyover) and incipient earthquakes on the subsatellite track, from the measured distributions of the electron density and temperature and the neutral particle temperature along the satellite orbit. To do this, the measured ionospheric plasma parameter distributions are connected to the coordinates on the subsatellite track.It is shown that local disturbances in the electron density and temperature and neutral particle temperature distributions in the satellite orbit in the ionosphere may serve as indicators of seismic activity on the subsatellite track. The epicenters of incipient earthquakes may be set off from other plasma parameter disturbance sources associated with seismic activity using information provided by special monitoring and survey centers that monitor the current seismic situation.     

Early analysis of precise orbit and clock offset determination for the satellites of the global BeiDou-3 system

Eight new-generation BeiDou satellites (BeiDou-3) have been launched into Medium Earth Orbit (MEO), allowing for global coverage since March 2018, and they are equipped with new hydrogen atomic clocks and updated rubidium clocks. Firstly, we analyzed the signals for the carrier-to-noise-density ratio (C/N0) and pseudorange multipath (MP) by using international GNSS (Global Navigation Satellite System) Monitoring and Assessment System (iGMAS) station data, and found that B1C has a lower C/N0, and B2a has the same level of C/N0 as the B1I and B3I signals. For pseudorange multipath, compared with the BeiDou-2 satellites, the obvious systematic variation of MP scatters related to the elevation angle is greatly improved for the BeiDou-3 and BeiDou-3e satellites signals. For the signals of the BeiDou-3 satellites, the order of the Root Mean Square (RMS) values of multipath and noise is B3I?<?B1I?<?B2a?<?B1C. Then, the comparison of the precise orbit determination and clock offset determination for the BeiDou-2, BeiDou-3, and BeiDou-3 experimental (BeiDou-3e) satellites was done by using 10 stations from iGMAS. The 3D precision of the 24?h orbit overlap is 24.55, 25.61, and 23.35?cm for the BeiDou-3, BeiDou-3e, and BeiDou-2 satellites, respectively. BeiDou-3 satellite has a comparable precision to that of the BeiDou-2 satellite. For the precision of clock offset estimation, the Standard Deviation (STD) of the BeiDou-3 MEO satellite is 0.350?ns, which is an improvement of 0.042?ns over that of the BeiDou-2 MEO satellite. The stabilities of the BeiDou-3 and BeiDou-3e onboard clocks are better than those of BeiDou-2 by factors of 2.84 and 1.61 at an averaging time of 1000 and 10,000?s, respectively.     

超低轨道卫星应用离子电推进技术方案

低地球轨道大气环境对诸如科学探测和对地观测卫星的阻尼作用十分明显，而且阻尼随太阳和地磁活动以及昼夜、季节交替变化范围宽。为了保证卫星轨道精度或飞行状态满足任务要求，需要利用推进系统对卫星受到的阻尼进行实时或间歇式补偿以实现轨道或飞行状态的保持。针对轨道高度220~268km的无拖曳飞行和轨道维持应用，基于卫星轨道阻尼变化和有效载荷指标要求分析，研究确定了离子电推进技术指标、推力调节方案、系统组成、推力控制方案和在轨应用策略，并对推力调节方案进行了试验验证。结果表明，与无拖曳飞行卫星任务匹配的离子电推进指标为推力调节范围1~20mN，推力分辨率优于12μN，与对地观测卫星轨道维持任务匹配的指标为推力调节范围1～25mN，推力分辨率100μN。研究提出的针对超低轨道卫星应用需求的高精度推力连续调节离子电推进技术方案，具有工程任务针对性和参考价值。     

Online scheduling of distributed Earth observation satellite system under rigid communication constraints

This paper focused on online scheduling of distributed Earth observation satellite system in a dynamic environment. The objective was to maximize the total profit of the overall system by efficiently coordinating the different satellites with stochastic arrival of urgent tasks, subject to rigid communication and observation time window constraints. We formulated this problem as a single-task, single-robot, time-extended assignment problem with intra-schedule dependency using the multi-robot task assignment taxonomy and formulated the sub-problem after releasing urgent tasks in a mixed-integer linear programming model. We first described the online scheduling algorithm for a single satellite, then we proposed the modified consensus-based bundle algorithm (m-CBBA) and modified asynchronous consensus-based bundle algorithm (m-ACBBA) with synchronous and asynchronous communication, respectively. Compared with initial versions of CBBA and ACBBA, the modified versions added the communication loop prediction phases to efficiently utilize scarce communication opportunities and reduce the communication requirements. Additionally, we introduce two contract net protocol (CNP)-based algorithms for comparison, respectively SingleItem-CNP-based (SI-CNP) algorithm and Batch-CNP-based (BA-CNP) algorithm. Computational experiments indicated that both the total profit and percentage of scheduled urgent tasks achieved by the m-ACBBA and m-CBBA algorithms were much higher than those achieved by both SI-CNP and BA-CNP. Additionally, the number of communications needed by either m-ACBBA or m-CBBA algorithm was lower than that by SI-CNP. When the communication cost in the system is high, the m-CBBA algorithm is preferred because it balances the profit and the required number of communications. When the communication cost is low, the m-ACBBA algorithm is preferred because it achieves high total profit and high percentage of scheduled urgent tasks.     

导航卫星中断频次的分析方法

面向导航卫星中断频次的定量分析需求,该文分析导航卫星中断产生的主要原因,给出中断频次分析的流程,并针对分析过程中的三个关键问题,研究提出具体实施方法,包括通过相关性分析快速定位底层中断事件,通过中断树建立指标分析模型,并融合在轨数据、地面试验数据快速预估得到底层功能异常率等。最后通过示例进一步说明中断频次分析过程。该文方法已应用于北斗导航卫星工程。     

An empirical solar radiation pressure model for satellites moving in the orbit-normal mode

We present a family of empirical solar radiation pressure (SRP) models suited for satellites orbiting the Earth in the orbit normal (ON) mode. The proposed ECOM-TB model describes the SRP accelerations in the so-called terminator coordinate system. The choice of the coordinate system and the SRP parametrization is based on theoretical assumptions and on simulation results with a QZS-1-like box-wing model, where the SRP accelerations acting on the solar panels and on the box are assessed separately. The new SRP model takes into account that in ON-mode the incident angle of the solar radiation on the solar panels is not constant like in the yaw-steering (YS) attitude mode. It depends on the elevation angle of the Sun above the satellite’s orbital plane. The resulting SRP vector acts, therefore, not only in the Sun-satellite direction, but has also a component normal to it. Both components are changing as a function of the incident angle. ECOM-TB has been used for precise orbit determination (POD) for QZS-1 and BeiDou2 (BDS2) satellites in medium (MEO) and inclined geosynchronous Earth orbits (IGSO) based on IGS MGEX data from 2014 and 2015. The resulting orbits have been validated with SLR, long-arc orbit fits, orbit misclosures, and by the satellite clock corrections based on the orbits. The validation results confirm that—compared to ECOM2—ECOM-TB significantly (factor 3–4) improves the POD of QZS-1 in ON-mode for orbits with different arc lengths (one, three, and five days). Moderate orbit improvements are achieved for BDS2 MEO satellites—especially if ECOM-TB is supported by pseudo-stochastic pulses (the model is then called ECOM-TBP). For BDS2 IGSOs, ECOM-TB with its 9 SRP parameters appears to be over-parameterized. For use with BDS2 IGSO spacecraft we therefore developed a minimized model version called ECOM-TBMP, which is based on the same axis decomposition as ECOM-TB, but has only 2 SRP parameters and is supported by pseudo-stochastic parameters, as well. This model shows a similar performance as ECOM-TB with short arcs, but an improved performance with (3-day) long-arcs. The new SRP models have been activated in CODE’s IGS MGEX solution in Summer 2018. Like the other ECOM models the ECOM-TB derivatives might be used together with an a priori model.