Interpolating solutions of the Poisson equation in the disk based on Radon projections

We consider an algebraic method for reconstruction of a function satisfying the Poisson equation with a polynomial right-hand side in the unit disk. The given data, besides the right-hand side, is assumed to be in the form of a finite number of values of Radon projections of the unknown function. We first homogenize the problem by finding a polynomial which satisfies the given Poisson equation. This leads to an interpolation problem for a harmonic function, which we solve in the space of harmonic polynomials using a previously established method. For the special case where the Radon projections are taken along chords that form a regular convex polygon, we extend the error estimates from the harmonic case to this Poisson problem. Finally we give some numerical examples.     

A new approach to constant term identities and Selberg-type integrals

Selberg-type integrals that can be turned into constant term identities for Laurent polynomials arise naturally in conjunction with random matrix models in statistical mechanics. Built on a recent idea of Karasev and Petrov we develop a general interpolation based method that is powerful enough to establish many such identities in a simple manner. The main consequence is the proof of a conjecture of Forrester related to the Calogero–Sutherland model. In fact we prove a more general theorem, which includes Aomoto's constant term identity at the same time. We also demonstrate the relevance of the method in additive combinatorics.     

Convergence of FFT‐based homogenization for strongly heterogeneous media

The FFT‐based homogenization method of Moulinec–Suquet has recently attracted attention because of its wide range of applicability and short computational time. In this article, we deduce an optimal a priori error estimate for the homogenization method of Moulinec–Suquet, which can be interpreted as a spectral collocation method. Such methods are well‐known to converge for sufficiently smooth coefficients. We extend this result to rough coefficients. More precisely, we prove convergence of the fields involved for Riemann‐integrable coercive coefficients without the need for an a priori regularization. We show that our L² estimates are optimal and extend to mildly nonlinear situations and L^p estimates for p in the vicinity of 2. The results carry over to the case of scalar elliptic and curl ? curl‐type equations, encountered, for instance, in stationary electromagnetism. Copyright © 2014 John Wiley & Sons, Ltd.     

Acceleration of self‐consistent field convergence in ab initio molecular dynamics simulation with multiconfigurational wave function

The Lagrange interpolation of molecular orbital (LIMO) method, which reduces the number of self‐consistent field iterations in ab initio molecular dynamics simulations with the Hartree–Fock method and the Kohn–Sham density functional theories, is extended to the theory of multiconfigurational wave functions. We examine two types of treatments for the active orbitals that are partially occupied. The first treatment, as denoted by LIMO(C), is a simple application of the conventional LIMO method to the union of the inactive core and the active orbitals. The second, as denoted by LIMO(S), separately treats the inactive core and the active orbitals. Numerical tests to compare the two treatments clarify that LIMO(S) is superior to LIMO(C). Further applications of LIMO(S) to various systems demonstrate its effectiveness and robustness. © 2014 Wiley Periodicals, Inc.     

Wavelet multiresolution interpolation Galerkin method for nonlinear boundary value problems with localized steep gradients

The wavelet multiresolution interpolation for continuous functions defined on a finite interval is developed in this study by using a simple alternative of transformation matrix. The wavelet multiresolution interpolation Galerkin method that applies this interpolation to represent the unknown function and nonlinear terms independently is proposed to solve the boundary value problems with the mixed Dirichlet-Robin boundary conditions and various nonlinearities, including transcendental ones, in which the discretization process is as simple as that in solving linear problems, and only common two-term connection coefficients are needed. All matrices are independent of unknown node values and lead to high efficiency in the calculation of the residual and Jacobian matrices needed in Newton’s method, which does not require numerical integration in the resulting nonlinear discrete system. The validity of the proposed method is examined through several nonlinear problems with interior or boundary layers. The results demonstrate that the proposed wavelet method shows excellent accuracy and stability against nonuniform grids, and high resolution of localized steep gradients can be achieved by using local refined multiresolution grids. In addition, Newton’s method converges rapidly in solving the nonlinear discrete system created by the proposed wavelet method, including the initial guess far from real solutions.

     

闭区间[0,1]上的分形插值函数的分数阶积分

介绍了分形插值函数和迭代函数系统以及v阶黎曼-刘准尔分数阶积分的概念及相关定理.利用这些概念及定理讨论了分形插值函数的分数阶积分在[0,1]上连续性及判定[0,1]上的分形插值函数的分数阶积分也是[0,1]上的分形插值函数,并给予了证明.     

一种基于谱线幅值归零的加窗插值FFT谐波检测方法

提出了一种基于谱线幅值归零的加窗插值FFT谐波检测方法用于分析复杂电网信号中的各次谐波.该检测方法将每一次在离散频谱中找到的最大谱线及其左右两侧的数根谱线幅值归零,直到各次谐波分析完毕.通过仿真实验分析了在基波频率稳定时对模拟信号的测量精度,研究了基波频率波动对谐波分析的影响.通过对EAST实验数据的分析,验证了该谐波...     

构造低次有理插值函数的一种方法

关于有理插值的算法已有很多[1,4,5],受二元多项式插值迭加算法[6]的启发,我们给出一种简便的求低次有理插值函数的方法,同时给出有理插值函数存在的充分条件,便于检验.所给方法具有可操作性和实际应用价值,且具有较好的灵活性.     

优化端点条件的平面二次均匀B样条插值曲线

在利用反求法构造B样条插值曲线时,往往需要选取端点条件。 因此,可对端点条件进行优化选取,使得构造的B样条插值曲线满足特定要求。提出了一种利用曲线内能极小选取平面二次均匀B样条插值曲线端点条件的算法。首先给出了二次均匀B样条插值曲线分控制顶点与首个控制顶点（即端点条件）的递推关系式;然后给出了利用曲线内能极小优化选取首个控制顶点的算法,证明了利用该算法构造的C¹连续二次均匀B样条插值曲线为保形插值,并通过数值算例证明了算法的有效性;最后,为便于实际应用,基于MATLAB平台设计了算法所对应的图形用户界面,用户通过简单的操作即可获得光顺的C¹连续二次均匀B样条保形插值曲线。     

有耗介质层上多导体传输线的电磁耦合时域分析方法

目前,针对空间电磁场作用有耗介质层上传输线的电磁耦合,仍缺乏有效的数值分析方法.因此,本文提出一种高效的时域混合算法,很好地解决了有耗介质层上传输线电磁耦合建模难的问题.首先,对经典传输线方程进行改进,推导了适用于有耗介质层上多导体传输线电磁耦合分析的修正传输线方程.然后,结合时域有限差分方法和相应插值技术,求解修正传输线方程,获得多导线及其端接负载上的电压和电流响应,并实现空间电磁场辐射与多导线瞬态响应的同步计算.最后,通过相应计算实例的数值模拟,与CST软件的仿真结果进行对比,验证了时域混合算法的正确性和高效性.