New reals: Can live with them,can live without them

We give a self‐contained proof of the preservation theorem for proper countable support iterations known as “tools‐preservation”, “Case A” or “first preservation theorem” in the literature. We do not assume that the forcings add reals. (© 2006 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

几乎各向同性的高维空间分数阶扩散方程的分块快速正则Hermite分裂预处理方法

一类空间分数阶扩散方程经过有限差分离散后所得到的离散线性方程组的系数矩阵是两个对角矩阵与Toeplitz型矩阵的乘积之和.在本文中,对于几乎各向同性的二维或三维空间分数阶扩散方程的离散线性方程组,采用预处理Krylov子空间迭代方法,我们利用其系数矩阵的特殊结构和具体性质构造了一类分块快速正则Hermite分裂预处理子.通过理论分析,我们证明了所对应的预处理矩阵的特征值大部分都聚集于1的附近.数值实验也表明,这类分块快速正则Hermite分裂预处理子可以明显地加快广义极小残量(GMRES)方法和稳定化的双共轭梯度(BiCGSTAB)方法等Krylov子空间迭代方法的收敛速度.     

M-矩阵平方根的一类迭代法

矩阵平方根在数学的许多应用中起着重要的作用.本文研究M-矩阵平方根的计算问题,提出一种计算正则M-矩阵平方根的迭代方法.首先将这个问题转化为M-矩阵代数Riccati方程,进而提出一种有效的方法来求解这个特殊的MARE.理论分析表明,该方法在一定条件下是收敛的.数值实验表明该方法是可行的,且优于二项式迭代法.     

A Procedure for Transforming Impedance Spectra for the Determination of Mechanism of Electrochemical Reactions

Transformation of impedance spectra into relaxation time spectra (RTS) is used for determining contributions of individual processes of the oxygen electroreduction reaction (OER) to the polarization resistance of the electrochemical cell. The transformation technique involves the solution of the convolution equation found with the aid of a modified Van Cittert iteration algorithm checked on model impedance spectra. The technique, when used to analyze impedance spectra of electrochemical cells air|Pt|YSZ|YSZ + Pt|air, shows that the conversion of a globular structure of the YSZ + Pt cermet layer to a columnar one is accompanied by a change of peak amplitudes in RTS. The revealed RTS dynamics when heated to 750°C is compared with peculiarities of individual processes in OER.     

高次非球面光路计算中的两个重要问题

提出在高次非球面的光路计算中运用变曲率的“辅助球面”概念,对某些文献在此方面存在的问题进行分析,提出迭代逼近精度的数值判据及浮动处理方法。     

Root determination by use of Padé approximants

In this paper we describe a new technique for generating iteration formulas — of arbitrary order — for determining a zero (assumed simple) of a functionf, assumed analytic in a region containing the zero. The 1/p Padé Approximant (p0) to the functiong(t)f(z) is formed wherez=w+t, using the Taylor series forf at the pointw, an approxination to the zero off. The value oft for which the 1/p Padé Approximant vanishes provides the basis of iteration formulas of orderp+2.Some known iteration formulas, e.g., Newton-Raphson's, Halley's and Kiss's of order of convergence two, three and four, are directly obtained by settingp=0,1 and 2, respectively.     

用方阵乘幂求和法求线性方程组的数值解

介绍了一种新型的,不同于传统的雅克比或高斯塞德尔迭代法的,求解线性方程组的方阵乘幂求和法,并引入了方阵意义上求积分的龙贝格法.该算法成立须以方阵A为实阵,非奇异且主对角元素占优.该法较雅克比或高斯塞德尔迭代的计算量小,特别有助于求解大型线性方程组的问题.     

Joint distribution of direction and length of typical I-segment in a homogeneous random planar tessellation stable under iteration

The paper deals with homogeneous random planar tessellations stable under iteration (random STIT tessellations). The length distribution of the typical I-segment is already known in the isotropic case [8]. In the present paper, the anisotropic case is treated. Then also the direction of the typical I-segment is of interest. The joint distribution of direction and length of the typical I-segment is evaluated. As a first step, the corresponding joint distribution for the so-called typical remaining I-segment is derived. Dedicated to the 80th birthday of Klaus Krickeberg     

A new subspace iteration method for the algebraic Riccati equation

We consider the numerical solution of the continuous algebraic Riccati equation A*X + XA ? XFX + G = 0, with F = F*,G = G* of low rank and A large and sparse. We develop an algorithm for the low‐rank approximation of X by means of an invariant subspace iteration on a function of the associated Hamiltonian matrix. We show that the sought‐after approximation can be obtained by a low‐rank update, in the style of the well known Alternating Direction Implicit (ADI) iteration for the linear equation, from which the new method inherits many algebraic properties. Moreover, we establish new insightful matrix relations with emerging projection‐type methods, which will help increase our understanding of this latter class of solution strategies. Copyright © 2014 John Wiley & Sons, Ltd.