Lagrange multiplier method used in BESⅢ kinematic fitting

In this paper,the kinematic fitting with the Lagrange multiplier method has been studied for BESⅢ experiment.First we introduce the Lagrange multiplier method and implement kinematic constraints.Then we present the performance of the kinematic fitting algorithm.With the kinematic fitting,we can improve the resolution of track parameters and reduce the background.     

Lagrange multiplier method used in BESⅢ kinematic fitting

In this paper,the kinematic fitting with the Lagrange multiplier method has been studied for BESⅢ experiment.First we introduce the Lagrange multiplier method and implement kinematic constraints.Then we present the performance of the kinematic fitting algorithm.With the kinematic fitting,we can improve the resolution of track parameters and reduce the background.     

Improved Faddeev-Jackiw quantization of the electromagnetic field and Lagrange multiplier fields

We use the improved Faddeev-Jackiw quantization method to quantize the electromagnetic field and its Lagrange multiplier fields. The method's comparison with the usual Faddeev-Jackiw method and the Dirac method is given. We show that this method is equivalent to the Dirac method and also retains all the merits of the usual Faddeev-Jackiw method. Moreover, it is simpler than the usual one if one needs to obtain new secondary constraints. Therefore, the improved Faddeev-Jackiw method is essential. Meanwhile, we find the new meaning of the Lagrange multipliers and explain the Faddeev-Jackiw generalized brackets concerning the Lagrange multipliers.     

Variational iteration method for solving time-fractional diffusion equations in porous the medium

<正>The variational iteration method is successfully extended to the case of solving fractional differential equations, and the Lagrange multiplier of the method is identified in a more accurate way.Some diffusion models with fractional derivatives are investigated analytically,and the results show the efficiency of the new Lagrange multiplier for fractional differential equations of arbitrary order.     

Estimating the bound for the generalized Lorenz system

A chaotic system is bounded, and its trajectory is confined to a certain region which is called the chaotic attractor. No matter how unstable the interior of the system is, the trajectory never exceeds the chaotic attractor. In the present paper, the sphere bound of the generalized Lorenz system is given, based on the Lyapunov function and the Lagrange multiplier method. Furthermore, we show the actual parameters and perform numerical simulations.     

Variational iteration method for solving compressible Euler equations

This paper applies the variational iteration method to obtain approximate analytic solutions of compressible Euler equations in gas dynamics.This method is based on the use of Lagrange multiplier for identification of optimal values of parameters in a functional.Using this method,a rapid convergent sequence is produced which converges to the exact solutions of the problem.Numerical results and comparison with other two numerical solutions verify that this method is very convenient and efficient.     

Decay vertex reconstruction and 3-dimensional lifetime determination at BESⅢ

This paper focuses mainly on the vertex reconstruction of resonance particles with a relatively long lifetime such as Ks^0, A, as well as on lifetime measurements using a 3-dimensional fit. The kinematic constraints between the production and decay vertices and the decay vertex fitting algorithm based on the least squares method are both presented. Reconstruction efficiencies including experimental resolutions are discussed. The results and systematic errors are calculated based on a Monte Carlo simulation.     

Kinematics of τ two-body decay near τ threshold at BESⅢ

The kinematic properties of two-body decay near τ threshold are studied according to the special capacity of the BEPC accelerator and the BESⅢ detector.Explicitly presented are the transformations of energy and momentum of hadronic particles between different reference frames,and the corresponding distributions.A brand new method is proposed to obtain the energy spread of the accelerator by fitting the energy distribution of hadron from τ semi-leptonic decays.     

Approximate analytic solutions for a generalized Hirota—Satsuma coupled KdV equation and a coupled mKdV equation

<正>This paper applies the variational iteration method to obtain approximate analytic solutions of a generalized Hirota-Satsuma coupled Korteweg-de Vries(KdV) equation and a coupled modified Korteweg-de Vries(mKdV) equation. This method provides a sequence of functions which converges to the exact solution of the problem and is based on the use of the Lagrange multiplier for the identification of optimal values of parameters in a functional.Some examples are given to demonstrate the reliability and convenience of the method and comparisons are made with the exact solutions.     

Kinematics of τ two-body decay near τ threshold at BESⅢ

The kinematic properties of two-body decay near τ threshold are studied according to the special capacity of the BEPC accelerator and the BESⅢ detector.Explicitly presented are the transformations of energy and momentum of hadronic particles between different reference frames,and the corresponding distributions.A brand new method is proposed to obtain the energy spread of the accelerator by fitting the energy distribution of hadron from τ semi-leptonic decays.     

A Method for Determining the Phase Transition Point of the Bose-Einstein Condensation by Judging the Fitting Errors

In this paper, We demonstrate an approach to determine the phase transition point and critical temperature of Bose-Einstein condensation. During the fitting sequence of time-of-flight images, for just one picture we take four kinds of different calculations instead of only one. By calculating the quantitative least-square errors, which have always been neglected before, we find out that this value can act as a criterion to judge the status of atom clouds. Using this criterion, we can not only discriminate the status around the phase transition point, but can also find the critical point precisely. Also with this method, we can achieve the totally automatical running of calculating programs without human＇s judgments.     

Polynomial curve fitting method for refraction-angle extraction in diffraction enhanced imaging

X-ray diffraction sorbing low-Z sample. How enhanced imaging （DEI） is applied to extract phase information from to inspect internal structures of weakly abraw images measured in different positions of rocking curve is the key problem of DEI. In this paper, we present an effective extraction method called polynomial curve fitting method, in order to extract accurate information angular in a fast speed. It is compared with the existing methods such as multiple-images statistical method and Gaussian curve fitting method. The experiments results on a plastic cylinder and a black ant at the Beijing Synchrotron Radiation Facility prove that the polynomial curve fitting method can obtain most approximate refraction-angle values and its computation speed is 10 times faster than the Gaussian curve fitting method.     

Characterizing a proton beam with two different methods in beam halo experiments

In beam halo experiments, it is very important to correctly characterize the RFQ output proton beam. In order to simulate the beam dynamics properly, we must first know the correct initial beam parameters. We have used two different methods, quadrupole scans and multi-wire scanners to determine the transverse phase-space properties of the proton beam. The experimental data were analyzed by fitting to the 3-D nonlinear simulation code IMPACT. For the quadrupole scan method, we found that the RMS beam radius and the measured beam-core profiles agreed very well with the simulations. For the multi-wire scanner method, we choose the case of a matched beam. By fitting the IMPACT simulation results to the measured data, we obtained the Courant-Snyder parameters and the emittance of the beam. The difference between the two methods is about eight percent, which is acceptable in our experiments.     

Characterizing a proton beam with two different methods in beam halo experiments

In beam halo experiments, it is very important to correctly characterize the RFQ output proton beam. In order to simulate the beam dynamics properly, we must first know the correct initial beam parameters. We have used two different methods, quadrupole scans and multi-wire scanners to determine the transverse phase-space properties of the proton beam. The experimental data were analyzed by fitting to the 3-D nonlinear simulation code IMPACT. For the quadrupole scan method, we found that the RMS beam radius and the measured beam-core profiles agreed very well with the simulations. For the multi-wire scanner method, we choose the case of a matched beam. By fitting the IMPACT simulation results to the measured data, we obtained the Courant-Snyder parameters and the emittance of the beam. The difference between the two methods is about eight percent, which is acceptable in our experiments.     

Methods for Thickness Determination of SiC Homoepilayers by Using Infrared Reflectance Spectroscopy

Infrared reflectance spectra have been widely used to measure layer thickness based either on calculation or on curve fitting, and two traditional methods for thickness determination have been studied. Considering the disadvantages of those two methods, we propose a new fitting model based on the fitting of the fringe order difference. In comparison with the measured results, the new fitting model shows its superiority not only for its stable and accurate results which have great agreement with the results from SEM, but also for its simple and quick fitting process.     

Non-deterministic fatigue life analysis using convex set models

The non-probabilistic approach to fatigue life analysis was studied using the convex models-interval, ellipsoidal and multiconvex models. The lower and upper bounds of the fatigue life were obtained by using the second-order Taylor series and Lagrange multiplier method. The solving process for derivatives of the implicit life function was presented. Moreover, a median ellipsoidal model was proposed which can take into account the sample blind zone and almost impossibility of concurrence of some small probability events. The Monte Carlo method for multi-convex model was presented, an important alternative when the analytical method does not work. A project example was given. The feasibility and rationality of the presented approach were verified. It is also revealed that the proposed method is conservative compared to the traditional probabilistic method, but it is a useful complement when it is difficult to obtain the accurate probability densities of parameters.     

A new complex variable element-free Galerkin method for two-dimensional potential problems

In this paper, based on the element-free Galerkin (EFG) method and the improved complex variable moving least- square (ICVMLS) approximation, a new meshless method, which is the improved complex variable element-free Galerkin (ICVEFG) method for two-dimensional potential problems, is presented. In the method, the integral weak form of control equations is employed, and the Lagrange multiplier is used to apply the essential boundary conditions. Then the corresponding formulas of the ICVEFG method for two-dimensional potential problems are obtained. Compared with the complex variable moving least-square (CVMLS) approximation proposed by Cheng, the functional in the ICVMLS approximation has an explicit physical meaning. Furthermore, the ICVEFG method has greater computational precision and efficiency. Three numerical examples are given to show the validity of the proposed method.     

Estimating the directions of arrival based on multi-subarray subspace fitting

A multi-subarray subspace fitting method which take mutual coupling among array elements into account to estimate the directions of arrival was presented. The mutual coupling matrix of uniform linear arrays is modeled with a banded symmetric Toeplitz matrix. According to the DOF (Degree of Freedom) of mutual coupling matrix, part of the array elements at the two sides are neglected. The theoretical characteristics of ESPRIT algorithm, one of the subspace fitting algorithms, are studied. For large N and uniform linear array, the MCLS ESPRIT estimation error is asymptotically jointly Gaussian distributed with zero means, and its covariance expression is obtained. It is known from the simulation that there exists a subarray with lowest estimation error for certain number of the array elements and the signal sources, and its performance of estimating the directions of arrival is close to the ideal situation when mutual coupling does not exist while its variance has some increase. The method has been applied to bathymetric sidescan sonar with high resolution, and good results have been obtained. At the cost of increasing the number of the array elements, the method can reduce the affection of the mutual coupling among array elements.     

An improved method of energy calibration for position-sensitive silicon detectors

Energy calibration of resistive charge division-based position-sensitive silicon detectors is achieved by parabolic fitting in the traditional method, where the systematic variations of vertex and curvature of the parabola with energy must be considered. In this paper we extend the traditional method in order to correct the fitting function, simplify the procedure of calibration and improve the experimental data quality. Instead of a parabolic function as used in the traditional method, a new function describing the relation of position and energy is introduced.The energy resolution of the 8.088 Me V α decay of213 Rn is determined to be about 87 ke V(FWHM), which is better than the result of the traditional method, 104 ke V(FWHM). The improved method can be applied to the energy calibration of resistive charge division-based position-sensitive silicon detectors with various performances.     

Measurement of Random Surface Parameters of Weak Scatterers using the Speckle Contrast Method

Theoretical analysis shows that the deviation roughness ωand the lateral correlation length ξ of a weak scattering object determine the different properties of the contrast of the speckles in the image plane of a 4f system.Experimentally,we have measured the data of the speckle contrast versus the radius R of the variable filtering aperture.By fitting the theoretical results to these data,we extract the parameters ω and ξ of the random surfaces with Gaussian correlation.This method can determine the two parameters simultaneously and independently,and pre-calibrations are not needed.