An algorithm for the computation of consistent initial values for differential–algebraic equations

We use boundary value methods to compute consistent initial values for fully implicit nonlinear differential-algebraic equations. The obtained algorithm uses variable order formulae and a deferred correction technique to evaluate the error. A rigorous theory is stated for nonlinear index 1, 2 and 3 DAEs of Hessenberg form. Numerical tests on classical index 1, 2 and 3 DAE problems are reported. This revised version was published online in June 2006 with corrections to the Cover Date.     

奇异类内离差矩阵条件下的Fisher最优判据

特征提取是模式分类与识别的重要环节，Fisher最优判据是特征提取的基本方法之一．本文提出了一种计算奇异类内离差矩阵条件下Fisher最优判据的新方法，并给出了计算步骤．     

Inexact Newton methods for the nonlinear complementarity problem

An exact Newton method for solving a nonlinear complementarity problem consists of solving a sequence of linear complementarity subproblems. For problems of large size, solving the subproblems exactly can be very expensive. In this paper we study inexact Newton methods for solving the nonlinear, complementarity problem. In such an inexact method, the subproblems are solved only up to a certain degree of accuracy. The necessary accuracies that are needed to preserve the nice features of the exact Newton method are established and analyzed. We also discuss some extensions as well as an application. This research was based on work supported by the National Science Foundation under grant ECS-8407240.     

鲤鱼体中鱼腥味物质的提取和鉴定

采用同时蒸馏萃取法提取鲤鱼鱼体中的挥发和半挥发性有机物，用GC／MS从提取物中分析鉴定出鱼腥味和疑是鱼腥味物质16种。包括醛、烯醛、酮和呋喃4类化合物。其中己醛、庚醛和2，4-二烯癸醛被确认为鱼腥味的化合物。研究发现，这些物质总量在鲤鱼体的鱼鳞、鱼鳃和鱼肉的分布呈现下降趋势，且其总量与鱼腥味强度之间具有可比较的对应关系。     

A recursive-faulting model of distributed damage in confined brittle materials

We develop a model of distributed damage in brittle materials deforming in triaxial compression based on the explicit construction of special microstructures obtained by recursive faulting. The model aims to predict the effective or macroscopic behavior of the material from its elastic and fracture properties; and to predict the microstructures underlying the microscopic behavior. The model accounts for the elasticity of the matrix, fault nucleation and the cohesive and frictional behavior of the faults. We analyze the resulting quasistatic boundary value problem and determine the relaxation of the potential energy, which describes the macroscopic material behavior averaged over all possible fine-scale structures. Finally, we present numerical calculations of the dynamic multi-axial compression experiments on sintered aluminum nitride of Chen and Ravichandran [1994. Dynamic compressive behavior of ceramics under lateral confinement. J. Phys. IV 4, 177-182; 1996a. Static and dynamic compressive behavior of aluminum nitride under moderate confinement. J. Am. Soc. Ceramics 79(3), 579-584; 1996b. An experimental technique for imposing dynamic multiaxial compression with mechanical confinement. Exp. Mech. 36(2), 155-158; 2000. Failure mode transition in ceramics under dynamic multiaxial compression. Int. J. Fracture 101, 141-159]. The model correctly predicts the general trends regarding the observed damage patterns; and the brittle-to-ductile transition resulting under increasing confinement.     

Global optimization by random perturbation of the gradient method with a fixed parameter

The paper deals with the global minimization of a differentiable cost function mapping a ball of a finite dimensional Euclidean space into an interval of real numbers. It is established that a suitable random perturbation of the gradient method with a fixed parameter generates a bounded minimizing sequence and leads to a global minimum: the perturbation avoids convergence to local minima. The stated results suggest an algorithm for the numerical approximation of global minima: experiments are performed for the problem of fitting a sum of exponentials to discrete data and to a nonlinear system involving about 5000 variables. The effect of the random perturbation is examined by comparison with the purely deterministic gradient method.     

Scaling by Binormalization

We present an iterative algorithm (BIN) for scaling all the rows and columns of a real symmetric matrix to unit 2-norm. We study the theoretical convergence properties and its relation to optimal conditioning. Numerical experiments show that BIN requires 2–4 matrix–vector multiplications to obtain an adequate scaling, and in many cases significantly reduces the condition number, more than other scaling algorithms. We present generalizations to complex, nonsymmetric and rectangular matrices.     

Exact solutions and geometric phase factor of time-dependent three-generator quantum systems

There exist a number of typical and interesting systems and/or models, which possess three-generator Lie-algebraic structure, in atomic physics, quantum optics, nuclear physics and laser physics. The well-known fact that all simple 3-generator algebras are either isomorphic to the algebra sl (2, C) or to one of its real forms enables us to treat these time-dependent quantum systems in a unified way. By making use of both the Lewis-Riesenfeld invariant theory and the invariant-related unitary transformation formulation, the present paper obtains exact solutions of the time-dependent Schr?dinger equations governing various three-generator Lie-algebraic quantum systems. For some quantum systems whose time-dependent Hamiltonians have no quasialgebraic structures, it is shown that the exact solutions can also be obtained by working in a sub-Hilbert-space corresponding to a particular eigenvalue of the conserved generator (i.e., the time-independent invariant that commutes with the time-dependent Hamiltonian). The topological property of geometric phase factors and its adiabatic limit in time-dependent systems is briefly discussed. Received 6 July 2002 / Received in final form 21 October 2002 Published online 11 February 2003     

Ion chromatography–inductively coupled plasma mass spectrometry determination of arsenic species in marine samples

Arsenic speciation analysis in marine samples was performed using ion chromatography (IC) with inductively coupled plasma mass spectrometry (ICP‐MS) detection. The separation of eight arsenic species, viz. arsenite, monomethyl arsonic acid, dimethylarsinic acid, arsenate, arsenobetaine, tetramethylarsine oxide, arsenocholine and tetramethylarsonium ion was achieved on a Dionex AS4A (weaker anion exchange column) by using a nitric acid pH gradient eluent (pH 3.3 to 1.3). The entire separation was accomplished in 12 min. The detection limits for the eight arsenic species by IC–ICP‐MS were in the range 0.03–1.6 µ g l^?1, based on 3σ of the blank response (n = 6). The repeatability and day‐to‐day reproducibility were calculated to be less than 10% (residual standard deviation) for all eight species. The method was validated by analyzing a certified reference material (DORM‐2, dogfish muscle) and then successfully applied to several marine samples, e.g. oyster, fish muscle, shrimp and marine algae. The low power microwave digestion was employed for the extraction of arsenic from seafood products. Copyright © 2004 John Wiley & Sons, Ltd.