Performance Analysis for Matched-Field Source Localization: Simulations and Experimental Results

Matched-field methods concern estimation of source locations and/or ocean environmental parameters by exploiting full wave modeling of acoustic waveguide propagation. Typical estimation performance demonstrates two fundamental limitations. First, sidelobe ambiguities dominate the estimation at low signal-to-noise ratio (SNR), leading to a threshold performance behavior. Second, most matched-field algorithms show a strong sensitivity to environmental/system mismatch, introducing biased estimates at high SNR. In this paper, some theoretical developments on matched-field performance analysis are summarized, including Bayesian performance bounds and probabilistic ambiguity analysis, both incorporating environmental/system uncertainty/mismatch. Performance analysis is then implemented for source localization in a typical shallow water environment chosen from the Shallow Water Evaluation Cell Experiments (SWellEX). The performance predictions describe the simulations of the maximum-likelihood estimator (MLE) well, including the mean-square error (MSE) in all SNR regions as well as the bias at high SNR. The threshold SNR and bias predictions are also validated through SWellEX experimental data processing. The results suggest the current environmental, acoustic, and statistical modeling has developed to such a level that the optimum theoretical matched-field performance can be achieved in a well-controlled experiment.     

基于二阶锥规划的空域维纳滤波

提出了基于二阶锥规划的空域维纳滤波。基于噪声与阵列流型不相关的特性,将阵列的方向向量作为期望通过二阶锥规划方法对阵列快拍进行维纳滤波。利用维纳滤波抑制噪声的能力,提高阵列的信号检测和方位估计性能。通过仿真研究,验证了单目标时,其低信噪比下的信号检测性能优于Bartlett波束形成,其方位估计的信噪比门限要低于Bartlett波束形成和MUSIC算法,因而是一种优良的波束形成器。该波束形成器用于多波束测深仪和侧扫声纳等海洋测绘设备的阵列,可以有效的提高设备的作用距离和测量精度。     

A comparison of estimators for the generalised Pareto distribution

The generalised Pareto distribution (GPD) is often used to model the distribution of storm peak wave heights exceeding a high threshold, from which return values can be calculated. There are large differences in the performance of various parameter and quantile estimators for the GPD. Commonly used estimation methods such as maximum likelihood or probability weighted moments are not optimal, especially for smaller sample sizes. The performance of several estimators for the GPD is compared by the Monte Carlo simulation and the implications for estimating return values of significant wave height are discussed. Of the estimators compared, the likelihood-moment (LM) estimator has close to the lowest bias and variance over a wide range of sample sizes and GPD shape parameters. The LM estimator always exists, is simple to compute and has a low sensitivity to choice of threshold. It is recommended that the LM estimator is used for calculating return values of significant wave height when the sample size is less than 500. For sample sizes above 500 the NEW estimator of Zhang and Stephens (2009) can give accurate results for low computational cost.     

Towed array shape estimation using frequency-wavenumber data

Towed array beamforming algorithms require accurate array shape information in order to perform properly. Very often, these algorithms assume the array is linear. Unfortunately, the mechanical forces on the array due to ship motion and sea dynamics can change the shape of the array, which degrades the performance of the beamforming algorithm. A data-driven approach to estimating the relative shape of a nominally linear array is presented. The algorithm is robust in that it optimally combines information contained in a wide band of frequencies and source bearings. At the heart of the algorithm is a maximum-likelihood (ML) estimation scheme. The Cramer-Rao lower bound is derived and compared to the performance of the ML estimator. The utility of the algorithm is verified using both simulated and actual towed array data experiments     

Channel Modeling and Threshold Signal Processing in Underwater Acoustics: An Analytical Overview

An overview of underwater acoustic channel modeling and threshold signal processing is presented, which emphasizes the inhomogeneous, random, and non-Ganssian nature of the generalized channel, combined with appropriate weak-signal detection and estimation. Principal attention is given to the formal structuring of the scattered and ambient acoustic noise fields, as well as that of the desired signal, including both fading and Doppler "smear" phenomena. The role of general receiving arrays is noted, as well as their impact on spatial and temporal signal processing and beam forming, as indicated by various performance measures in detection and estimation. The emphasis here is on limiting optimum threshold systems, with some attention to suboptimum cases. Specific first-order probability density functions (pdf's) for the non-Ganssian components of typical underwater acoustic noise environments are included along with their field covariances. Several examples incorporating these pdf's are given, to illustrate the applications and general methods involved. The fundamental role of the detector structure in determining the associated optimum estimators is noted: the estimators arc specific linear or nonlinear functionals of the original optimum detector algorithm, depending on the criterion (i.e., minimization of the chosen error or cost function) selected. Results for both coherent and incoherent modes of reception are presented, reflecting the fact that frequently signal epoch is not known initially at the receiver. To supplement the general discussion, a selected list of references is included, to provide direct access to specific detailed problems, techniques, and results, for which the present paper is only a guide.     

基于高光谱遥感的渤海海冰厚度半经验模型

Sea ice thickness is one of the most important input parameters for the prevention and mitigation of sea ice disasters and the prediction of local sea environments and climates. Estimating the sea ice thickness is currently the most important issue in the study of sea ice remote sensing. With the Bohai Sea as the study area, a semiempirical model of the sea ice thickness(SEMSIT) that can be used to estimate the thickness of first-year ice based on existing water depth estimation models and hyperspectral remote sensing data according to an optical radiative transfer process in sea ice is proposed. In the model, the absorption and scattering properties of sea ice in different bands(spectral dimension information) are utilized. An integrated attenuation coefficient at the pixel level is estimated using the height of the reflectance peak at 1 088 nm. In addition, the surface reflectance of sea ice at the pixel level is estimated using the 1 550–1 750 nm band reflectance. The model is used to estimate the sea ice thickness with Hyperion images. The first validation results suggest that the proposed model and parameterization scheme can effectively reduce the estimation error associated with the sea ice thickness that is caused by temporal and spatial heterogeneities in the integrated attenuation coefficient and sea ice surface. A practical semi-empirical model and parameterization scheme that may be feasible for the sea ice thickness estimation using hyperspectral remote sensing data are potentially provided.     

Prediction of short-term distributions of load extremes of offshore wind turbines

This paper proposes a new methodology to select an optimal threshold level to be used in the peak over threshold (POT) method for the prediction of short-term distributions of load extremes of offshore wind turbines. Such an optimal threshold level is found based on the estimation of the variance-to-mean ratio for the occurrence of peak values, which characterizes the Poisson assumption. A generalized Pareto distribution is then fitted to the extracted peaks over the optimal threshold level and the distribution parameters are estimated by the method of the maximum spacing estimation. This methodology is applied to estimate the short-term distributions of load extremes of the blade bending moment and the tower base bending moment at the mudline of a monopile-supported 5MW offshore wind turbine as an example. The accuracy of the POT method using the optimal threshold level is shown to be better, in terms of the distribution fitting, than that of the POT methods using empirical threshold levels. The comparisons among the short-term extreme response values predicted by using the POT method with the optimal threshold levels and with the empirical threshold levels and by using direct simulation results further substantiate the validity of the proposed new methodology.     

Transient Hydroelastic Response of VLFS by FEM with Impedance Boundary Conditions in Time Domain

A time-dependent finite element method (FEM) is developed to analyze the transient hydroelastie responses of very large floating structures (VLFS) subjected to dynamic loads. The hydrodynamic problem is formulated based on the linear theory of fluid and the structural response is analyzed based on the thin plate theory. The FEM truncates the unbounded fluid domain by introducing an artificial boundary surface, thus defining a finite computational domain. At this boundary surface an impedance boundary conditions are applied so that no wave reflections occur. In the proposed scheme, all of the procedures are processed directly in time domain, which is efficient for nonlinear analyses of structure floating on unbounded fluid. Numerical results indicate acceptable accuracy of the proposed method.     

地波雷达目标检测跟踪联合处理粒子滤波方法

针对固定粒子数PF-TBD算法计算量大、复杂环境下地波雷达海上船只目标检测与跟踪性能不佳的问题,本文将粒子滤波方法应用于地波雷达船只目标检测与跟踪中,提出了基于自适应粒子滤波的地波雷达目标检测与跟踪联合处理方法。该方法结合地波雷达回波谱中目标展宽特性,充分利用了地波雷达回波谱中面目标的粒子权重信息来设置粒子自适应采样策略,提高了目标检测和跟踪联合处理的效果。通过地波雷达实测数据的目标跟踪结果及与同步AIS信息的比对分析,结果表明:提出的检测跟踪联合处理方法在对低信噪比、快速机动等复杂环境下的多目标跟踪时,可提高目标整体跟踪性能。     

航海雷达图像噪声抑制方法研究

针对航海雷达图像上常见的同频干扰噪声,提出了一种有效的抑制方法。该方法采用阈值分割技术对噪声进行检测,然后对检测出的噪声点进行插值修复。使用X波段航海雷达图像对该方法进行检验。采用信噪比、偏差指数和相关系数等指标进行定量评价。结果表明,提出的新方法既能显著地抑制航海雷达图像的同频干扰噪声,又能较好地保持原始图像的信息,具有良好的降噪效果。它的综合性能明显优于传统的均值滤波和中值滤波方法。     

海流剖面仪信号处理系统中最小二乘改进算法的应用研究

ADCP技术是海洋测流领域的热门技术,而频谱估计方法的研究是该技术的核心。文章介绍一种最小二乘的改进算法(ELMS),利用该算法可以很好的提取出淹没在白噪声下的正弦频谱。和传统的快速傅立叶变换方法(FFT)相比,该算法极大的压缩了数据量的要求,而且在信噪比为0DB的时候,仍能够较准确的估计出信号的谱峰位置。实验结果辩明该算法具有良好的鲁棒性。     

Performance comparison of high resolution bearing estimationalgorithms using simulated and sea test data

The performance of both the Capon and the MUSIC high resolution bearing estimation algorithms is investigated using both simulated data and sea test data collected with an experimental planar array. The major problem with these estimators is their sensitivity to both system errors and deviations from the assumed noise model. To alleviate this problem, two methods for preprocessing the data before they are input into the high-resolution algorithm are investigated: beam space and sector focused stability. The performance of both high-resolution estimators is examined, using both types of preprocessing, and the results are compared with those for the standard element-space (ES) techniques, assuming both finite cross-spectral-matrix (CSM) averaging errors and weakening target strengths. For the Capon estimator the performance is only superior to the standard element space technique when the CSM is calculated using a small number of averages. For the MUSIC estimator, both preprocessing techniques give clearly superior results over standard space techniques, with the SFS preprocessor performing the best     

基于广义极值分布的设计波高推算

简介了广义极值分布函数及其3种参数估计方法,包括极大似然(ML)、线性矩(LM)和间隔最大积(MPS)估计的计算方法。使用广义极值分布函数推算了北部湾涠洲岛海域3个波向的年波高极值序列设计波高,并与Weibull分布、Gumbel分布和皮尔逊Ⅲ型分布的推算结果加以对比。分析表明,涠洲岛海域极值波高服从于广义极值Ⅲ型分布,拟合优度检验结果表明广义极值分布能更好地拟合极值波高;MPS方法是一种优良的参数估计法,推算的设计波高可作为海岸环境工程设计的首要参考值。     

A closed-form solution to bearings-only target motion analysis

Bearings-only target motion analysis is a nonlinear state estimation problem in which the noise corrupted angle of arrival measurements of an emitted signal are used to obtain estimates of the source's range, bearing, course, and speed. The estimation process is complicated by unusual observability properties that render the quality of the estimate highly dependent on both the measurement noise levels and the source-observer geometry. Solutions that use recursive Kalman filtering approach or batch-style algorithms have been reported. The nonlinear batch style estimators for this process require iterative solution methods and under certain scenarios can be sensitive to initial conditions. Pseudolinear solutions that alleviate some of the difficulties with the iterative batch algorithms have been proposed. Although early versions of the pseudolinear filter suffered from biased estimates, subsequent improvements appear to have reduced the bias problem. This paper discusses a new pseudolinear solution based on the observable parameters from individual data segments defined by periods of constant observer velocity (termed “legs”). This solution is a true closed-form solution to the bearings-only target motion analysis problem. Although theoretically interesting, the technique does suffer under conditions of poor observability. A practical pseudolinear estimate, that does not suffer from the same observability problems, is developed and related to the first solution. Algorithm performance results, obtained from computer simulation, are presented. For the scenarios examined, the technique provides good state estimates under conditions of high observability. As observability conditions deteriorate, the solution does develop biases. However, it may still be useful for initializing an iterative nonlinear batch-style estimation algorithm     

Adaptive output feedback control based on DRFNN for AUV

The tracking control problem of AUV in six degrees-of-freedom (DOF) is addressed in this paper. In general, the velocities of the vehicles are very difficult to be accurately measured, which causes full state feedback scheme to be not feasible. Hence, an adaptive output feedback controller based on dynamic recurrent fuzzy neural network (DRFNN) is proposed, in which the location information is only needed for controller design. The DRFNN is used to online estimate the dynamic uncertain nonlinear mapping. Compared to the conventional neural network, DRFNN can clearly improve the tracking performance of AUV due to its less inputs and stronger memory features. The restricting condition for the estimation of the external disturbances and network's approximation errors, which is often given in the existing literatures, is broken in this paper. The stability analysis is given by Lyapunov theorem. Simulations illustrate the effectiveness of the proposed control scheme.     

Robust high-resolution direction-of-arrival estimation via signaleigenvector domain

A robust high-resolution direction-of-arrival (DOA) estimation approach for coherent/noncoherent sources is presented. The approach is based on the fact that the signal eigenvectors of the covariance matrix are a linear combination of the direction vectors that contain the DOA information. By applying a high-resolution frequency estimation algorithm to an element sequence from a combination of the signal eigenvectors, the approach achieves better performance at low SNR than the conventional methods. It is shown that the improvement in performance increases with the number of snapshots. For example, the resolution improvement of the proposed signal eigenvector domain approach over spatial-smoothed minimum-norm is about 2.5 dB and 7 dB for 20 and 100 snapshots, respectively     

The estimation of the shape of an array using a hidden Markov model

A hidden Markov model (HMM) technique for the estimation of the shape of a towed array is presented. It is assumed that there is a far-field source radiating sound containing possibly weak spectral lines. The technique uses either the Fourier coefficients at a given frequency computed from a single time block or the maximal eigenvector of a sample spectral covariance matrix. The technique is illustrated using several simulations. The results of these simulations indicate that the HMM technique yields shape and bearing estimates more accurate than those provided by a maximum-likelihood array shape estimation technique     

Laboratory investigation of damping of gravity-capillary waves on the surface of turbulized liquid

Investigation of damping of gravity-capillary waves (GCWs) in the presence of turbulence is a classical hydrodynamic problem which has important geophysical applications, one of which is related with the problem of forming a radar and optical image of a ship wake on wavy water surface. In this work a new method for the laboratory study of surface wave damping in turbulized liquid is described and the results are presented. The damping of standing GCWs by turbulence on the water surface in a tank mounted on a vibration table is studied. GCWs and turbulence are excited using a two-frequency mode of vibration table oscillations. A high-frequency small amplitude signal is used for parametric GCW excitation; a low-frequency large amplitude signal is used for generating turbulence due to water flowing through a fixed perforated grid submerged into the tank. The coefficient of GCW damping is determined by measured threshold of parametric excitation of the waves; turbulence characteristics are determined by the PIV and PTV techniques. Dependences of GCW damping coefficients on their frequency at different turbulence intensities are obtained, estimates for turbulent viscosity are presented, and a comparison with empirical models proposed earlier is performed.     

Passive sonar detection and localization by matched velocityfiltering

This paper presents a new bearings-only method of detecting and tracking low signal-to-noise ratio (SNR) wideband targets on a constant course and velocity trajectory. A track-before-detect strategy based on matched velocity filtering is adopted using spatial images constructed from a sequence of power bearing map (PBM) estimates accumulated during a track. To lower the threshold SNR for detection, a discrete bank of matched velocity filters integrates the PBM images over a range of hypothesized trajectories, such an approach eliminates the need to estimate the number of targets since signal detection is determined by comparing the output of each matched filter (MF) to a decision threshold. The distribution of the MF output is derived based on a single point target in diffuse noise assumption. Receiver operating characteristic curves show a definite detection gain under low SNR conditions for matched velocity filtering (track-before-detect) over detection from a single PBM     

Four-element CODAR beam forming

An HF radar called the Coastal Ocean Dynamics Applications Radar (CODAR) is presently being used in several forms to measure ocean surface parameters. The original version was developed by the National Oceanic and Atmospheric Administration (NOAA) and utilizes a four-element receive array. The array consists of four equally spaced elements arranged on a circle with a radius of 0.2151 wavelengths (at 25.4 MHz). It was designed to measure ocean currents using a direction-finding technique based on an extension to a simple two-element interferometer. The problem of determining the bearing of a radiating source can be readily shown to be equivalent to that incurred in spectral estimation. In an attempt to improve upon the processing of existing data, modern nonlinear spectral estimation techniques are applied in a beam-forming bearing estimation procedure and compared against several direction-finding algorithms. Enhancement of bearing estimators via analysis of the eigenstructure of a spatial correlation matrix is included. Antenna response patterns are calculated and used to investigate properties of direction-finding algorithms. Simulated data are used for a comparison of direction finding and beam forming. The asymmetrical bias of each method is investigated to determine its effect on the error in estimating the angle of arrival of a radar target.