The main aim of this paper is to study the local anisotropic interpolation error estimates. We show that the interpolation of a nonconforming element satisfy the anisotropic property for both the second and fourth ord...The main aim of this paper is to study the local anisotropic interpolation error estimates. We show that the interpolation of a nonconforming element satisfy the anisotropic property for both the second and fourth order problems.展开更多
In this paper, a-posteriori error estimators are proposed for the Legendre spectral Galerkin method for two-point boundary value problems. The key idea is to postprocess the Galerkin approximation, and the analysis sh...In this paper, a-posteriori error estimators are proposed for the Legendre spectral Galerkin method for two-point boundary value problems. The key idea is to postprocess the Galerkin approximation, and the analysis shows that the postproeess improves the order of convergence. Consequently, we obtain asymptotically exact aposteriori error estimators based on the postprocessing results. Numerical examples are included to illustrate the theoretical analysis.展开更多
We present several numerical methods and establish their error estimates for the discretization of the nonlinear Dirac equation (NLDE) in the nonrelativistic limit regime, involving a small dimensionless parameter 0...We present several numerical methods and establish their error estimates for the discretization of the nonlinear Dirac equation (NLDE) in the nonrelativistic limit regime, involving a small dimensionless parameter 0 〈 ε〈〈1 which is inversely proportional to the speed of light. In this limit regime, the solution is highly oscillatory in time, i.e., there are propagating waves with wavelength O( ε^2) and O(1) in time and space, respectively. We begin with the conservative Crank-Nicolson finite difference (CNFD) method and establish rigorously its error estimate which depends explicitly on the mesh size h and time step τ- as well as the small parameter 0 〈 ε≤1 Based on the error bound, in order to obtain 'correct' numerical solutions in the nonrelativistic limit regime, i.e., 0 〈 ε≤1 , the CNFD method requests the ε-scalability: τ- = O(ε3) and h = O(√ε). Then we propose and analyze two numerical methods for the discretization of NLDE by using the Fourier spectral discretization for spatial derivatives combined with the exponential wave integrator and time- splitting technique for temporal derivatives, respectively. Rigorous error bounds for the two numerical methods show that their ε-scalability is improved to τ = O(ε2) and h = O(1) when 0 〈 ε 〈〈 1. Extensive numerical results are reported to confirm our error estimates.展开更多
Based on the work of Xu and Zhou(2000),this paper makes a further discussion on conforming finite elements approximation for Steklov eigenvalue problems,and proves a local a priori error estimate and a new local a pos...Based on the work of Xu and Zhou(2000),this paper makes a further discussion on conforming finite elements approximation for Steklov eigenvalue problems,and proves a local a priori error estimate and a new local a posteriori error estimate in ||·||1,Ω0 norm for conforming elements eigenfunction,which has not been studied in existing literatures.展开更多
In this paper,a residual type of a posteriori error estimator for the general second order elliptic eigenpair approximation by the mixed finite element method is derived and analyzed,based on a type of superconvergenc...In this paper,a residual type of a posteriori error estimator for the general second order elliptic eigenpair approximation by the mixed finite element method is derived and analyzed,based on a type of superconvergence result of the eigenfunction approximation.Its efficiency and reliability are proved by both theoretical analysis and numerical experiments.展开更多
Consider a semiparametric regression model with linear time series errors Y_k= x′ _kβ + g(t_k) + ε_k, 1 ≤ k ≤ n, where Y_k's are responses, x_k =(x_(k1),x_(k2),···,x_(kp))′ and t_k ∈ T is con...Consider a semiparametric regression model with linear time series errors Y_k= x′ _kβ + g(t_k) + ε_k, 1 ≤ k ≤ n, where Y_k's are responses, x_k =(x_(k1),x_(k2),···,x_(kp))′ and t_k ∈ T is contained in R are fixed design points, β =(β_1,β_2,···,β_p)′ is an unknown parameter vector, g(·) is an unknown bounded real-valuedfunction defined on a compact subset T of the real line R, and ε_k is a linear process given byε_k = ∑ from j=0 to ∞ of ψ_je_(k-j), ψ_0=1, where ∑ from j=0 to ∞ of |ψ_j| < ∞, and e_j,j=0, +-1, +-2,···, ard i.i.d. random variables. In this paper we establish the asymptoticnormality of the least squares estimator of β, a smooth estimator of g(·), and estimators of theautocovariance and autocorrelation functions of the linear process ε_k.展开更多
This work is concerned with time stepping finite element methods for abstract second order evolution problems. We derive optimal order a posteriori error estimates and a posteriori nodal superconvergence error estimat...This work is concerned with time stepping finite element methods for abstract second order evolution problems. We derive optimal order a posteriori error estimates and a posteriori nodal superconvergence error estimates using the energy approach and the duality argument. With the help of the a posteriori error estimator developed in this work, we will further propose an adaptive time stepping strategy. A number of numerical experiments are performed to illustrate the reliability and efficiency of the a posteriori error estimates and to assess the effectiveness of the proposed adaptive time stepping method.展开更多
This paper investigates the finite element approximation of a class of parameter estimation problems which is the form of performance as the optimal control problems governed by bilinear parabolic equations,where the ...This paper investigates the finite element approximation of a class of parameter estimation problems which is the form of performance as the optimal control problems governed by bilinear parabolic equations,where the state and co-state are discretized by piecewise linear functions and control is approximated by piecewise constant functions.The authors derive some a priori error estimates for both the control and state approximations.Finally,the numerical experiments verify the theoretical results.展开更多
This paper considers the Legendre Galerkin spectral approximation for the unconstralnea optimal control problems. The authors derive a posteriori error estimate for the spectral approximation scheme of optimal control...This paper considers the Legendre Galerkin spectral approximation for the unconstralnea optimal control problems. The authors derive a posteriori error estimate for the spectral approximation scheme of optimal control problem. By choosing the appropriate basis functions, the stiff matrix of the discretization equations is sparse. And the authors use the Fast Legendre Transform to improve the efficiency of this method. Two numerical experiments demonstrating our theoretical results are presented.展开更多
In this paper,optimize-then-discretize,variational discretization and the finite volume method are applied to solve the distributed optimal control problems governed by a second order hyperbolic equation.A semi-discre...In this paper,optimize-then-discretize,variational discretization and the finite volume method are applied to solve the distributed optimal control problems governed by a second order hyperbolic equation.A semi-discrete optimal system is obtained.We prove the existence and uniqueness of the solution to the semidiscrete optimal system and obtain the optimal order error estimates in L ∞(J;L 2)-and L ∞(J;H 1)-norm.Numerical experiments are presented to test these theoretical results.展开更多
The authors are concerned with a zero-flux type initial boundary value problem for scalar conservation laws.Firstly,a kinetic formulation of entropy solutions is established.Secondly,by using the kinetic formulation a...The authors are concerned with a zero-flux type initial boundary value problem for scalar conservation laws.Firstly,a kinetic formulation of entropy solutions is established.Secondly,by using the kinetic formulation and kinetic techniques,the uniqueness of entropy solutions is obtained.Finally,the parabolic approximation is studied and an error estimate of order η 1/3 between the entropy solution and the viscous approximate solutions is established by using kinetic techniques,where η is the size of artificial viscosity.展开更多
The authors consider the partially linear model relating a response Y to predictors (x, T) with a mean function x^Tβ0 + g(T) when the x's are measured with an additive error. The estimators of parameter β0 are...The authors consider the partially linear model relating a response Y to predictors (x, T) with a mean function x^Tβ0 + g(T) when the x's are measured with an additive error. The estimators of parameter β0 are derived by using the nearest neighbor-generalized randomly weighted least absolute deviation (LAD for short) method. The resulting estimator of the unknown vector 30 is shown to be consistent and asymptotically normal. In addition, the results facilitate the construction of confidence regions and the hypothesis testing for the unknown parameters. Extensive simulations are reported, showing that the proposed method works well in practical settings. The proposed methods are also applied to a data set from the study of an AIDS clinical trial group.展开更多
文摘The main aim of this paper is to study the local anisotropic interpolation error estimates. We show that the interpolation of a nonconforming element satisfy the anisotropic property for both the second and fourth order problems.
基金supported partially by the innovation fund of Shanghai Normal Universitysupported partially by NSERC of Canada under Grant OGP0046726.
文摘In this paper, a-posteriori error estimators are proposed for the Legendre spectral Galerkin method for two-point boundary value problems. The key idea is to postprocess the Galerkin approximation, and the analysis shows that the postproeess improves the order of convergence. Consequently, we obtain asymptotically exact aposteriori error estimators based on the postprocessing results. Numerical examples are included to illustrate the theoretical analysis.
基金supported by the Ministry of Education of Singapore(Grant No.R146-000-196-112)National Natural Science Foundation of China(Grant No.91430103)
文摘We present several numerical methods and establish their error estimates for the discretization of the nonlinear Dirac equation (NLDE) in the nonrelativistic limit regime, involving a small dimensionless parameter 0 〈 ε〈〈1 which is inversely proportional to the speed of light. In this limit regime, the solution is highly oscillatory in time, i.e., there are propagating waves with wavelength O( ε^2) and O(1) in time and space, respectively. We begin with the conservative Crank-Nicolson finite difference (CNFD) method and establish rigorously its error estimate which depends explicitly on the mesh size h and time step τ- as well as the small parameter 0 〈 ε≤1 Based on the error bound, in order to obtain 'correct' numerical solutions in the nonrelativistic limit regime, i.e., 0 〈 ε≤1 , the CNFD method requests the ε-scalability: τ- = O(ε3) and h = O(√ε). Then we propose and analyze two numerical methods for the discretization of NLDE by using the Fourier spectral discretization for spatial derivatives combined with the exponential wave integrator and time- splitting technique for temporal derivatives, respectively. Rigorous error bounds for the two numerical methods show that their ε-scalability is improved to τ = O(ε2) and h = O(1) when 0 〈 ε 〈〈 1. Extensive numerical results are reported to confirm our error estimates.
基金supported by National Natural Science Foundation of China(Grant Nos.11201093 and 11161012)
文摘Based on the work of Xu and Zhou(2000),this paper makes a further discussion on conforming finite elements approximation for Steklov eigenvalue problems,and proves a local a priori error estimate and a new local a posteriori error estimate in ||·||1,Ω0 norm for conforming elements eigenfunction,which has not been studied in existing literatures.
基金supported by National Natural Science Foundation of China(Grant Nos.11001259,11031006,11071265,11201501 and 91230110)National Basic Research Program of China(973 Project)(Grant No. 2011CB309703)+3 种基金International S&T Cooperation Program of China(Grant No. 2010DFR00700)Croucher Foundation of Hong Kong Baptist Universitythe National Center for Mathematics and Interdisciplinary Science,CAS,the President Foundation of AMSS-CASthe Fundamental Research Funds for the CentralUniversities(Grant No. 2012121003)
文摘In this paper,a residual type of a posteriori error estimator for the general second order elliptic eigenpair approximation by the mixed finite element method is derived and analyzed,based on a type of superconvergence result of the eigenfunction approximation.Its efficiency and reliability are proved by both theoretical analysis and numerical experiments.
基金CHEN Min's work is supported by Grant No. 70221001 and No. 70331001 from NNSFC and Grant No. KZCX2-SW-118 from CAS.
文摘Consider a semiparametric regression model with linear time series errors Y_k= x′ _kβ + g(t_k) + ε_k, 1 ≤ k ≤ n, where Y_k's are responses, x_k =(x_(k1),x_(k2),···,x_(kp))′ and t_k ∈ T is contained in R are fixed design points, β =(β_1,β_2,···,β_p)′ is an unknown parameter vector, g(·) is an unknown bounded real-valuedfunction defined on a compact subset T of the real line R, and ε_k is a linear process given byε_k = ∑ from j=0 to ∞ of ψ_je_(k-j), ψ_0=1, where ∑ from j=0 to ∞ of |ψ_j| < ∞, and e_j,j=0, +-1, +-2,···, ard i.i.d. random variables. In this paper we establish the asymptoticnormality of the least squares estimator of β, a smooth estimator of g(·), and estimators of theautocovariance and autocorrelation functions of the linear process ε_k.
基金supported by National Natural Science Foundation of China(Grant Nos.1117121911161130004 and 11101199)+1 种基金E-Institutes of Shanghai Municipal Education Commission(Grant No.E03004)Program for New Century Excellent Talents in Fujian Province University(Grant No.JA12260)
文摘This work is concerned with time stepping finite element methods for abstract second order evolution problems. We derive optimal order a posteriori error estimates and a posteriori nodal superconvergence error estimates using the energy approach and the duality argument. With the help of the a posteriori error estimator developed in this work, we will further propose an adaptive time stepping strategy. A number of numerical experiments are performed to illustrate the reliability and efficiency of the a posteriori error estimates and to assess the effectiveness of the proposed adaptive time stepping method.
基金supported by the National Natural Science Foundation of China under Grant Nos.11101025,11071080,11171113the National Natural Science Foundation of China under Grant No.11126279+1 种基金the Fundamental Research Funds for the Central Universitiesthe Youth Foundation of Tianyuan Mathematics
文摘This paper investigates the finite element approximation of a class of parameter estimation problems which is the form of performance as the optimal control problems governed by bilinear parabolic equations,where the state and co-state are discretized by piecewise linear functions and control is approximated by piecewise constant functions.The authors derive some a priori error estimates for both the control and state approximations.Finally,the numerical experiments verify the theoretical results.
基金supported by the Foundation for Talent Introduction of Guangdong Provincial University,Guangdong Province Universities and Colleges Pearl River Scholar Funded Scheme(2008)the National Natural Science Foundation of China under Grant No.10971074
文摘This paper considers the Legendre Galerkin spectral approximation for the unconstralnea optimal control problems. The authors derive a posteriori error estimate for the spectral approximation scheme of optimal control problem. By choosing the appropriate basis functions, the stiff matrix of the discretization equations is sparse. And the authors use the Fast Legendre Transform to improve the efficiency of this method. Two numerical experiments demonstrating our theoretical results are presented.
基金supported by National Natural Science Foundation of China(Grant Nos.11261011,11271145 and 11031006)Foundation of Guizhou Science and Technology Department(Grant No.[2011]2098)+2 种基金Foundation for Talent Introduction of Guangdong Provincial UniversitySpecialized Research Fund for the Doctoral Program of Higher Education(Grant No. 20114407110009)the Project of Department of Education of Guangdong Province(Grant No. 2012KJCX0036)
文摘In this paper,optimize-then-discretize,variational discretization and the finite volume method are applied to solve the distributed optimal control problems governed by a second order hyperbolic equation.A semi-discrete optimal system is obtained.We prove the existence and uniqueness of the solution to the semidiscrete optimal system and obtain the optimal order error estimates in L ∞(J;L 2)-and L ∞(J;H 1)-norm.Numerical experiments are presented to test these theoretical results.
基金supported by the National Natural Science Foundation of China (No. 10971135)the Program for New Century Excellent Talents of the Ministry of Education of China (No. NCET-07-0546)+2 种基金the University Young Teacher Sciences Foundation of Anhui Province (No. 2010SQRL145)Shanghai Jiao Tong University Innovation Fund for Postgraduates (No. AE071202)the Quality Project Fund of Fuyang Teachers College (No. 2010JPKC07)
文摘The authors are concerned with a zero-flux type initial boundary value problem for scalar conservation laws.Firstly,a kinetic formulation of entropy solutions is established.Secondly,by using the kinetic formulation and kinetic techniques,the uniqueness of entropy solutions is obtained.Finally,the parabolic approximation is studied and an error estimate of order η 1/3 between the entropy solution and the viscous approximate solutions is established by using kinetic techniques,where η is the size of artificial viscosity.
文摘The authors consider the partially linear model relating a response Y to predictors (x, T) with a mean function x^Tβ0 + g(T) when the x's are measured with an additive error. The estimators of parameter β0 are derived by using the nearest neighbor-generalized randomly weighted least absolute deviation (LAD for short) method. The resulting estimator of the unknown vector 30 is shown to be consistent and asymptotically normal. In addition, the results facilitate the construction of confidence regions and the hypothesis testing for the unknown parameters. Extensive simulations are reported, showing that the proposed method works well in practical settings. The proposed methods are also applied to a data set from the study of an AIDS clinical trial group.