The phonon density of states of YBa2Cu3O6.92 is obtained with the method of Fourier decon-volution and Tikhonov-s regularization technique from the experimental specific heat data. The result is compared with that of ...The phonon density of states of YBa2Cu3O6.92 is obtained with the method of Fourier decon-volution and Tikhonov-s regularization technique from the experimental specific heat data. The result is compared with that of neutron inelastic scattering.展开更多
The r\|matrices and classical Poisson structures are constructed for x\| and t n\|constrained flows of the modified Jaulent\|Miodek (MJM) hierarchy.The Lax matrix is used to study the separation of variables method f...The r\|matrices and classical Poisson structures are constructed for x\| and t n\|constrained flows of the modified Jaulent\|Miodek (MJM) hierarchy.The Lax matrix is used to study the separation of variables method for these constrained flows. The Jacobi inversion problem for the MJM equation is obtained through the factorization of the MJM equation and the separability of the constrained flows. This is analogous to separation of variables for solving the MJM equation.展开更多
By using Lax representation, we study the separation of variables for x- and tn-finitedimensional integrable Hamiltonian system (FDIHS) obtained from the factorization of AKNShierarchy. Then the separability of X- and...By using Lax representation, we study the separation of variables for x- and tn-finitedimensional integrable Hamiltonian system (FDIHS) obtained from the factorization of AKNShierarchy. Then the separability of X- and tn-FDIHS and the factorization of AKNS hierarchy give rise to the Jacobi inversion problem for soliton equations in AKNS hierarchy. By a standard Jacobi inversion technique, the soliton equations can be solved in terms of Riemann theta function.展开更多
Solutions of inverse problems are required in various fields of science and engineering. The concept of network inversion has been studied as a neural-network-based solution to inverse problems. In general, inverse pr...Solutions of inverse problems are required in various fields of science and engineering. The concept of network inversion has been studied as a neural-network-based solution to inverse problems. In general, inverse problems are not limited to a real-valued area. Recently, complex-valued neural networks have been actively studied in the field of neural networks. As an extension of network inversion to complex numbers, a complex-valued network inversion has been proposed. Moreover, inverse problems for estimating the parameters of distributed generation systems such as distributed energy plants or smart grids from observed electric circuit data have been studied in the field of natural energy. These emphasize the need to handle complex numbers in an alternating current (AC) circuit. In this paper, the authors propose an application of the complex-valued network inversion to the inverse estimation of a distributed generation. Further, the authors confirm the effectiveness of the complex-valued network inversion on the basis of simulation results.展开更多
In this article,we consider the diffusion equation with multi-term time-fractional derivatives.We first derive,by a subordination principle for the solution,that the solution is positive when the initial value is non-...In this article,we consider the diffusion equation with multi-term time-fractional derivatives.We first derive,by a subordination principle for the solution,that the solution is positive when the initial value is non-negative.As an application,we prove the uniqueness of solution to an inverse problem of determination of the temporally varying source term by integral type information in a subdomain.Finally,several numerical experiments are presented to show the accuracy and efficiency of the algorithm.展开更多
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.展开更多
The presence of non-gray radiative properties in a reheating furnace’s medium that absorbs,emits,and involves non-gray creates more complex radiative heat transfer problems.Furthermore,it adds difficulty to solving t...The presence of non-gray radiative properties in a reheating furnace’s medium that absorbs,emits,and involves non-gray creates more complex radiative heat transfer problems.Furthermore,it adds difficulty to solving the coupled conduction,convection,and radiation problem,leading to suboptimal efficiency that fails to meet real-time control demands.To overcome this difficulty,comparable gray radiative properties of non-gray media are proposed and estimated by solving an inverse problem.However,the required iteration numbers by using a least-squares method are too many and resulted in a very low inverse efficiency.It is necessary to present an efficient method for the equivalence.The Levenberg-Marquardt algorithm is utilized to solve the inverse problem of coupled heat transfer,and the gray-equivalent radiative characteristics are successfully recovered.It is our intention that the issue of low inverse efficiency,which has been observed when the least-squares method is employed,will be resolved.To enhance the performance of the Levenberg-Marquardt algorithm,a modification is implemented for determining the damping factor.Detailed investigations are also conducted to evaluate its accuracy,stability of convergence,efficiency,and robustness of the algorithm.Subsequently,a comparison is made between the results achieved using each method.展开更多
In this paper, the inverse spectral problem of Sturm-Liouville operator with boundary conditions and jump conditions dependent on the spectral parameter is investigated. Firstly, the self-adjointness of the problem an...In this paper, the inverse spectral problem of Sturm-Liouville operator with boundary conditions and jump conditions dependent on the spectral parameter is investigated. Firstly, the self-adjointness of the problem and the eigenvalue properties are given, then the asymptotic formulas of eigenvalues and eigenfunctions are presented. Finally, the uniqueness theorems of the corresponding inverse problems are given by Weyl function theory and inverse spectral data approach.展开更多
Electrical impedance tomography (EIT) aims to reconstruct the conductivity distribution using the boundary measured voltage potential. Traditional regularization based method would suffer from error propagation due to...Electrical impedance tomography (EIT) aims to reconstruct the conductivity distribution using the boundary measured voltage potential. Traditional regularization based method would suffer from error propagation due to the iteration process. The statistical inverse problem method uses statistical inference to estimate unknown parameters. In this article, we develop a nonlinear weighted anisotropic total variation (NWATV) prior density function based on the recently proposed NWATV regularization method. We calculate the corresponding posterior density function, i.e., the solution of the EIT inverse problem in the statistical sense, via a modified Markov chain Monte Carlo (MCMC) sampling. We do numerical experiment to validate the proposed approach.展开更多
Because of the constraint mode of the inversion objective function in the traditional resistivity-inversion method of electromagnetic-propagation resistivity logging while drilling(EPR-LWD),obvious differences appear ...Because of the constraint mode of the inversion objective function in the traditional resistivity-inversion method of electromagnetic-propagation resistivity logging while drilling(EPR-LWD),obvious differences appear in the radial and vertical investigation characteristics between the amplitude-ratio and phase difference,which affect the practical application of EPR-LWD data.In this paper,according to the EPR-LWD data,a self-adaptive constraint resistivity-inversion method,which adopts a self-adaptive constraint weighted expression in the objective function to balance the contributions of the phase difference and amplitude attenuation,is proposed.A particle swarm optimization algorithm is also introduced to eliminate the dependence of the accuracy and convergence on the initial value of the inversion.According to the inversion results of multiple classical formation models for EPR-LWD,the differences between the adaptive constraint inversion-resistivity logs with the traditional amplitude-ratio and the phase difference of the resistivity logs are discussed in detail.The results demonstrate that the adaptive resistivity logs take into account the advantages of the amplitude-ratio logs in the radial investigation and phase difference logs in the vertical resolution.Further,it is superior in thin-layer identification and invasion-effect appraisal compared with the single-amplitude-ratio and phase difference logs.The inversion results can provide a theoretical reference for research on the resistivity-inversion method of electromagnetic wave LWD.展开更多
A class of matrix inverse problems minimizing ‖A-‖ F on the linear manifold l A={A∈R n×m |‖AX-B‖ F=min} is considered. The perturbation analysis of the solution to these problems is carried out. Th...A class of matrix inverse problems minimizing ‖A-‖ F on the linear manifold l A={A∈R n×m |‖AX-B‖ F=min} is considered. The perturbation analysis of the solution to these problems is carried out. The perturbation upper bounds of the solution are given for both the consistent and inconsistent cases. The obtained preturbation upper bounds are with respect to the distance from the perturbed solution to the unperturbed manifold.展开更多
We present a method for solving partial differential equations using artificial neural networks and an adaptive collocation strategy.In this procedure,a coarse grid of training points is used at the initial training s...We present a method for solving partial differential equations using artificial neural networks and an adaptive collocation strategy.In this procedure,a coarse grid of training points is used at the initial training stages,while more points are added at later stages based on the value of the residual at a larger set of evaluation points.This method increases the robustness of the neural network approximation and can result in significant computational savings,particularly when the solution is non-smooth.Numerical results are presented for benchmark problems for scalar-valued PDEs,namely Poisson and Helmholtz equations,as well as for an inverse acoustics problem.展开更多
Model error is one of the key factors restricting the accuracy of numerical weather prediction (NWP). Considering the continuous evolution of the atmosphere, the observed data (ignoring the measurement error) can ...Model error is one of the key factors restricting the accuracy of numerical weather prediction (NWP). Considering the continuous evolution of the atmosphere, the observed data (ignoring the measurement error) can be viewed as a series of solutions of an accurate model governing the actual atmosphere. Model error is represented as an unknown term in the accurate model, thus NWP can be considered as an inverse problem to uncover the unknown error term. The inverse problem models can absorb long periods of observed data to generate model error correction procedures. They thus resolve the deficiency and faultiness of the NWP schemes employing only the initial-time data. In this study we construct two inverse problem models to estimate and extrapolate the time-varying and spatial-varying model errors in both the historical and forecast periods by using recent observations and analogue phenomena of the atmosphere. Numerical experiment on Burgers' equation has illustrated the substantial forecast improvement using inverse problem algorithms. The proposed inverse problem methods of suppressing NWP errors will be useful in future high accuracy applications of NWP.展开更多
Errors inevitably exist in numerical weather prediction (NWP) due to imperfect numeric and physical parameterizations. To eliminate these errors, by considering NWP as an inverse problem, an unknown term in the pred...Errors inevitably exist in numerical weather prediction (NWP) due to imperfect numeric and physical parameterizations. To eliminate these errors, by considering NWP as an inverse problem, an unknown term in the prediction equations can be estimated inversely by using the past data, which are presumed to represent the imperfection of the NWP model (model error, denoted as ME). In this first paper of a two-part series, an iteration method for obtaining the MEs in past intervals is presented, and the results from testing its convergence in idealized experiments are reported. Moreover, two batches of iteration tests were applied in the global forecast system of the Global and Regional Assimilation and Prediction System (GRAPES-GFS) for July-August 2009 and January-February 2010. The datasets associated with the initial conditions and sea surface temperature (SST) were both based on NCEP (National Centers for Environmental Prediction) FNL (final) data. The results showed that 6th h forecast errors were reduced to 10% of their original value after a 20-step iteration. Then, off-line forecast error corrections were estimated linearly based on the 2-month mean MEs and compared with forecast errors. The estimated error corrections agreed well with the forecast errors, but the linear growth rate of the estimation was steeper than the forecast error. The advantage of this iteration method is that the MEs can provide the foundation for online correction. A larger proportion of the forecast errors can be expected to be canceled out by properly introducing the model error correction into GRAPES-GFS.展开更多
Wavelet network, a class of neural network consisting of wavelets, is proposed to solve the inverse kinematics problem in robotic manipulator. A wavelet network suitable for dealing with multi-input and multi-output s...Wavelet network, a class of neural network consisting of wavelets, is proposed to solve the inverse kinematics problem in robotic manipulator. A wavelet network suitable for dealing with multi-input and multi-output system is constructed. The network is optimized by reducing the number of wavelets handling large dimension problem according to the sample data. The algorithms for sparseness analysis of input data and fitting wavelets to the output data with orthogonal method are introduced. Then Levenberg-Marquardt algorithm is used to train the network. Simulation results showed that this method is capable of solving the inverse kinematics problem for PUMA560.展开更多
This paper presents an inverse problem in analytical dynamics. The inverse problem is to construct the Lagrangian when the integrals of a system are given. Firstly, the differential equations are obtained by using the...This paper presents an inverse problem in analytical dynamics. The inverse problem is to construct the Lagrangian when the integrals of a system are given. Firstly, the differential equations are obtained by using the time derivative of the integrals. Secondly, the differential equations can be written in the Lagrange equations under certain conditions and the Lagrangian can be obtained. Finally, two examples are given to illustrate the application of the result.展开更多
A novel method based on the relevance vector machine(RVM) for the inverse scattering problem is presented in this paper.The nonlinearity and the ill-posedness inherent in this problem are simultaneously considered.T...A novel method based on the relevance vector machine(RVM) for the inverse scattering problem is presented in this paper.The nonlinearity and the ill-posedness inherent in this problem are simultaneously considered.The nonlinearity is embodied in the relation between the scattered field and the target property,which can be obtained through the RVM training process.Besides,rather than utilizing regularization,the ill-posed nature of the inversion is naturally accounted for because the RVM can produce a probabilistic output.Simulation results reveal that the proposed RVM-based approach can provide comparative performances in terms of accuracy,convergence,robustness,generalization,and improved performance in terms of sparse property in comparison with the support vector machine(SVM) based approach.展开更多
In the present work, we investigate the inverse problem of reconstructing the parameter of an integro-differential parabolic equation, which comes from pollution problems in porous media, when the final observation is...In the present work, we investigate the inverse problem of reconstructing the parameter of an integro-differential parabolic equation, which comes from pollution problems in porous media, when the final observation is given. We use the optimal control framework to establish both the existence and necessary condition of the minimizer for the cost func- tional. Furthermore, we prove the stability and the local uniqueness of the minimizer. Some numerical results will be presented and discussed.展开更多
The spectral distribution exp( ), where {} are the eigenvalues of the negative Laplacian -△=- in the (x^1,x^2)-plane, is studied for a variety of domains, where -∞< t <∞ and i=(1/2)(-1) . The dependence of (...The spectral distribution exp( ), where {} are the eigenvalues of the negative Laplacian -△=- in the (x^1,x^2)-plane, is studied for a variety of domains, where -∞< t <∞ and i=(1/2)(-1) . The dependence of (t)on the connectivity of a domain and the boundary conditions are analyzed. Particular attention is given to a general bounded domain Ω in R^2 with a smooth boundary Ω, where a finite number of piecewise smooth Dirichlet, Neumann and Robin boundary conditions on the piecewise smooth parts Γj(j = 1,……,n) of Ω are considered such that Some geometrical properties of Ω(e.g., the area of Ω, the total lengths of the boundary, the curvature of its boundary, etc.) are determined, from the asymptotic expansions of (t) for |t| → 0.展开更多
To determine the wall thickness, conductivity and permeability of a ferromagnetic plate, an inverse problem is established with measured values and calculated values of time-domain induced voltage in pulsed eddy curre...To determine the wall thickness, conductivity and permeability of a ferromagnetic plate, an inverse problem is established with measured values and calculated values of time-domain induced voltage in pulsed eddy current testing on the plate. From time-domain analytical expressions of the partial derivatives of induced voltage with respect to parameters,it is deduced that the partial derivatives are approximately linearly dependent. Then the constraints of these parameters are obtained by solving a partial linear differential equation. It is indicated that only the product of conductivity and wall thickness, and the product of relative permeability and wall thickness can be determined accurately through the inverse problem with time-domain induced voltage. In the practical testing, supposing the conductivity of the ferromagnetic plate under test is a fixed value, and then the relative variation of wall thickness between two testing points can be calculated via the ratio of the corresponding inversion results of the product of conductivity and wall thickness. Finally, this method for wall thickness measurement is verified by the experiment results of a carbon steel plate.展开更多
文摘The phonon density of states of YBa2Cu3O6.92 is obtained with the method of Fourier decon-volution and Tikhonov-s regularization technique from the experimental specific heat data. The result is compared with that of neutron inelastic scattering.
基金Supported by the National Basic Research Project forNonlinear Sciences and the Doctorate DissertationFoundation of Tsinghua University
文摘The r\|matrices and classical Poisson structures are constructed for x\| and t n\|constrained flows of the modified Jaulent\|Miodek (MJM) hierarchy.The Lax matrix is used to study the separation of variables method for these constrained flows. The Jacobi inversion problem for the MJM equation is obtained through the factorization of the MJM equation and the separability of the constrained flows. This is analogous to separation of variables for solving the MJM equation.
文摘By using Lax representation, we study the separation of variables for x- and tn-finitedimensional integrable Hamiltonian system (FDIHS) obtained from the factorization of AKNShierarchy. Then the separability of X- and tn-FDIHS and the factorization of AKNS hierarchy give rise to the Jacobi inversion problem for soliton equations in AKNS hierarchy. By a standard Jacobi inversion technique, the soliton equations can be solved in terms of Riemann theta function.
文摘Solutions of inverse problems are required in various fields of science and engineering. The concept of network inversion has been studied as a neural-network-based solution to inverse problems. In general, inverse problems are not limited to a real-valued area. Recently, complex-valued neural networks have been actively studied in the field of neural networks. As an extension of network inversion to complex numbers, a complex-valued network inversion has been proposed. Moreover, inverse problems for estimating the parameters of distributed generation systems such as distributed energy plants or smart grids from observed electric circuit data have been studied in the field of natural energy. These emphasize the need to handle complex numbers in an alternating current (AC) circuit. In this paper, the authors propose an application of the complex-valued network inversion to the inverse estimation of a distributed generation. Further, the authors confirm the effectiveness of the complex-valued network inversion on the basis of simulation results.
基金supported by National Natural Science Foundation of China(12271277)the Open Research Fund of Key Laboratory of Nonlinear Analysis&Applications(Central China Normal University),Ministry of Education,China.
文摘In this article,we consider the diffusion equation with multi-term time-fractional derivatives.We first derive,by a subordination principle for the solution,that the solution is positive when the initial value is non-negative.As an application,we prove the uniqueness of solution to an inverse problem of determination of the temporally varying source term by integral type information in a subdomain.Finally,several numerical experiments are presented to show the accuracy and efficiency of the algorithm.
基金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 the Na⁃tional Natural Science Foundation of China(No.12172078)the Fundamental Research Funds for the Central Univer⁃sities(No.DUT24MS007).
文摘The presence of non-gray radiative properties in a reheating furnace’s medium that absorbs,emits,and involves non-gray creates more complex radiative heat transfer problems.Furthermore,it adds difficulty to solving the coupled conduction,convection,and radiation problem,leading to suboptimal efficiency that fails to meet real-time control demands.To overcome this difficulty,comparable gray radiative properties of non-gray media are proposed and estimated by solving an inverse problem.However,the required iteration numbers by using a least-squares method are too many and resulted in a very low inverse efficiency.It is necessary to present an efficient method for the equivalence.The Levenberg-Marquardt algorithm is utilized to solve the inverse problem of coupled heat transfer,and the gray-equivalent radiative characteristics are successfully recovered.It is our intention that the issue of low inverse efficiency,which has been observed when the least-squares method is employed,will be resolved.To enhance the performance of the Levenberg-Marquardt algorithm,a modification is implemented for determining the damping factor.Detailed investigations are also conducted to evaluate its accuracy,stability of convergence,efficiency,and robustness of the algorithm.Subsequently,a comparison is made between the results achieved using each method.
文摘In this paper, the inverse spectral problem of Sturm-Liouville operator with boundary conditions and jump conditions dependent on the spectral parameter is investigated. Firstly, the self-adjointness of the problem and the eigenvalue properties are given, then the asymptotic formulas of eigenvalues and eigenfunctions are presented. Finally, the uniqueness theorems of the corresponding inverse problems are given by Weyl function theory and inverse spectral data approach.
文摘Electrical impedance tomography (EIT) aims to reconstruct the conductivity distribution using the boundary measured voltage potential. Traditional regularization based method would suffer from error propagation due to the iteration process. The statistical inverse problem method uses statistical inference to estimate unknown parameters. In this article, we develop a nonlinear weighted anisotropic total variation (NWATV) prior density function based on the recently proposed NWATV regularization method. We calculate the corresponding posterior density function, i.e., the solution of the EIT inverse problem in the statistical sense, via a modified Markov chain Monte Carlo (MCMC) sampling. We do numerical experiment to validate the proposed approach.
基金supported by the Foundation of Key Laboratory of Exploration Technology for Oil and Gas Resources of the Ministry of Education, Yangtze University, Wuhan (No. K201812)the Foundation of State Key Laboratory of Petroleum Resources and Prospecting, China University of Petroleum, Beijing (No. PRP/open-1704)the Foundation of Education of Hubei Province, China (No. Q20171304)。
文摘Because of the constraint mode of the inversion objective function in the traditional resistivity-inversion method of electromagnetic-propagation resistivity logging while drilling(EPR-LWD),obvious differences appear in the radial and vertical investigation characteristics between the amplitude-ratio and phase difference,which affect the practical application of EPR-LWD data.In this paper,according to the EPR-LWD data,a self-adaptive constraint resistivity-inversion method,which adopts a self-adaptive constraint weighted expression in the objective function to balance the contributions of the phase difference and amplitude attenuation,is proposed.A particle swarm optimization algorithm is also introduced to eliminate the dependence of the accuracy and convergence on the initial value of the inversion.According to the inversion results of multiple classical formation models for EPR-LWD,the differences between the adaptive constraint inversion-resistivity logs with the traditional amplitude-ratio and the phase difference of the resistivity logs are discussed in detail.The results demonstrate that the adaptive resistivity logs take into account the advantages of the amplitude-ratio logs in the radial investigation and phase difference logs in the vertical resolution.Further,it is superior in thin-layer identification and invasion-effect appraisal compared with the single-amplitude-ratio and phase difference logs.The inversion results can provide a theoretical reference for research on the resistivity-inversion method of electromagnetic wave LWD.
文摘A class of matrix inverse problems minimizing ‖A-‖ F on the linear manifold l A={A∈R n×m |‖AX-B‖ F=min} is considered. The perturbation analysis of the solution to these problems is carried out. The perturbation upper bounds of the solution are given for both the consistent and inconsistent cases. The obtained preturbation upper bounds are with respect to the distance from the perturbed solution to the unperturbed manifold.
文摘We present a method for solving partial differential equations using artificial neural networks and an adaptive collocation strategy.In this procedure,a coarse grid of training points is used at the initial training stages,while more points are added at later stages based on the value of the residual at a larger set of evaluation points.This method increases the robustness of the neural network approximation and can result in significant computational savings,particularly when the solution is non-smooth.Numerical results are presented for benchmark problems for scalar-valued PDEs,namely Poisson and Helmholtz equations,as well as for an inverse acoustics problem.
基金Project supported by the Special Scientific Research Project for Public Interest(Grant No.GYHY201206009)the Fundamental Research Funds for the Central Universities,China(Grant Nos.lzujbky-2012-13 and lzujbky-2013-11)the National Basic Research Program of China(Grant Nos.2012CB955902 and 2013CB430204)
文摘Model error is one of the key factors restricting the accuracy of numerical weather prediction (NWP). Considering the continuous evolution of the atmosphere, the observed data (ignoring the measurement error) can be viewed as a series of solutions of an accurate model governing the actual atmosphere. Model error is represented as an unknown term in the accurate model, thus NWP can be considered as an inverse problem to uncover the unknown error term. The inverse problem models can absorb long periods of observed data to generate model error correction procedures. They thus resolve the deficiency and faultiness of the NWP schemes employing only the initial-time data. In this study we construct two inverse problem models to estimate and extrapolate the time-varying and spatial-varying model errors in both the historical and forecast periods by using recent observations and analogue phenomena of the atmosphere. Numerical experiment on Burgers' equation has illustrated the substantial forecast improvement using inverse problem algorithms. The proposed inverse problem methods of suppressing NWP errors will be useful in future high accuracy applications of NWP.
基金funded by the National Natural Science Foundation Science Fund for Youth (Grant No.41405095)the Key Projects in the National Science and Technology Pillar Program during the Twelfth Fiveyear Plan Period (Grant No.2012BAC22B02)the National Natural Science Foundation Science Fund for Creative Research Groups (Grant No.41221064)
文摘Errors inevitably exist in numerical weather prediction (NWP) due to imperfect numeric and physical parameterizations. To eliminate these errors, by considering NWP as an inverse problem, an unknown term in the prediction equations can be estimated inversely by using the past data, which are presumed to represent the imperfection of the NWP model (model error, denoted as ME). In this first paper of a two-part series, an iteration method for obtaining the MEs in past intervals is presented, and the results from testing its convergence in idealized experiments are reported. Moreover, two batches of iteration tests were applied in the global forecast system of the Global and Regional Assimilation and Prediction System (GRAPES-GFS) for July-August 2009 and January-February 2010. The datasets associated with the initial conditions and sea surface temperature (SST) were both based on NCEP (National Centers for Environmental Prediction) FNL (final) data. The results showed that 6th h forecast errors were reduced to 10% of their original value after a 20-step iteration. Then, off-line forecast error corrections were estimated linearly based on the 2-month mean MEs and compared with forecast errors. The estimated error corrections agreed well with the forecast errors, but the linear growth rate of the estimation was steeper than the forecast error. The advantage of this iteration method is that the MEs can provide the foundation for online correction. A larger proportion of the forecast errors can be expected to be canceled out by properly introducing the model error correction into GRAPES-GFS.
文摘Wavelet network, a class of neural network consisting of wavelets, is proposed to solve the inverse kinematics problem in robotic manipulator. A wavelet network suitable for dealing with multi-input and multi-output system is constructed. The network is optimized by reducing the number of wavelets handling large dimension problem according to the sample data. The algorithms for sparseness analysis of input data and fitting wavelets to the output data with orthogonal method are introduced. Then Levenberg-Marquardt algorithm is used to train the network. Simulation results showed that this method is capable of solving the inverse kinematics problem for PUMA560.
基金Project supported by the National Natural Science Foundation of China (Grant Nos 10272021, 10572021) and the Doctoral Programme Foundation of Institution of Higher Education of China (Grant No 20040007022).
文摘This paper presents an inverse problem in analytical dynamics. The inverse problem is to construct the Lagrangian when the integrals of a system are given. Firstly, the differential equations are obtained by using the time derivative of the integrals. Secondly, the differential equations can be written in the Lagrange equations under certain conditions and the Lagrangian can be obtained. Finally, two examples are given to illustrate the application of the result.
基金Project supported by the National Natural Science Foundation of China (Grant No. 61071022)the Graduate Student Research and Innovation Program of Jiangsu Province,China (Grant No. CXZZ11-0381)
文摘A novel method based on the relevance vector machine(RVM) for the inverse scattering problem is presented in this paper.The nonlinearity and the ill-posedness inherent in this problem are simultaneously considered.The nonlinearity is embodied in the relation between the scattered field and the target property,which can be obtained through the RVM training process.Besides,rather than utilizing regularization,the ill-posed nature of the inversion is naturally accounted for because the RVM can produce a probabilistic output.Simulation results reveal that the proposed RVM-based approach can provide comparative performances in terms of accuracy,convergence,robustness,generalization,and improved performance in terms of sparse property in comparison with the support vector machine(SVM) based approach.
基金supported in part by the CNRST Morocco,the Volkswagen Foundation:Grant number I/79315Hydromed project
文摘In the present work, we investigate the inverse problem of reconstructing the parameter of an integro-differential parabolic equation, which comes from pollution problems in porous media, when the final observation is given. We use the optimal control framework to establish both the existence and necessary condition of the minimizer for the cost func- tional. Furthermore, we prove the stability and the local uniqueness of the minimizer. Some numerical results will be presented and discussed.
文摘The spectral distribution exp( ), where {} are the eigenvalues of the negative Laplacian -△=- in the (x^1,x^2)-plane, is studied for a variety of domains, where -∞< t <∞ and i=(1/2)(-1) . The dependence of (t)on the connectivity of a domain and the boundary conditions are analyzed. Particular attention is given to a general bounded domain Ω in R^2 with a smooth boundary Ω, where a finite number of piecewise smooth Dirichlet, Neumann and Robin boundary conditions on the piecewise smooth parts Γj(j = 1,……,n) of Ω are considered such that Some geometrical properties of Ω(e.g., the area of Ω, the total lengths of the boundary, the curvature of its boundary, etc.) are determined, from the asymptotic expansions of (t) for |t| → 0.
基金supported by the National Defense Basic Technology Research Program of China(Grant No.Z132013T001)
文摘To determine the wall thickness, conductivity and permeability of a ferromagnetic plate, an inverse problem is established with measured values and calculated values of time-domain induced voltage in pulsed eddy current testing on the plate. From time-domain analytical expressions of the partial derivatives of induced voltage with respect to parameters,it is deduced that the partial derivatives are approximately linearly dependent. Then the constraints of these parameters are obtained by solving a partial linear differential equation. It is indicated that only the product of conductivity and wall thickness, and the product of relative permeability and wall thickness can be determined accurately through the inverse problem with time-domain induced voltage. In the practical testing, supposing the conductivity of the ferromagnetic plate under test is a fixed value, and then the relative variation of wall thickness between two testing points can be calculated via the ratio of the corresponding inversion results of the product of conductivity and wall thickness. Finally, this method for wall thickness measurement is verified by the experiment results of a carbon steel plate.