带小参数的反应—扩散方程组的数值方法

本文讨论带小参数的反应—扩散方程组的数值方法.由于边界层效应,使得这类问题的数值求解十分困难.我们根据奇异摄动理论和Green函数方法建立起一种适合求解这类问题的差分格式.在文中,我们引入了可行等距度α,并证明了若a≥2则格式在l₁(m)意义下一致收敛且收敛阶为O(h+△t).     

Adomian分解法求解非线性分数阶Fredholm积分微分方程

为了求解一类非线性分数阶Fredholm积分微分方程的数值解,本文将Adomian分解法(Adomian Decomposition Method,ADM)引入到非线性分数阶Fredholm微积分方程的求解中.将ADM多项式与分数阶积分定义有效结合,得到Adomian级数解.通过收敛性分析证明所得的级数解收敛于精确解,给出最大绝对截断误差.并结合实例,证明方法的有效性和实用性.     

富里叶级数的收敛速度

对正弦和余弦富立叶级数,通过合并相邻同号项,使其重排成交错级数.讨论了重排形成的交错级数的敛散性.指出根据自变量x的不同取值,该交错级数可能是单调递减或周期递减的级数.按照莱布尼茨判定法提出了不同精度要求的级数项数的计算公式.选取一到三阶收敛的富立叶级数计算了不同比值精度及差值精度要求的级数项数.计算表明,在x的取值为2π的等分点时,富立叶级数的部分和随项数的增加单调地逼近其收敛值.在x的取值为其它点时,富立叶级数的部分和随项数的增加围绕收敛值上下变动,周期地逼近其收敛值.低收敛阶富立叶级数的收敛速度较慢.要达到0.01%的精度,一收敛阶富立叶级数需要数万项,二收敛阶富立叶级数也需要数百项.在不同计算点处,要达到相同的计算精度,需要的级数项数差别较大.     

一类分数阶Langevin方程block-by-block算法的数值分析

分数阶Langevin方程有重要的科学意义和工程应用价值,基于经典block-by-block算法,求解了一类含有Caputo导数的分数阶Langevin方程的数值解.Block-by-block算法通过引入二次Lagrange基函数插值,构造出逐块收敛的非线性方程组,通过在每一块耦合求得分数阶Langevin方程的数值解.在0<α<1条件下,应用随机Taylor展开证明block-by-block算法是3+α阶收敛的,数值试验表明在不同α和时间步长h取值下,block-by-block算法具有稳定性和收敛性,克服了现有方法求解分数阶Langevin方程速度慢精度低的缺点,表明block-by-block算法求解分数阶Langevin方程是高效的.     

非线性随机分数阶延迟积分微分方程Euler-Maruyama方法的强收敛性

本文研究了一类新的模型问题:非线性随机分数阶延迟积分微分方程.当方程中的漂移项和扩散项满足全局Lipschitz条件和线性增长条件时,基于压缩映射原理给出了该方程解存在唯一的充分条件.由于理论求解的困难,构造了一种数值方法(Euler-Maruyama方法),并证得强收敛阶为α-1/2,α∈(1/2,1].最后通过数值试验,验证了这一理论结果.     

P-adic数域Qp上的级数理论

主要研究了P-adic数域上无穷级数的收敛问题,给出了P-adic数域上数项级数,函数项级数收敛的充要条件,并作了完整的证明.得到了P-adic数域上级数收敛比实数域上级数收敛更容易判断的结论.     

机械求积法及其外推算法求解Helmholtz非线性边界积分方程

当Helmholtz微分方程转化为非线性边界积分方程后,可以利用机械求积法求得近似解,此方法具有较高的收敛精度阶O(h3)和较低的计算复杂度.构造机械求积法时,一个非线性方程系统通过离散非线性积分方程得到.此外,每个矩阵元素的值都不需要计算任何奇异积分.根据渐近紧理论和Stepleman定理,整个系统的稳定性和收敛性得到了证明.利用h3-Richardson外推算法,收敛精度阶可以提高到O(h5).为了求解非线性方程组,利用Ostrowski不动点定理研究了Newton的解的收敛性.几个算例从数值上说明了本算法的有效性.     

用格子Boltzmann方法模拟非线性热传导方程

对具有非线性源项和非线性扩散项的热传导方程建立格子Boltzmann求解模型.在演化方程中增加了两个关于源项分布函数的微分算子,对演化方程实施Chapman-Enskog展开.通过对演化方程的进一步改进,恢复出具有高阶截断误差的宏观方程.对不同参数选取下的非线性热传导方程进行了数值模拟,数值解与精确解吻合得很好.该模型也可以用于同类型的其他偏微分方程的数值计算中.     

求解一维对流方程的高精度紧致差分格式

本文提出数值求解一维对流方程的一种两层隐式紧致差分格式,采用泰勒级数展开法以及对截断误差余项中的三阶导数进行修正的方法对时间和空间导数进行离散.格式的截断误差为O(τ~4+τ~2h~2+h~4),即该格式在时间和空间上均可达到四阶精度.利用von Neumann方法分析得到该格式是无条件稳定的.通过数值实验验证了本文格式的精确性和稳定性.     

时间延迟扩散-波动分数阶微分方程有限差分方法

本文提出求解时间延迟扩散-波动分数阶微分方程有限差分方法，方程中对时间的一阶导函数用α阶（0 < α < 1） Caputo分数阶导数代替.文章中利用Lubich线性多步法对分数阶微分进行差分离散，且文章利用分段区间证明该方法是稳定的，且利用数值实验加以验证.     

Asymptotic and Partial Asymptotic Hankel Operators on $$H^2(\mathbb{D}^n)$$ H 2 ( D n ) )

In this paper, we generalize the concept of asymptotic Hankel operators on H²(D) to the Hardy space H²(Dⁿ) (over polydisk) in terms of asymptotic Hankel and partial asymptotic Hankel operators and investigate some properties in case of its weak and strong convergence. Meanwhile, we introduce i^th-partial Hankel operators on H²(Dⁿ) and obtain a characterization of its compactness for n > 1. Our main results include the containment of Toeplitz algebra in the collection of all strong partial asymptotic Hankel operators on H²(Dⁿ). It is also shown that a Toeplitz operator with symbol φ is asymptotic Hankel if and only if φ is holomorphic function in L^∞(Tⁿ).     

Numerical differentiation for two-dimensional scattered data

In this paper, we propose a regularization method for numerical differentiation of two-dimensional mildly scattered input data. A regularized solution is constructed based on the Green's function. The existence and uniqueness of the regularized solution are proved and the convergence estimates are provided under a simple choice of regularization parameter. Numerical results show that our method is quite effective. One of the advantages for our proposed method is that the basis functions are independent of input data.     

旋波媒质中谱域电并矢Green函数的分解

介绍了一种推导无耗、互易和无界旋波媒质中谱域并矢Green函数表达式的新方法·　这种方法以Hemholtz定理以及并矢Diracδ函数的无散和无旋分解为基础,首先将电矢量的并矢Green函数方程分解成无散电矢量的并矢Green函数方程和无旋电矢量的并矢Green函数方程,然后经Fourier变换导出了旋波媒质中谱域电并矢Green函数的无散分解表达式和无旋分解表达式·　用这种方法推导旋波媒质中并矢Green函数就可以避免必须用波场的分解法和并矢Green函数的本征函数展开法·     

Ewald's Method Revisited: Rapidly Convergent Series Representations of Certain Green's Functions

We propose a justification of Ewald's method of obtaining rapidly convergent series for the Green's function of the 3-dimensional Helmholtz equation. Our point of view enables us to extend the method to Green's functions for the Helmholtz equation in certain domains of ^d with quite general boundary conditions.     

ON THE SHARING VALUES OF ALGEBROID FUNCTIONS AND THEIR DERIVATIVES

In this paper,we define the shared value of an algebroid function and its derivative on its Riemann surface.By considering the relationship between the shared values and the branch points of algebroid functions and their derivatives,we obtain some uniqueness theorems of algebroid functions sharing values with their derivatives,which extend 3 IM shared values theorem of nonconstant meromorphic functions and their derivatives obtained by Mues-Steinmetz and Gundersen.     

基于增广Lagrange函数的RQP方法

Recursive quadratic programming is a family of techniques developd by Bartholomew-Biggs and other authors for solving nonlinear programming problems.This paperdescribes a new method for constrained optimization which obtains its search di-rections from a quadratic programming subproblem based on the well-known aug-mented Lagrangian function.It avoids the penalty parameter to tend to infinity.We employ the Fletcher‘s exact penalty function as a merit function and the use of an approximate directional derivative of the function that avoids the need toevaluate the second order derivatives of the problem functions.We prove that thealgorithm possesses global and superlinear convergence properties.At the sametime, numerical results are reported.     

Carathéodory-Toeplitz and Nehari problems for matrix valued almost periodic functions

In this paper the positive and strictly contractive extension problems for almost periodic matrix functions are treated. We present necessary and sufficient conditions for the existence of extensions in terms of Toeplitz and Hankel operators on Besicovitch spaces and Lebesgue spaces. Furthermore, when a solution exists a special extension (the band extension) is constructed which enjoys a maximum entropy property. A linear fractional parameterization of the set of all extensions is also provided. The techniques used in the proofs include factorizations of matrix valued almost periodic functions and a general algebraic scheme called the band method.

     

一个扩散问题的自然边界元法与有限元法组合

本文讨论由Helmholtz方程描述的扩散问题的自然边界元法与有限元法的组合．取一个圆作为公共边界，用Fourier展开建立边界积分方程，将无界区域上的问题化为有界区域上的非局部边值问题．在变分方程中公共边界上的未知量只包含函数本身而不包含其法向导数，从而减少了未知数的数目，并且边界元剐度矩阵只有极少量不同的元素，有利于数值计算．这种组台方法优越于建立在直接边界元法基础上的组合方法．文中证明了变分解的唯一性，数值解的收敛性和误差估计．最后讨论了数值技术并给出一个算倒．     

铁磁链Landau-Lifshitz方程的显式差分法

正如在研究流体动力学时,Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程起着十分重要的作用一样,在对于非平衡态磁学的研究中,描述连续铁磁体自旋场发展过程的 Ｌａｎｄａｎ－Ｌｉｆｓｈｉｔｚ方程［１］起着十分重要的作用［２］．一九九三年,美国和印度签署了一个大约 ２８０万美元的合作研究计划,在三年的时间里,对Ｌａｎｄａｕ－Ｌｉｆｓｈｉｔｚ铁磁链方程进行研究．在无阻尼的情况下,它为一完全可积的孤立子系统［３,４,５］。很多物理学家研究了它的孤立子解的存在性、逆散射方法以及相互碰撞［３,４,５］．关于解的存在性, Ａｌｏｎ…     

Treatment of singularities in the method of fundamental solutions for two-dimensional Helmholtz-type equations

We investigate a meshless method for the accurate and non-oscillatory solution of problems associated with two-dimensional Helmholtz-type equations in the presence of boundary singularities. The governing equation and boundary conditions are approximated by the method of fundamental solutions (MFS). It is well known that the existence of boundary singularities affects adversely the accuracy and convergence of standard numerical methods. The solutions to such problems and/or their corresponding derivatives may have unbounded values in the vicinity of the singularity. This difficulty is overcome by subtracting from the original MFS solution the corresponding singular functions, without an appreciable increase in the computational effort and at the same time keeping the same MFS approximation. Four examples for both the Helmholtz and the modified Helmholtz equations are carefully investigated and the numerical results presented show an excellent performance of the approach developed.