The method of extracting Green's function between stations from cross correlation has proven to be effective theoretically and experimentally. It has been widely applied to surface wave tomography of the crust and up...The method of extracting Green's function between stations from cross correlation has proven to be effective theoretically and experimentally. It has been widely applied to surface wave tomography of the crust and upmost mantle. However, there are still controversies about why this method works. Snieder employed stationary phase approximation in evaluating contribution to cross correlation function from scatterers in the whole space, and concluded that it is the constructive interference of waves emitted by the scatterers near the receiver line that leads to the emergence of Green's function. His derivation demonstrates that cross correlation function is just the convolution of noise power spectrum and the Green's function. However, his derivation ignores influence from the two stationary points at infinities, therefore it may fail when attenuation is absent. In order to obtain accurate noise-correlation function due to scatters over the whole space, we compute the total contribution with numerical integration in polar coordinates. Our numerical computation of cross correlation function indicates that the incomplete stationary phase approximation introduces remarkable errors to the cross correlation function, in both amplitude and phase, when the frequency is low with reasonable quality factor Q. Our results argue that the dis- tance between stations has to be beyond several wavelengths in order to reduce the influence of this inaccuracy on the applications of ambient noise method, and only the station pairs whose distances are above several (〉5) wavelengths can be used.展开更多
The structure of any a.s. self-similar set K(w) generated by a class of random elements {gn,wσ} taking values in the space of contractive operators is given and the approximation of K(w) by the fixed points {Pn,wσ} ...The structure of any a.s. self-similar set K(w) generated by a class of random elements {gn,wσ} taking values in the space of contractive operators is given and the approximation of K(w) by the fixed points {Pn,wσ} of {gn,ow} is obtained. It is useful to generate the fractal in computer.展开更多
We study iterative processes of stochastic approximation for finding fixed points of weakly contractive and nonexpansive operators in Hilbert spaces under the condition that operators are given with random errors. We ...We study iterative processes of stochastic approximation for finding fixed points of weakly contractive and nonexpansive operators in Hilbert spaces under the condition that operators are given with random errors. We prove mean square convergence and convergence almost sure (a.s.) of iterative approximations and establish both asymptotic and nonasymptotic estimates of the convergence rate in degenerate and non-degenerate cases. Previously the stochastic approximation algorithms were studied mainly for optimization problems.展开更多
In this paper, we consider a Lorentz space with a mixed norm of periodic functions of many variables. We obtain the exact estimation of the best M-term approximations of Nikol'skii's and Besov's classes in the Lore...In this paper, we consider a Lorentz space with a mixed norm of periodic functions of many variables. We obtain the exact estimation of the best M-term approximations of Nikol'skii's and Besov's classes in the Lorentz space with the mixed norm.展开更多
Recently, we developed the projective truncation approximation for the equation of motion of two-time Green's functions(Fan et al., Phys. Rev. B 97, 165140(2018)). In that approximation, the precision of results d...Recently, we developed the projective truncation approximation for the equation of motion of two-time Green's functions(Fan et al., Phys. Rev. B 97, 165140(2018)). In that approximation, the precision of results depends on the selection of operator basis. Here, for three successively larger operator bases, we calculate the local static averages and the impurity density of states of the single-band Anderson impurity model. The results converge systematically towards those of numerical renormalization group as the basis size is enlarged. We also propose a quantitative gauge of the truncation error within this method and demonstrate its usefulness using the Hubbard-I basis. We thus confirm that the projective truncation approximation is a method of controllable precision for quantum many-body systems.展开更多
In this paper, we apply the two-time Green's function method, and provide a simple way to study themagnetic properties of one-dimensional spin-(S, s) Heisenberg ferromagnets.The magnetic susceptibility and correla...In this paper, we apply the two-time Green's function method, and provide a simple way to study themagnetic properties of one-dimensional spin-(S, s) Heisenberg ferromagnets.The magnetic susceptibility and correlationfunctions are obtained by using the Tyablikov decoupling approximation.Our results show that the magnetic susceptibilityand correlation length are a monotonically decreasing function of temperature regardless of the mixed spins.It isfound that in the case of S = s, our results of one-dimensional mixed-spin model is reduced to be those of the isotropicferromagnetic Heisenberg chain in the whole temperature region.Our results for the susceptibility are in agreement withthose obtained by other theoretical approaches.展开更多
In this paper, Lorentz space of functions of several variables and Besov's class are considered. We establish an exact approximation order of Besov's class by partial sums of Fourier's series for multiple trigonome...In this paper, Lorentz space of functions of several variables and Besov's class are considered. We establish an exact approximation order of Besov's class by partial sums of Fourier's series for multiple trigonometric system.展开更多
There exist many iterative methods for computing the maximum likelihood estimator but most of them suffer from one or several drawbacks such as the need to inverse a Hessian matrix and the need to find good initial ap...There exist many iterative methods for computing the maximum likelihood estimator but most of them suffer from one or several drawbacks such as the need to inverse a Hessian matrix and the need to find good initial approximations of the parameters that are unknown in practice. In this paper, we present an estimation method without matrix inversion based on a linear approximation of the likelihood equations in a neighborhood of the constrained maximum likelihood estimator. We obtain closed-form approximations of solutions and standard errors. Then, we propose an iterative algorithm which cycles through the components of the vector parameter and updates one component at a time. The initial solution, which is necessary to start the iterative procedure, is automated. The proposed algorithm is compared to some of the best iterative optimization algorithms available on R and MATLAB software through a simulation study and applied to the statistical analysis of a road safety measure.展开更多
In this work we slwly linear polynomial operators preserving some consecutive i-convexities and leaving in-verant the polynomtals up to a certain degree. First we study the existence of an incompatibility between the ...In this work we slwly linear polynomial operators preserving some consecutive i-convexities and leaving in-verant the polynomtals up to a certain degree. First we study the existence of an incompatibility between the conservation of cenain i-cotivexities and the invariance of a space of polynomials. Interpolation properties are obtained and a theorem by Berens and DcVore about the Bernstein's operator ts extended. Finally, from these results a genera'ized Bernstein's operator is obtained.展开更多
This paper presents a design method of H<sub>2</sub> and H<sub>∞</sub>-feedback control loop for nonlinear smooth gene networks that are in control affine form. Formulaic solution methodology ...This paper presents a design method of H<sub>2</sub> and H<sub>∞</sub>-feedback control loop for nonlinear smooth gene networks that are in control affine form. Formulaic solution methodology for solving the nonlinear partial differential equations, namely the Hamilton-Jacobi-Bellman and Hamilton-Jacobi-Isaacs equations through successive Galerkin’s approximation is implemented and the results are compared. Throughout the implementation, there were several caveats that need to be further resolved for practical applications in general cases. Such issues and the clarification of causes are mathematically established and reviewed.展开更多
In 1805, Thomas Young was the first to propose an equation(Young's equation) to predict the value of the equilibrium contact angle of a liquid on a solid. On the basis of our predecessors, we further clarify that ...In 1805, Thomas Young was the first to propose an equation(Young's equation) to predict the value of the equilibrium contact angle of a liquid on a solid. On the basis of our predecessors, we further clarify that the contact angle in Young's equation refers to the super-nano contact angle. Whether the equation is applicable to nanoscale systems remains an open question. Zhu et al. [College Phys. 4 7(1985)] obtained the most simple and convenient approximate formula, known as the Zhu–Qian approximate formula of Young's equation. Here, using molecular dynamics simulation, we test its applicability for nanodrops. Molecular dynamics simulations are performed on argon liquid cylinders placed on a solid surface under a temperature of 90 K, using Lennard–Jones potentials for the interaction between liquid molecules and between a liquid molecule and a solid molecule with the variable coefficient of strength a. Eight values of a between 0.650 and 0.825 are used. By comparison of the super-nano contact angles obtained from molecular dynamics simulation and the Zhu–Qian approximate formula of Young's equation, we find that it is qualitatively applicable for nanoscale systems.展开更多
Aggregate nearest neighbor(ANN) search retrieves for two spatial datasets T and Q, segment(s) of one or more trajectories from the set T having minimum aggregate distance to points in Q. When interacting with large am...Aggregate nearest neighbor(ANN) search retrieves for two spatial datasets T and Q, segment(s) of one or more trajectories from the set T having minimum aggregate distance to points in Q. When interacting with large amounts of trajectories, this process would be very time-consuming due to consecutive page loads. An approximate method for finding segments with minimum aggregate distance is proposed which can improve the response time. In order to index large volumes of trajectories, scalable and efficient trajectory index(SETI) structure is used. But some refinements are provided to temporal index of SETI to improve the performance of proposed method. The experiments were performed with different number of query points and percentages of dataset. It is shown that proposed method besides having an acceptable precision, can reduce the computation time significantly. It is also shown that the main fraction of search time among load time, ANN and computing convex and centroid, is related to ANN.展开更多
The private quantum channel (PQC) maps any quantum state to the maximally mixed state for the discrete as well as the bosonic Gaussian quantum systems, and it has fundamental meaning on the quantum cryptographic tasks...The private quantum channel (PQC) maps any quantum state to the maximally mixed state for the discrete as well as the bosonic Gaussian quantum systems, and it has fundamental meaning on the quantum cryptographic tasks and the quantum channel capacity problems. In this paper, we primally introduce a notion of approximate private quantum channel (<em>ε</em>-PQC) on <em>fermionic</em> Gaussian systems (<em>i.e.</em>, <em>ε</em>-FPQC), and construct its explicit form of the fermionic (Gaussian) private quantum channel. First of all, we suggest a general structure for <em>ε</em>-FPQC on the fermionic Gaussian systems with respect to the Schatten <em>p</em>-norm class, and then we give an explicit proof of the statement in the trace norm case. In addition, we study that the cardinality of a set of fermionic unitary operators agrees on the <em>ε</em>-FPQC condition in the trace norm case. This result may give birth to intuition on the construction of emerging fermionic Gaussian quantum communication or computing systems.展开更多
The bounds on the discrepancy of approximate solutions constructed by Gedunov's scheme to IVP of isentropic equations of gas dynamics are obtained, Three well-knowu results obtained by Lax for shock waves with sma...The bounds on the discrepancy of approximate solutions constructed by Gedunov's scheme to IVP of isentropic equations of gas dynamics are obtained, Three well-knowu results obtained by Lax for shock waves with small jumps for general quasilinear hyperbolic systems of conservation laws are extended to shock waves for isentropic equations of gas dynamics in a bounded invariant region with ρ=0 as one of boundries of the region. Two counterexamples are given to show that two iuequalities given by Godunov do not hold for all rational numbers γ∈(1, 3]. It seems that the approach by Godunov to obtain the forementioned bounds may not be possible.展开更多
Dykstra’s alternating projection algorithm was proposed to treat the problem of finding the projection of a given point onto the intersection of some closed convex sets. In this paper, we first apply Dykstra’s alter...Dykstra’s alternating projection algorithm was proposed to treat the problem of finding the projection of a given point onto the intersection of some closed convex sets. In this paper, we first apply Dykstra’s alternating projection algorithm to compute the optimal approximate symmetric positive semidefinite solution of the matrix equations AXB = E, CXD = F. If we choose the initial iterative matrix X<sub>0</sub> = 0, the least Frobenius norm symmetric positive semidefinite solution of these matrix equations is obtained. A numerical example shows that the new algorithm is feasible and effective.展开更多
基金supported by the National Natural Science Foundation of China (No. 40674027)CAS outstanding 100 research program,MOST program 2007FY220100
文摘The method of extracting Green's function between stations from cross correlation has proven to be effective theoretically and experimentally. It has been widely applied to surface wave tomography of the crust and upmost mantle. However, there are still controversies about why this method works. Snieder employed stationary phase approximation in evaluating contribution to cross correlation function from scatterers in the whole space, and concluded that it is the constructive interference of waves emitted by the scatterers near the receiver line that leads to the emergence of Green's function. His derivation demonstrates that cross correlation function is just the convolution of noise power spectrum and the Green's function. However, his derivation ignores influence from the two stationary points at infinities, therefore it may fail when attenuation is absent. In order to obtain accurate noise-correlation function due to scatters over the whole space, we compute the total contribution with numerical integration in polar coordinates. Our numerical computation of cross correlation function indicates that the incomplete stationary phase approximation introduces remarkable errors to the cross correlation function, in both amplitude and phase, when the frequency is low with reasonable quality factor Q. Our results argue that the dis- tance between stations has to be beyond several wavelengths in order to reduce the influence of this inaccuracy on the applications of ambient noise method, and only the station pairs whose distances are above several (〉5) wavelengths can be used.
基金Supported by NNSF of China and the Foundation of Wuhan University
文摘The structure of any a.s. self-similar set K(w) generated by a class of random elements {gn,wσ} taking values in the space of contractive operators is given and the approximation of K(w) by the fixed points {Pn,wσ} of {gn,ow} is obtained. It is useful to generate the fractal in computer.
文摘We study iterative processes of stochastic approximation for finding fixed points of weakly contractive and nonexpansive operators in Hilbert spaces under the condition that operators are given with random errors. We prove mean square convergence and convergence almost sure (a.s.) of iterative approximations and establish both asymptotic and nonasymptotic estimates of the convergence rate in degenerate and non-degenerate cases. Previously the stochastic approximation algorithms were studied mainly for optimization problems.
基金supported by the Ministry of Education and Science of Republic Kazakhstan(Grant No.5129/GF4)partially by the Russian Academic Excellence Project(agreement between the Ministry of Education and Science of the Russian Federation and Ural Federal University No.02.A03.21.006 of August 27,2013)
文摘In this paper, we consider a Lorentz space with a mixed norm of periodic functions of many variables. We obtain the exact estimation of the best M-term approximations of Nikol'skii's and Besov's classes in the Lorentz space with the mixed norm.
基金Project supported by the National Key Basic Research Program of China(Grant No.2012CB921704)the National Natural Science Foundation of China(Grant No.11374362)+1 种基金the Fundamental Research Funds for the Central Universitiesthe Research Funds of Renmin University of China(Grant No.15XNLQ03)
文摘Recently, we developed the projective truncation approximation for the equation of motion of two-time Green's functions(Fan et al., Phys. Rev. B 97, 165140(2018)). In that approximation, the precision of results depends on the selection of operator basis. Here, for three successively larger operator bases, we calculate the local static averages and the impurity density of states of the single-band Anderson impurity model. The results converge systematically towards those of numerical renormalization group as the basis size is enlarged. We also propose a quantitative gauge of the truncation error within this method and demonstrate its usefulness using the Hubbard-I basis. We thus confirm that the projective truncation approximation is a method of controllable precision for quantum many-body systems.
基金Supported by the Natural Science Foundation of Guangdong Province under Grant No.8151009001000055
文摘In this paper, we apply the two-time Green's function method, and provide a simple way to study themagnetic properties of one-dimensional spin-(S, s) Heisenberg ferromagnets.The magnetic susceptibility and correlationfunctions are obtained by using the Tyablikov decoupling approximation.Our results show that the magnetic susceptibilityand correlation length are a monotonically decreasing function of temperature regardless of the mixed spins.It isfound that in the case of S = s, our results of one-dimensional mixed-spin model is reduced to be those of the isotropicferromagnetic Heisenberg chain in the whole temperature region.Our results for the susceptibility are in agreement withthose obtained by other theoretical approaches.
文摘In this paper, Lorentz space of functions of several variables and Besov's class are considered. We establish an exact approximation order of Besov's class by partial sums of Fourier's series for multiple trigonometric system.
文摘There exist many iterative methods for computing the maximum likelihood estimator but most of them suffer from one or several drawbacks such as the need to inverse a Hessian matrix and the need to find good initial approximations of the parameters that are unknown in practice. In this paper, we present an estimation method without matrix inversion based on a linear approximation of the likelihood equations in a neighborhood of the constrained maximum likelihood estimator. We obtain closed-form approximations of solutions and standard errors. Then, we propose an iterative algorithm which cycles through the components of the vector parameter and updates one component at a time. The initial solution, which is necessary to start the iterative procedure, is automated. The proposed algorithm is compared to some of the best iterative optimization algorithms available on R and MATLAB software through a simulation study and applied to the statistical analysis of a road safety measure.
基金This work was supported by Junta de Andalucia. Grupo de investigacion Matematica Aplioada. Codao 1107
文摘In this work we slwly linear polynomial operators preserving some consecutive i-convexities and leaving in-verant the polynomtals up to a certain degree. First we study the existence of an incompatibility between the conservation of cenain i-cotivexities and the invariance of a space of polynomials. Interpolation properties are obtained and a theorem by Berens and DcVore about the Bernstein's operator ts extended. Finally, from these results a genera'ized Bernstein's operator is obtained.
文摘This paper presents a design method of H<sub>2</sub> and H<sub>∞</sub>-feedback control loop for nonlinear smooth gene networks that are in control affine form. Formulaic solution methodology for solving the nonlinear partial differential equations, namely the Hamilton-Jacobi-Bellman and Hamilton-Jacobi-Isaacs equations through successive Galerkin’s approximation is implemented and the results are compared. Throughout the implementation, there were several caveats that need to be further resolved for practical applications in general cases. Such issues and the clarification of causes are mathematically established and reviewed.
基金Project supported by the National Natural Science Foundation of China(Grant No.11072242)the Key Scientific Studies Program of Hebei Province Higher Education Institute,China(Grant No.ZD2018301)Cangzhou National Science Foundation,China(Grant No.177000001)
文摘In 1805, Thomas Young was the first to propose an equation(Young's equation) to predict the value of the equilibrium contact angle of a liquid on a solid. On the basis of our predecessors, we further clarify that the contact angle in Young's equation refers to the super-nano contact angle. Whether the equation is applicable to nanoscale systems remains an open question. Zhu et al. [College Phys. 4 7(1985)] obtained the most simple and convenient approximate formula, known as the Zhu–Qian approximate formula of Young's equation. Here, using molecular dynamics simulation, we test its applicability for nanodrops. Molecular dynamics simulations are performed on argon liquid cylinders placed on a solid surface under a temperature of 90 K, using Lennard–Jones potentials for the interaction between liquid molecules and between a liquid molecule and a solid molecule with the variable coefficient of strength a. Eight values of a between 0.650 and 0.825 are used. By comparison of the super-nano contact angles obtained from molecular dynamics simulation and the Zhu–Qian approximate formula of Young's equation, we find that it is qualitatively applicable for nanoscale systems.
文摘Aggregate nearest neighbor(ANN) search retrieves for two spatial datasets T and Q, segment(s) of one or more trajectories from the set T having minimum aggregate distance to points in Q. When interacting with large amounts of trajectories, this process would be very time-consuming due to consecutive page loads. An approximate method for finding segments with minimum aggregate distance is proposed which can improve the response time. In order to index large volumes of trajectories, scalable and efficient trajectory index(SETI) structure is used. But some refinements are provided to temporal index of SETI to improve the performance of proposed method. The experiments were performed with different number of query points and percentages of dataset. It is shown that proposed method besides having an acceptable precision, can reduce the computation time significantly. It is also shown that the main fraction of search time among load time, ANN and computing convex and centroid, is related to ANN.
文摘The private quantum channel (PQC) maps any quantum state to the maximally mixed state for the discrete as well as the bosonic Gaussian quantum systems, and it has fundamental meaning on the quantum cryptographic tasks and the quantum channel capacity problems. In this paper, we primally introduce a notion of approximate private quantum channel (<em>ε</em>-PQC) on <em>fermionic</em> Gaussian systems (<em>i.e.</em>, <em>ε</em>-FPQC), and construct its explicit form of the fermionic (Gaussian) private quantum channel. First of all, we suggest a general structure for <em>ε</em>-FPQC on the fermionic Gaussian systems with respect to the Schatten <em>p</em>-norm class, and then we give an explicit proof of the statement in the trace norm case. In addition, we study that the cardinality of a set of fermionic unitary operators agrees on the <em>ε</em>-FPQC condition in the trace norm case. This result may give birth to intuition on the construction of emerging fermionic Gaussian quantum communication or computing systems.
基金Project supported by National Natural Science Foundation of China.
文摘The bounds on the discrepancy of approximate solutions constructed by Gedunov's scheme to IVP of isentropic equations of gas dynamics are obtained, Three well-knowu results obtained by Lax for shock waves with small jumps for general quasilinear hyperbolic systems of conservation laws are extended to shock waves for isentropic equations of gas dynamics in a bounded invariant region with ρ=0 as one of boundries of the region. Two counterexamples are given to show that two iuequalities given by Godunov do not hold for all rational numbers γ∈(1, 3]. It seems that the approach by Godunov to obtain the forementioned bounds may not be possible.
文摘Dykstra’s alternating projection algorithm was proposed to treat the problem of finding the projection of a given point onto the intersection of some closed convex sets. In this paper, we first apply Dykstra’s alternating projection algorithm to compute the optimal approximate symmetric positive semidefinite solution of the matrix equations AXB = E, CXD = F. If we choose the initial iterative matrix X<sub>0</sub> = 0, the least Frobenius norm symmetric positive semidefinite solution of these matrix equations is obtained. A numerical example shows that the new algorithm is feasible and effective.