带插值的自适应网格重构算法求解带移动热源的反应扩散方程的理论分析

本文研究一个带插值的网格重构算法求解一类带移动热源的反应扩散方程. 算法包括两步: 第一步是用旧时间网层上的计算解计算新时间层上的空间网格; 第二步是使用有限差分方法在新时间层 空间网格上离散方程, 并且将旧时间层上计算解的插值作为初始值. 对于时间, 我们获得了一阶收敛结果. 对于空间, 我们证明了使用线性插值算法的一阶收敛性和使用二次插值算法的二阶收敛性. 数值例子肯定了本文的理论结果.     

闭凸集约束下线性矩阵方程求解的松弛交替投影算法

研究线性矩阵方程AXB=C在闭凸集合R约束下的数值迭代解法.所考虑的闭凸集合R为(1)有界矩阵集合,(2)Q-正定矩阵集合和(3)矩阵不等式解集合.构造松弛交替投影算法求解上述问题,并用算子理论证明了由该算法生成的序列具有弱收敛性.给出了矩阵方程AXB=C求对称非负解和对称半正定解的数值算例,大量数值实验验证了该算法的可行性和高效性,并说明该算法与交替投影算法和谱投影梯度算法比较在迭代效率上的明显优势.     

高波数波动问题的多水平方法

本文主要讨论求解高波数Helmholtz方程的多水平方法,主要回顾了一些具有代表性的多重网格方法.如Erlangga等人的shifted Laplacian预处理的多重网格法;Elman等提出的修正的多重网格方法;以及我们的基于连续内罚有限元(CIP-FEM)离散代数系统的多水平算法.最后还介绍了求解高波数时谐Maxwell方程的CIP-FEM离散代数系统的多水平算法.     

用自适应MG法,FMG法求解Helmholtz方程

本文针对Ｖ循环、Ｗ循环和多重网格法中最优光滑次数及循环体个数难以确定的缺点，以Ｈｅｌｍｈｏｌｔｚ方程为例给出自适应的多重网格算法和自适应的完全多重网格算法。     

带移动奇异源热方程的移动网格方法

研究了一个带若干奇异源热方程的数值求解,其源的移动由一个常微分方程描述.基于移动观察区域和区域分解思想提出了一个移动网格预估校正算法.网格方程可自然的通过并行高效求解,算法避免了跳跃信息[u]的计算而使物理方程的离散格式变得非常简单,且仍保持了空间上的二阶收敛性.数值例子验证了算法的收敛性和高效性,并模拟了非线性源函数带来的爆破现象.     

一类加速的模系对称超松弛迭代方法定价双资产美式期权

构造和研究了一类加速的模系对称超松弛迭代方法,用来求解由双资产美式期权定价模型离散出来的线性互补问题.理论分析给出该算法的收敛性条件.数值实验表明,该方法对于求解双资产美式期权定价模型是有效的,并且优于经典的模系超松弛迭代方法和模系对称超松弛迭代方法.     

混合方程的多重网格算法

本文给出了混合方程的多重网格算法以及收敛性的证明。为简单起见,我们用Jacobi迭代作为光滑过程。     

非凸非光滑不可分离优化的线性对称邻近ADMM收敛性分析

交替方向乘子法(ADMM)是一种求解可分离优化问题的简单有效的方法,相关研究已经较为完善.然而,当目标函数存在耦合项时,对ADMM算法收敛性的研究还处于初期.文章针对非凸非光滑不可分离优化问题,基于对称交替方向乘子法(SADMM),结合线性化技术,提出了一种新的线性对称邻近ADMM.在一定的假设条件下,证明了算法生成的序列有界并收敛至增广拉格朗日函数的稳定点.其次,当辅助函数满足Kurdyka-Lojasiewicz性质时,证明了算法的强收敛性.最后,数值实验的结果表明了算法的有效性.     

非协调Wilson有限元的多重网格方法

§1.引言 非协调Wilson有限元[1—3]对解弹性力学方程有实用价值,在工程上有用。本文分析Wilson元的多重网格法,给出用多重网格方法求得的近似解按L~2模和能量模的最佳收敛阶误差估计。对于W-循环,可以证明其计算量与离散空间的维数为同一量级O(N_k)。 考虑二阶椭圆Dirchlet边值问题:     

大Reynolds数Navier-Stokes方程带回溯技巧的两水平方法

本文考虑了一种求解大Reynolds数定常Navier-Stokes方程带回溯(backtracking)技巧的两水平有限元方法.其基本思想是,首先在一粗网格上求解带有亚格子模型稳定项的Navier-Stokes方程,然后在细网格上求解一个亚格子模型稳定化的线性Newton问题,最后又回到粗网格上求解线性化的校正问题.通过适当的稳定化参数和粗细网格尺寸的选取,本文的算法能取得最优渐近收敛阶.数值实验检验了理论分析的正确性和算法的有效性.     

用常微分方程模型分析预防和隔离措施对SARS发病率的影响

通过建立常微分方程模型 ,分析了预防和隔离措施对 SARS发病率的影响 ,并把计算结果与实际统计数据进行了比较 ,结果表明 ,及时高效的预防和隔离措施能够有效地控制 SARS的传播 .     

On the distributions of the ratios of the extreme roots to the trace of the Wishart matrix

In this paper, the authors cosider the derivation of the exact distributions of the ratios of the extreme roots to the trace of the Wishart matrix. Also, exact percentage points of these distributions are given and their applications are discussed.     

On the dependence of the growth rate on the length of the defining relator

Let { } be a sequence of finitely presented groups with generating setA={a1, …, am}, and letRk be the symmetrized set of words over the alphabetA∪A−1 obtained from the defining words and their inverses by all cyclic shifts. We shall assume that the words inRk are cyclically irreducible, and their lengths tend to ∞ ask increases. In the paper, it is proved that ifRk satisfies the small cancellation conditionC'(1/6) and the number of relators increases not very rapidly with increasingk, then the growth rate ψ(Gk) tends to 2m−1 ask→∞. Translated fromMatematicheskie Zametki, Vol. 65, No. 4, pp. 611–617, April, 1999.     

On the stability of the solution of the thermistor problem

The stability of the stationary solution of the thermistor problem 1s proved using a Liapunov functional for a class of physically relevant electrical conductivity.     

宏观因素影响下的系统中元件重要性研究

为研究复杂系统在工作环境中其组成元件对系统安全运行的重要性,将汪培庄先生的因素空间理论与笔者提出的空间事故树理论相结合,构造了一套元件重要性研究方法.构建系统T={U,C,D},将元件作为研究对象集合U,系统工作的宏观环境作为因素集C,元件重要性排序集作为D.对宏观环境中的工作时间a1和温度a_2进行划分形成不同的状态区域S_q,计算在S_q中元件xj的失效权重γ(AS_q(x_j))和在S_q中系统T的失效权重δ(AS_q(T))),从而得到x_j在S状态下的等效失效权重Z(AS_q(x_j)),研究状态S_q下的原件重要性排序D_η,及元件x_j失效性对a_1及a_2的敏感性.使用一个实际的电气系统维修情况统计资料,使用上述方法进行了研究,结果表明:不同工作环境下元件对系统的重要程度是不同的.元件对温度和使用时间是敏感的,并得到了在10³0°且5030°且50⁷5d环境下工作系统可靠性是最高的结论.在给定工作环境下,重要性大的元件多储备,重要性小的元件少储备,以满足系统维修需要,并指导实际工程.     

一般正则泛函的W-极小的正则性

本文讨论了一般的正则泛函：Ｆ（ｕ；Ω）＝∫_{Ωｆ（ｘ，ｕ，Ｄｕ）ｄｘ的局部Ｗ－极小的Ｃ^１，α正则性．获得了处理Ｗ－极小ｕ的Ｈ?ｌｄｅｒ连续的指数估计．}     

On the Decomposition of the Riesz Operator and the Expansion of the Riesz Semigroup

For a Riesz operator T on a reflexive Banach space X with nonzero eigenvalues denote by E(λ_i; T) the eigen-projection corresponding to an eigenvalue λ_i. In this paper we will show that if the operator sequence is uniformly bounded, then the Riesz operator T can be decomposed into the sum of two operators T_p and T_r: T = T_p + T_r, where T_p is the weak limit of T_n and T_r is quasi-nilpotent. The result is used to obtain an expansion of a Riesz semigroup T(t) for t ≥ τ. As an application, we consider the solution of transport equation on a bounded convex body.     

On the relation between the A-polynomial and the Jones polynomial

This paper shows that the noncommutative generalization of the A-polynomial of a knot, defined using Kauffman bracket skein modules, together with finitely many colored Jones polynomials, determines the remaining colored Jones polynomials of the knot. It also shows that under certain conditions, satisfied for example by the unknot and the trefoil knot, the noncommutative generalization of the A-polynomial determines all colored Jones polynomials of the knot.

     

On the rates of the other law of the logarithm

Let X, X1 , X2 , . . . be i.i.d. random variables, and set Sn = X1 +···+Xn , Mn = maxk≤n |Sk|, n ≥1. Let an = o( (n)(1/2)/logn). By using the strong approximation, we prove that, if EX = 0, VarX = σ2 0 and E|X| 2+ε ∞ for some ε 0, then for any r 1, lim ε1/(r-1)(1/2) [ε-2-(r-1)]∞∑n=1 nr-2 P{Mn ≤εσ (π2n/(8log n))(1/2) + an } = 4/π . We also show that the widest a n is o( n(1/2)/logn).     

关于循环子半群的结构与数量问题及拟环的特征与结构

彻底解决了所有循环半群及其子群的结构和数量问题,并讨论了拟群分解问题,同时,对群论基本定理作了部分推广,并给出了定理的另一部分不可推广的反例,最后,建立了一类特殊环－拟环。