解初值问题y″=g(x，y)的含参数块方法

1引言对于二阶常微分方程的初值问题y″=g(x,y),y(x_0)=y_0,y′(x_0)=y_0′,x_0(?)x(?)T(1)的数值解法的研究引起人们的广泛兴趣．对于直接积分(1),自从1976年J．D．Lambert和I．A．Waston提出二阶P-稳定方法和1978年G．Dahlquist证明P-稳定常系数线性多步方法的最高相容阶不超过2的重要结论以来,截止目前,已积累了许多高于2阶的P-稳定方法．例如,修正的Numerov方法,混合法(特殊形式RK的方法),多导法,Obrechkoff方法,显式RKN方法,单隐方法和对角隐式RKN方法等(顺便指出,文献[5,16]中所说的高阶方法的相容阶均不超过4)．所有这些方法,有些相     

Unconventional antiferromagnetic correlations of the doped Haldane gapsystem Y 2 BaNi 1 - x Zn x O 5

We make a new proposal to describe the very low temperature susceptibility of the doped Haldane gap compound Y₂BaNi_1-xZn_xO₅. We propose a new mean field model relevant for this compound. The ground state of this mean field model is unconventional because antiferromagnetism coexists with random dimers. We present new susceptibility experiments at very low temperature. We obtain a Curie-Weiss susceptibility χ( T ) ∼ C /(Θ + T ) as expected for antiferromagnetic correlations but we do not obtain a direct signature of antiferromagnetic long range order. We explain how to obtain the “impurity” susceptibility ( T ) by subtracting the Haldane gap contribution to the total susceptibility. In the temperature range [1 K, 300 K] the experimental data are well fitted by T ( T ) = C _imp 1 + T _imp / T . In the temperature range [100 mK, 1 K] the experimental data are well fitted by T ( T ) = A ln( T / T _c ), where T _c increases with x. This fit suggests the existence of a finite Néel temperature which is however too small to be probed directly in our experiments. We also obtain a maximum in the temperature dependence of the ac-susceptibility ( T ) which suggests the existence of antiferromagnetic correlations at very low temperature. Received 17 July 2001     

三级近似束流光学计算程序

为了精确计算束流在离子光学系统中的传输,用Visual FORTRAN 6.5语言编写了一个计算程序,长约13000行. 此程序可以计算由三圆筒单透镜、三膜片单透镜、双元筒透镜、均匀场静电加速管、磁四极透镜、六极磁铁、静电四极透镜、偏转磁铁、螺线管透镜、ExB~正交电磁场分析器、静电偏转器、漂浮管、QWR(Quarter Wave Resonators)和SLR(Split Loop Resonators)射频加速元件等元件任意组成的离子光学系统. 粒子轨迹的计算可精确到三级近似. 粒子的分布类型也可以有多种选择. 程序具有最优化计算功能,即可以自动调整元件的参数,以实现所需要的光学条件. 各元件之后的横向和纵向相图以及系统的束流包络线以图形方式显示在屏幕上.     

Nondegeneracy and uniqueness of positive solutions for Robin problem of second order ordinary differential equations and its applications

This article studies positive solutions of Robin problem for semi-linear second order ordinary differential equations. Nondegeneracy and uniqueness results are proven for homogeneous differential equations. Necessary and sufficient conditions for the existence of one or two positive solutions for inhomogeneous differential equations or differential equations with concave-convex nonlinearities are obtained by making use of the nondegeneracy and uniqueness results for positive solutions of homogeneous differential equations.     

Operator stable laws: Series representations and domains of normal attraction

LetX,X ₁,X ₂,... be i.i.d. random vectors in ^d. The limit laws that can arise by suitable affine normalizations of the partial sums,S _n=X ₁+...+X _n, are calledoperator-stable laws. These laws are a natural extension to ^d of the stable laws on. Thegeneralized domain of attraction of [GDOA()] is comprised of all random vectorsX whose partial sums can be affinely normalized to converge to . If the linear part of the affine transformation is restricted to take the formn ^–B for some exponent operatorB naturally associated to thenX is in thegeneralized domain of normal attraction of [GDONA()]. This paper extends the theory of operator-stable laws and their domains of attraction and normal attraction.     

The upper and lower second order directional derivatives of a sup-type function

The purpose of this paper is to give a formula for expressing the second order directional derivatives of the sup-type functionS(x) = sup{f(x, t); t T} in terms of the first and second derivatives off(x, t), whereT is a compact set in a metric space and we assume thatf, f/x and ² f/x ² are continuous on ⁿ × T. We will give a geometrical meaning of the formula. We will moreover give a sufficient condition forS(x) to be directionally twice differentiable.     

Minimum concave cost production system: A further generalization of multi-echelon model

We will consider a concave minimization problem associated with a series production system in which raw material is processed inm consecutive facilities. The products at some facility are either sent to the next facility or stocked in the warehouse. The amount of demand for the final products during periodi, i = 1,,n, are known in advance. Our problem is to minimize the sum of processing, holding and backlogging cost, all of which are assumed to be concave.The origin of this model is the classical economic lot size problem of Wagner and Whitin and was extensively studied by Zangwill. This model is very important from the theoretical as well as practical point of view and this is one of the very rare instances in which polynomial time algorithm has been constructed for concave minimization problems.The purpose of this paper is to extend the model further to the situation in which time lag is associated with processing at each facility. We will propose an efficient O(n ⁴ m) algorithm for this class of problems.     

On the spectrum of the Heisenberg Hamiltonian

The quantum, antiferromagnetic, spin-1/2 Heisenberg Hamiltonian on thed-dimensional cubic lattice ^d is considered for any dimensiond. First the anisotropic case is considered for small transversal coupling and a convergent expansion is given for a family of eigenprojections which is complete in all finite-volume truncations. Then the general case is considered, for which an upper bound to the ground-state energy is given which is optimal for strong enough anisotropy. This bound is expressed through a functional involving the statistical expectation value at finite temperature of a certain correlation function of an Ising model defined on the lattice ^d itself.     

纳米碳酸锌的制备和热分解动力学参数的测定

Nanosized zinc carbonate was prepared in the droplet of nano-reactor based on microemusion. The influences of the concentration of the surfactant and reactant on the diameter of nanoparticles were studied and the kinetic pa-rameters of thermal decomposition reaction were determined. It can be shown from experimental results， the di-ameter of the droplet in the microemulsion which have dissolved reactant zinc nitrate and sal volatile is between 4.9nm and 7.7nm， 6.2nm and 12.4nm， respectively. Nanosized zinc carbonate prepared by the method of mi-croemusion has shown good dispersion， narrow distribution and light agglomeration. The particle size of nanosized zinc carbonate is between 5nm and 40nm， and its kinetic parameters of thermal decomposition reaction-activation energy E and reaction order n are 110kJ·mol^-1 and 0.9， respectively.     

减压渣油供氢剂减黏裂化研究

减黏裂化作为一种不生成焦炭的热加工方法主要涉及两类化学反应，裂解反应和缩合反应。裂解反应使减压渣油的平均分子量及其胶团直径变小，改善渣油的倾点和黏度；缩合反应使减压渣油及其中间产物中的芳烃、烯烃等缩合成大分子量产物，产生新的胶核，甚至生成焦炭。常规减黏裂化过程中，这两类化学反应主要与反应温度和反应时间有关。在热作用下渣油中的大分子裂解为H／C原子比相对较高的饱和烃自由基和H／C原子比相对较低的芳香性自由基。后者还可失去自身的环烷氢而使其H／C原子比进一步降低。如果体系中有足量的氢化芳烃存在，那么氢化芳烃分子上的活泼氢就可能转移到芳香性自由基的单电子位，将自由基封闭，从而阻止芳香性自由基之间的相互缩合，抑制使分子长大的缩合反应。如果体系中的氢化芳烃含量不足。则热反应体系中就没有足够的活泼氢将芳香性自由基封闭，芳香性自由基缩合反应的几率较高，造成过早地生焦。这种来自渣油分子自身的活泼氢的化学行为与渣油热反应生焦诱导期有一定关系。