Strategies for in-orbit calibration of drag-free control systems

Drag-Free Satellites (DFS) are a class of scientific satellite missions designed for research on fundamental physics as well as geodesy. They consist, basically, of a small inner satellite (test mass) located in a cavity inside a larger satellite, the normal one. The Drag-Free Attitude Control System (DFACS) is the most complex technology on-board these satellites. This key technology allows the residual accelerations on experiments on board the satellites to be significantly reduced. In order to achieve this very low disturbance environment (for some missions <10⁻¹⁴ g) the drag-free control system has to be optimized. This optimization process is required because of uncertainties in system parameters that demand a robustness of the control system. This paper will present approaches for in-orbit calibration of drag-free control systems. The discussion includes modeling, with scale factors and cross couplings, possible excitation signals, comparison of different parameter identification/estimation methods as well as simulation results.     

GNSS广域差分改正数估计方法研究

随着我国北斗三号基本系统的正式运行,基于地面监测站的广域差分增强系统成为进一步提升卫星导航定位精度的手段之一。在码噪声多径误差修正(CNMC)的基础上,使用等效钟差方法实现GNSS卫星轨道与钟差误差的解耦,并依据卫星轨道运动的动力学特性,引入希尔差分方程描述卫星轨道误差变化,实现对轨道误差的实时卡尔曼滤波估计。基于GPS实测数据,对改正前后的用户等效测距误差(UERE)和定位精度进行了对比研究。实验结果表明,采用该方法,UERE标准差由改正前的0.456 m减小至0.227 m,降幅达到50.22%;整体水平定位误差(95%置信区间)由0.981 m减小至0.782 m,垂直定位误差(95%置信区间)由1.991 m减小至1.131 m,分别提升了20.29%和43.19%,差分改正效果明显。     

一种频域插值变化迭代频率估计方法

针对现有频率估计算法存在的复杂度高、频率估计能力弱、估计结果均方差大等缺点,在固定迭代AM（Aboutanios—Mulgrew）无偏频率估计算法基础上,提出一种频域插值变化迭代频率估计算法,推导了不同迭代参数实现无偏估计的充分条件,证明了有偏估计时本算法的收敛性和偏离度,通过设置不同迭代参数,可以实现无偏或有偏估计。仿真分析表明：当具有较高信噪比时,在整个频率估计范围内,该方法均方误差接近CRLB（Cramer-RaoLowerBound,克拉美一罗下限）;当FFT（FastFourierTransform,快速傅里叶变换）粗估计残余频率接近0．5时,该方法的均方误差优于CRLB,为CRLB的96％。     

A novel multistage estimation of signal parameters

A multistage estimation scheme is presented for estimating the parameters of a received carrier signal possibly phase-modulated by unknown data and experiencing very high Doppler, Doppler rate, etc. Such a situation arises, for example, in the case of the Global Positioning Systems (GPS). In the proposed scheme, the first-stage estimator operates as a coarse estimator of the frequency and its derivatives, resulting in higher RMS estimation errors but with a relatively small probability of the frequency estimation error exceeding one-half of the sampling frequency (an event termed cycle slip). The second stage of the estimator operates on the error signal available from the first stage, refining the overall estimates, and in the process also reduces the number of cycle slips. The first-stage algorithm is a modified least-squares algorithm operating on the differential signal model and referred to as differential least squares (DLS). The second-stage algorithm is an extended Kalman filter, which yields the estimate of the phase as well as refining the frequency estimate. A major advantage of the proposed algorithm is a reduction in the threshold for the received carrier power-to-noise power spectral density ratio (CNR) as compared with the threshold achievable by either of the algorithms alone     

附加双差模糊度网解约束的GNSS相位延迟估计

针对现有相位延迟估计方法未利用双差模糊度网解约束的特点,提出了一种附加双差模糊度网解约束的相位延迟估计方法。首先以整周双差模糊度网解为真值更新精密单点定位(PPP)估计的非差浮点模糊度,再通过网平差法估计未校正相位延迟(UPD)。实验结果表明,PPP双差模糊度与双差模糊度网解一致性好,其宽、窄巷中误差分别为0.07周和0.11周,误差大于1周的PPP双差模糊度主要出现在卫星初升阶段。约束改变了非差浮点模糊度,从而改变了参与UPD网平差的测站,使14%新升起卫星的窄巷UPD在短时间内较无约束有大于0.2周的差异。附加约束/无约束的全天星间单差UPD差异在大于99.9%的置信水平下满足零均值假设,表明整周双差模糊度网解约束方法与无约束方法估计UPD相关产品具有等效性。研究结果可为参考网非差模糊度解算和PPP与网络实时动态定位(RTK)的融合工作提供参考。     

基于逆系统方法的多输入多输出非高斯驱动信号生成

多输入多输出(MIMO)振动环境试验是振动试验的新方法,该方法比传统的单输入单输出(SISO)方法更真实地模拟了结构在实际工作下的振动环境。传统随机振动试验的频域方法所模拟的是平稳高斯信号,但是实际振动环境往往是超高斯的,且超高斯激励和高斯激励对结构的损害差异很大,因此对MIMO超高斯振动环境试验的研究十分重要。驱动信号的产生是振动环境试验的关键。提出一种基于逆系统的MIMO超高斯驱动信号生成方法。首先,根据给定的响应参考功率谱密度(PSD)和参考峭度值采用相位调节法生成参考响应超高斯伪随机信号;其次,利用逆系统方法生成超高斯伪随机驱动信号;最后,对伪随机驱动信号进行时域随机化处理得到超高斯真随机驱动信号。对一根两输入两输出的悬臂梁模型进行了仿真验证,结果表明,采用本文方法所得到的驱动信号加载在梁上后,其输出响应谱与给定参考谱之间的误差完全满足工程中常用的±3 dB的要求,而且峭度也与给定的参考值十分接近,证实了所提方法的有效性。     

HOI延迟对LEO卫星简化动力学POD的影响

通常使用无电离层(IF)线性组合(LC)消除低地球轨道(LEO)卫星简化动力学精密定轨(POD)一阶电离层延迟误差,忽略了高阶电离层(HOI)延迟误差。随着LEO卫星POD技术的发展,计算不同轨道高度的HOI延迟并探索其变化已成为进一步提高POD精度的重要手段。首先,使用国际参考电离层-2016(IRI-2016)和国际地磁参考场第13代(IGRF-13)模型,计算电离层穿刺点(IPP)位置和地磁场强度。其次,使用平滑星载GNSS数据计算电离层斜路径总电子含量(STEC)。然后,分别计算GOCE、GRACE-A和SWARM-A/B卫星的二阶和三阶电离层延迟。最后,评估了HOI延迟对LEO卫星重叠轨道分析、卫星激光测距(SLR)检核和精密科学轨道(PSO)比较结果的影响。实验结果表明：HOI延迟对LEO卫星简化动力学POD的影响大约在毫米至厘米的数量级上;HOI延迟对LEO卫星简化动力学POD外符合精度的影响分别达到0.92,0.22,0.21和0.18 mm;随着LEO卫星轨道高度的增加,HOI延迟对LEO卫星简化动力学POD的影响减小。     

一种基于轨道根数约束的最优制导方法

针对航天器空间变轨的制导问题,研究了一种基于轨道根数约束的最优制导方法。在地心惯性坐标系下直接建立航天器的最优控制模型,给出了位置速度表达式和协态变量初值之间的关系;进一步,在不约束真近点角的前提下,推导了5个轨道根数的约束方程,并通过对轨道根数和最优控制理论中协态方程的特性分析,获得了另外两个约束方程。协态变量初值可直接通过求解7个完整约束方程组获得,进而得到最优推力方向。仿真验证了所提制导方法的有效性。     

编队飞行小卫星空间状态测量方法研究

介绍了编队飞行小卫星特点，提出了小卫星空间状态测量方案。着重论述了星间相对状态自主测量方法，计算了伪码和载波相位跟踪环路、时钟基准频差以及测量设备延迟带来的星间相对距离测量误差，讨论了降低误差的方法。在对误差项进行综合分析后，得出相对距离测量误差可达到厘米级的结论。     

Efficient Utilization of Orbit/Frequency for Satellite Broadcasting

In this paper some basic concepts involved in efficient orbit/spectrum utilization for broadcast satellites are discussed. A theoretical analysis leads to the determination of interference zones and minimum number of frequency channels required for homogeneous systems operating in the 800-MHz and 12-GHz bands. A crossed-beam arrangement of satellite systems is proposed, in which the beams from satellites to service areas cross each other. This arrangement has a more efficient orbit spectrum utilization compared to the identical-longitude approach, in which the subsatellite point has the same longitude as the center of the service area, Based upon practical considerations, the elemental number of frequency channels required for world-wide coverage has been calculated in relation to satellite orbital position, for single-beam and multibeam satellite systems, using the identical-longitude and crossed -beam arrangements.     

不同基准钟对BDS实时钟差估计的影响分析CSCD

BDS原始观测量是站星之间的相对时间延迟,在实时卫星钟差估计过程中需要引入一个基准钟,求解该基准约束下的钟差产品。基于两种不同的基准约束条件实时估计BDS卫星钟差,并从实时钟差的估计精度、钟差频率特性(频率漂移率、频率准确度、频率稳定度)及钟差预报精度等方面分析其对BDS实时钟差估计的影响。算例结果表明,在两种不同的基准模式下,估计得到的BDS实时卫星钟差性能基本一致,在实际使用中可根据情况采用不同的基准钟进行钟差估计,为BDS实时卫星钟差估计时基准钟的选择提供了参考。     

一种稳健的利用分布式小卫星获取宽域、高分辨SAR图像的方法

 传统星载合成孔径雷达(SAR)受最小天线面积条件的限制,其横向分辨率与测绘带宽度相矛盾,采用编队飞行的分布式小卫星SAR体制,可解决这一矛盾,即能够同时获得大测绘带、高分辨率的地面场景的SAR图像。主要研究如何在误差环境下抑制由小天线引入的多普勒模糊的信号处理方法。该方法利用子空间投影技术对导向矢量进行精确估计,在小卫星存在基线及其他误差时仍然能够稳健地恢复高分辨的二维SAR图像。     

Adaptive robust cubature Kalman filtering for satellite attitude estimation

This paper is concerned with the adaptive robust cubature Kalman filtering problem for the case that the dynamics model error and the measurement model error exist simultaneously in the satellite attitude estimation system. By using Hubel-based robust filtering methodology to correct the measurement covariance formulation of cubature Kalman filter, the proposed filtering algorithm could effectively suppress the measurement model error. To further enhance this effect and reduce the impact of the dynamics model error, two different adaptively robust filtering algorithms, one with the optimal adaptive factor based on the estimated covariance matrix of the predicted residuals and the other with multiple fading factors based on strong tracking algorithm, are developed and applied for the satellite attitude estimation. The quaternion is employed to represent the global attitude parameter, and three-dimensional generalized Rodrigues parameters are introduced to define the local attitude error. A multiplicative quaternion error is derived from the local attitude error to maintain quaternion normalization constraint in the filter. Simulation results indicate that the proposed novel algorithm could exhibit higher accuracy and faster convergence compared with the multiplicative extended Kalman filter, the unscented quaternion estimator, and the adaptive robust unscented Kalman filter.     

Minimum Number of Satellites for Three-Dimensional Continuous Worldwide Coverage

Global positioning by means of satellites requires simultaneous observation by at least four satellites. The problem is to determine the minimum number of satellites and the corresponding orbital geometry necessary to satisfy this requirement on a continuous basis. To model the problem, a fixed number of users are assumed uniformly distributed in a known manner over the surface of the earth, and the satellites are restricted to exist in either three or four orbital planes. However, the orbit radius and inclination angle are left as variables. Under these assumptions, and starting with a small number of satellites which will be increased afterwards, an algorithm is developed to determine the visibility of satellites at each surface location. In this way it is possible to specify the minimum number of satellites needed by any desired orbital geometry. It is found that the number of satellites required for three-dimensional continuous worldwide coverage decreases as the orbit radius is increased. There appears to be no general trend regarding the effect of the inclination angle on the minimum number of satellites.     

低照度小样本限制下的失效卫星相对位姿估计与优化

低照度图像信息受损严重,会导致失效卫星的位姿估计精度和鲁棒性降低。基于此,提出了无监督生成式对抗网络低照度图像增强模型。生成器以U-Net网络为基础,并设计密集残差连接结构。判别器设计为全局-局部的双辨别器结构,由传统的单一标量扩展为多标量判别。在小样本的条件下,基于进化训练与并行训练方式改进基于SinGAN的数据增广方法。最后,在基于ORB-SLAM位姿初始化的基础上,建立特征信息的局部地图,克服位姿估计对参考帧的依赖;通过关键帧ROI的稠密匹配,建立关于平面法向量和单目相机安装高度的非线性优化模型求解尺度因子;通过闭环检测后的相似性变换,构建关键帧集合的联合位姿图优化方程,实现对位姿矩阵的全局校正。实验结果表明：测量稳定后,低照度图像的姿态角误差最大值为4°,而图像增强后的姿态角误差最大为0.5°;对于以角速度20°/s运动的失效卫星旋转5周,相对静止下的跟踪测量为5周,1 m水平方向机动下的跟踪测量为4.5周。可以满足失效卫星相对姿态测量的任务需求。     

非理想均匀圆阵的波达方向和多普勒频率联合估计

针对均匀圆阵存在一般阵列误差 (如阵元的幅相误差和安装位置误差等 )的情况 ,提出了多个信号的波达方向和多普勒频率估计方法。直接利用均匀圆阵的阵列流形 ,采用波达矩阵法估计各个信号的多普勒频率。由一般阵列误差的统计特性构造加权矩阵 ,采用加权总体最小二乘法估计各个信号的波达方向。此方法具有鲁棒性强等特点。计算机仿真证明了此方法的有效性     

A Weighted Inner Product Estimator in the Geometric Algebra of Euclidean 3-Space for Source Localization Using an EM Vector-sensor

In this paper, the source localization by utilizing the measurements of a single electromagnetic (EM) vector-sensor is investigated in the framework of the geometric algebra of Euclidean 3-space. In order to describe the orthogonality among the electric and magnetic measurements, two multivectors of the geometric algebra of Euclidean 3-space (G3) are used to model the outputs of a spatially collocated EM vector-sensor. Two estimators for the wave propagation vector estimation are then formulated by the inner product between a vector and a bivector in the G3. Since the information used by the two estimators is different, a weighted inner product estimator is then proposed to fuse the two estimators together in the sense of the minimum mean square error (MMSE). Analytical results show that the statistical performances of the weighted inner product estimator are always better than its traditional cross product counterpart. The efficacy of the weighted inner product estimator and the correctness of the analytical predictions are demonstrated by simulation results.     

基于相似变换的涡扇发动机智能动态推力估计方法

智能推力估计面临飞行包线大、工作状态多变带来的数据采集和处理问题,获得的训练数据难以覆盖整个飞行包线的各种过渡工作状态,为此本文提出一种基于相似变换的推力估计数据处理方法。通过机理分析选择推力估计器输入,以相似变换对推力估计的输入和输出数据进行处理,并设计了基于输入延迟的深层动态神经网络来实现动态推力估计。非训练数据区域的动态仿真结果表明,相似变换后,深层动态神经网络的最大推力估计误差降低了62.20%,平均误差降低了43.50%;未进行相似变换时,相比深层静态神经网络,深层动态神经网络的最大推力估计误差降低了43.42%,平均误差降低了2.35%,仿真结果表明了本文所提出的数据处理方法和动态推力估计结构有效性。     

Identification of reduced-order model for an aeroelastic system from flutter test data

Recently, flutter active control using linear parameter varying (LPV) framework has attracted a lot of attention. LPV control synthesis usually generates controllers that are at least of the same order as the aeroelastic models. Therefore, the reduced-order model is required by syn-thesis for avoidance of large computation cost and high-order controller. This paper proposes a new procedure for generation of accurate reduced-order linear time-invariant (LTI) models by using sys-tem identification from flutter testing data. The proposed approach is in two steps. The well-known poly-reference least squares complex frequency (p-LSCF) algorithm is firstly employed for modal parameter identification from frequency response measurement. After parameter identification, the dominant physical modes are determined by clear stabilization diagrams and clustering tech-nique. In the second step, with prior knowledge of physical poles, the improved frequency-domain maximum likelihood (ML) estimator is presented for building accurate reduced-order model. Before ML estimation, an improved subspace identification considering the poles constraint is also proposed for initializing the iterative procedure. Finally, the performance of the proposed procedure is validated by real flight flutter test data.     

导航辅助下基于LAS码的天基多目标测控

中继卫星的天基测控是解决中低轨道航天器大范围、长时间、多目标测控的有效途径,但用户星到中继卫星的时延误差、多普勒频移校正误差等参数会对多目标测控带来干扰。针对该问题,首先分析了定位误差引起的时延误差、多普勒频移校正误差等因素,设计了天基多目标前向链路遥控SMA(S-band Multiple Access,S频段多址)信号形式和反向遥测链路SMA信号形式,选择与时延校正误差相匹配的LAS(Large Area Synchronous,大区域同步)码作为反向遥测信息的扩频码,建立了导航辅助的终端时延和频率预校正方案模型,可以有效消除多用户干扰。仿真表明,当Eb/N0≥10.5dB时,前向遥控信息误码率pe≤1×10~(-6),反向遥测信息误码率pe≤1×10~(-6),使用LAS码比Gold序列约有2dB性能改善,为基于我国天链卫星的中低轨道卫星稳定运行、载人航天交会对接以及后续空间站建设等任务的测控提供重要参考。