In this paper,the super-inverse iterative method is proposed to compute the accurateand complete eigen-solutions for anti-plane cracks/notches with multi-materials,arbitrary opening an-gles and various surface conditi...In this paper,the super-inverse iterative method is proposed to compute the accurateand complete eigen-solutions for anti-plane cracks/notches with multi-materials,arbitrary opening an-gles and various surface conditions.Taking the advantage of the knowledge of the variation forms ofthe eigen-functions,a series of numerical techniques are proposed to simplify the computation andspeed up the convergence rate of the inverse iteration.A number of numerical examples are given todemonstrate the excellent accuracy,efficiency and reliability of the proposed approach.展开更多
In this paper, an algorithm based on a shifted inverse power iteration for computing generalized eigenvalues with corresponding eigenvectors of a large scale sparse symmetric positive definite matrix pencil is present...In this paper, an algorithm based on a shifted inverse power iteration for computing generalized eigenvalues with corresponding eigenvectors of a large scale sparse symmetric positive definite matrix pencil is presented. It converges globally with a cubic asymptotic convergence rate, preserves sparsity of the original matrices and is fully parallelizable. The algebraic multilevel itera-tion method (AMLI) is used to improve the efficiency when symmetric positive definite linear equa-tions need to be solved.展开更多
Explicit Exact and Approximate Inverse Preconditioners for solving complex linear systems are introduced. A class of general iterative methods of second order is presented and the selection of iterative parameters is ...Explicit Exact and Approximate Inverse Preconditioners for solving complex linear systems are introduced. A class of general iterative methods of second order is presented and the selection of iterative parameters is discussed. The second order iterative methods behave quite similar to first order methods and the development of efficient preconditioners for solving the original linear system is a decisive factor for making the second order iterative methods superior to the first order iterative methods. Adaptive preconditioned Conjugate Gradient methods using explicit approximate preconditioners for solving efficiently large sparse systems of algebraic equations are also presented. The generalized Approximate Inverse Matrix techniques can be efficiently used in conjunction with explicit iterative schemes leading to effective composite semi-direct solution methods for solving large linear systems of algebraic equations.展开更多
An inverse problem for identification of the coefficient in heat-conduction equation is considered. After reducing the problem to a nonlinear ill-posed operator equation, Newton type iterative methods are considered. ...An inverse problem for identification of the coefficient in heat-conduction equation is considered. After reducing the problem to a nonlinear ill-posed operator equation, Newton type iterative methods are considered. The implicit iterative method is applied to the linearized Newton equation, and the key step in the process is that a new reasonable a posteriori stopping rule for the inner iteration is presented. Numerical experiments for the new method as well as for Tikhonov method and Bakushikskii method are given, and these results show the obvious advantages of the new method over the other ones.展开更多
This paper proposes a sensitivity analysis method for engineering parameters using interval analyses.This method substantially extends the application of interval analysis method.In this scheme,parameter intervals and...This paper proposes a sensitivity analysis method for engineering parameters using interval analyses.This method substantially extends the application of interval analysis method.In this scheme,parameter intervals and decision-making target intervals are determined using the interval analysis method.As an example,an inverse analysis method for uncertainty is presented.The intervals of unknown parameters can be obtained by sampling measured data.Even for limited measured data,robust results can also be obtained with the inverse analysis method,which can be intuitively evaluated by the uncertainty expressed in terms of an interval.For complex nonlinear problems,an iteratively optimized inverse analysis model is proposed.In a given set of loose parameter intervals,all the unknown parameter intervals that satisfy the measured information can be obtained by an iteratively optimized inverse analysis model.The influences of measured precisions and the number of parameters on the results of the inverse analysis are evaluated.Finally,the uniqueness of the interval inverse analysis method is discussed.展开更多
This paper studies to numerical solutions of an inverse heat conduction problem.The effect of algorithms based on the Newton-Tikhonov method and the Newton-implicit iterative method is investigated,and then several mo...This paper studies to numerical solutions of an inverse heat conduction problem.The effect of algorithms based on the Newton-Tikhonov method and the Newton-implicit iterative method is investigated,and then several modifications are presented.Numerical examples show the modified algorithms always work and can greatly reduce the computational costs.展开更多
This paper investigates the numerical solution of the uncertain inverse heat conduction problem. Uncertainties present in the system parameters are modelled through triangular convex normalized fuzzy sets. In the solu...This paper investigates the numerical solution of the uncertain inverse heat conduction problem. Uncertainties present in the system parameters are modelled through triangular convex normalized fuzzy sets. In the solution process, double parametric forms of fuzzy numbers are used with the variational iteration method(VIM). This problem first computes the uncertain temperature distribution in the domain. Next, when the uncertain temperature measurements in the domain are known, the functions describing the uncertain temperature and heat flux on the boundary are reconstructed. Related example problems are solved using the present procedure. We have also compared the present results with those in [Inf.Sci.(2008) 178 1917] along with homotopy perturbation method(HPM) and [Int. Commun. Heat Mass Transfer(2012) 3930] in the special cases to demonstrate the validity and applicability.展开更多
Microscopic traffic simulations are useful for solving various trafficrelated problems,e.g.traffic jams and accidents,local and global environmental and energy problems,maintaining mobility in aging societies,and evac...Microscopic traffic simulations are useful for solving various trafficrelated problems,e.g.traffic jams and accidents,local and global environmental and energy problems,maintaining mobility in aging societies,and evacuation planning for natural as well as man-made disasters.The origin-destination(OD)matrix is often used as the input to represent traffic demands into traffic simulators.In this study,we propose an indirect method for estimating the OD matrix using a traffic simulator as an internal model.The proposed method is designed to output results that are consistent with the input of the simulator.The method consists of the following steps:(1)calculating link traffic volume from the OD matrix,and(2)updating the matrix.The estimated matrix is updated iteratively until it converges to a predefined tolerance level.Numerical experiments are then conducted using the proposed method on a grid network and on a representation of an actual road network.Finally,we discuss the characteristics of the proposed method and the non-negative constraint for the traffic volume.展开更多
In this research article,we interrogate two new modifications in inverse Weierstrass iterative method for estimating all roots of non-linear equation simultaneously.These modifications enables us to accelerate the con...In this research article,we interrogate two new modifications in inverse Weierstrass iterative method for estimating all roots of non-linear equation simultaneously.These modifications enables us to accelerate the convergence order of inverse Weierstrass method from 2 to 3.Convergence analysis proves that the orders of convergence of the two newly constructed inverse methods are 3.Using computer algebra system Mathematica,we find the lower bound of the convergence order and verify it theoretically.Dynamical planes of the inverse simultaneous methods and classical iterative methods are generated using MATLAB(R2011b),to present the global convergence properties of inverse simultaneous iterative methods as compared to classical methods.Some non-linear models are taken from Physics,Chemistry and engineering to demonstrate the performance and efficiency of the newly constructed methods.Computational CPU time,and residual graphs of the methods are provided to present the dominance behavior of our newly constructed methods as compared to existing inverse and classical simultaneous iterative methods in the literature.展开更多
A class of general inverse matrix techniques based on adaptive algorithmic modelling methodologies is derived yielding iterative methods for solving unsymmetric linear systems of irregular structure arising in complex...A class of general inverse matrix techniques based on adaptive algorithmic modelling methodologies is derived yielding iterative methods for solving unsymmetric linear systems of irregular structure arising in complex computational problems in three space dimensions. The proposed class of approximate inverse is chosen as the basis to yield systems on which classic and preconditioned iterative methods are explicitly applied. Optimized versions of the proposed approximate inverse are presented using special storage (k-sweep) techniques leading to economical forms of the approximate inverses. Application of the adaptive algorithmic methodologies on a characteristic nonlinear boundary value problem is discussed and numerical results are given.展开更多
We discuss the incomplete semi-iterative method (ISIM) for an approximate solution of a linear fixed point equations x=Tx+c with a bounded linear operator T acting on a complex Banach space X such that its resolvent h...We discuss the incomplete semi-iterative method (ISIM) for an approximate solution of a linear fixed point equations x=Tx+c with a bounded linear operator T acting on a complex Banach space X such that its resolvent has a pole of order k at the point 1. Sufficient conditions for the convergence of ISIM to a solution of x=Tx+c, where c belongs to the range space of R(I-T) k, are established. We show that the ISIM has an attractive feature that it is usually convergent even when the spectral radius of the operator T is greater than 1 and Ind 1T≥1. Applications in finite Markov chain is considered and illustrative examples are reported, showing the convergence rate of the ISIM is very high.展开更多
Since China power grids have a hierarchical architecture in operation and man- agement, centralized computation patterns are difficult to meet the demands of small-signal-stability analysis of the bulk interconnected ...Since China power grids have a hierarchical architecture in operation and man- agement, centralized computation patterns are difficult to meet the demands of small-signal-stability analysis of the bulk interconnected power systems. A dis- tributed eigenvalue algorithm derived from the inverse iteration method is pro- posed. It can not only obtain eigenvalues and eigenvectors from power system state matrix but also provide participation factors of all generators. In the comput- ing process, the algorithm only requires exchanging data of boundary nodes and a small amount of other information of different regions. Therefore, it is very suitable to be deployed in a WAN (wide area network) based distributed environment. The algorithm has been tested on an IEEE39 system.展开更多
The inverse synthetic aperture radar(ISAR) imaging can be converted into a sparse reconstruction problem and solved by the l_1-norm minimization algorithm. The basis matrix in sparse ISAR imaging is usually characteri...The inverse synthetic aperture radar(ISAR) imaging can be converted into a sparse reconstruction problem and solved by the l_1-norm minimization algorithm. The basis matrix in sparse ISAR imaging is usually characterized by the unknown rotation rate of a moving target, thus the rotation rate and the sparse signal should be jointly estimated. Especially due to the imperfect coarse motion compensation, we consider the phase error correction problem in the context of the sparse signal reconstruction. To address this issue, we propose an iterative reweighted method,which jointly estimates the rotation rate, corrects the phase error and reconstructs a high resolution ISAR image. The proposed method gives a gradual and interweaved iterative process to refine the unknown parameters to achieve the best sparse representation for the ISAR signals. Particularly, in ISAR image reconstruction,the l_1-norm minimization algorithm is sensitive to user parameters.Setting these user parameters are not trivial and the reconstruction performance depends significantly on their choices. Then, we consider an expansion-compression variance-component(ExCoV) based method, which is automatic and demands no prior knowledge about signal-sparsity or measurement-noise levels. Both numerical and electromagnetic data experiments are implemented to show the effectiveness of the proposed method. It is shown that the proposed method can estimate the rotation rate and correct the phase errors simultaneously, and its superior performance is proved in terms of high resolution ISAR image.展开更多
A brief review of the works of the author and his co-authors on the application of nonlinear analysis, numerical and analytical methods for solving the nonlinear inverse problems (synthesis problems) for optimizing th...A brief review of the works of the author and his co-authors on the application of nonlinear analysis, numerical and analytical methods for solving the nonlinear inverse problems (synthesis problems) for optimizing the different types of radiating systems, is presented in the paper. The synthesis problems are formulated in variational statements and further they are reduced to research and numerical solution of nonlinear integral equations of Hammerstein type. The existence theorems are proof, the investigation methods of nonuniqueness problem of solutions and numerical algorithms of finding the optimal solutions are proved.展开更多
In the paper we present a new method to invert the interior structure in the basement or ancient hidden hill by us-ing magnetotelluric (MT) data with seismic data constraint. We first obtain the thickness and resistiv...In the paper we present a new method to invert the interior structure in the basement or ancient hidden hill by us-ing magnetotelluric (MT) data with seismic data constraint. We first obtain the thickness and resistivity of each layer above the basement or buried hill by the inversion of seismic and log data and create a geoelectrical model for the layers above the basement or hidden hill. Then with the reference to the inversion of 1D MT data, a geoelectrical model for the layers below the basement or hidden hill is created. On the basis of the above initial model, we present an effective and practical forward method, i.e., a model-matched approach to conduct forward inversion arithmetic. Finally, by the method of conjugate gradient iteration, a forward and backward iterative cal-culation is made. Taking No. 618 profile of Shengli Oil Field as an example, we have found out that the tectonic information that is unreflective in the seismic data below the basement is better reflected in the inversion result.展开更多
基金The project is supported by the Natural Science Foundation of China.
文摘In this paper,the super-inverse iterative method is proposed to compute the accurateand complete eigen-solutions for anti-plane cracks/notches with multi-materials,arbitrary opening an-gles and various surface conditions.Taking the advantage of the knowledge of the variation forms ofthe eigen-functions,a series of numerical techniques are proposed to simplify the computation andspeed up the convergence rate of the inverse iteration.A number of numerical examples are given todemonstrate the excellent accuracy,efficiency and reliability of the proposed approach.
文摘In this paper, an algorithm based on a shifted inverse power iteration for computing generalized eigenvalues with corresponding eigenvectors of a large scale sparse symmetric positive definite matrix pencil is presented. It converges globally with a cubic asymptotic convergence rate, preserves sparsity of the original matrices and is fully parallelizable. The algebraic multilevel itera-tion method (AMLI) is used to improve the efficiency when symmetric positive definite linear equa-tions need to be solved.
文摘Explicit Exact and Approximate Inverse Preconditioners for solving complex linear systems are introduced. A class of general iterative methods of second order is presented and the selection of iterative parameters is discussed. The second order iterative methods behave quite similar to first order methods and the development of efficient preconditioners for solving the original linear system is a decisive factor for making the second order iterative methods superior to the first order iterative methods. Adaptive preconditioned Conjugate Gradient methods using explicit approximate preconditioners for solving efficiently large sparse systems of algebraic equations are also presented. The generalized Approximate Inverse Matrix techniques can be efficiently used in conjunction with explicit iterative schemes leading to effective composite semi-direct solution methods for solving large linear systems of algebraic equations.
文摘An inverse problem for identification of the coefficient in heat-conduction equation is considered. After reducing the problem to a nonlinear ill-posed operator equation, Newton type iterative methods are considered. The implicit iterative method is applied to the linearized Newton equation, and the key step in the process is that a new reasonable a posteriori stopping rule for the inner iteration is presented. Numerical experiments for the new method as well as for Tikhonov method and Bakushikskii method are given, and these results show the obvious advantages of the new method over the other ones.
基金Supported by the National Natural Science Foundation of China(50978083)the Fundamental Research Funds for the Central Universities(2010B02814)
文摘This paper proposes a sensitivity analysis method for engineering parameters using interval analyses.This method substantially extends the application of interval analysis method.In this scheme,parameter intervals and decision-making target intervals are determined using the interval analysis method.As an example,an inverse analysis method for uncertainty is presented.The intervals of unknown parameters can be obtained by sampling measured data.Even for limited measured data,robust results can also be obtained with the inverse analysis method,which can be intuitively evaluated by the uncertainty expressed in terms of an interval.For complex nonlinear problems,an iteratively optimized inverse analysis model is proposed.In a given set of loose parameter intervals,all the unknown parameter intervals that satisfy the measured information can be obtained by an iteratively optimized inverse analysis model.The influences of measured precisions and the number of parameters on the results of the inverse analysis are evaluated.Finally,the uniqueness of the interval inverse analysis method is discussed.
基金Project supported by the Key Disciplines of Shanghai Municipality (Grant No.S30104)the Shanghai Leading Academic Discipline Project (Grant No.J50101)
文摘This paper studies to numerical solutions of an inverse heat conduction problem.The effect of algorithms based on the Newton-Tikhonov method and the Newton-implicit iterative method is investigated,and then several modifications are presented.Numerical examples show the modified algorithms always work and can greatly reduce the computational costs.
基金the UGC, Government of India, for financial support under the Rajiv Gandhi National Fellowship (RGNF)
文摘This paper investigates the numerical solution of the uncertain inverse heat conduction problem. Uncertainties present in the system parameters are modelled through triangular convex normalized fuzzy sets. In the solution process, double parametric forms of fuzzy numbers are used with the variational iteration method(VIM). This problem first computes the uncertain temperature distribution in the domain. Next, when the uncertain temperature measurements in the domain are known, the functions describing the uncertain temperature and heat flux on the boundary are reconstructed. Related example problems are solved using the present procedure. We have also compared the present results with those in [Inf.Sci.(2008) 178 1917] along with homotopy perturbation method(HPM) and [Int. Commun. Heat Mass Transfer(2012) 3930] in the special cases to demonstrate the validity and applicability.
文摘Microscopic traffic simulations are useful for solving various trafficrelated problems,e.g.traffic jams and accidents,local and global environmental and energy problems,maintaining mobility in aging societies,and evacuation planning for natural as well as man-made disasters.The origin-destination(OD)matrix is often used as the input to represent traffic demands into traffic simulators.In this study,we propose an indirect method for estimating the OD matrix using a traffic simulator as an internal model.The proposed method is designed to output results that are consistent with the input of the simulator.The method consists of the following steps:(1)calculating link traffic volume from the OD matrix,and(2)updating the matrix.The estimated matrix is updated iteratively until it converges to a predefined tolerance level.Numerical experiments are then conducted using the proposed method on a grid network and on a representation of an actual road network.Finally,we discuss the characteristics of the proposed method and the non-negative constraint for the traffic volume.
文摘In this research article,we interrogate two new modifications in inverse Weierstrass iterative method for estimating all roots of non-linear equation simultaneously.These modifications enables us to accelerate the convergence order of inverse Weierstrass method from 2 to 3.Convergence analysis proves that the orders of convergence of the two newly constructed inverse methods are 3.Using computer algebra system Mathematica,we find the lower bound of the convergence order and verify it theoretically.Dynamical planes of the inverse simultaneous methods and classical iterative methods are generated using MATLAB(R2011b),to present the global convergence properties of inverse simultaneous iterative methods as compared to classical methods.Some non-linear models are taken from Physics,Chemistry and engineering to demonstrate the performance and efficiency of the newly constructed methods.Computational CPU time,and residual graphs of the methods are provided to present the dominance behavior of our newly constructed methods as compared to existing inverse and classical simultaneous iterative methods in the literature.
文摘A class of general inverse matrix techniques based on adaptive algorithmic modelling methodologies is derived yielding iterative methods for solving unsymmetric linear systems of irregular structure arising in complex computational problems in three space dimensions. The proposed class of approximate inverse is chosen as the basis to yield systems on which classic and preconditioned iterative methods are explicitly applied. Optimized versions of the proposed approximate inverse are presented using special storage (k-sweep) techniques leading to economical forms of the approximate inverses. Application of the adaptive algorithmic methodologies on a characteristic nonlinear boundary value problem is discussed and numerical results are given.
基金Project1 990 1 0 0 6 supported by National Natural Science Foundation of China,Doctoral Foundation of China,Chi-na Scholarship council and Laboratory of Computational Physics in Beijing of Chinathe second author is also supportedby the State Major Key
文摘We discuss the incomplete semi-iterative method (ISIM) for an approximate solution of a linear fixed point equations x=Tx+c with a bounded linear operator T acting on a complex Banach space X such that its resolvent has a pole of order k at the point 1. Sufficient conditions for the convergence of ISIM to a solution of x=Tx+c, where c belongs to the range space of R(I-T) k, are established. We show that the ISIM has an attractive feature that it is usually convergent even when the spectral radius of the operator T is greater than 1 and Ind 1T≥1. Applications in finite Markov chain is considered and illustrative examples are reported, showing the convergence rate of the ISIM is very high.
基金Supported by the National Basic Research Program of China (Grant No. 2004CB217903)
文摘Since China power grids have a hierarchical architecture in operation and man- agement, centralized computation patterns are difficult to meet the demands of small-signal-stability analysis of the bulk interconnected power systems. A dis- tributed eigenvalue algorithm derived from the inverse iteration method is pro- posed. It can not only obtain eigenvalues and eigenvectors from power system state matrix but also provide participation factors of all generators. In the comput- ing process, the algorithm only requires exchanging data of boundary nodes and a small amount of other information of different regions. Therefore, it is very suitable to be deployed in a WAN (wide area network) based distributed environment. The algorithm has been tested on an IEEE39 system.
基金supported by the National Natural Science Foundation(61302148)
文摘The inverse synthetic aperture radar(ISAR) imaging can be converted into a sparse reconstruction problem and solved by the l_1-norm minimization algorithm. The basis matrix in sparse ISAR imaging is usually characterized by the unknown rotation rate of a moving target, thus the rotation rate and the sparse signal should be jointly estimated. Especially due to the imperfect coarse motion compensation, we consider the phase error correction problem in the context of the sparse signal reconstruction. To address this issue, we propose an iterative reweighted method,which jointly estimates the rotation rate, corrects the phase error and reconstructs a high resolution ISAR image. The proposed method gives a gradual and interweaved iterative process to refine the unknown parameters to achieve the best sparse representation for the ISAR signals. Particularly, in ISAR image reconstruction,the l_1-norm minimization algorithm is sensitive to user parameters.Setting these user parameters are not trivial and the reconstruction performance depends significantly on their choices. Then, we consider an expansion-compression variance-component(ExCoV) based method, which is automatic and demands no prior knowledge about signal-sparsity or measurement-noise levels. Both numerical and electromagnetic data experiments are implemented to show the effectiveness of the proposed method. It is shown that the proposed method can estimate the rotation rate and correct the phase errors simultaneously, and its superior performance is proved in terms of high resolution ISAR image.
文摘A brief review of the works of the author and his co-authors on the application of nonlinear analysis, numerical and analytical methods for solving the nonlinear inverse problems (synthesis problems) for optimizing the different types of radiating systems, is presented in the paper. The synthesis problems are formulated in variational statements and further they are reduced to research and numerical solution of nonlinear integral equations of Hammerstein type. The existence theorems are proof, the investigation methods of nonuniqueness problem of solutions and numerical algorithms of finding the optimal solutions are proved.
文摘In the paper we present a new method to invert the interior structure in the basement or ancient hidden hill by us-ing magnetotelluric (MT) data with seismic data constraint. We first obtain the thickness and resistivity of each layer above the basement or buried hill by the inversion of seismic and log data and create a geoelectrical model for the layers above the basement or hidden hill. Then with the reference to the inversion of 1D MT data, a geoelectrical model for the layers below the basement or hidden hill is created. On the basis of the above initial model, we present an effective and practical forward method, i.e., a model-matched approach to conduct forward inversion arithmetic. Finally, by the method of conjugate gradient iteration, a forward and backward iterative cal-culation is made. Taking No. 618 profile of Shengli Oil Field as an example, we have found out that the tectonic information that is unreflective in the seismic data below the basement is better reflected in the inversion result.