Toeplitz矩阵的循环延拓及其分析与应用

In this paper, we have presented the combined preconditioner which is derived from k =±-1~(1/2) circulant extensions of the real symmetric positive-definite Toeplitz matrices, proved it with great efficiency and stability and shown that it is easy to make error analysis and to remove the boundary effect with the combined preconditioner. This paper has also presented the methods for the direct and inverse computation of the real Toeplitz sets of equations and discussed many problems correspondingly, especially replaced the Toeplitz matrices with the combined preconditoners for analysis. The paper has also discussed the spectral analysis and boundary effect. Finally, as an application in geophysics, the paper makes some discussion about the squared root of a real matrix which comes from the Laplace algorithm.

Microbiologic Diagnosis of Prosthetic Shoulder Infection by Use of Implant Sonication

We recently described a sonication technique for the diagnosis of prosthetic knee and hip infections. We compared periprosthetic tissue culture to implant sonication followed by sonicate fluid culture for the diagnosis of prosthetic shoulder infection. One hundred thirty-six patients undergoing arthroplasty revision or resection were studied; 33 had definite prosthetic shoulder infections and 2 had probable prosthetic shoulder infections. Sonicate fluid culture was more sensitive than periprosthetic tissue culture for the detection of definite prosthetic shoulder infection (66.7 and 54.5%, respectively; P = 0.046). The specificities were similar (98.0% and 95.1%, respectively; P = 0.26). Propionibacterium acnes was the commonest species detected among culture-positive definite prosthetic shoulder infection cases by periprosthetic tissue culture (38.9%) and sonicate fluid culture (40.9%). All subjects from whom P. acnes was isolated from sonicate fluid were male. We conclude that sonicate fluid culture is useful for the diagnosis of prosthetic shoulder infection.

Sparse kernel density estimation technique based on zero-norm constraint

A sparse kernel density estimator is derived based on the zero-norm constraint, in which the zero-norm of the kernel weights is incorporated to enhance model sparsity. The classical Parzen window estimate is adopted as the desired response for density estimation, and an approximate function of the zero-norm is used for achieving mathemtical tractability and algorithmic efficiency. Under the mild condition of the positive definite design matrix, the kernel weights of the proposed density estimator based on the zero-norm approximation can be obtained using the multiplicative nonnegative quadratic programming algorithm. Using the -optimality based selection algorithm as the preprocessing to select a small significant subset design matrix, the proposed zero-norm based approach offers an effective means for constructing very sparse kernel density estimates with excellent generalisation performance.

Simplified recursive identifier for ARMA processes

Symmetrical behaviour of the covariance matrix and the positive-definite criterion are used to simplify identification of single-input/single-output systems using recursive least squares. Simulation results are obtained and these are compared with ordinary recursive least squares. The adaptive nature of the identifier is verified by varying the system parameters on convergence.

Boundary integral methods for singularly perturbed boundary value problems

In this paper we consider boundary integral methods applied to boundary value problems for the positive definite Helmholtz-type problem -DeltaU + alpha U-2 = 0 in a bounded or unbounded domain, with the parameter alpha real and possibly large. Applications arise in the implementation of space-time boundary integral methods for the heat equation, where alpha is proportional to 1/root deltat, and deltat is the time step. The corresponding layer potentials arising from this problem depend nonlinearly on the parameter alpha and have kernels which become highly peaked as alpha --> infinity, causing standard discretization schemes to fail. We propose a new collocation method with a robust convergence rate as alpha --> infinity. Numerical experiments on a model problem verify the theoretical results.

Nonlinear symmetric stability of planetary atmospheres

The energy–Casimir method is applied to the problem of symmetric stability in the context of a compressible, hydrostatic planetary atmosphere with a general equation of state. Formal stability criteria for symmetric disturbances to a zonally symmetric baroclinic flow are obtained. In the special case of a perfect gas the results of Stevens (1983) are recovered. Finite-amplitude stability conditions are also obtained that provide an upper bound on a certain positive-definite measure of disturbance amplitude.

A unified theory of available potential energy

Traditional derivations of available potential energy, in a variety of contexts, involve combining some form of mass conservation together with energy conservation. This raises the questions of why such constructions are required in the first place, and whether there is some general method of deriving the available potential energy for an arbitrary fluid system. By appealing to the underlying Hamiltonian structure of geophysical fluid dynamics, it becomes clear why energy conservation is not enough, and why other conservation laws such as mass conservation need to be incorporated in order to construct an invariant, known as the pseudoenergy, that is a positive‐definite functional of disturbance quantities. The available potential energy is just the non‐kinetic part of the pseudoenergy, the construction of which follows a well defined algorithm. Two notable features of the available potential energy defined thereby are first, that it is a locally defined quantity, and second, that it is inherently definable at finite amplitude (though one may of course always take the small‐amplitude limit if this is appropriate). The general theory is made concrete by systematic derivations of available potential energy in a number of different contexts. All the well known expressions are recovered, and some new expressions are obtained. The possibility of generalizing the concept of available potential energy to dynamically stable basic flows (as opposed to statically stable basic states) is also discussed.

Nonlinear wave-activity conservation laws and Hamiltonian structure for the two-dimensional anelastic equations

Exact, finite-amplitude, local wave-activity conservation laws are derived for disturbances to steady flows in the context of the two-dimensional anelastic equations. The conservation laws are expressed entirely in terms of Eulerian quantities, and have the property that, in the limit of a small-amplitude, slowly varying, monochromatic wave train, the wave-activity density A and flux F, when averaged over phase, satisfy F = cgA where cg is the group velocity of the waves. For nonparallel steady flows, the only conserved wave activity is a form of disturbance pseudoenergy; when the steady flow is parallel, there is in addition a conservation law for the disturbance pseudomomentum. The above results are obtained not only for isentropic background states (which give the so-called “deep form” of the anelastic equations), but also for arbitrary background potential-temperature profiles θ0(z) so long as the variation in θ0(z) over the depth of the fluid is small compared with θ0 itself. The Hamiltonian structure of the equations is established in both cases, and its symmetry properties discussed. An expression for available potential energy is also derived that, for the case of a stably stratified background state (i.e., dθ0/dz > 0), is locally positive definite; the expression is valid for fully three-dimensional flow. The counterparts to results for the two-dimensional Boussinesq equations are also noted.

On phantom thermodynamics

The thermodynamic properties of dark energy fluids described by an equation of state parameter omega = p/rho are rediscussed in the context of FRW type geometries. Contrarily to previous claims, it is argued here that the phantom regime omega < -1 is not physically possible since that both the temperature and the entropy of every physical fluids must be always positive definite. This means that one cannot appeal to negative temperature in order to save the phantom dark energy hypothesis as has been recently done in the literature. Such a result remains true as long as the chemical potential is zero. However, if the phantom fluid is endowed with a non-null chemical potential, the phantom field hypothesis becomes thermodynamically consistent, that is, there are macroscopic equilibrium states with T > 0 and S > 0 in the course of the Universe expansion. (C) 2008 Elsevier B.V. All rights reserved.

Newton`s iterates can converge to non-stationary points

In this note we discuss the convergence of Newton`s method for minimization. We present examples in which the Newton iterates satisfy the Wolfe conditions and the Hessian is positive definite at each step and yet the iterates converge to a non-stationary point. These examples answer a question posed by Fletcher in his 1987 book Practical methods of optimization.

Noether and topological currents equivalence and soliton/particle correspondence in affine sl(2)((1)) Toda theory coupled to matter

A submodel of the so-called conformal affine Toda model coupled to the matter field (CATM) is defined such that its real Lagrangian has a positive-definite kinetic term for the Toda field and a usual kinetic term for the (Dirac) spinor field. After spontaneously broken the conformal symmetry by means of BRST analysis, we end up with an effective theory, the off-critical affine Toda model coupled to the matter (ATM). It is shown that the ATM model inherits the remarkable properties of the general CATM model such as the soliton solutions, the particle/soliton correspondence and the equivalence between the Noether and topological currents. The classical solitonic spectrum of the ATM model is also discussed. (C) 2001 Elsevier B.V. B.V. All rights reserved.

Zero sets of bivariate Hermite polynomials

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Analysis and numerical solution of the Peterlin viscoelastic model

Liquids and gasses form a vital part of nature. Many of these are complex fluids with non-Newtonian behaviour. We introduce a mathematical model describing the unsteady motion of an incompressible polymeric fluid. Each polymer molecule is treated as two beads connected by a spring. For the nonlinear spring force it is not possible to obtain a closed system of equations, unless we approximate the force law. The Peterlin approximation replaces the length of the spring by the length of the average spring. Consequently, the macroscopic dumbbell-based model for dilute polymer solutions is obtained. The model consists of the conservation of mass and momentum and time evolution of the symmetric positive definite conformation tensor, where the diffusive effects are taken into account. In two space dimensions we prove global in time existence of weak solutions. Assuming more regular data we show higher regularity and consequently uniqueness of the weak solution. For the Oseen-type Peterlin model we propose a linear pressure-stabilized characteristics finite element scheme. We derive the corresponding error estimates and we prove, for linear finite elements, the optimal first order accuracy. Theoretical error of the pressure-stabilized characteristic finite element scheme is confirmed by a series of numerical experiments.

Hyalotekite, (Ba,Pb,K)(4)(Ca,Y)(2)Si-8(B,Be)(2)(Si,B)(2)O28F, a Tectosilicate Related to Scapolite: New Structure Refinement, Phase Transitions and a Short-Range Ordered 3B Superstructure

Hyalotekite, a framework silicate of composition (Ba,Pb,K)(4)(Ca,Y)(2)Si-8(B,Be)(2) (Si,B)(2)O28F, is found in relatively high-temperature(greater than or equal to 500 degrees C) Mn skarns at Langban, Sweden, and peralkaline pegmatites at Dara-i-Pioz, Tajikistan. A new paragenesis at Dara-i-Pioz is pegmatite consisting of the Ba borosilicates leucosphenite and tienshanite, as well as caesium kupletskite, aegirine, pyrochlore, microcline and quartz. Hyalotekite has been partially replaced by barylite and danburite. This hyalotekite contains 1.29-1.78 wt.% Y2O3, equivalent to 0.172-0.238 Y pfu or 8-11% Y on the Ca site; its Pb/(Pb+Ba) ratio ranges 0.36-0.44. Electron microprobe F contents of Langban and Dara-i-Pioz hyalotekite range 1.04-1.45 wt.%, consistent with full occupancy of the F site. A new refinement of the structure factor data used in the original structural determination of a Langban hyalotekite resulted in a structural formula, (Pb1.96Ba1.86K0.18)Ca-2(B1.76Be0.24)(Si1.56B0.44)Si8O28F, consistent with chemical data and all cations with positive-definite thermal parameters, although with a slight excess of positive charge (+57.14 as opposed to the ideal +57.00). An unusual feature of the hyalotekite framework is that 4 of 28 oxygens are non-bridging; by merging these 4 oxygens into two, the framework topology of scapolite is obtained. The triclinic symmetry of hyalotekite observed at room temperature is obtained from a hypothetical tetragonal parent structure via a sequence of displacive phase transitions. Some of these transitions are associated with cation ordering, either Pb-Ba ordering in the large cation sites, or B-Be and Si-B ordering on tetrahedral sites. Others are largely displacive but affect the coordination of the large cations (Pb, Ba, K, Ca). High-resolution electron microscopy suggests that the undulatory extinction characteristic of hyalotekite is due to a fine mosaic microstructure. This suggests that at least one of these transitions occurs in nature during cooling, and that it is first order with a large volume change. A diffuse superstructure observed by electron diffraction implies the existence of a further stage of short-range cation ordering which probably involves both (Pb,K)-Ba and (BeSi,BB)-BSi.

Sequential and parallel synchronous alternating iterative methods

The so-called parallel multisplitting nonstationary iterative Model A was introduced by Bru, Elsner, and Neumann [Linear Algebra and its Applications 103:175-192 (1988)] for solving a nonsingular linear system Ax = b using a weak nonnegative multisplitting of the first type. In this paper new results are introduced when A is a monotone matrix using a weak nonnegative multisplitting of the second type and when A is a symmetric positive definite matrix using a P -regular multisplitting. Also, nonstationary alternating iterative methods are studied. Finally, combining Model A and alternating iterative methods, two new models of parallel multisplitting nonstationary iterations are introduced. When matrix A is monotone and the multisplittings are weak nonnegative of the first or of the second type, both models lead to convergent schemes. Also, when matrix A is symmetric positive definite and the multisplittings are P -regular, the schemes are also convergent.