GNSS clock corrections densification at SHAO: from 5 min to 30 s

High frequency multi-GNSS zero-difference applications like Precise Orbit Determination(POD)for Low Earth Orbiters(LEO)and high frequency kinematic positioning require corresponding high-rate GNSS clock corrections.The determination of the GNSS clocks in the orbit determination process is time consuming,especially in the combined GPS/GLONASS processing.At present,a large number of IGS Analysis Centers(AC)provide clock corrections in 5-min sampling and only a few ACs provide clocks in 30-s sampling for both GPS and GLONASS.In this paper,an efficient epoch-difference GNSS clock determination algorithm is adopted based on the algorithm used by the Center for Orbit Determination in Europe(CODE).The clock determination procedure of the GNSS Analysis Center at Shanghai Astronomical Observatory(SHAO)and the algorithm is described in detail.It is shown that the approach greatly speeds up the processing,and the densified 30-s clocks have the same quality as the 5-min clocks estimated based on a zero-difference solution.Comparing the densified 30-s GNSS clocks provided by SHAO with that of IGS and its ACs,results show that our 30-s GNSS clocks are of the same quality as that of IGS.Allan deviation also gives the same conclusion.Further validation of the SHAO 30-s clock product is performed in kinematic PPP and LEO POD.Results indicate that the positions have the same accuracy when using SHAO 30-s GNSS clocks or IGS(and its AC)finals.The robustness of the algorithm and processing approach ensure its extension to provide clocks in 5-s or even higher frequencies.The implementation of the new approach is simple and it could be delivered as a black-box to the current scientific software packages.     

基于伪相关函数的多级电平编码符号信号通用无模糊跟踪方法

针对新一代全球导航卫星系统(GNSS)中多级电平编码符号(MCS)信号存在的跟踪模糊问题,本文提出了一种通用的MCS信号无模糊跟踪方法.首先推导了不同MCS信号互相关函数的统一表达式,并给出了伪相关函数的定义;然后深入分析了实现无模糊跟踪需要满足的约束条件,推导了两路参考信号的通用构造方法以及相互之间的关系,为具体MCS信号的求解提供了极大的便利;进而给出了利用本文方法的GNSS接收机码跟踪环路模型.作为MCS信号的特例,分别讨论了本文方法在四种二进制偏移载波信号跟踪中的应用.仿真结果表明,本文方法能够有效解决MCS信号的跟踪模糊问题,具有良好的性能和广阔的应用前景.     

基于软件接收机的卫星导航信号模拟器性能指标测试

针对卫星导航信号模拟器的性能指标测试评估问题,提出了基于软件接收机的性能指标测试方法;通过分析模拟器性能,包括对动态性能、相位噪声、通道间时延一致性和静态定位精度的综合分析,给出了相应的测试方法和测试流程,并结合卫星导航信号模拟器产生的实际信号对各项性能指标的测试结果进行分析;经测试结果表明,软件接收机可以有效地对卫星导航信号模拟器各项性能指标进行测试,测试结果符合理论预期,满足模拟器性能指标测试评估的需求。     

一种新的卫星钟差Kalman滤波噪声协方差估计方法

采用Kalman滤波方法进行钟差参数计算和预报时, 需确定Kalman滤波噪声协方差矩阵. 针对这一问题, 提出了一种新的卫星钟差Kalman滤波噪声协方差估计方法, 通过建立新息的相关函数序列与未知的噪声参数间的线性函数模型, 采用最小二乘法进行噪声参数估计. 采用精密钟差数据进行钟差参数估计和预报分析, 结果表明, 该方法具有较好的收敛性, 并与顾及随机噪声模型的开窗分类因子自适应抗差估计方法进行对比分析, 验证了新方法的正确性和有效性.     

一种新型GNSS/MEMS紧组合导航系统设计

针对探空火箭在高动态与非线性环境下高精度导航问题,基于DSP6747+FPGA双核构架,设计了一种低成本的新型GNSS/MEMS一体化紧组合导航系统。同时,在直接法滤波的基础上采用了改进的自适应渐消扩展卡尔曼粒子滤波(PF-AFEKF)算法,解决了粒子退化问题并充分利用了当前观测信息以及自适应权重因子抑制了滤波发散,实现了高动态、非线性状态下的GNSS/MEMS深度信息融合。最后,进行了地面跑车试验,试验结果证明了所设计的低成本GNSS/MEMS紧组合导航系统具有较高的导航精度,容错能力强。     

Image quality assessment in multiband DCT domain based on SSIM

The performance of image quality assessment method based on SSIM (structural similarity) is better than the PSNR (peak signal to noise ratio), but the assessment effects of SSIM is poor for seriously blurred image, therefore, the model that combined HVS (human visual sensitivity) and SSIM was established. The basic idea is based on the human eye's sensitivity to different frequency distortion image, the image is two-dimensional discrete cosine transform frequency component into low, mid, high-frequency component, to obtain the frequency component of light, contrast and structural information, using Pearson coefficient for weight and sum processing to the sub-image according to frequency bands of different sensitive degree, finally, get the sharpness of the image. Through nonlinear regression analysis of objective assessment and DMOS, experiments showed that this method was closer to human perception than SSIM and GSSIM for serious blurred distortion image. At the same time, compared to conventional algorithm MAE (mean absolute error), MSE (mean square error) and PSNR, this model was more consistent with human visual characteristics.     

总能系统全工况特性数学模拟的研究

总能系统分析研究的主要方法是数学模拟：数学模型和模拟分析。随着热力系统复杂化、大型化，传统的方法面临着新的挑战，而基于通用性组合思想的模块化建模成为当前系统数学模拟发展的新方向[1,4]。但模块化建模的研究仅在化工部门过程系统方面进展较快，而对总能系统来说，还刚刚起步。本文研究总能系统（包括核能发电系统）全工况特性一维稳态模型的建模和模拟分析的新构思。1总能系统模块化建模模块化建模利用预先建立的典型单元模块模型，通过一定手段实现一定功能的组合，从而构造出拥有大量部件的系统模型。其主要内容有：大系统分解…     

流体力学拉氏守恒滑移线算法设计

本文在Wilkins滑移线算法的基础上, 设计了从区网格边人工黏性对主点施加算法. 结合拉氏“相容性格式”, 通过定义“接触力”和“接触力做功”, 利用局部修正内能的方法设计了总能量守恒的流体力学拉氏格式. 该修正方法可保证对称性、守恒性; 提高数值模拟分辨率.     

清华物理系量子信息研究的新进展

在清华大学物理系成立60周年之际,我们对近年来清华大学物理系量子信息研究的主要进展情况作一介绍,包括量子搜索算法研究,核磁共振量子计算的实验研究,量子通讯的理论与实验研究.在量子搜索算法研究方面,我们提出了量子搜索算法的相位匹配,纠正了当时的一种错误观点,并且提出了一种成功率为100%的量子搜索算法,改进了Grover算法;在核磁共振量子计算实验方面,我们实现了2到7个量子比特的多种量子算法的实验演示;在量子通讯方面,我们提出了分布式传输的量子通讯的思想,应用于量子密钥分配、量子秘密共享、量子直接安全通讯等方面,构造了多个量子通讯的理论方案.在实验室,我们实现了2米距离的空间量子密码通讯的演示实验.     

Computational leakage: Grover's algorithm with imperfections

We study the effects of dissipation or leakage on the time evolution of Grover's algorithm for a quantum computer. We introduce an effective two-level model with dissipation and randomness (imperfections), which is based upon the idea that ideal Grover's algorithm operates in a 2-dimensional Hilbert space. The simulation results of this model and Grover's algorithm with imperfections are compared, and it is found that they are in good agreement for appropriately tuned parameters. It turns out that the main features of Grover's algorithm with imperfections can be understood in terms of two basic mechanisms, namely, a diffusion of probability density into the full Hilbert space and a stochastic rotation within the original 2-dimensional Hilbert space. Received 12 August 2002 / Received in final form 14 October 2002 Published online 4 February 2003     

Global modeling 2~(nd)-order ionospheric delay and its effects on GNSS precise positioning

Ionospheric delay is one of the major error sources in GNSS navigation and positioning.Nowadays,the dual-frequency technique is the most widely used in ionospheric refraction correction.However,dual-frequency measurements can only eliminate the first-order term of ionospheric delay,while the effect of the second-order term on GNSS observations may be several centimeters.In this paper,two models,the International Reference Ionosphere (IRI) 2007 and International Geomagnetic Reference Field (IGRF) 11 are used...     

Bayesian fusion of GNSS,ITS-G5 and IR–UWB data for robust cooperative vehicular localization

In the automotive domain, Cooperative Localization (CLoc) is a new promising paradigm that aims at outperforming conventional Global Navigation Satellite Systems (GNSS) in terms of positioning accuracy, robustness, and service continuity, by relying on Vehicle-to-Vehicle (V2V) communications and hybrid data fusion. However, the growing number and the variety of the sensors aboard vehicles raise unprecedented challenges, especially in the context of distributed fusion approaches. This paper thus compares parametric and nonparametric Bayesian data fusion engines (e.g., based on cooperative variants of the Extended Kalman Filter (EKF) and Particle Filter (PF), respectively), while validating a CLoc scheme suitable to Vehicular Ad Hoc Networks (VANETs). More particularly, absolute position information from both onboard GNSS receiver and ITS-G5 V2V messages, as well as relative distance measurements based on the Impulse Radio–Ultra-Wideband (IR–UWB) technology, are combined into a single location solution that is hopefully more robust and more accurate than that of standalone GNSS. First, we investigate V2V ranging accuracy on a highway under real mobility conditions. In the same environment, we then provide offline validations of CLoc positioning, confirming significant performance gains through cooperation over conventional GNSS, even in case of poor initialization. In this specific context, the PF solution is thus shown to yield even better accuracy in comparison with EKF, thanks to its fine robustness against faced non-linear dynamics and non-Gaussian noise processes. Finally, we illustrate the resilience of the proposed solution under temporary GNSS denial.     

Generalized Correntropy Criterion-Based Performance Assessment for Non-Gaussian Stochastic Systems

Control loop performance assessment (CPA) is essential in the operation of industrial systems. In this paper, the shortcomings of existing performance assessment methods and indicators are summarized firstly, and a novel evaluation method based on generalized correntropy criterion (GCC) is proposed to evaluate the performance of non-Gaussian stochastic systems. This criterion could characterize the statistical properties of non-Gaussian random variables more fully, so it can be directly used as the assessment index. When the expected output of the given system is unknown, generalized correntropy is used to describe the similarity of two random variables in the joint space neighborhood controlled and take it as the criterion function of the identification algorithms. To estimate the performance benchmark more quickly and accurately, a hybrid-EDA (H-EDA) combined with the idea of “wading across the stream algorithm” is proposed to obtain the system parameters and disturbance noise PDF. Through the simulation of a single loop feedback control system under different noise disturbances, the effectiveness of the improved algorithm and new indexes are verified.     

一种新的适用于水声信道的常模类盲均衡算法

为了缓解常数模盲均衡算法收敛速度缓慢、稳态误差大的问题,考虑到初始化权值对CMA算法的重要影响,本文利用遗传算法对CMA算法进行了有效改进,引入了小样本重用的思想,给出了一种新的适用于水声信道的常模类盲均衡算法。计算机仿真研究证明,该算法不仅大大加快了收敛速度,而且有效地降低了稳态误差。     

基于加权实数编码遗传算法的超材料优化设计

针对超材料优化设计参数多、权重各异的特点,提出一种基于加权实数编码的遗传算法,该算法采用等位基因或双倍基因实现对遗传基因的加权编码,是二进制编码中加权思想的推广,与普通遗传算法相比,加权实数编码遗传算法加入了人工选择因素,既可以加快种群收敛速度,又可以提高算法的求解质量,尤其适用于种群规模较大的遗传优化计算,文中以一个超材料吸波体的优化设计为例,对算法进行了验证。     

Impact of GNSS radio occultation bending angle data assimilation in YH4DVAR system

Global Navigation Satellite System （GNSS） radio occultation measurements have been assimilated into the four- dimensional variational data assimilation system （YH4DVAR） using a one-dimensional bending angle operator （GBAO） as a new type of observation. For the sake of verifying the impact of GNSS radio occultation （RO） measurements to the data assimilation system, three experiments have been conducted. The statistical results of the analysis error experiment and forecast skill experiment show that the GNSS RO measurements have an impact on the analysis system. The typhoon forecast experiment shows the impact on the important weather process. They all have a positive impact on the weather forecast. Lastly, we look forward to future work on the observation system simulation experiment （OSSE） to investigate the impact of GNSS RO measurements as a function of observation number, which is an effective method to estimate the saturation of the observation number.     

Variational Quantum Algorithms for the Steady States of Open Quantum Systems

The solutions of the problems related to open quantum systems have attracted considerable interest.We propose a variational quantum algorithm to find the steady state of open quantum systems.In this algorithm,we employ parameterized quantum circuits to prepare the purification of the steady state and define the cost function based on the Lindblad master equation,which can be efficiently evaluated with quantum circuits.We then optimize the parameters of the quantum circuit to find the steady state.Numerical simulations are performed on the one-dimensional transverse field Ising model with dissipative channels.The result shows that the fidelity between the optimal mixed state and the true steady state is over 99%.This algorithm is derived from the natural idea of expressing mixed states with purification and it provides a reference for the study of open quantum systems.     

一类基于ENO插值的守恒重映算法

在大变形流体力学问题的数值模拟中,常常会涉及到计算网格的重分.基于不同网格的物理量传递便是所谓的重映技术.基于ENO插值的思想,发展了一类适用于任意网格的ENO守恒重映算法,并给出了数值结果.     

一种气象观测数据求导的新方法

气象观测数据不可避免地带有观测误差,而对带有误差的离散数据求导是一个不适定的反问题.为了解决这一长期困扰气象工作者的问题,利用吉洪诺夫正则化思想,提出了矩形区域内观测数据的二维一阶偏导数重构算法.通过系列模拟观测数据试验检验了该算法的性能,表明该算法有效且计算精度高.另外,应用该算法对气象观测资料进行客观分析,结果表明：该算法作为一种新的客观分析方法是可行的,并且还能提高小尺度天气系统的识别能力. 关键词： 数值微分 吉洪诺夫正则化 客观分析     

Parallel re-initialization of level set functions on distributed unstructured tetrahedral grids

Level set functions are employed to track interfaces in various application areas including simulation of two-phase flows and image segmentation. Often, a re-initializing algorithm is incorporated to transform a numerically instable level set function to a signed distance function. In this note, we present a parallel algorithm for re-initializing level set functions on unstructured, three-dimensional tetrahedral grids. The main idea behind this new domain decomposition approach is to combine a parallel brute-force re-initializing algorithm with an efficient way to compute distances between the interface and grid points. Time complexity and error analysis of the algorithm are investigated. Detailed numerical experiments demonstrate the accuracy and scalability on up to 128 processes.