Variable-step multilevel preconditioning methods,I: Self-adjoint and positive definite elliptic problems

When solving large size systems of equations by preconditioned iterative solution methods, one normally uses a fixed preconditioner which may be defined by some eigenvalue information, such as in a Chebyshev iteration method. In many problems, however, it may be more effective to use variable preconditioners, in particular when the eigenvalue information is not available. In the present paper, a recursive way of constructing variable-step of, in general, nonlinear multilevel preconditioners for selfadjoint and coercive second-order elliptic problems, discretized by the finite element method is proposed. The preconditioner is constructed recursively from the coarsest to finer and finer levels. Each preconditioning step requires only block-diagonal solvers at all levels except at every k₀, k₀ ≥ 1 level where we perform a sufficient number ν, ν ≥ 1 of GCG-type variable-step iterations that involve the use again of a variable-step preconditioning for that level. It turns out that for any sufficiently large value of k₀ and, asymptotically, for ν sufficiently large, but not too large, the method has both an optimal rate of convergence and an optimal order of computational complexity, both for two and three space dimensional problem domains. The method requires no parameter estimates and the convergence results do not depend on the regularity of the elliptic problem.     

Fastiterative methods for least squares estimations

Least squares estimations have been used extensively in many applications, e.g. system identification and signal prediction. When the stochastic process is stationary, the least squares estimators can be found by solving a Toeplitz or near-Toeplitz matrix system depending on the knowledge of the data statistics. In this paper, we employ the preconditioned conjugate gradient method with circulant preconditioners to solve such systems. Our proposed circulant preconditioners are derived from the spectral property of the given stationary process. In the case where the spectral density functions() of the process is known, we prove that ifs() is a positive continuous function, then the spectrum of the preconditioned system will be clustered around 1 and the method converges superlinearly. However, if the statistics of the process is unknown, then we prove that with probability 1, the spectrum of the preconditioned system is still clustered around 1 provided that large data samples are taken. For finite impulse response (FIR) system identification problems, our numerical results show that annth order least squares estimator can usually be obtained inO(n logn) operations whenO(n) data samples are used. Finally, we remark that our algorithm can be modified to suit the applications of recursive least squares computations with the proper use of sliding window method arising in signal processing applications.Research supported in part by HKRGC grant no. 221600070, ONR contract no. N00014-90-J-1695 and DOE grant no. DE-FG03-87ER25037.     

Fourier acceleration of iterative processes in disordered systems

Technical details are given on how to use Fourier acceleration with iterative processes such as relaxation and conjugate gradient methods. These methods are often used to solve large linear systems of equations, but become hopelessly slow very rapidly as the size of the set of equations to be solved increases. Fourier acceleration is a method designed to alleviate these problems and result in a very fast algorithm. The method is explained for the Jacobi relaxation and conjugate gradient methods and is applied to two models: the random resistor network and the random central-force network. In the first model, acceleration works very well; in the second, little is gained. We discuss reasons for this. We also include a discussion of stopping criteria.     

共轭线性对称性及其对PT-对称量子理论的应用

传统量子系统的哈密顿是自伴算子,哈密顿的自伴性不仅保证系统遵循酉演化和保持概率守恒,而且也保证了它自身具有实的能量本征值,这类系统称为自伴量子系统.然而,确实存在一些物理系统(如PT-对称量子系统),其哈密顿不是自伴的,这类系统称为非自伴量子系统.为了深入研究PT-对称量子系统,并考虑到算子PT的共轭线性性,首先讨论了共轭线性算子的一些性质,包括它们的矩阵表示和谱结构等;其次,分别研究了具有共轭线性对称性和完整共轭线性对称性的线性算子,通过它们的矩阵表示,给出了共轭线性对称性和完整共轭线性对称性的等价刻画;作为应用,得到了关于PT-对称及完整PT-对称算子的一些有趣性质,并通过一些具体例子,说明了完整PT-对称性对张量积运算不具有封闭性,同时说明了完整PT-对称性既不是哈密顿算子在某个正定内积下自伴的充分条件,也不是必要条件.     

A Two-Level Preconditioned Conjugate-Gradient Method in Distorted and Structured Grids

In this paper, we propose a new two-level preconditioned C-G method which uses the quadratic smoothing and the linear correction in distorted but topologically structured grid. The CPU time of this method is less than that of the multigrid preconditioned C-G method (MGCG) using the quadratic element, but their accuracy is almost the same. Numerical experiments and eigenvalue analysis are given and the results show that the proposed two-level preconditioned method is efficient.     

无直透光和共轭像的数字全息

根据数字全息的可计算特性，提出用图像处理的方法直接对数字全息图在空域进行预处理，使其再现时不出现直透光和共轭像，大幅度提高再现象的质量。并进行了实验验证。     

Identification of source terms in wave equation with dynamic boundary conditions

This paper studies an inverse hyperbolic problem for the wave equation with dynamic boundary conditions. It consists of determining some forcing terms from the final overdetermination of the displacement. First, the Fréchet differentiability of the Tikhonov functional is studied, and a gradient formula is obtained via the solution of an associated adjoint problem. Then, the Lipschitz continuity of the gradient is proved. Furthermore, the existence and the uniqueness for the minimization problem are discussed. Finally, some numerical experiments for the reconstruction of an internal wave force are implemented via a conjugate gradient algorithm.     

基于功能核酸的水中汞离子荧光检测方法

二价汞离子作为一种重要的环境污染物,一直以来受到国内外广泛关注。基于T-T错配的Hg2+检测极大依赖于汞选择性寡聚核苷酸（MSO）的设计,利用SYBR Green I对目前所报道的汞离子探针进行了优化,研究了若干种MSO探针与Hg2+的结合响应,在对探针二级结构进行分析和讨论的基础上,提出最优的富T探针序列,并由此建立了一种基于SYBR Green I的水中汞离子快速、便捷的荧光检测方法。最终测得3种实际水样的加标回收率在82.8％~101.8％之间,相对标准偏差小于15%,表明该方法受环境基质的影响较小,可应用于实际水样中的汞离子检测。     

一类可以产生独立同分布密钥流的混沌系统

构造了一类与帐篷映射拓扑同构的混沌系统, 并根据拓扑共轭变换关系给出了此类混沌系统产生独立、均匀分布密钥流序列的采样规则. 理论证明和数值模拟, 均验证了结论的有效性. 本文为产生独立同分布密钥流提供了更多的非线性系统选择. 实验结果证明利用本文定理产生的密钥流能够通过美国信息技术管理改革法案的随机数检测标准(FIPS PUB 140-2)和美国国家标准与技术研究院安全检测标准(NIST SP800-22), 符合密钥流的选取标准. 关键词： 独立同分布 混沌系统 帐篷映射 拓扑同构     

六聚吡咯/丝组缀合物结构的NMR和 ESI-MS/MS表征

利用片断合成法合成了具有DNA切割能力的六聚吡咯/丝组缀合物, 检测了该化合物的¹H、¹³C NMR 和ESI-MS/ MS 图谱, 确证了该化合物的结构, 通过¹H-¹H COSY, HSQC,HMBC等2D NMR 技术对其¹H 和¹³C NMR 数据进行了归属和解析, 并探讨了其ESI-MS/ MS 质谱裂解规律.