一类常微分方程组的紧支撑样条小波方法探讨

应用三次紧支撑样条小波插值函数得到了求一类常微分方程组数值解的隐式公式,并求得到其局部截断误差为O(-h5).在此基础上给出1个显式校正求解公式,并讨论得到其局部截断误差为O(-h4).     

Boundedness and Dissipativity of Truncated Rotations on Uniform Planar Lattices

The finiteness of computer arithmetic can lead to some dramatic differences between the behaviour of a continuous dynamical system and a computer simulation. A thorough rigorous theoretical analysis of what may or what does happen is usually extremely difficult and to date little has been done even in relatively simple contexts. The comparative behaviour of a rotation mapping in the plane and on a uniform lattice in the plane is one such example. Simulations show that the rounding operator applied to a planar rotation mapping more or less preserves the qualitative behaviour of the original mapping, whereas the application of the truncation operator to a planar rotation can lead to quite different dynamical features. In this paper a theoretical justification of the properties of the planar rotation mappings under truncation to a, uniform integer lattice is provided, in particular properties of boundedness and dissipativity are investigated.     

Initial data truncation for multivariate output of discrete-event simulations using the Kalman filter

Data truncation is a commonly accepted method of dealing with initialization bias in discrete-event simulation. An algorithm for determining the appropriate initial-data truncation point for multivariate output is proposed. The technique entails averaging across independent replications and estimating a steady-state output model in a state-space framework. A Bayesian technique called Multiple Model Adaptive Estimation (MMAE) is applied to compute a time varying estimate of the output's steady-state mean vector. This MMAE implementation features the use, in paralle, of a bank of Kalman filters. Each filter is constructed under a different assumption concerning the output's steady-state mean vector. One of the filters assumes that the steady-state mean vector is accurately reflected by an estimate, called the assumed steady-state mean vector, taken from the last half of the simulation data. As the filters process the output through the effective transient, this particular filter becomes more likely (in a Bayesian sense) to be the best filter to represent the data and the MMAE mean estimator is influenced increasingly towards the assumed steady-state mean vector. The estimated truncation point is selected when a norm of the MMAE mean vector estimate is within a small tolerance of the assumed steady-state mean vector. A Monte Carlo analysis using data from simulations of open and closed queueing models is used to evaluate the technique. The evaluation criteria include the ability to construct accurate and reliable confidence regions for the mean response vector based on the truncated sequences.     

Search of truncation of (N−1) electron basis containing full connected triple excitations in computing main and satellite ionization potentials via Fock‐space coupled cluster approach

A valence‐universal multireference coupled cluster (VUMRCC) theory, realized via the eigenvalue independent partitioning (EIP) route, has been implemented with full inclusion of triples excitations for computing and analyzing the entire main and several satellite peaks in the ionization potential spectra of several molecules. The EIP‐VUMRCC method, unlike the traditional VUMRCC theory, allows divergence‐free homing‐in to satellite roots which would otherwise have been plagued by intruders, and is thus numerically more robust to obtain more efficient and dependable computational schemes allowing more extensive use of the approach. The computed ionization potentials (IPs) as a result of truncation of the (N−1) electron basis manifold involving virtual functions such as 2h‐p and 3h‐2p by different energy thresholds varying from 5 to 15 a.u. with 1 a.u. intervals as well as thresholds such as 20, 25, and 30 a.u. have been carefully looked into. Cutoff at around 25 a.u. turns out to be an optimal threshold. Molecules such as C₂H₄ and C₂H₂ (X = D,T), and N₂ and CO (X = D,T,Q) with Dunning's cc‐pVXZ bases have been investigated to determine all main and 2h‐p shake‐up and 3h‐2p double shake‐up satellite IPs. We believe that the present work will pave the way to a wider application of the method by providing main and satellite IPs for some problematic N‐electron closed shell systems. © 2013 Wiley Periodicals, Inc.     

R3中一类Kirchhoff型方程变号解的存在性及集中性

本文研究了一类Kirchhoff型方程。利用极大极小原理及惩罚函数方法,证明了上述方程变号解的存在性及集中性,我们的结果推广了文献[4]的结果。     

Distributed and multicore QuBiLS-MIDAS software v2.0: Computing chiral,fuzzy, weighted and truncated geometrical molecular descriptors based on tensor algebra

Advances to the distributed, multi-core and fully cross-platform QuBiLS-MIDAS software v2.0 ( http://tomocomd.com/qubils-midas ) are reported in this article since the v1.0 release. The QuBiLS-MIDAS software is the only one that computes atom-pair and alignment-free geometrical MDs (3D-MDs) from several distance metrics other than the Euclidean distance, as well as alignment-free 3D-MDs that codify structural information regarding the relations among three and four atoms of a molecule. The most recent features added to the QuBiLS-MIDAS software v2.0 are related (a) to the calculation of atomic weightings from indices based on the vertex-degree invariant (e.g., Alikhanidi index); (b) to consider central chirality during the molecular encoding; (c) to use measures based on clustering methods and statistical functions to codify structural information among more than two atoms; (d) to the use of a novel method based on fuzzy membership functions to spherically truncate inter-atomic relations; and (e) to the use of weighted and fuzzy aggregation operators to compute global 3D-MDs according to the importance and/or interrelation of the atoms of a molecule during the molecular encoding. Moreover, a novel module to compute QuBiLS-MIDAS 3D-MDs from their headings was also developed. This module can be used either by the graphical user interface or by means of the software library. By using the library, both the predictive models built with the QuBiLS-MIDAS 3D-MDs and the QuBiLS-MIDAS 3D-MDs calculation can be embedded in other tools. A set of predefined QuBiLS-MIDAS 3D-MDs with high information content and low redundancy on a set comprised of 20,469 compounds is also provided to be employed in further cheminformatics tasks. This set of predefined 3D-MDs evidenced better performance than all the universe of Dragon (v5.5) and PaDEL 0D-to-3D MDs in variability studies, whereas a linear independence study proved that these QuBiLS-MIDAS 3D-MDs codify chemical information orthogonal to the Dragon 0D-to-3D MDs. This set of predefined 3D-MDs would be periodically updated as long as new results be achieved. In general, this report highlights our continued efforts to provide a better tool for a most suitable characterization of compounds, and in this way, to contribute to obtaining better outcomes in future applications.     

Recent advances in surface x-ray diffraction and the potential for determining structure-sensitivity relations in single-crystal electrocatalysis

Determining how electrode structure governs the performance of an electrocatalyst requires techniques capable of probing structure at the atomic scale, often in situ and operando. In recent years, there have been numerous advances in the main experimental techniques for determining the structure of the electrochemical interface. In situ/operando synchrotron surface x-ray diffraction measurements are key to investigate the atomic structure of the electrode surfaces as well as understand the structure-reactivity relations in electrocatalysis. Here we discuss some recent improvements that have taken place in surface x-ray diffraction and how we expect them to lead to an enhanced understanding of electrocatalysis.     

A study of properties of adaptive quadratic grids for three-dimensional near-surface field computations

Curved geometries and the corresponding near-surface fields typically require a large number of linear computational elements. High-order numerical solvers have been primarily used with low-order meshes. There is a need for curved, high-order computational elements. Typical near-surface meshes consist of hexahedral and/or prismatic elements. The present work studies the employment of quadratic meshes that are relatively coarse for field simulations. Directionally quadratic high-order elements are proposed for the near-surface field regions. The quadratic meshes are compared with the conventional low-order ones in terms of accuracy and efficiency. The cases considered include closed surface volume calculations, as well as computation of gradients of several analytic fields. A special method of adaptive local quadratic meshes is proposed and evaluated. Truncation error analysis for quadratic grids yields comparison with the conventional linear hexahedral/prismatic meshes, which are subject to typical distortions such as stretching, skewness, and torsion.     

四阶奇异摄动边值问题在自适应网格上的一致收敛分析

we study a difference scheme for the fourth-order singular pertur-bation differential equation on the Bakhvalov-Shishkin grid by Green‘‘s function.The method is shown to be uniformly convergent with respect to the perturbation parameter,of order N^-2 in the maxmum norm on Bakhvalov-Shishkin meshes.Numerical results support our theoretical results.     

Non-linear elliptic problems having natural growth and <Emphasis Type="Italic">L</Emphasis><Superscript>1</Superscript> data in Orlicz spaces

An existence result for a strongly non-linear elliptic equation with natural growth condition on the non-linearity and L¹-data is proved. Mathematics Subject Classification (2000) 35J25, 35J60