基于COS-2S模型的ADCP海浪方向谱估计

美国RD公司声学多普勒流速剖面仪(ADCP)使用的波向估计算法是迭代的最大似然算法(IMLM),IMLM可以准确估计出主浪向,但波向估计结果在主浪向以外其他方向上存在由噪声带来的能量分布,使得估计结果变得不准确。设计了基于COS-2S模型的波向估计方法,模型中只含有主波向,算法在准确估计主波向同时有效去除其他方向上噪声带来的能量分布。将一种混合遗传算法(HGA)代替遗传算法(GA)应用到波浪模型的参数求解中,HGA相对GA收敛更快,可有效找出全局最优解。仿真分析和实测数据表明,HGA和IMLM都能准确估计出波浪的主浪向,但HGA可有效去除主浪向以外其他方向上噪声带来的能量分布,HGA比IMLM更加适用于海洋各种噪声存在条件下的测量环境。     

多波束测深异常的两种趋势面检测算法比较

利用三次样条插值算法模拟海底地形曲面,并加入高斯白噪声和不同数量的异常值作为多波束测深数据模拟值。分别基于最小二乘估计和高崩溃污染率抗差估计两种算法建立趋势面模型,通过各自的异常值标定准则对模拟数据进行测深异常检测,比较和分析了两种算法的处理结果,并得出相应结论。最后,利用上述两算法对多波束实测数据进行处理,结果表明,经高崩溃污染率抗差趋势面异常值检测后的数据能够更为准确地反映海底的真实情况。     

基于BP神经网络的航磁补偿方法研究

水中未爆弹危害极大，航空磁探因探测效率高的优势，常用于水中未爆弹的探测。飞机干扰磁场的存在限制了航空磁探的发展。传统的 T-L 模型待求参数间具有很强的复共线性，难以满足高精度磁补偿的要求，而神经网络算法具有容错率高、求解精度高的特点。首先使用 BP 神经网络建立干扰磁场的数学模型； 随后，通过仿真生成干扰磁场以及未爆弹目标磁场信号对算法进行验证；最后，利用四旋翼无人机平台进行目标探测试验，补偿改善比超过 20。试验结果表明，算法可以用于提高平台对水下目标探测的精度。     

变系数分位数回归模型的贝叶斯推断

主要考察一种基于非对称Laplace分布的变系数分位数回归模型,并利用可逆跳的MCMC算法来估计了变系数函数。这种估计方法充分考虑了每一个变系数函数的差异性,允许不同的系数函数有不同的节点个数和位置。本文最后以具有典型性的波士顿房价作为变系数分位数回归的贝叶斯估计的例子,说明了该模型在现实问题中应用的合理性。     

基于复contourlet域隐马尔科夫树模型的海面溢油合成孔径雷达图像相干斑抑制

海面溢油SAR图像中的相干斑噪声严重影响了后续的图像分割、特征提取和分类.为了更有效地抑制海面溢油SAR图像相干斑,文中提出了一种基于复contourlet域隐马尔科夫树模型的海面溢油SAR图像相干斑抑制方法.首先对观测图像取对数并进行复contourlet变换;然后在复contourlet域中用隐马尔科夫树模型对相邻尺度间的带通方向子带系数进行建模,并依据贝叶斯最小均方误差准则估计无噪系数;最后进行逆复contourlet变换和指数变换,得到相干斑抑制后的图像.大量实验结果表明,与Lee、Kuan、Frost及Gamma Map等4种经典滤波方法以及小波域和contourlet域隐马尔科夫树模型方法相比,文中方法从主观视觉和客观定量评价两方面来看综合性能更为优越,是一种行之有效的SAR遥感图像海面溢油检测的预处理方法.     

基于FVCOM模式钱塘江河口涌潮三维数值模拟研究

应用FVCOM模型进行钱塘江涌潮的二维和三维数值模拟,并使用2007年10月大潮期间的现场观测结果进行模型验证。结果表明,平面二维模拟结果与实测值吻合较好,与基于Godunov算法求解钱塘江涌潮二维特征的结果基本一致。三维计算结果再现了流速在垂向上的差异,以及在落潮转向涨潮时刻,涨潮先从底部开始,垂向上从底层向表层逐步从落潮流转为涨潮流。敏感性分析显示模型网格大小和曼宁系数的选取对模拟结果有较大的影响。     

改进的基于非高斯双参数模型的图像滤波算法

以描述小波系数分布规律的非高斯双参数模型为基础,介绍和分析了BiShrink滤波算法,指出了其存在的不足。提出的改进算法采用冗余小波变换替代正交小波变换,将子带内小波系数的局部相关性纳入滤波过程,给出了局部自适应的阈值估计策略,再通过双参数联合收缩函数达到系数收缩的目的。实验结果表明,改进算法同时兼顾了子带内小波系数之间的相关性和尺度间系数的传播特性,在有效滤除噪声的同时,较好地保持了图像的细节信息。     

吉尼系数的计算及正态分布的吉尼系数

本文给出了吉尼系数的一般计算公式 ,推导出正态分布的吉尼系数渐进无偏估计 ,并以实例对比了传统的吉尼系数的算法——万分法与本文中提出的估计方法。     

病态半参数回归模型的岭估计求解

针对半参数回归模型求解过程可能出现的法方程病态问题，提出了用岭估计原则改进半参数模型的求解。通过模拟算例将岭估计解法和其他方法进行了比较，结果表明，岭估计解法能较好地解决半参数回归模型求解过程中的病态问题。     

正交解调及复自相关频率估计在海流测量中的应用

正交解调技术在使用复自相关算法进行频率估计中具有非常重要的作用,简要介绍了正交解调技术和复自相关算法在声学多普勒海流计中的应用.     

基于自适应抗差最小二乘的DEM数据粗差剔除

数字高程模型(DEM)生产中最关键的环节是质量控制,而粗差是影响DEM质量的重要因素之一。简单介绍了自适应最小二乘估计的原理,并将其用于DEM趋势面拟合中,结果表明该方法具有良好的抗粗差的能力,用含粗差的数据拟合出的曲面与实际地形比较符合。     

基于VMD和WinLSP相结合的GNSS-R海平面高度估测模型

全球导航卫星系统反射（Global Navigation Satellite System-Reflectometry,GNSS-R）技术是一种新兴的监测海平面高度变化的技术。本文依据GNSS-R技术中的信噪比分析法的原理,通过分析其分离趋势项和提取振荡频率的过程,建立了新的估测模型以提高反演精度。针对传统模型存在的信号分离不佳的问题,本文提出使用变分模态分解（Variational Mode Decomposition,VMD）算法替换传统的最小二乘拟合法（Least Squares Fitting, LSF）进行趋势项分量的分离。在此基础上,本文引入基于凯塞窗函数改进的LSP（Lomb-Scargle Periodogram）频谱分析法（记为WinLSP）来减弱因频谱泄露带来的反演误差。在瑞典翁萨拉的GTGU站和美国阿拉斯加州的SC02站开展的海平面高度反演实验结果表明,本文建立的估测模型相比于传统模型具有更高的反演精度。基于VMD+WinLSP估测模型得到的GTGU站反演结果的均方根误差（RMSE）、相关系数和反演点数分别为4.70 cm、0.98和5 647。与传统的LSF+LSP估测模型相比,反演精度和GNSS数据利用率分别提高了约29.7%和15.0%。SC02站的RMSE、相关系数和反演点数分别为14.34 cm、0.99和1 785,反演精度和GNSS数据利用率分别提高了约12.3%和9.4%。     

Estimation of Rapidly Time-Varying Sparse Channels

The estimation of sparse shallow-water acoustic communication channels and the impact of estimation performance on the equalization of phase coherent communication signals are investigated. Given sufficiently wide transmission bandwidth, the impulse response of the shallow-water acoustic channel is often sparse as the multipath arrivals become resolvable. In the presence of significant surface waves, the multipath arrivals associated with surface scattering fluctuate rapidly over time, in the sense that the complex gain, the arrival time, and the Dopplers of each arrival all change dynamically. A sparse channel estimation technique is developed based on the delay-Doppler-spread function representation of the channel. The delay-Doppler-spread function may be considered as a first-order approximation to the rapidly time-varying channel in which each channel component is associated with Doppler shifts that are assumed constant over an averaging interval. The sparse structure of the delay-Doppler-spread function is then exploited by sequentially choosing the dominant components that minimize a least squares error. The advantage of this approach is that it captures both the channel structure as well as its dynamics without the need of explicit dynamic channel modeling. As the symbols are populated with the sample Dopplers, the increase in complexity depends on the channel Doppler spread and can be significant for a severely Doppler-spread channel. Comparison is made between nonsparse recursive least squares (RLS) channel estimation, sparse channel impulse response estimation, and estimation using the proposed approach. The results are demonstrated using experimental data. In training mode, the proposed approach shows a 3-dB reduction in signal prediction error. In decision-directed mode, it improves significantly the robustness of the performance of the channel-estimate-based equalizer against rapid channel fluctuations.     

岭估计在多项式曲面拟合GPS高程中的应用

探讨了岭估计方法在GPS高程拟合中的应用,并给出了岭估计中确定岭参数k值的一种新方法——方差扩大因子法。通过某地区GPS水准网数据进行实验验证,结果表明,与传统多项式曲面拟合模型比较,岭估计方法具有拟合精度高、适用复杂的地形等特点。     

THE COMPUTATION OF OPTIMUM OCEAN SHIP ROUTING

A computational method for the ocean ship routing based on optimization theory is presented in this paper. In the computational programme "breaking circle method" is used, by which a principal contradictory line can be obtained in the operations research. The automatic selection of calculating network and the initial values have been realized. The quantitative estimation for the change of navigating speed under the different conditions of ocean wave and current has been done. The objectivly quantitative computational method is provided for ocean ship routing and forecast service following ship tracks. The calculating examples and practical experiments of ships have proved that this computational method is highly effective.     

A synthetic aperture radar sea surface distribution estimation by n-order Bézier curve and its application in ship detection

To dates,most ship detection approaches for single-pol synthetic aperture radar(SAR) imagery try to ensure a constant false-alarm rate(CFAR).A high performance ship detector relies on two key components:an accurate estimation to a sea surface distribution and a fine designed CFAR algorithm.First,a novel nonparametric sea surface distribution estimation method is developed based on n-order Bézier curve.To estimate the sea surface distribution using n-order Bézier curve,an explicit analytical solution is derived based on a least square optimization,and the optimal selection also is presented to two essential parameters,the order n of Bézier curve and the number m of sample points.Next,to validate the ship detection performance of the estimated sea surface distribution,the estimated sea surface distribution by n-order Bézier curve is combined with a cell averaging CFAR(CA-CFAR).To eliminate the possible interfering ship targets in background window,an improved automatic censoring method is applied.Comprehensive experiments prove that in terms of sea surface estimation performance,the proposed method is as good as a traditional nonparametric Parzen window kernel method,and in most cases,outperforms two widely used parametric methods,K and G0 models.In terms of computation speed,a major advantage of the proposed estimation method is the time consuming only depended on the number m of sample points while independent of imagery size,which makes it can achieve a significant speed improvement to the Parzen window kernel method,and in some cases,it is even faster than two parametric methods.In terms of ship detection performance,the experiments show that the ship detector which constructed by the proposed sea surface distribution model and the given CA-CFAR algorithm has wide adaptability to different SAR sensors,resolutions and sea surface homogeneities and obtains a leading performance on the test dataset.     

海上风的边界层模式

由于航海、海上开发作业等对海洋上风和海浪的预报提出越来越高的要求,而海浪、风暴潮等海洋水文要素的数值计算和预报,迫切需要解决海洋上风场的精确计算。但是,复杂的海面结构,大气稳定度的影响以及风、浪之间动量的交换等,使海上风的理论计算遇到很多困难,至今大部分工作是依靠统计方法。利用天气预报的形势场计算地转风或梯度风,以及它与海面摩擦、大气稳定度的经验订正关系。     

An adjoint method for the assimilation of statistical characteristics into eddy-resolving ocean models

The study investigates perspectives of the parameter estimation problem with the adjoint method in eddy-resolving models. Sensitivity to initial conditions resulting from the chaotic nature of this type of model limits the direct application of the adjoint method by predictability. Prolonging the period of assimilation is accompanied by the appearance of an increasing number of secondary minima of the cost function that prevents the convergence of this method. In the framework of the Lorenz model it is shown that averaged quantities are suitable for describing invariant properties, and that secondary minima are for this type of data transformed into stochastic deviations. An adjoint method suitable for the assimilation of statistical characteristics of data and applicable on time scales beyond the predictability limit is presented. The approach assumes a greater predictability for averaged quantities. The adjoint to a prognostic model for statistical moments is employed for calculating cost function gradients that ignore the fine structure resulting from secondary minima. Coarse resolution versions of eddy-resolving models are used for this purpose. Identical twin experiments are performed with a quasigeostrophic model to evaluate the performance and limitations of this approach in improving models by estimating parameters. The wind stress curl is estimated from a simulated mean stream function. A very simple parameterization scheme for the assimilation of second-order moments is shown to permit the estimation of gradients that perform efficiently in minimizing cost functions.     

改进的单纯形法及其在非线性参数估计中的应用

采用了改进的单纯形法，该算法利用遗传算法的变异操作来增强单纯形法的全局搜索能力，既兼顾了全局收敛性，又具有高效的局部搜索能力。并通过实例验证了该算法的优越性，以及应用于测量上的非线性最小二乘估计的可行性。     

Viscous free surface effect on coastal embankment hydrodynamics

A finite-differnece method was used to calculate the nonlinear hydrodynamic pressures acting on the coastal embankment faces by seismic-wave actions. The nonlinearity of free surface flow, convective acceleration, viscosity and surface tension of fluid are included in the analysis. The kinematic and dynamic free surface boundary conditions are employed for calculating the horizontal fluid velocity, pressure at the free surface and the surface profile of the fluid. The time-dependent water surface is transformed to the horizontal plane, and the flow field is mapped onto a rectangular, making it convenient to model the complex sea bottom geometry and the wavy water surface by the finite-difference method. Fully nonlinear and weakly nonlinear dynamic free surface conditions are used and compared. The effects of surface tension of fluid are also discussed. The nonslip boundary condition is applied on the most part of the interface between fluid and solid face, except the region near the intersection between free surface and wall face. The numerical results are presented for various water depths and ground motion intensities, and their associate viscous effects on coastal embankment hydrodynamics are discussed.