二维保单调保守恒插值算子

基于一个一维保单调保守恒插值算子,利用不完全双二次插值提出一个二维保单调保守恒插值算子.从插值逼近角度,通过几个数值实验验证该插值算子有效.用得到的二维插值算子作为结构网格自适应加密(structured adaptive mesh refinement,SAMR)算法中的细化插值算子,求解几个二维Euler方程数值例子,结果表明,提出的二维插值算子有效.     

SAMR网格上扩散方程有限体格式的逼近性与两层网格算法

针对结构自适应加密网格(SAMR)上扩散方程的求解,分析几种有限体格式的逼近性,同时设计和分析一种两层网格算法.首先,讨论一种常见的守恒型有限体格式,并给出网格加密区域和细化/粗化插值算子的条件;接着,通过在粗细界面附近引入辅助三角形单元,消除粗细界面处的非协调单元,设计了一种保对称有限体元(SFVE)格式,分析表明,该格式具有更好的逼近性,且对网格加密区域和插值算子的限制更弱;最后,为SFVE格式构造一种两层网格(TL)算法,理论分析和数值实验表明该算法的一致收敛性.     

遗传算法的几点改进及其在脉冲整形试验中的应用

为了解决脉冲整形实验中经常碰到的遗传算法收敛速度慢,早熟等问题,我们对传统的遗传算法进行了几点改进,例如：将两个个体间的欧几里得距离作为判断是否进行交叉操作的判据之一,而不再仅仅依靠个体的适应度值（fitness）,这样能有效地保持种群的基因多样性,提高交叉算子的效率;第二,引入多个交叉算子共同作用于种群. 由于算子的组合效应,共同作用产生的子代适应度值要优于任何一个算子单独作用时产生的子代适应度值. 因而可以产生更大的探索范围,防止算法收敛在某个局部最优解;第三,为了提高收敛速度,我们提出一种新的插值方式：非线性插值,即依据频谱的强度大小决定插值点的密度. 我们初步将此改进算法应用到飞秒整形光路输出光的相位补偿实验中,得到了比较令人满意的结果.      

求解非线性偏微分方程的自适应小波精细积分法

以Burgers方程为例,提出了一种求解偏微分方程的自适应多层插值小波配置法,通过引入一种具有插值特性的拟Shannon小波并利用插值小波理论构造了多层自适应插值小波算子,从而在空间实现了偏微分方程的自适应离散.另外,精细时程积分方法和外推法的引入不但有助于提高求解速度和数值结果的精度,而且使时间积分步长的选取具有了自适应性.     

梯度算子选择对基于梯度的亚像素位移算法的影响

基于梯度的亚像素位移算法中不同的灰度梯度计算方法对计算结果的影响进行了研究,列出了几种常用灰度梯度计算方法,并用对相邻5像素点的离散灰度进行三次样条插值的方法推导出两种灰度梯度算子。然后对计算机模拟生成的散斑图像对用各种梯度算子计算得出的亚像素位移与预先设定的真实值做了比较,给出了相应的均值误差和均值相对误差。最后的真实试件的刚体平移和单向拉伸实验结果表明,Barron算子是所列的几种梯度算法中最为精确、稳定的梯度算子。     

基于L1范数和正交梯度算子的超分辨率重建

针对超分辨率图像重建的病态问题,设计了一种新的自适应超分辨率图像序列重建算法。该算法在L1范数重建框架下,利用金字塔算法与Lucas-Kanade算法相结合的方法实现图像配准,获得亚像素的运动估计;通过引入移位算子给出了基于正交梯度算子的正则项的实现方法,并从自适应的角度选择正则化参数,最后通过最速下降法求解模型的目标泛函最小值。结果表明:对于模拟实验和真实序列实验,该方法相比于样条插值算法、Tikhonov正则化算法、双边全变差重建算法都有一定的优势,能够取得更好的复原效果,并且由于正则项较为简单,重建所需时间相对减少。     

求解静电场问题的新瀑布型代数二重网格法

针对静电场离散化代数方程组,通过选用一种快速的粗化算法,提出一种新的插值算子,构造一种新的瀑布型代数二重网格法。数值实验表明新算法大大减低了计算时间。     

基于Shape Context和尺度不变特征变换的多模图像自动配准方法

提出了基于修正的尺度不变特征变换（SIFT）特征提取和Shape Context特征描述算子相结合的多模图像自动配准算法,该算法利用修正的SIFT算法提取多模图像中的特征点,然后采用Shape Context算子描述特征点,利用特征点周围区域边缘点的梯度方向形成特征向量。采用欧氏距离作为匹配标准对多模图像中特征点进行初始匹配,然后通过RANSAC算法消除误匹配的特征点对,并采用最小二乘法计算仿射变换参数,最后通过仿射变换和双线性插值实现图像配准。对红外图像和可见光图像的配准实验结果表明了本算法的有效性和稳定性。     

利用完全复数极限学习机增强近场声全息空间分辨率

针对稀疏测量阵列条件下近场声全息重建结果空间分辨率不足的问题,提出了一种基于完全复数极限学习机的全息声压插值方法。该方法首先将已测量的全息面复声压和对应的测点坐标组成训练样本输入完全复数极限学习机,接着把插值点的坐标代入训练好的极限学习机,得到相应位置的复声压,实现全息数据的插值。利用插值后的全息数据进行重建,并与不做插值处理的重建结果和传统插值处理后的重建结果比较。仿真和实验结果均表明:与不做插值相比,该方法在不增加传声器的条件下显著提高了重建结果的空间分辨率。与基于支持向量机或传统极限学习机的插值方法相比,该方法速度更快,插值后重建结果精度更高。同时,通过添加噪声干扰验证了该方法的稳健性。      

求解静电场偏微分方程的新代数多重网格法

多重网格法是求解偏微分方程大规模离散化方程最有效的方法,针对静电场偏微分方程,讨论一致线性有限元剖分下的拉格朗日有限元方程的代数多重网格法,给出了一种新的粗化算法和构造插值算子的途径。数值实验表明,新的代数多重网格法的有效性。     

On an interpolation model for the transition operator for Markov and non-Markov processes

A phenomenological interpolation model for the transition operator of a stationary Markov process is shown to be equivalent to the simplest difference approximation in the master equation for the conditional density. Comparison with the formal solution of the Fokker-Planck equation yields a criterion for the choice of the correlation time in the approximate solution. The interpolation model is shown to be form-invariant under an iteration-cum-rescaling scheme. Next, we go beyond Markov processes to find the effective time-development operator (the counterpart of the conditional density) in the following very general situation: the stochastic interruption of the systematic evolution of a variable by an arbitrary stationary sequence of randomizing pulses. Continuous-time random walk theory with a distinct first-waiting-time distribution is used, along with the interpolation model for the transition operator, to obtain the solution. Convenient closed-form expressions for the ‘averaged’ time-development operator and the autocorrelation function are presented in various special cases. These include (i) no systematic evolution, but correlated pulses; (ii) systematic evolution interrupted by an uncorrelated (Poisson) sequence of pulses.     

Sixth order compact scheme combined with multigrid method and extrapolation technique for 2D poisson equation

We develop a sixth order finite difference discretization strategy to solve the two dimensional Poisson equation, which is based on the fourth order compact discretization, multigrid method, Richardson extrapolation technique, and an operator based interpolation scheme. We use multigrid V-Cycle procedure to build our multiscale multigrid algorithm, which is similar to the full multigrid method (FMG). The multigrid computation yields fourth order accurate solution on both the fine grid and the coarse grid. A sixth order accurate coarse grid solution is computed by using the Richardson extrapolation technique. Then we apply our operator based interpolation scheme to compute sixth order accurate solution on the fine grid. Numerical experiments are conducted to show the solution accuracy and the computational efficiency of our new method, compared to Sun–Zhang’s sixth order Richardson extrapolation compact (REC) discretization strategy using Alternating Direction Implicit (ADI) method and the standard fourth order compact difference (FOC) scheme using a multigrid method.     

基于CCD的无接触式的物体测量

基于CCD的物体测量平台，首先对物体图像进行预处理，继而利用Canny算子提取物体边缘，再通过本文改进的Hough变换实现对被测物体边界直线检测和重构，最后通过多项式插值算法实现物体像素尺寸到物体实际尺寸的转换．实验证明，此方法能够实现无接触式物体尺寸较高精度的测量．     

Generalized form of the direct and adjoint radiative transfer equations

The main goal of this paper is to give a rigorous derivation of the generalized form of the direct (also referenced as forward) and adjoint radiative transfer equations. The obtained expressions coincide with expressions derived by Ustinov [Adjoint sensitivity analysis of radiative transfer equation: temperature and gas mixing ratio weighting functions for remote sensing of scattering atmospheres in thermal IR. JQSRT 2001;68:195-211]. However, in contrast to [Ustinov EA. Adjoint sensitivity analysis of radiative transfer equation: temperature and gas mixing ratio weighting functions for remote sensing of scattering atmospheres in thermal IR. JQSRT 2001;68:195-211] we formulate the generalized form of the direct radiative transfer operator fully independent from its adjoint. To illustrate the application of the derived adjoint radiative transfer operator we consider the angular interpolation problem in the framework of the discrete ordinate method widely used to solve the radiative transfer equation. It is shown that under certain conditions the usage of the solution of the adjoint radiative transfer equation for the angular interpolation of the intensity can be computationally more efficient than the commonly used source function integration technique.     

Interpolation theory and the Wigner-Yanase-Dyson-Lieb concavity

The Wigner-Yanase-Dyson-Lieb concavity is naturally captured in the frame of interpolation theory. Among other results, a certain generalization (involving operator monotone functions) of this concavity in the context of general von Neumann algebras is obtained. Also, a close relationship between the above subjects and F. Hansen's inequality is clarified. All results are proved by using simple variational expressions of involved quantities.Supported in part by the National Science Foundation, grant number MCS-8102158     

Coefficient Jump-Independent Approximation of the Conforming and Nonconforming Finite Element Solutions

A counterexample is constructed. It confirms that the error of conforming finite element solution is proportional to the coefficient jump, when solving interface elliptic equations. The Scott-Zhang operator is applied to a nonconforming finite element. It is shown that the nonconforming finite element provides the optimal order approximation in interpolation, in $L^2$-projection, and in solving elliptic differential equation, independent of the coefficient jump in the elliptic differential equation. Numerical tests confirm the theoretical finding.     

Ising Correlations and Elliptic Determinants

Correlation functions of the two-dimensional Ising model on the periodic lattice can be expressed in terms of form factors—matrix elements of the spin operator in the basis of common eigenstates of the transfer matrix and translation operator. Free-fermion structure of the model implies that any multiparticle form factor is given by the pfaffian of a matrix constructed from the two-particle ones. Crossed two-particle form factors can be obtained by inverting a block of the matrix of linear transformation induced on fermions by the spin conjugation. We show that the corresponding matrix is of elliptic Cauchy type and use this observation to solve the inversion problem explicitly. Non-crossed two-particle form factors are then obtained using theta functional interpolation formulas. This gives a new simple proof of the factorized formulas for periodic Ising form factors, conjectured by A. Bugrij and one of the authors.