扩展卡尔曼滤波和不敏卡尔曼滤波在实时雷达回波反演大气波导中的应用

为了改善雷达回波反演大气波导(RFC)方面存在的单时次、单方位角反演的问题,提出利用扩展卡尔曼滤波和不敏卡尔曼滤波的反演算法对大气波导结构的多方位角实时跟踪反演. 在卡尔曼滤波方法中分别给出大气波导结构的参数化方程、观测方程、滤波算法的状态转移方程,最后导出滤波反演算法的迭代求解流程. 在大气波导结构不随时间变化和随时间变化的两种条件下,对扩展卡尔曼滤波和不敏卡尔曼滤波算法进行数值实验. 实验结果表明,不敏卡尔曼滤波更适用于RFC这高度非线性反演问题,它可能今后为大气波导结构多方位角实时跟踪反演的业务化运行提供理论基础与技术保证. 关键词： 大气波导 雷达回波 扩展卡尔曼滤波 不敏卡尔曼滤波     

雷达回波资料反演海洋波导中观测值权重的确定

提出采用一种与电磁波传播损耗相类似的指数衰减形式来描绘观测值权重随距离的变化,结合标准的统计反演算法来测试不同的观测值权值对反演速度与精度的影响.结果表明,指数衰减形式能较好地描绘观测值权值随距离的变化情况;其中衰减指数的选取是关键.合理地选取衰减指数可较好地改进反演的速度与精度.对于不同频率的电磁波,最优衰减指数是各不相同的. 关键词： 海洋波导 观测值权重 衰减指数 反演算法     

Monitoring of the ducting by a ground-based GPS receiver

In this paper, we describe the estimation of low-altitude refractivity structure from simulation and real ground-based GPS delays. The vertical structure of the refractive environment is modeled using three parameters, i.e., duct height, duct thickness, and duct slope. The refractivity model is implemented with a priori constraints on the duct height, thickness, and strength, which might be derived from soundings or numerical weather-prediction models. A ray propagation model maps the refractivity structure into a replica field. Replica fields are compared with the simulation observed data using a squared-error objective function. A global search for the three environmental parameters is performed using genetic algorithm. The inversion is assessed by comparing the refractivity profiles from the radiosondes to those estimated. This technique could provide near-real-time estimation of ducting effect. The results suggest that ground-based GPS provides significant atmospheric refractivity information, despite certain fundamental limitations of ground-based measurements. Radiosondes typically are launched just a few times daily. Consequently, estimates of temporally and spatially varying refractivity that assimilate GPS delays could substantially improve over-estimates using radiosonde data alone.     

Refractivity estimation from radar sea clutter

This paper addresses the problem of estimating lower atmospheric refractivity under the nonstandard propagation conditions frequently encountered in low altitude maritime radar applications. The vertical structure of the refractive environment is modeled by using a five-parameter model, and the horizontal structure is modeled as range-independent. The electromagnetic propagation in the troposphere is simulated by using a split-step fast Fourier transform based on parabolic approximation to the wave equation. A global search marked as a modified genetic algorithm (MGA) for the 5 environmental parameters is performed by using a genetic algorithm (GA) integrated with a simulated annealing technique. The retrieved results from simulated runs demonstrate the ability of this method to make atmospheric refractivity estimations. A comparison with the classical GA and the Bayesian Markov Chain Monte Carlo (Bayesian-MCMC) technique shows that the MGA can not only shorten the inverse time but also improve the inverse precision. For real data cases, the inversion values do not match the reference data very well. The inverted profile, however, can be used to synoptically describe the real refractive structure.     

Inversion for atmosphere duct parameters using real radar sea clutter

This paper addresses the problem of estimating the lower atmospheric refractivity (M profile) under nonstandard propagation conditions frequently encountered in low altitude maritime radar applications. The vertical structure of the refractive environment is modeled using five parameters and the horizontal structure is modeled using five parameters. The refractivity model is implemented with and without a priori constraint on the duct strength as might be derived from soundings or numerical weather-prediction models. An electromagnetic propagation model maps the refractivity structure into a replica field. Replica fields are compared with the observed clutter using a squared-error objective function. A global search for the 10 environmental parameters is performed using genetic algorithms. The inversion algorithm is implemented on the basis of S-band radar sea-clutter data from Wallops Island, Virginia (SPANDAR). Reference data are from range-dependent refractivity profiles obtained with a helicopter. The inversion is assessed (ⅰ) by comparing the propagation predicted from the radar-inferred refractivity profiles with that from the helicopter profiles, (ⅱ) by comparing the refractivity parameters from the helicopter soundings with those estimated. This technique could provide near-real-time estimation of ducting effects.     

利用粒子滤波从雷达回波实时跟踪反演大气波导

粒子滤波(particle filter,PF)是利用蒙特卡洛仿真方法处理递推估计问题的非线性滤波算法,这种方法不受模型线性和高斯假设的约束,是处理非线性非高斯动态系统状态估计的有效算法,适用于雷达回波反演大气波导(RFC)这类非线性非高斯问题.文中分别介绍了PF的基本思想和具体算法实现步骤,最后导出PF反演算法的迭代求解格式.数值试验结果表明,与扩展卡尔曼滤波(extended kalman filter,EKF)和不敏卡尔曼滤波(unscented kalman filter,UKF)相比,PF更适用于RFC这类高度非线性反演问题,可有效提高反演结果的稳定性和精度.     

一种全球定位系统接收机单点定位的新算法及数值试验

借助于Bancroft方法快速决定全球定位系统中接收机的大概位置，可作为观测方程线性化的初值.采用最小二乘法来求解方程时，观测数据存在较大误差(即出现粗差)的情况下，最小二乘解会出现不稳定性，即问题呈不适定性.为了克服解的不适定性，引入了正则化方法，采用正则化方法及正则化参数最优选择来提高抗粗差能力.实验结果表明，利用这种定位方法的抗粗差能力有一定的改进，这对实时、快速定位有重要意义. 关键词： 全球定位系统 Bancroft方法 最小二乘法 正则化     

基于COSMIC掩星资料的全球第二对流层顶详细特征

使用2006年12月–2008年11月COSMIC (constellation observing system for meteorology, ionosphere and climate)掩星湿廓线资料对第二对流层顶的全球分布特征进行统计, 对比三个站点的无线电探空仪和COSMIC的对流层顶资料, 研究结果表明: 1) 第二对流层顶的出现频率在副热带急流区较高,冬季在北半球为50%–70%, 在南半球为20%–40%; 2) 赤道带的第二对流层顶出现频率约为20%–26%, 与越赤道急流和对流层顶上的毛卷云有关; 3) 在副热带急流区, 第一对流层顶的温度基本高于第二对流层顶; 4) 在热带, 对流层顶厚度和第二对流层顶出现频率随纬度减小; 热带以外, 对流层顶厚度随纬度增加, 在冬半球60^o有最大值7–8 km; 5) 单站点对流层顶的日变化剧烈程度与COSMIC和探空仪的对流层顶高度偏差正相关．     

Ray tracing/correlation approach to estimation of surface-based duct parameters from radar clutter

This paper describes a technique to estimate surface-based duct parameters by using a simple ray tracing/correlation method. The approach is novel in that it incorporates the Spearman rank-order correlation scheme between the observed surface clutter and the surface ray density for a given propagation path. The simulation results and the real data results both demonstrate the ability of this method to estimate surface-based duct parameters. Compared with the results obtained by a modified genetic algorithm combined with the parabolic wave equation, the results retrieved from the ray tracing/correlation scheme show a minor reduction in accuracy but a great improvement on computation time. Therefore the ray tracing/correlation method might be used as a precursor to more sophisticated and slower techniques, such as genetic algorithm and particle filters, by narrowing the parameter search space and providing a comprehensive and more efficient estimation algorithm.     

变分伴随正则化方法从雷达回波反演海洋波导 Ⅰ.理论推导部分

针对传统统计反演算法在雷达回波反演海洋波导(RFC)方面计算量过大的问题,提出一种变分伴随正则化物理反演算法.在变分伴随方法中分别导出切线性模式、伴随方程及伴随边界条件、伴随方程求解表达式、以及泛函梯度的数学表达形式;讨论了伴随方法求泛函梯度对复方程如何协调的问题.考虑到雷达电磁波传播的特征,选择了适当的正则化项来解决反演中的不适定问题.最后给出反演算法实施中的迭代格式.