Optimal variable shape parameter for multiquadric based RBF-FD method

In this follow up paper to our previous study in Bayona et al. (2011) [2], we present a new technique to compute the solution of PDEs with the multiquadric based RBF finite difference method (RBF-FD) using an optimal node dependent variable value of the shape parameter. This optimal value is chosen so that, to leading order, the local approximation error of the RBF-FD formulas is zero. In our previous paper (Bayona et al., 2011) [2] we considered the case of an optimal (constant) value of the shape parameter for all the nodes. Our new results show that, if one allows the shape parameter to be different at each grid point of the domain, one may obtain very significant accuracy improvements with a simple and inexpensive numerical technique. We analyze the same examples studied in Bayona et al. (2011) [2], both with structured and unstructured grids, and compare our new results with those obtained previously. We also find that, if there are a significant number of nodes for which no optimal value of the shape parameter exists, then the improvement in accuracy deteriorates significantly. In those cases, we use generalized multiquadrics as RBFs and choose the exponent of the multiquadric at each node to assure the existence of an optimal variable shape parameter.     

RBF-FD formulas and convergence properties

The local RBF is becoming increasingly popular as an alternative to the global version that suffers from ill-conditioning. In this paper, we study analytically the convergence behavior of the local RBF method as a function of the number of nodes employed in the scheme, the nodal distance, and the shape parameter. We derive exact formulas for the first and second derivatives in one dimension, and for the Laplacian in two dimensions. Using these formulas we compute Taylor expansions for the error. From this analysis, we find that there is an optimal value of the shape parameter for which the error is minimum. This optimal parameter is independent of the nodal distance. Our theoretical results are corroborated by numerical experiments.     

A numerical method for one-dimensional nonlinear sine-Gordon equation using multiquadric quasi-interpolation

In this paper, we use a univariate multiquadric quasi-interpolation scheme to solve the one-dimensional nonlinear sine-Gordon equation that is related to many physical phenomena. We obtain a numerical scheme by using the derivative of the quasi-interpolation to approximate the spatial derivative and a difference scheme to approximate the temporal derivative. The advantage of the obtained scheme is that the algorithm is very simple so that it is very easy to implement. The results of numerical experiments are presented and compared with analytical solutions to confirm the good accuracy of the presented scheme.     

Optimal switching policy for performance enhancement of distributed parameter systems based on event-driven control

This paper aims to improve the performance of a class of distributed parameter systems for the optimal switching of actuators and controllers based on event-driven control. It is assumed that in the available multiple actuators, only one actuator can receive the control signal and be activated over an unfixed time interval, and the other actuators keep dormant.After incorporating a state observer into the event generator, the event-driven control loop and the minimum inter-event time are ultimately bounded. Based on the event-driven state feedback control, the time intervals of unfixed length can be obtained. The optimal switching policy is based on finite horizon linear quadratic optimal control at the beginning of each time subinterval. A simulation example demonstrate the effectiveness of the proposed policy.     

Group-theoretical method for physical property tensors of quasicrystals

In addition to the phonon variable there is the phason variable in hydrodynamics for quasicrystals. These two kinds of hydrodynamic variables have different transformation properties. The phonon variable transforms under the vector representation, whereas the phason variable transforms under another related representation. Thus, a basis (or a set of basis functions) in the representation space should include such two kinds of variables. This makes it more difficult to determine the physical property tensors of quasicrystals. In this paper the group-theoretical method is given to determine the physical property tensors of quasicrystals. As an illustration of this method we calculate the third-order elasticity tensors of quasicrystals with five-fold symmetry by means of basis functions. It follows that the linear phonon elasticity is isotropic, but the nonlinear phonon elasticity is anisotropic for pentagonal quasicrystals. Meanwhile, the basis functions are constructed for all noncrystallographic point groups of quasicrystals.     

Application of quantitative structure-activity relationship to the determination of binding constant based on fluorescence quenching

Quantitative structure-activity relationship (QSAR) model was used to predict and explain binding constant (log K) determined by fluorescence quenching. This method allowed us to predict binding constants of a variety of compounds with human serum albumin (HSA) based on their structures alone. Stepwise multiple linear regression (MLR) and nonlinear radial basis function neural network (RBFNN) were performed to build the models. The statistical parameters provided by the MLR model (R²=0.8521, RMS=0.2678) indicated satisfactory stability and predictive ability while the RBFNN predictive ability is somewhat superior (R²=0.9245, RMS=0.1736). The proposed models were used to predict the binding constants of two bioactive components in traditional Chinese medicines (isoimperatorin and chrysophanol) whose experimental results were obtained in our laboratory and the predicted results were in good agreement with the experimental results. This QSAR approach can contribute to a better understanding of structural factors of the compounds responsible for drug-protein interactions, and can be useful in predicting the binding constants of other compounds.     

求解Kuramoto-Sivashinsky方程的平移基无单元Galerkin方法

Kuramoto-Sivashinsky方程是一种可以描述复杂混沌现象的高阶非线性演化方程.方程中高阶导数项的存在,使得传统无单元Galerkin方法采用高次多项式基函数构造形函数时,形函数违背了一致性条件.因此,本文提出了一种采用平移多项式基函数的无单元Galerkin方法.与传统无单元Galerkin方法相比,该方法在方程离散时依然采用Galerkin进行离散,但形函数的构造采用了基于平移多项式基函数的移动最小二乘近似.通过对具有行波解和混沌现象的Kuramoto-Sivashinsky方程的数值模拟,验证了本文方法的有效性.     

超声波法测量大气折射率结构常数

 用AMK-02型超声波大气参数综合测量仪同时测量了常规气象参数和大气折射率结构常数C²_n,分析并找出了该仪器给出的C²_n在转换时刻与热丝温度脉动仪的测量结果差异较大的原因,是由于超声波法测量温度脉动精度太低,只能达到0.01 K。同时建立了另一种适用于超声波仪计算C²_n的方法,即利用该仪器准确的常规气象参数计算C²_n,其计算结果与热丝温度脉动仪的测量结果吻合较好。     

一种基于参数标定的投影条纹周期简易校正方法

基于三角测量原理的条纹投影轮廓测量系统中,倾斜投影到参考平面上的条纹将产生周期展宽现象,引起相位失真甚至影响测量精度.论文以条纹位置为控制变量推导出条纹周期校正的线性数学模型,通过简便的标定获得模型参数,由此反算出新的待投影条纹,并在参考平面上获得周期分布的投影条纹.实验结果表明校正后的条纹周期变化范围在±0.1像素内...     

势问题的径向基函数——局部边界积分方程方法

将基于径向基函数构造的具有插值特性的近似函数和局部边界积分方程方法相结合，建立了求解势问题的径向基函数——局部边界积分方程方法，推导了相应离散方程.与其他边界积分方程的无网格方法相比，本文方法具有数值实现过程简单、计算量小、精度高的优点，并可直接施加边界条件.最后通过算例说明了该方法的有效性. 关键词： 径向基函数 无网格方法 局部边界积分方程 势问题     

Variational piecewise constant level set methods for shape optimization of a two-density drum

We apply the piecewise constant level set method to a class of eigenvalue related two-phase shape optimization problems. Based on the augmented Lagrangian method and the Lagrange multiplier approach, we propose three effective variational methods for the constrained optimization problem. The corresponding gradient-type algorithms are detailed. The first Uzawa-type algorithm having applied to shape optimization in the literature is proven to be effective for our model, but it lacks stability and accuracy in satisfying the geometry constraint during the iteration. The two other novel algorithms we propose can overcome this limitation and satisfy the geometry constraint very accurately at each iteration. Moreover, they are both highly initial independent and more robust than the first algorithm. Without penalty parameters, the last projection Lagrangian algorithm has less severe restriction on the time step than the first two algorithms. Numerical results for various instances are presented and compared with those obtained by level set methods. The comparisons show effectiveness, efficiency and robustness of our methods. We expect our promising algorithms to be applied to other shape optimization and multiphase problems.     

Spline based least squares integration for two-dimensional shape or wavefront reconstruction

In this work, we present a novel method to handle two-dimensional shape or wavefront reconstruction from its slopes. The proposed integration method employs splines to fit the measured slope data with piecewise polynomials and uses the analytical polynomial functions to represent the height changes in a lateral spacing with the pre-determined spline coefficients. The linear least squares method is applied to estimate the height or wavefront as a final result. Numerical simulations verify that the proposed method has less algorithm errors than two other existing methods used for comparison. Especially at the boundaries, the proposed method has better performance. The noise influence is studied by adding white Gaussian noise to the slope data. Experimental data from phase measuring deflectometry are tested to demonstrate the feasibility of the new method in a practical measurement.     

The EPS method: A new method for constructing pseudospectral derivative operators

We develop a new type of derivative matrix for pseudospectral methods. The norm of these matrices grows at the optimal rate O(N²) for N-by-N matrices, in contrast to standard pseudospectral constructions that result in O(N⁴) growth of the norm. The smaller norm has a big advantage when using the derivative matrix for solving time dependent problems such as wave propagation. The construction is based on representing the derivative operator as an integral kernel, and does not rely on the interpolating polynomials. In particular, we construct second derivative matrices that incorporate Dirichlet or Neumann boundary conditions on an interval and on the disk, but the method can be used to construct a wide variety of commonly used operators for solving PDEs and integral equations. The construction can be used with any quadrature, including traditional Gauss–Legendre quadratures, but we have found that by using quadratures based on prolate spheroidal wave functions, we can achieve a near optimal sampling rate close to two points per wavelength, even for non-periodic problems. We provide numerical results for the new construction and demonstrate that the construction achieves similar or better accuracy than traditional pseudospectral derivative matrices, while resulting in a norm that is orders of magnitude smaller than the standard construction. To demonstrate the advantage of the new construction, we apply the method for solving the wave equation in constant and discontinuous media and for solving PDEs on the unit disk. We also present two compression algorithms for applying the derivative matrices in O(N log N) operations.     

An improved method for shape measurement using two-dimensional digital image correlation

Shape measurement is a significant application of digital image correlation (DIC). An improved method that combines a rotatable plane mirror is proposed to measure the shape of an immovable object. In this method, two images, one before and the other after rotating the plane mirror, are obtained and then in-plane translation which related to the shape of the detected object can be calculated by the use of two-dimensional digital image correlation (2D DIC). The relationship between the in-plane translation and the shape of the object is described. Experimental results show that the proposed method is feasible for shape and distance measurement with high accuracy.     

Photoacoustic method for the measurement of the nonlinearity parameter of liquids

I.IntroductionThenonlinearityparameterB/Aisanimportantacousticalparameteroffluids.Hitherto,manytheoreticalandexperimentalstudiesforthedeterminahonofthenon1inearityparamctersof1iquidshavebeenrcported.ThemcthodsfordetCrminationofthenonlinearityparamenterofthe1iquidcanbeseparatedintotwocategories:thermodynamicmethod[llandfinite-amp1itudeacousticwavcmethod.The1attercancseparatedfurtherintoharmonicwavemcthodt2]andpulsemethod.Karabutovetal.[31andBozkhovetaI.[4lusedthepulsedlasertogencratetheplane…     

一种基于激光辐照热效应的薄膜参数反演方法

本文研究并建立了一种基于激光辐照热效应的薄膜参数反演方法.首先给出激光辐照薄膜产生温升问题的热传导理论模型,并利用拉普拉斯变换得到了膜层和基底温度场的解析解;然后以膜层和基底的导热系数为反演参数,基于非线性共轭梯度算法给出反演基本原理及流程,并推导得到了反演过程中灵敏度系数的解析表达式;以aluminum,silver,copper和gold四种金属薄膜为例,通过与有限元法的计算结果对比验证了温度场解析解的正确性;最后结合四种金属薄膜进行了参数反演,通过考察分析不同随机噪声等条件下的参数反演结果,验证了本文方法在薄膜参数反演精度与反演效率等方面的有效性.反演结果显示:本文方法具有较高的反演精度和效率,在迭代截止误差为10~(-7)时只需用少于20次迭代就能收敛;在测量数据中加入的随机噪声越小,反演的迭代收敛次数就越少,即使是在迭代初值与反演结果相差较大时,用包含5%随机噪声的测量数据反演也能快速收敛.本文提出的薄膜参数反演方法不仅适用于反演导热系数,也可扩展用于反演膜层反射系数或吸收率等参数,具有一定的适用性.本文方法对于激光加工或激光损伤过程中的参数反演及优化具有一定的指导意义.     

一种基于混沌系统部分序列参数辨识的混沌保密通信方法

本文提出一种混沌保密通信方法,即混沌系统的部分序列用于混沌系统参数辨识其他序列用于通信保密.利用混沌蚁群优化算法对部分序列混沌系统进行参数辨识,以达到了解混沌系统全部信息的目的.在参数辨识过程中引入参数空间和蚁群空间,通过空间变换函数使参数空间与蚁群空间之间相互变换.文中使用Lorenz系统进行数值试验,其结果验证混沌系统部分序列参数辨识及混沌保密通信的可行性.     

A 3D shape retrieval method for orthogonal fringe projection based on a combination of variational image decomposition and variational mode decomposition

The orthogonal fringe projection technique has as wide as long practical application nowadays. In this paper, we propose a 3D shape retrieval method for orthogonal composite fringe projection based on a combination of variational image decomposition (VID) and variational mode decomposition (VMD). We propose a new image decomposition model to extract the orthogonal fringe. Then we introduce the VMD method to separate the horizontal and vertical fringe from the orthogonal fringe. Lastly, the 3D shape information is obtained by the differential 3D shape retrieval method (D3D). We test the proposed method on a simulated pattern and two actual objects with edges or abrupt changes in height, and compare with the recent, related and advanced differential 3D shape retrieval method (D3D) in terms of both quantitative evaluation and visual quality. The experimental results have demonstrated the validity of the proposed method.     

三维复杂温度场高精度声学测量方法

声学温度场检测技术通过多路径声波传播时间数据,反演被测区域的温度分布。提供了一种高精度的三维复杂温度场的声学测量方法。首先从射线声学角度给出了三维非均匀温度场中声波传播路径的数学模型。在此基础上,将三维温度场的重建问题转化为声波传播路径的求解和温度场的反演问题,建立了基于多项式修正径向基函数(RBF-PR)和改进的Tikhonov正则化三维温度场重建算法。采用两种典型的炉膛三维温度场模型,在信噪比SNR=35 dB下进行了数值模拟,分析了声波传播路径在非均匀温度场中的弯曲特性、算法的重建质量和抗噪性,同时进行了实际炉膛内二维温度场的重建。结果表明了提出的考虑声线弯曲的温度场重建算法具有精度高,抗噪性强、适用性好的特点,为实现高精度的复杂温度场的声学测量提供了有效方法。     

一种基于遗传算法的混沌系统参数估计方法

通过构造一个适当的适应度函数,将混沌系统的参数估计问题转化为一个参数的寻优问题,然后利用遗传算法的全局优化搜索能力对其进行求解.以典型的Lorenz混沌系统为例进行了数值模拟.实际数值模拟表明,使用这种方法可以有效地对混沌系统的参数进行估计 关键词： 混沌系统 参数估计 遗传算法