Error Covariance Analysis of a Karhunen-Loeve Random Field Estimator

An error covariance analysis of a two-dimensional Karhunen-Loeve random field estimator (KLE) for gridded gravity data is presented without actually using the data. Attention is focused particularly on the north-south deflection component of the disturbance vector so that the rms value of its residual error may not exceed 3 percent of the signal rms for grid spacing of not less than 5 nm, when signal-to-noise ratio in the data is varied within reasonable limits. To achieve this rather stringent goal, a discrete KLE derived entirely from a two-dimensional grid configuration is needed. The KLE for geodesy applications was developed initially by Bose. Using continuous measurements and the orthogonality relations of the Karhunen-Loeve (KL) eigenvectors in a continuous domain, continuous KL gain coefficients ?mnj and consequently the continuous KL gain vector Kmn were obtained; thus a satisfactory continuous KLE was achieved. However the discrete version of the KLE was then given by keeping the old (continuous) ?mnj, and simply computing the KL eigenvectors on the grid coordinates to evaluate the discrete KL gain vector Kmn. Clearly this is not a satisfactory discrete KLE since , ?mnj are not constants but rather depend on the measurement noise variances. So if the measurements are now going to be discrete, one must rederive ?mnj for this situation.     

Error Covariance Analysis of a Karhunen-Loeve Collocation Estimator

An error covariance analysis of a two-dimensional gravity compensation technique (KLC) employing a Karhunen-Loeve gravity disturbance model and the linear least-square collocation algorithm for its estimation is presented, without actually using any data. Its performance is compared with another gravity compensation technique (KLE), whose error covariance analysis was previously presented by Gupta. From the mismodeling analysis, KLC appears to be superior to KLE.     

A Modified Karhunen-Loeve Model of Gravity Disturbance Vector

A modified anisotropic Karhunen-Loeve (KL) model for the gravity disturbing potential is proposed. This model, unlike the previous KL model of Bose, does not vanish on the boundaries of a local region. Laplace's equation for the gravity disturbing potential and flatearth Vening-Meinesz relationships of the gravity disturbance vectorcomponents (?,?,?g) when the model is isotropic, are exploited to determine the unknown parameters of the model. The resulting model, in addition to having random boundary values, has three more parameters than those found in Bose's KL model. This should provide extra flexibility in fitting the model to the noisy data available in a local finite region.     

一种基于火箭橇轨道的动基座重力 测量系统评价方法

传统的重力仪精度评价方法为重力仪测量值与静态定点值或飞行状态下的重力线网格参考值进行比较,根据差值或相对异常值来评价重力仪的精度指标,但其缺点是飞行轨迹重复性差.火箭橇轨道是进行火箭橇测试的重要组成部分,它严格按照设计建设,长达数公里,拥有完整已知的地理数据信息,且沿着火箭橇轨道布有众多的绝对重力基准点.提出一种基于火箭橇轨道的重力测量系统评价方法,通过火箭橇轨道上的众多绝对重力基准点和重力仪测量得到的重力值比较,评价重力仪的精度指标性能.其优点是该方法能够连续、快速、可重复评价重力测量系统的精度指标,为重力仪研制提供精确的评价手段和鉴定方法.     

跨音速喷管流场波系分布的数值仿真

应用分离系数矩阵方法对跨音速喷管流场进行数值仿真, 同时计算喷管内马赫波分布。采用5151个网格点的计算结果清晰反映出喷管内膨胀波和压缩波的生成和反射。数值仿真与实验结果在宽广马赫数范围内吻合一致, 并揭示了波系分区分布, 反映出波系影响下跨音速喷管流场参数分布的规律性。      

基于径向基函数与混合背景网格的动态网格变形方法

发展了一种基于径向基函数(RBF)插值和混合背景网格映射的动态网格变形方法。混合背景网格由运动物面附近的各向异性网格和偏离运动物面的各向同性网格组成。物面变形或运动后,包含在各向异性背景网格中的所有物面网格点的位移由指定的物面运动规律得到,而各向同性背景网格点的位移通过径向基函数插值获得。然后变形后的计算网格利用计算网格与混合背景网格之间固定的映射关系代数插值得到。各向同性网格距离物面一定的距离,消除了物面保形对径向基函数插值的要求,从而可以减少径向基函数基点的数量,以提高插值的效率和可靠性。最后,通过二维和三维算例对该方法的变形网格质量、参数影响、网格拓扑依赖性、复杂外形变形能力及网格变形效率进行了测试。测试结果表明,该方法能够高效生成大变形下高质量的变形网格。     

重力和重力梯度基准图构建的多参量联合网格化

高精度、高分辨率的重力及重力梯度基准图是决定潜艇水下辅助导航定位精度的关键因素。我国海洋重力测量目前的主要比例尺为1:1000000，在局部区域可达1:500000，测线间隔分布相对稀疏。传统的网格化插值技术在远离测点位置时容易产生虚假异常，或异常特征发生偏移。构建了全张量重力梯度数据的6个梯度分量的联合网格化方法,利用各个梯度分量与引力位在波数域中的关系重构了引力位，实现了重力和重力梯度数据的再计算，从而实现了网格处理。通过模型数据和实测数据，验证了该方法在插值异常分辨率和位置准确性上具有优势。     

Fast approximate maximum likelihood period estimation from incomplete timing data

 To estimate the period of a periodic point process from noisy and incomplete observations, the classical periodogram algorithm is modified. The original periodogram algorithm yields an estimate by performing grid search of the peak of a spectrum, which is equivalent to the periodogram of the periodic point process, thus its performance is found to be sensitive to the chosen grid spacing. This paper derives a novel grid spacing formula, after finding a lower bound of the width of the spectral mainlobe. By employing this formula, the proposed new estimator can determine an appropriate grid spacing adaptively, and is able to yield approximate maximum likelihood estimate (MLE) with a computational complexity of O(n²). Experimental results prove that the proposed estimator can achieve better trade-off between statistical accuracy and complexity, as compared to existing methods. Simulations also show that the derived grid spacing formula is also applicable to other estimators that operate similarly by grid search.     

基于径向基函数改进的Delaunay图映射动网格方法

为了提高大变形下动网格生成的效率和质量,基于径向基函数插值技术对原始Delaunay图映射动网格方法进行了改进。首先通过带多项式的径向基函数插值方法计算出背景网格远场边界点的位移;然后将背景网格位移插值到计算网格;最后利用衰减函数将计算网格远场位移衰减为零。基于矩形旋转网格变形实例比较了改进方法与原始Delaunay图映射动网格方法之间的差别,并研究了径向基基点数目和衰减函数类型对变形网格质量的影响。矩形旋转网格变形实例说明该方法能够精确恢复出Delaunay背景网格的旋转特性。NACA 0012翼型、NLR 7301两段翼和M6机翼网格变形算例进一步证明,通过添加控制点,该方法能够不重构背景网格实现大变形下高质量动网格的生成。     

Precision large scale air traffic surveillance using IMM/assignmentestimators

We present the development and implementation of a multisensor-multitarget tracking algorithm for large scale air traffic surveillance based on interacting multiple model (IMM) state estimation combined with a 2-dimensional assignment for data association. The algorithm can be used to track a large number of targets from measurements obtained with a large number of radars. The use of the algorithm is illustrated on measurements obtained from 5 FAA radars, which are asynchronous, heterogeneous, and geographically distributed over a large area. Both secondary radar data (beacon returns from cooperative targets) as well as primary radar data (skin returns from noncooperative targets) are used. The target IDs from the beacon returns are not used in the data association. The surveillance region includes about 800 targets that exhibit different types of motion. The performance of an IMM estimator with linear motion models is compared with that of the Kalman filter (KF). A number of performance measures that can be used on real data without knowledge of the ground truth are presented for this purpose. It is shown that the IMM estimator performs better than the KF. The advantage of fusing multisensor data is quantified. It is also shown that the computational requirements in the multisensor case are lower than in single sensor case, Finally, an IMM estimator with a nonlinear motion model (coordinated turn) is shown to further improve the performance during the maneuvering periods over the IMM with linear models     

基于STL文件的笛卡尔网格生成方法研究

采用从STL文件中直接获取模型几何数据的方法,利用ADT数据结构快速剔除相同顶点,提出一种新的“推进查找”方法替代以往网格单元是否在物面内部的判定方法,避免了对物面的自封闭要求,无需进行几何重构,直接完成网格生成,更进一步提高了自适应笛卡尔网格生成速度和自动化程度,为实现CAD与CFD的之间的数据对接提供了一种新的方式。     

Delaunay非结构网格生成之布点技术

Ｄｅｌａｕｎａｙ 非结构网格生成技术是 Ｃ Ｆ Ｄ 领域广为应用的有效技术之一。如何布置流场中的计算点是 Ｄｅｌａｕｎａｙ 三角剖分的一个关键问题。本文以二维非结构网格生成为例， 介绍三种布点策略， 即重心布点、外接圆心布点和 Ｖｏｒｏｎｏｉ 边布点， 算例表明三种方法均可生成满意的非结构网格， 其中 Ｖｏｒｏｎｏｉ 边布点方法生成的网格质量更好一些。三种方法均可方便地推广到三维问题。     

Ground target tracking with variable structure IMM estimator

In this paper we present the design of a Variable Structure Interacting Multiple Model (VS-IMM) estimator for tracking groups of ground targets on constrained paths using Moving Target Indicator (MTI) reports obtained from an airborne sensor. The targets are moving along a highway, with varying obscuration due to changing terrain conditions. In addition, the roads can branch, merge or cross-the scenario represents target convoys along a realistic road network with junctions, changing terrains, etc. Some of the targets may also move in an open field. This constrained motion estimation problem is handled using an IMM estimator with varying mode sets depending on the topography, The number of models in the IMM estimator, their types and their parameters are modified adaptively, in real-time, based on the estimated position of the target and the corresponding road/visibility conditions. This topography-based variable structure mechanism eliminates the need for carrying all the possible models throughout the entire tracking period as in the standard IMM estimator, significantly improving performance and reducing computational load. Data association is handled using an assignment algorithm. The estimator is designed to handle a very large number of ground targets simultaneously. A simulated scenario consisting of over one hundred targets is used to illustrate the selection of design parameters and the operation of the tracker. Performance measures are presented to contrast the benefits of the VS-IMM estimator over the Kalman filter and the standard IMM estimator, The VS-IMM estimator is then combined with multidimensional assignment to gain “time-depth.” The additional benefit of using higher dimensional assignment algorithms for data association is also evaluated     

用非结构动网格技术计算确定再入飞行器配平攻角

提出了用数值计算直接确定飞行器配平攻角的自由振动法.该方法采用非结构动网格技术,使飞行器在气动力作用下绕通过质心的轴线自由振动,最终衰减到配平攻角.构造了一种MUSCL类型有限体积格式,求解考虑了动网格效应的三维非定常Euler方程,采用弹簧近似技术实现非结构网格的变形运动.作为例证,计算确定了载人飞船返回舱的配平攻角,计算结果与实验结果接近.该方法外形适应能力强,效率较高,具有较好的应用前景.     

Determination of optimal samples for robot calibration based on error similarity

Industrial robots are used for automatic drilling and riveting.The absolute position accuracy of an industrial robot is one of the key performance indexes in aircraft assembly,and can be improved through error compensation to meet aircraft assembly requirements.The achievable accuracy and the difficulty of accuracy compensation implementation are closely related to the choice of sampling points.Therefore,based on the error similarity error compensation method,a method for choosing sampling points on a uniform grid is proposed.A simulation is conducted to analyze the influence of the sample point locations on error compensation.In addition,the grid steps of the sampling points are optimized using a statistical analysis method.The method is used to generate grids and optimize the grid steps of a Kuka KR-210 robot.The experimental results show that the method for planning sampling data can be used to effectively optimize the sampling grid.After error compensation,the position accuracy of the robot meets the position accuracy requirements.     

飞机起落架滑行载荷识别

 在飞机滑行方向跑道路面不平度为平稳随机函数的假设下,采用滑行过程分段等速化、起落架动态特性分段线性化方法,把原非线性非平稳随机振动问题,简化为分段线性平稳随机振动问题。在频域中采用频响函数描述飞机动态特性,计及全机刚体平动和俯仰模态及结构前ｎ阶弹性模态,给出了从飞机重心处加速度观测值识别起落架接地点载荷自谱密度函数的方法。起落架当量线性参数采用迭代求解。算例表明,本方法和程序可行。     

Eruptive prominences as sources of magnetic clouds in the solar wind

Large amounts of coronal material are propelled outward into interplanetary space by Coronal Mass Ejections (CMEs). Thus one might expect to find evidence for expanding flux ropes in the solar wind as well. To prove this assumption magnetic clouds were analyzed and correlated with H-observations of disappearing filaments. When clouds were found to be directly associated with a disappearing filament, the magnetic structure of the cloud was compared with that of the associated filament. Additionally the expansion of magnetic clouds was examined over a wide range of the heliosphere and compared with the expansion observed for erupting prominences.     

Orbital reference frame estimation with power spectral density constraints for drag-free satellites

The drag-free satellites are widely used in the field of fundamental science as they enable the high-precision measurement in pure gravity fields. This paper investigates the estimation of local orbital reference frame (LORF) for drag-free satellites. An approach, taking account of the combi-nation of the minimum estimation error and power spectral density (PSD) constraint in frequency domain, is proposed. Firstly, the relationship between eigenvalues of estimator and transfer func-tion is built to analyze the suppression and amplification effect on input signals and obtain the eigenvalue range. Secondly, an optimization model for state estimator design with minimum estima-tion error in time domain and PSD constraint in frequency domain is established. It is solved by the sequential quadratic programming (SQP) algorithm. Finally, the orbital reference frame estimation of low-earth-orbit satellite is taken as an example, and the estimator of minimum variance with PSD constraint is designed and analyzed using the method proposed in this paper.     

二维轴流压气机转—静非定常Euler方程数值计算

将Euler方程写在包括转子和静子的绝对坐标系下,通过在转—静之间设计一个可以随时间扭动的平行四边形网格计算了二维轴流压气机转—静非定常流场。由于整个方程统一写在动网格下,转—静之间没有使用插值。      

A Kalman Filter Based Tracking Scheme with Input Estimation

Beginning with the derivation of a least squares estimator that yields an estimate of the acceleration input vector, this paper first develops a detector for sensing target maneuvers and then develops the combination of the estimator, detector, and a "simple" Kalman filter to form a tracker for maneuvering targets. Finally, some simulation results are presented. A relationship between the actual residuals, assuming target maneuvers, and the theoretical residuals of the "simple" Kalman filter that assumes no maneuvers, is first formulated. The estimator then computes a constant acceleration input vector that best fits that relationship. The result is a least squares estimator of the input vector which can be used to update the "simple" Kalman filter. Since typical targets spend considerable periods of time in the constant course and speed mode, a detector is used to guard against automatic updating of the "simple" Kalman filter. A maneuver is declared, and updating performed, only if the norm of the estimated input vector exceeds a threshold. The tracking sclheme is easy to implement and its capability is illustrated in three tracking examples.