经验正交函数与广义数值环境模式重构声速剖面的比较

本文利用Argo资料对经验正交函数(EOF)和GDEM模式拟合声速剖面的效果进行了比较。结果表明,两种方法都能较好地重构局地声速剖面。EOF方法的拟合精度与选取的EOF阶数有关,阶数越高,精度越高;GDEM模式的拟合精度介于EOF的前5阶到前10阶之间。EOF易受到样本数量和采样深度的限制,难以提供完整的声速剖面;而GDEM模式能够将剖面扩展到更深的深度,与仅有上层海洋采样数据的声速剖面相比,结合深海水团的信息可使拟合结果达到更高的精度。     

浅水声速剖面用经验正交函数(EOF)表示的可行性研究

本文研究利用经验正交函数（Ｅ０Ｆ）及少量参数以近似表示实测浅水声速剖面的方法,研究了样本数,不同浅水域的样杜和算法的关系,并得出结论：一定区域的浅水声速剖面,以经验正交函数（Ｅ０Ｆ）近似表示,仍可以达到较好的精度,这对于声场层析及被动定位总理２有着很重要的意义。     

浅水声速剖面用经验正交函数(EOF)表示的可行性研究

本文研究利用经验正交函数（ＥＯＦ）及少量参数以近似表示实测浅水声速剖面的方法,研究了 样本数,不同浅水域的样本和算法的关系,并得出结论：一定区域的浅水声速剖面,以经验正交函数 （ＥＯＦ）近似表示,仍可以达到较好的精度。这对于声场层析及被动定位问题有着很重要的意义。     

激光雷达信号的可变间隔阈值经验模式分解去噪法

针对采用经验模式分解直接阈值(EMD-DT)和经验模式分解间隔阈值(EMD-IT)在激光雷达回波信号的去噪应用中会产生的模态混叠现象,采用一种可变间隔阈值的经验模式分解(EMD-SIT)的去噪方法。首先,对信号进行经验模式分解。然后,采用过零率方法将分解出的含有噪声的固有模态函数分离。最后,应用过零点阈值,设立一个新的可变阈值,将EMD-IT和EMD-DT有效融合对信号进行去噪。通过与多种阈值的仿真对比以及激光雷达的回波信号去噪实验,结果表明该方法可以有效地去除噪声,抑制模态混叠,较EMD-IT和EMD-DT更具有优越性,因此有着很好的应用前景。     

基于希尔伯特-黄变换方法的同轴粒子全息图中粒子轴向位置提取方法

在颗粒的数字全息中,传统方法应用数值再现反演计算得到颗粒的相关信息,其中再现的判焦过程中存在繁琐费时的缺点.本文提出了一种基于希尔伯特-黄变换方法的同轴粒子全息图分析方法.将粒子同轴全息图沿中心往外得到的单个粒子径向强度分布作为初始信号,根据希尔伯特-黄变换方法中的经验模态分解首先将信号分解成几个本征模态函数,通过分析第一本征模态函数的希尔伯特谱,由拉依达法则剔除部分奇异点后做最小二乘线性回归分析,得到粒子的空间位置.该方法不需要对全息图进行重建,由信号自身构建基函数,有很强的适应性.理论模拟和实验证明,该方法计算速度快,准确度高,有望应用于全息图的实时在线分析.     

A robust algorithm for high-resolution dynamic MRI based on the partially separable functions model

A recently developed partially separable functions (PSF) model can be used to generate high-resolution dynamic magnetic resonance imaging (MRI). However, this method could not robustly reconstruct high-quality MR images because the estimation of the PSF parameters is often interfered by the noise of the sampled MR data. To improve the robustness of MRI reconstruction using the PSF model, we proposed a new algorithm to estimate the PSF parameters by jointly using robust principal component analysis and modified truncated singular value decomposition regularization methods, instead of using the least square fitting method in the original PSF model. The experiment results of in vivo cardiac MRI demonstrated that the proposed algorithm can robustly reconstruct dynamic MR images with higher signal-to-noise ratio and clearer anatomical structures in comparison with the previous PSF model.     

利用经验模态分解的高频水声信号滤波方法*

研究了一种高频水声信号的滤波问题,提出了一种改进的经验模态分解加小波阈值滤波方法。首先对信号进行带通滤波处理,再进行经验模态分解,将分解得到的各个模态转换为频域信号,采用小波软阈值方法对这些频域信号进行滤波,最后对信号进行重构,并转换为时域信号。经数值仿真与试验数据验证表明此方法是可行有效的,与原基于经验模态分解的小波阈值滤波方法相比,本方法滤波效果较好:对不同输入信噪比的仿真信号进行滤波后,本方法的输出信噪比最大提高17.41 d B,滤波后所得信号与加噪前纯信号的相关系数最大提高0.90;对实验数据进行滤波后,不同时间段信号的相关系数最大提高0.62。     

Statistical techniques analysis of SST and SLP in the Persian Gulf

Thirty-four years of data (1967-2000) are used to investigate the variability pattern relevant to air-sea interaction in the Persian Gulf. The patterns are derived using statistical techniques, such as empirical orthogonal function (EOF) and singular value decomposition (SVD). Statistical analysis methods are applied to determine the coupled modes of variability of monthly sea surface temperature (SST) and sea level pressure (SLP). The significant of the air-sea interaction is found by a strong resemblance between EOF and SVD eigenvectors and expansion coefficients of SST and SLP. We find that the four leading EOF patterns of SST together account for 99.8% of the total monthly SST variance and 94.4% of the SLP variance. The zero contour in the first SST EOF identified the front which separates the Persian Gulf cyclonic gyres.The SVD modes provide more information on the coupling between the fields than the modes obtained by EOF methods. Lagged correlation analysis between SVD1(SLP) and SVD1(SST) indicates that the coupling is strongest when SLP leads SST by −12, −6, 6 and 12 months. Therefore, the first mode of the SVD analysis seems to depict an air-to-sea forcing, in which the sea response to the atmospheric changes appears with an semiannual and interannual time lag.The two leading SVD modes of variability of the coupled SST and SLP fields account for 99.6% of the total variance. The main patterns of both variables of variability of both variables independently provide considerable information on the coupling, but only one of the two variables dominates each of the two first coupled modes.The first coupled mode of variability between the SST and atmospheric pressure can be described as a strengthening and weakening of the cyclonic gyres, which seems to force fluctuations in a north-south dipole structure in the SST by Ekman upwelling which is a wind-related process. The atmospheric forcing of the SST changes is detectable in the sea with a lag of 1 and 6 months.     

Orthogonal-triangular decomposition ghost imaging

Ghost imaging (GI) offers great potential with respect to conventional imaging techniques. However, there are still some obstacles for reconstructing images with high quality, especially in the case that the orthogonal measurement matrix is impossible to construct. In this paper, we propose a new scheme based on the orthogonal-triangular (QR) decomposition, named QR decomposition ghost imaging (QRGI) to reconstruct a better image with good quality. In the scheme, we can change the randomly non-orthogonal measurement matrix into orthonormal matrix by performing QR decomposition in two cases. (1) When the random measurement matrix is square, it can be firstly decomposed into an orthogonal matrix $\bm Q$ and an upper triangular matrix $\bm R$. Then let the off-diagonal values of $\bm R$ equal to 0.0, the diagonal elements of $\bm R$ equal to a constant $k$, where $k$ is the average of all values of the main diagonal, so the resulting measurement matrix can be obtained. (2) When the random measurement matrix is with full rank, we firstly compute its transpose, and followed with above QR operation. Finally, the image of the object can be reconstructed by correlating the new measurement matrix and corresponding bucket values. Both experimental and simulation results verify the feasibility of the proposed QRGI scheme. Moreover, the results also show that the proposed QRGI scheme could improve the imaging quality comparing to traditional GI (TGI) and differential GI (DGI). Besides, in comparison with the singular value decomposition ghost imaging (SVDGI), the imaging quality and the reconstruction time by using QRGI are similar to those by using SVDGI, while the computing time (the time consuming on the light patterns computation) is substantially shortened.     

Wavefront reconstruction algorithm for wavefront sensing based on binary aberration modes

We propose a new algorithm for wavefront sensing based on binary intensity modulation. The algorithm is based on the fact that a wavefront can be expended with a series of orthogonal and binary functions, the Walsh series. We use a spatial light modulator(SLM) to produce different binary-intensity-modulation patterns which are the simple linear transformation of the Walsh series. The optical fields under different binary-intensity-modulation patterns are detected with a photodiode.The relationships between the incident wavefront modulated with the patterns and their optical fields are built to determinate the coefficients of the Walsh series. More detailed and strict relationship equations are established with the algorithm by adding new modulation patterns according to the properties of the Walsh functions. An exact value can be acquired by solving the equations. Finally, with the help of phase unwrapping and smoothing, the wavefront can be reconstructed. The advantage of the algorithm is providing an analytical solution for the coefficients of the Walsh series to reconstruct the wavefront. The simulation experiments are presented and the effectiveness of the algorithm is demonstrated.     

基于遗传优化的西太平洋副高异常活动年份的动力预报模型反演

基于2000–2010年NECR/NECP的500 hPa位势场资料, 用EOF时空分解方法和动力模型重构思想, 通过遗传算法的全局优化搜索和并行计算途径, 开展了500 hPa位势场动力预报模型反演, 建立了刻画副高活动的非线性动力预报模型, 实现了副高活动的中长期预报. 模型预报试验表明, 该模型对副高的中长期活动, 尤其是异常活动能够较好地描述和预报, 进而为副高等复杂天气系统的预报探索了新的方法思路. 关键词： 副热带高压 500 hPa位势场 动力模型反演 遗传算法     

一种基于时空变换和压缩感知的磁共振螺旋采样的图像重建算法

螺旋采样磁共振快速成像在功能性成像、并行成像和动态成像等领域发挥着越来越重要的作用.螺旋采样图像重建的传统算法是用核函数将螺旋状分布的k空间数据插值到均匀网格中,再利用傅里叶变换和最小二乘法进行重建.但是基于网格化的算法对核函数过于依赖,在网格化过程中产生难以避免的误差.该文提出了基于时空变换和压缩感知的l1范数的最优化模型和重建算法.时空变换矩阵描述了空间上的磁共振图像与采集到的时域信号间的关系,使得算法直接使用采集到的数据作为保真约束项,避免了网格化过程产生的误差.此外,基于图像处理单元的并行计算被用来提高时空变换矩阵的运算速度,使得算法具有较强的应用价值.     

印度洋中北部声速剖面结构的时空变化及其物理机理研究

海洋的声速结构对水下声传播有重要影响,在印度洋中北部复杂多变的海洋物理和水文环境中, 获取声速剖面的时空统计分布规律对水下目标探测和水下声通信有重要意义. 由于垂直梯度法在声速结构分析中的局限性及其在印度洋中北部海域的适用性问题, 采用多元统计分析中的最优分割法对声速跃层进行分析,并应用最近10年的地转海洋学实时观测阵 数据对印度洋中北部海域声速剖面的特征量进行了计算,获得了声速跃层的垂直结构特征和时空变化规律; 还利用经验正交函数(EOF)表示方法,分析了印度洋中北部声速剖面拟合精度随EOF阶次的分布特点. 根据印度洋的海洋物理特征,揭示了声速剖面特征量时空演变的内在物理机理.研究结果表明: 最优分割法是适合印度洋声速结构的跃层判断方法,并提出了相应的判断准则参数; 声速剖面拟合精度随EOF阶次变化的区域性分布特征较明显,其季节性变化较小; 印度洋中北部的深海声道轴只在5°S以南明显存在,在15°S---25°S 附近海域存在三个跃层;印度洋中北部声速剖面结构可分为单跃层、双跃层Ⅰ型、 双跃层Ⅱ型和三跃层四种类型以及春夏秋冬四个季节模态. 声速剖面的分析结果对于水声传播和声纳系统的使用具有一定参考意义.     

基于改进变分模态分解的北极海域声速剖面分类*

应用支持向量机对北极声速剖面进行分类,特征量提取是关键。该文采用一种基于经验模态分解的改进变分模态分解算法,以准确提取声速剖面特征量。算法首先对声速剖面信号进行经验模态分解,依据最大类间方差原则划分各分量边际谱主频带,以相似度作为最小分解层数判断标准,获得最小分解层数,进行变分模态分解。对北极区海水声速实测数据(信号)处理表明,该方法可有效提取信号经验模态分解各分量的希尔伯特边际谱特征,进行支持向量机分类,实现对北极海域声速剖面的分类识别,解决以往人工分类耗时久的问题。     

Tracking-positioning of sound speed profiles and moving acoustic source in shallow water

An ensemble Kalman filter(EnKF) approach is proposed to perform sequential tracking of water column sound speed profile(SSP) using a moving acoustic source. First,the SSPs are discretized in depth and range, and are expressed by the empirical orthogonal functions(EOFs). Second, the acoustic source state information and the first three orders of EOF coefficients are expressed as the state variable, and the acoustic field information received by the vertical line array are the measured values. Successively, the state variables and measured values are used to establish the state-measure model. Last, the EnKF is utilized to track the state variables. The simulation results show that the root mean square error of SSP and the absolute error of source are all small, and thus the acoustic source tracking-positioning has high accuracy. Moreover, increasing the number of sample collection, the signal-to-noise ratio and the number of receiving elements can improve the tracking-positioning results. The method is verified using the experimental data of the East China Sea.     

浅海声速剖面与移动声源的跟踪定位

在水平非均匀分布的浅海环境中,针对移动声源跟踪时,声速剖面的变化会对声场产生影响,提出了一种利用集合卡尔曼滤波算法的声速剖面跟踪反演和移动声源跟踪定位的方法。首先,将声速剖面进行距离和深度的参数化表示,从而将对声速剖面的跟踪转化为对声速剖面前3阶经验正交函数系数的跟踪;其次,通过将声源状态信息和声速剖面信息表示为状态变量,而将垂直线列阵接收到的声场信息作为测量值建立状态-测量模型,然后利用集合卡尔曼滤波方法对模型状态变量进行跟踪。仿真结果得出:声速剖面跟踪反演的均方根误差和移动声源跟踪定位的绝对误差都非常小,对声源的跟踪定位精度很高。并且通过增加集合样本数、增加接收信号信噪比以及增加接收阵元数目都可以提高跟踪定位结果精度。最后,利用东海实验数据对本方法进行了验证。      

Near optimal performance joint semi-blind channel estimation and data detection techniques for Alamouti coded single-carrier (SC) MIMO communication systems

In this paper, two novel joint semi-blind channel estimation and data detection techniques are proposed and investigated for Alamouti coded single-carrier (SC) multiple-input multiple-output (MIMO) communication system using Rayleigh flat fading channel model. In the first novel semi-blind technique, blind channel estimation can be performed by using singular value decomposition (SVD) of received output autocorrelation matrix and training based channel estimation for orthogonal training symbols can be performed by using orthogonal pilot maximum likelihood (OPML) algorithm. Further using, that semi-blind channel estimate and received output, data detection is performed by using Maximum likelihood (ML) detection. Finally we derived new training symbols from error covariance matrix of estimated data and known orthogonal training symbols, which further applied to OPML algorithm for final channel estimate. In the second novel semi-blind technique, blind channel estimation can be performed by using matrix triangularization based on householder QR decomposition (H-QRD) of received output autocorrelation matrix instead of SVD decomposition. Other steps are same as the first novel technique to calculate data detection and final channel estimation. Simulation results are presented under 2-PSK, 4-PSK, 8-PSK and 16-QAM data modulation schemes using 2 transmitters and different combinations of receiver antennas to investigate the performances of novel techniques compare to conventional whitening rotation (WR) and rotation optimization maximum likelihood (ROML) based semi-blind channel estimation techniques. Result demonstrates that novel techniques outperform others by achieving near optimal performance.     

基于分类和曲线拟合的干涉超光谱图像压缩

根据干涉超光谱图像的特点,提出了一种基于图像分类与曲线拟合的干涉超光谱图像数据分解算法,结合内嵌比特平面编码技术实现干涉超光谱图像的压缩。与JPEG2000一样,该算法实现了有损、无损压缩的兼容。将干涉超光谱图像数据分为主干涉区域与非主干涉区域两类,针对主干涉区域提出了一种相似匹配算法,而对非主干涉区域采用经验模式分解和二次曲线拟合方法进行数据分析,两种分析算法结合起来能够有效地对谱线数据进行分解,从而有利于取得更好的压缩效果。仿真结果表明,提出的算法可以使无损压缩的输出码率降低0.2~0.4bit/pixel,而近无损、限失真压缩的重建图像质量相应提高。     

稀疏正交联合约束多通道非负矩阵分解声信号分离算法

针对复杂环境下多通道声信号分离问题,提出稀疏正交联合约束多通道非负矩阵分解声信号分离方法。首先设计基于多通道扩展坂仓斋藤(Itakura-Saito,IS)散度的稀疏正交联合约束项构造代价函数,给出信号稀疏和信号正交约束辅助函数,实现代价函数最小化求解。然后通过迭代更新规则设计,得到稀疏正交优化的多通道非负矩阵分解基矩阵和系数矩阵,讨论了稀疏正交约束对基矩阵和系数矩阵稀疏性与连续性影响。最后基于多通道信号空间特性,进行了非负矩阵分解基聚类以获得多通道非负矩阵分解声信号的分离结果。双通道音频数据与四通道声学目标分离实验数据测试表明,对音频数据,所提算法在性能指标信号失真比(SDR)上提高了0.84 dB,对于直升机声源数据,所提算法在SDR上提高了4.53 dB。