In the past decade,notable progress has been achieved in the development of the generalized finite difference method(GFDM).The underlying principle of GFDM involves dividing the domain into multiple sub-domains.Within...In the past decade,notable progress has been achieved in the development of the generalized finite difference method(GFDM).The underlying principle of GFDM involves dividing the domain into multiple sub-domains.Within each sub-domain,explicit formulas for the necessary partial derivatives of the partial differential equations(PDEs)can be obtained through the application of Taylor series expansion and moving-least square approximation methods.Consequently,the method generates a sparse coefficient matrix,exhibiting a banded structure,making it highly advantageous for large-scale engineering computations.In this study,we present the application of the GFDM to numerically solve inverse Cauchy problems in two-and three-dimensional piezoelectric structures.Through our preliminary numerical experiments,we demonstrate that the proposed GFDMapproach shows great promise for accurately simulating coupled electroelastic equations in inverse problems,even with 3%errors added to the input data.展开更多
We develop the generalized conditional symmetry (GCS) approach to solve the problem of dimensional reduction of Cauchy problems for the KdV-type equations. We characterize these equations that admit certain higheror...We develop the generalized conditional symmetry (GCS) approach to solve the problem of dimensional reduction of Cauchy problems for the KdV-type equations. We characterize these equations that admit certain higherorder GCSs and show the main reduction procedure by some examples. The obtained reductions cannot be derived within the framework of the standard Lie approach.展开更多
This paper is devoted to studying symmetry reduction of Cauchy problems for the fourth-order quasi-linear parabolic equations that admit certain generalized conditional symmetries (GCSs). Complete group classificati...This paper is devoted to studying symmetry reduction of Cauchy problems for the fourth-order quasi-linear parabolic equations that admit certain generalized conditional symmetries (GCSs). Complete group classification results are presented, and some examples are given to show the main reduction procedure.展开更多
We exploit higher-order conditional symmetry to reduce initial-value problems for evolution equations toCauchy problems for systems of ordinary differential equations (ODEs).We classify a class of fourth-order evoluti...We exploit higher-order conditional symmetry to reduce initial-value problems for evolution equations toCauchy problems for systems of ordinary differential equations (ODEs).We classify a class of fourth-order evolutionequations which admit certain higher-order generalized conditional symmetries (GCSs) and give some examples to showthe main reduction procedure.These reductions cannot be derived within the framework of the standard Lie approach,which hints that the technique presented here is something essential for the dimensional reduction of evolu tion equations.展开更多
L^p- L^q decay estimate of solution to Cauchy problem of a linear thermoviscoelastic system is studied. By using a diagonalization argument of frequency analysis, the coupled system will be decoupled micrologically. T...L^p- L^q decay estimate of solution to Cauchy problem of a linear thermoviscoelastic system is studied. By using a diagonalization argument of frequency analysis, the coupled system will be decoupled micrologically. Then with the help of the information of characteristic roots for the coefficient matrix of the system, L^p- L^q decay estimate of parabolic type of solution to the Cauchy problem is obtained.展开更多
This article is concerned with the ill-posed Cauchy problem associated with a densely defined linear operator A in a Banach space. A family of weak regularizing operators is introduced. If the spectrum of A is contain...This article is concerned with the ill-posed Cauchy problem associated with a densely defined linear operator A in a Banach space. A family of weak regularizing operators is introduced. If the spectrum of A is contained in a sector of right-half complex plane and its resolvent is polynomially bounded, the weak regularization for such ill-posed Cauchy problem can be shown by using the quasi-reversibilky method and regularized semigroups. Finally, an example is given.展开更多
In this paper, α-times integrated C-regularized cosine functions and mild α-times integrated C-existence families of second order are introduced. Equivalences are proved among α-times integrated C-regularized cosin...In this paper, α-times integrated C-regularized cosine functions and mild α-times integrated C-existence families of second order are introduced. Equivalences are proved among α-times integrated C-regularized cosine function for a linear operator A, C-wellposed of (α+1)-times abstract Cauchy problem and mild a -times integrated C-existence family of second order for A when the commutable condition is satisfied. In addition, if A = C-1AC, they are also equivalent to A generating the α -times integrated C-regularized cosine function. The characterization of an exponentially bounded mild α -times integrated C-existence family of second order is given out in terms of a Laplace transform.展开更多
In this paper, we study the regularity of mild solution for the following fractional abstract Cauchy problem Dt αu(t)=Au(t)+f(t), t ∈ (0,T] u(0)= x0 on a Banach space X with order α ∈ (0,1), where the fractional d...In this paper, we study the regularity of mild solution for the following fractional abstract Cauchy problem Dt αu(t)=Au(t)+f(t), t ∈ (0,T] u(0)= x0 on a Banach space X with order α ∈ (0,1), where the fractional derivative is understood in the sense of Caputo fractional derivatives. We show that if A generates an analytic α-times resolvent family on X and f ∈ Lp ([0,T];X) for some p > 1/α, then the mild solution to the above equation is in Cα-1/p[ò,T] for every ò > 0. Moreover, if f is H?lder continuous, then so are the Dt αu(t) and Au(t).展开更多
It is well known that the Cauchy problem for Laplace equations is an ill-posed problem in Hadamard’s sense.Small deviations in Cauchy data may lead to large errors in the solutions.It is observed that if a bound is i...It is well known that the Cauchy problem for Laplace equations is an ill-posed problem in Hadamard’s sense.Small deviations in Cauchy data may lead to large errors in the solutions.It is observed that if a bound is imposed on the solution,there exists a conditional stability estimate.This gives a reasonable way to construct stable algorithms.However,it is impossible to have good results at all points in the domain.Although numerical methods for Cauchy problems for Laplace equations have been widely studied for quite a long time,there are still some unclear points,for example,how to evaluate the numerical solutions,which means whether they can approximate the Cauchy data well and keep the bound of the solution,and at which points the numerical results are reliable?In this paper,the authors will prove the conditional stability estimate which is quantitatively related to harmonic measures.The harmonic measure can be used as an indicate function to pointwisely evaluate the numerical result,which further enables us to find a reliable subdomain where the local convergence rate is higher than a certain order.展开更多
This paper is devoted to the Cauchy problem for the generalized damped Boussinesq equation with a nonlinear source term in the natural energy space.With the help of linear time-space estimates,we establish the local e...This paper is devoted to the Cauchy problem for the generalized damped Boussinesq equation with a nonlinear source term in the natural energy space.With the help of linear time-space estimates,we establish the local existence and uniqueness of solutions by means of the contraction mapping principle.The global existence and blow-up of the solutions at both subcritical and critical initial energy levels are obtained.Moreover,we construct the sufficient conditions of finite time blow-up of the solutions with arbitrary positive initial energy.展开更多
We consider the higher-order Cauchy problem (ACP_n) x^(n)(t)=sum from i=0 to n-1 B_ix^(i)(t)_1x^(i)(0)=x_i for 0≤i≤n-1,where B_i(0≤i≤n-1) are closed linear operators on a Banach space X such that D=∩ i=0 n-1 D(B_...We consider the higher-order Cauchy problem (ACP_n) x^(n)(t)=sum from i=0 to n-1 B_ix^(i)(t)_1x^(i)(0)=x_i for 0≤i≤n-1,where B_i(0≤i≤n-1) are closed linear operators on a Banach space X such that D=∩ i=0 n-1 D(B_i)is dense in X. It is well known that the solvability and the well-posedness of (ACP_n)were studied only in some special cases, such as D(B_(n-1))?D(B_i) for 0≤i≤n-2 by F. Neu-brander and a factoring case by J. T. Sandefur. In this paper, by using some new results ofvector valued Laplace transforms given by W. Arenddt, we obtain some characterizations ofthe solvability and some sufficiency conditions of the well-posedness for general (ACP_n),which generalize F. Neubrander's results and the famous results for (ACP_1)展开更多
In this paper, on the one hand, we take the conventional quasi-reversibility method to obtain the error estimates of approximate solutions of the Cauchy problems for parabolic equations in a sub-domain of QT with stro...In this paper, on the one hand, we take the conventional quasi-reversibility method to obtain the error estimates of approximate solutions of the Cauchy problems for parabolic equations in a sub-domain of QT with strong restrictions to the measured boundary data. On the other hand, weakening the conditions on the measured data, then combining the duality method in optimization with the quasi-reversibility method, we solve the Cauchy problems for parabolic equations in the presence of noisy data. Using this method, we can get the proper regularization parameter ε that we need in the quasi-reversibility method and obtain the convergence rate of approximate solutions as the noise of amplitude δ tends to zero.展开更多
In this article,a level-set approach for solving nonlinear elliptic Cauchy problems with piecewise constant solutions is proposed,which allows the definition of a Tikhonov functional on a space of level-set functions....In this article,a level-set approach for solving nonlinear elliptic Cauchy problems with piecewise constant solutions is proposed,which allows the definition of a Tikhonov functional on a space of level-set functions.We provide convergence analysis for the Tikhonov approach,including stability and convergence results.Moreover,a numerical investigation of the proposed Tikhonov regularization method is presented.Newton-type methods are used for the solution of the optimality systems,which can be interpreted as stabilized versions of algorithms in a previous work and yield a substantial improvement in performance.The whole approach is focused on three dimensional models,better suited for real life applications.展开更多
If the second order problem u(t) + Bu(t) + Au(t) = f(t), u(0) =u(0) = 0 has L^p-maximal regularity, 1 〈 p 〈 ∞, the analyticity of the corresponding propagator of the sine type is shown by obtaining th...If the second order problem u(t) + Bu(t) + Au(t) = f(t), u(0) =u(0) = 0 has L^p-maximal regularity, 1 〈 p 〈 ∞, the analyticity of the corresponding propagator of the sine type is shown by obtaining the estimates of ‖λ(λ^2 + λB + A)^-1‖ and ‖B(λ^2 + λB + A)^-1‖ for λ∈ C with Reλ 〉 ω, where the constant ω≥ 0.展开更多
The authors prove well posedness in Gevrey classes of Cauchy problem for nonlinear hyper- bolic equations of constant multiplicity with Holder dependence on the time variable.
We consider the growth rate and quenching rate of the following problem with singular nonlinearityfor some positive constants b:, b2 (see Theorem 3.3 for the parametersfor some constantsHence, the solution (u, v) ...We consider the growth rate and quenching rate of the following problem with singular nonlinearityfor some positive constants b:, b2 (see Theorem 3.3 for the parametersfor some constantsHence, the solution (u, v) quenches at the originx = 0 at the same time '1' (see Theorem 4.3). We also tind various other conditions tor the solution to quench in a finite time and obtain the corresponding decay rate of the solution near the quenching time.展开更多
In this paper,the notion of C-semigroup of continuous module homomorphisms on a complete random normal(RN)module is introduced and investigated.The existence and uniqueness of solution to the Cauchy problem with respe...In this paper,the notion of C-semigroup of continuous module homomorphisms on a complete random normal(RN)module is introduced and investigated.The existence and uniqueness of solution to the Cauchy problem with respect to exponentially bounded C-semigroups of continuous module homomorphisms in a complete RN module are established.展开更多
In this paper,we address 3D inverse Cauchy issues of highly nonlinear elliptic equations in large cuboids by utilizing the new 3D homogenization functions of different orders to adapt all the specified boundary data.W...In this paper,we address 3D inverse Cauchy issues of highly nonlinear elliptic equations in large cuboids by utilizing the new 3D homogenization functions of different orders to adapt all the specified boundary data.We also add the average classification as an approximate solution to the nonlinear operator part,without requiring to cope with nonlinear equations to resolve the weighting coefficients because these constructions are owned many conditions about the true solution.The unknown boundary conditions and the result can be retrieved straightway by coping with a small-scale linear system when the outcome is described by a new 3D homogenization function,which is right to find the numerical solutions with the errors smaller than the level of noise being put on the over-specified Neumann conditions on the bottom of the cuboid.Besides,note that the new homogenization functions method(HFM)does not require dealing with the regularization and highly nonlinear equations.The robustness and accuracy of the HFM are verified by comparing the recovered results of several numerical experiments to the exact solutions in the entire region,even though a very large level of noise 50%is imposed on the over specified Neumann conditions.The numerical errors of our scheme are in the order of O(10^(−1))-O(10^(−4)).展开更多
Integrated semigroups which are not exponentially bounded have been studied. The integrated solutions of inhomogeneous and nonlinear abstract Caucby problems are discussed. The reults extends and improves the theorems...Integrated semigroups which are not exponentially bounded have been studied. The integrated solutions of inhomogeneous and nonlinear abstract Caucby problems are discussed. The reults extends and improves the theorems in [1].展开更多
The simplified linear model of Grad-Shafranov (GS) reconstruction can be reformulated into an inverse boundary value problem of Laplace's equation. Therefore, in this paper we focus on the method of solving the inv...The simplified linear model of Grad-Shafranov (GS) reconstruction can be reformulated into an inverse boundary value problem of Laplace's equation. Therefore, in this paper we focus on the method of solving the inverse boundary value problem of Laplace's equation. In the first place, the variational regularization method is used to deal with the ill- posedness of the Cauchy problem for Laplace's equation. Then, the 'L-Curve' principle is suggested to be adopted in choosing the optimal regularization parameter. Finally, a numerical experiment is implemented with a section of Neumann and Dirichlet boundary conditions with observation errors. The results well converge to the exact solution of the problem, which proves the efficiency and robustness of the proposed method. When the order of observation error δ is 10-1, the order of the approximate result error can reach 10-3.展开更多
基金the Natural Science Foundation of Shandong Province of China(Grant No.ZR2022YQ06)the Development Plan of Youth Innovation Team in Colleges and Universities of Shandong Province(Grant No.2022KJ140)the Key Laboratory ofRoad Construction Technology and Equipment(Chang’an University,No.300102253502).
文摘In the past decade,notable progress has been achieved in the development of the generalized finite difference method(GFDM).The underlying principle of GFDM involves dividing the domain into multiple sub-domains.Within each sub-domain,explicit formulas for the necessary partial derivatives of the partial differential equations(PDEs)can be obtained through the application of Taylor series expansion and moving-least square approximation methods.Consequently,the method generates a sparse coefficient matrix,exhibiting a banded structure,making it highly advantageous for large-scale engineering computations.In this study,we present the application of the GFDM to numerically solve inverse Cauchy problems in two-and three-dimensional piezoelectric structures.Through our preliminary numerical experiments,we demonstrate that the proposed GFDMapproach shows great promise for accurately simulating coupled electroelastic equations in inverse problems,even with 3%errors added to the input data.
基金supported by National Natural Science Foundation of China under Grant Nos.10447007 and 10671156the Natural Science Foundation of Shaanxi Province of China under Grant No.2005A13
文摘We develop the generalized conditional symmetry (GCS) approach to solve the problem of dimensional reduction of Cauchy problems for the KdV-type equations. We characterize these equations that admit certain higherorder GCSs and show the main reduction procedure by some examples. The obtained reductions cannot be derived within the framework of the standard Lie approach.
基金Supported by the National Natural Science Foundation of China under Grant No.10671156the Natural Science Foundation of Shaanxi Province of China under Grant No.SJ08A05
文摘This paper is devoted to studying symmetry reduction of Cauchy problems for the fourth-order quasi-linear parabolic equations that admit certain generalized conditional symmetries (GCSs). Complete group classification results are presented, and some examples are given to show the main reduction procedure.
基金National Natural Science Foundation of China under Grant Nos.10447007 and 10671156the Natural Science Foundation of Shaanxi Province of China under Grant No.2005A13
文摘We exploit higher-order conditional symmetry to reduce initial-value problems for evolution equations toCauchy problems for systems of ordinary differential equations (ODEs).We classify a class of fourth-order evolutionequations which admit certain higher-order generalized conditional symmetries (GCSs) and give some examples to showthe main reduction procedure.These reductions cannot be derived within the framework of the standard Lie approach,which hints that the technique presented here is something essential for the dimensional reduction of evolu tion equations.
基金supported by the National Natural Science Foundation of China (10771055)HNSF(07JJ3007)
文摘L^p- L^q decay estimate of solution to Cauchy problem of a linear thermoviscoelastic system is studied. By using a diagonalization argument of frequency analysis, the coupled system will be decoupled micrologically. Then with the help of the information of characteristic roots for the coefficient matrix of the system, L^p- L^q decay estimate of parabolic type of solution to the Cauchy problem is obtained.
基金This project was supported by TRAPOYT, the Key Project of Chinese Ministry of Education(104126) the NNSF of China(10371046)
文摘This article is concerned with the ill-posed Cauchy problem associated with a densely defined linear operator A in a Banach space. A family of weak regularizing operators is introduced. If the spectrum of A is contained in a sector of right-half complex plane and its resolvent is polynomially bounded, the weak regularization for such ill-posed Cauchy problem can be shown by using the quasi-reversibilky method and regularized semigroups. Finally, an example is given.
基金This project is supported by the Natural Science Foundation of China and Science Development Foundation of the Colleges and University of Shanghai.
文摘In this paper, α-times integrated C-regularized cosine functions and mild α-times integrated C-existence families of second order are introduced. Equivalences are proved among α-times integrated C-regularized cosine function for a linear operator A, C-wellposed of (α+1)-times abstract Cauchy problem and mild a -times integrated C-existence family of second order for A when the commutable condition is satisfied. In addition, if A = C-1AC, they are also equivalent to A generating the α -times integrated C-regularized cosine function. The characterization of an exponentially bounded mild α -times integrated C-existence family of second order is given out in terms of a Laplace transform.
文摘In this paper, we study the regularity of mild solution for the following fractional abstract Cauchy problem Dt αu(t)=Au(t)+f(t), t ∈ (0,T] u(0)= x0 on a Banach space X with order α ∈ (0,1), where the fractional derivative is understood in the sense of Caputo fractional derivatives. We show that if A generates an analytic α-times resolvent family on X and f ∈ Lp ([0,T];X) for some p > 1/α, then the mild solution to the above equation is in Cα-1/p[ò,T] for every ò > 0. Moreover, if f is H?lder continuous, then so are the Dt αu(t) and Au(t).
基金suported by the National Natural Science Foundation of China(Nos.11971121,12201386,12241103)Grant-in-Aid for Scientific Research(A)20H00117 of Japan Society for the Promotion of Science.
文摘It is well known that the Cauchy problem for Laplace equations is an ill-posed problem in Hadamard’s sense.Small deviations in Cauchy data may lead to large errors in the solutions.It is observed that if a bound is imposed on the solution,there exists a conditional stability estimate.This gives a reasonable way to construct stable algorithms.However,it is impossible to have good results at all points in the domain.Although numerical methods for Cauchy problems for Laplace equations have been widely studied for quite a long time,there are still some unclear points,for example,how to evaluate the numerical solutions,which means whether they can approximate the Cauchy data well and keep the bound of the solution,and at which points the numerical results are reliable?In this paper,the authors will prove the conditional stability estimate which is quantitatively related to harmonic measures.The harmonic measure can be used as an indicate function to pointwisely evaluate the numerical result,which further enables us to find a reliable subdomain where the local convergence rate is higher than a certain order.
基金supported by the National Natural Science Foundation of China(12301272)the Natural Science Foundation of Henan(202300410109)the Cultivation Programme for Young Backbone Teachers in Henan University of Technology,and the Innovative Funds Plan of Henan University of Technology(2020ZKCJ09).
文摘This paper is devoted to the Cauchy problem for the generalized damped Boussinesq equation with a nonlinear source term in the natural energy space.With the help of linear time-space estimates,we establish the local existence and uniqueness of solutions by means of the contraction mapping principle.The global existence and blow-up of the solutions at both subcritical and critical initial energy levels are obtained.Moreover,we construct the sufficient conditions of finite time blow-up of the solutions with arbitrary positive initial energy.
文摘We consider the higher-order Cauchy problem (ACP_n) x^(n)(t)=sum from i=0 to n-1 B_ix^(i)(t)_1x^(i)(0)=x_i for 0≤i≤n-1,where B_i(0≤i≤n-1) are closed linear operators on a Banach space X such that D=∩ i=0 n-1 D(B_i)is dense in X. It is well known that the solvability and the well-posedness of (ACP_n)were studied only in some special cases, such as D(B_(n-1))?D(B_i) for 0≤i≤n-2 by F. Neu-brander and a factoring case by J. T. Sandefur. In this paper, by using some new results ofvector valued Laplace transforms given by W. Arenddt, we obtain some characterizations ofthe solvability and some sufficiency conditions of the well-posedness for general (ACP_n),which generalize F. Neubrander's results and the famous results for (ACP_1)
基金supported by National Natural Science Foundation of China (Grant No.11226166)Scientific Research Fund of Hu'nan Provincial Education Department (Grant No.11C0052)
文摘In this paper, on the one hand, we take the conventional quasi-reversibility method to obtain the error estimates of approximate solutions of the Cauchy problems for parabolic equations in a sub-domain of QT with strong restrictions to the measured boundary data. On the other hand, weakening the conditions on the measured data, then combining the duality method in optimization with the quasi-reversibility method, we solve the Cauchy problems for parabolic equations in the presence of noisy data. Using this method, we can get the proper regularization parameter ε that we need in the quasi-reversibility method and obtain the convergence rate of approximate solutions as the noise of amplitude δ tends to zero.
文摘In this article,a level-set approach for solving nonlinear elliptic Cauchy problems with piecewise constant solutions is proposed,which allows the definition of a Tikhonov functional on a space of level-set functions.We provide convergence analysis for the Tikhonov approach,including stability and convergence results.Moreover,a numerical investigation of the proposed Tikhonov regularization method is presented.Newton-type methods are used for the solution of the optimality systems,which can be interpreted as stabilized versions of algorithms in a previous work and yield a substantial improvement in performance.The whole approach is focused on three dimensional models,better suited for real life applications.
基金Supported by National Natural Science Foundation of China (Grant No. 10672062)
文摘If the second order problem u(t) + Bu(t) + Au(t) = f(t), u(0) =u(0) = 0 has L^p-maximal regularity, 1 〈 p 〈 ∞, the analyticity of the corresponding propagator of the sine type is shown by obtaining the estimates of ‖λ(λ^2 + λB + A)^-1‖ and ‖B(λ^2 + λB + A)^-1‖ for λ∈ C with Reλ 〉 ω, where the constant ω≥ 0.
文摘The authors prove well posedness in Gevrey classes of Cauchy problem for nonlinear hyper- bolic equations of constant multiplicity with Holder dependence on the time variable.
基金supported by NSFC(11201380)the Fundamental Research Funds for the Central Universities(XDJK2012B007)+1 种基金Doctor Fund of Southwest University(SWU111021)Educational Fund of Southwest University(2010JY053)
文摘We consider the growth rate and quenching rate of the following problem with singular nonlinearityfor some positive constants b:, b2 (see Theorem 3.3 for the parametersfor some constantsHence, the solution (u, v) quenches at the originx = 0 at the same time '1' (see Theorem 4.3). We also tind various other conditions tor the solution to quench in a finite time and obtain the corresponding decay rate of the solution near the quenching time.
文摘In this paper,the notion of C-semigroup of continuous module homomorphisms on a complete random normal(RN)module is introduced and investigated.The existence and uniqueness of solution to the Cauchy problem with respect to exponentially bounded C-semigroups of continuous module homomorphisms in a complete RN module are established.
基金This work was financially supported by the National United University[Grant Numbers T110M20600].
文摘In this paper,we address 3D inverse Cauchy issues of highly nonlinear elliptic equations in large cuboids by utilizing the new 3D homogenization functions of different orders to adapt all the specified boundary data.We also add the average classification as an approximate solution to the nonlinear operator part,without requiring to cope with nonlinear equations to resolve the weighting coefficients because these constructions are owned many conditions about the true solution.The unknown boundary conditions and the result can be retrieved straightway by coping with a small-scale linear system when the outcome is described by a new 3D homogenization function,which is right to find the numerical solutions with the errors smaller than the level of noise being put on the over-specified Neumann conditions on the bottom of the cuboid.Besides,note that the new homogenization functions method(HFM)does not require dealing with the regularization and highly nonlinear equations.The robustness and accuracy of the HFM are verified by comparing the recovered results of several numerical experiments to the exact solutions in the entire region,even though a very large level of noise 50%is imposed on the over specified Neumann conditions.The numerical errors of our scheme are in the order of O(10^(−1))-O(10^(−4)).
文摘Integrated semigroups which are not exponentially bounded have been studied. The integrated solutions of inhomogeneous and nonlinear abstract Caucby problems are discussed. The reults extends and improves the theorems in [1].
基金Project supported by the National Natural Science Foundation of China(Grant No.41175025)
文摘The simplified linear model of Grad-Shafranov (GS) reconstruction can be reformulated into an inverse boundary value problem of Laplace's equation. Therefore, in this paper we focus on the method of solving the inverse boundary value problem of Laplace's equation. In the first place, the variational regularization method is used to deal with the ill- posedness of the Cauchy problem for Laplace's equation. Then, the 'L-Curve' principle is suggested to be adopted in choosing the optimal regularization parameter. Finally, a numerical experiment is implemented with a section of Neumann and Dirichlet boundary conditions with observation errors. The results well converge to the exact solution of the problem, which proves the efficiency and robustness of the proposed method. When the order of observation error δ is 10-1, the order of the approximate result error can reach 10-3.