The purpose of this paper is to show the preconditioned BMinPert algorithm and analyse the practical implementation. Then a posteriori backward error for BGMRES is given. Furthermore, we discuss their applications in ...The purpose of this paper is to show the preconditioned BMinPert algorithm and analyse the practical implementation. Then a posteriori backward error for BGMRES is given. Furthermore, we discuss their applications in color image restoration. The key differences between BMinPert and other methods such as BFGMRES-S(m, p<sub>f</sub>), GsGMRES and BGMRES are illustrated with numerical experiments which expound the advantages of BMinPert in the presence of sensitive data with ill-condition problems.展开更多
In this paper, a new kind of light beam called off-axial elliptical cosine-Gaussian beam (ECosGBs) is defined by using the tensor method. An analytical propagation expression for the ECosGBs passing through axially ...In this paper, a new kind of light beam called off-axial elliptical cosine-Gaussian beam (ECosGBs) is defined by using the tensor method. An analytical propagation expression for the ECosGBs passing through axially nonsymmetrical optical systems is derived by using vector integration. The intensity distributions of ECosGBs on the input plane, on the output plane with the equivalent Fresnel number being equal to 0.1 and on the focal plane are respectively illustrated for the propagation properties. The results indicate that an ECosGB is eventually transformed into an elliptical cosh- Gaussian beam. In other words, ECosGBs and cosh-Gaussian beams act in a reciprocal manner after propagation.展开更多
This paper gives the truncated version of the generalized minimum backward error algorithm(GMBACK)—the incomplete generalized minimum backward perturbation algorithm(IGMBACK)for large nonsymmetric linear systems.It i...This paper gives the truncated version of the generalized minimum backward error algorithm(GMBACK)—the incomplete generalized minimum backward perturbation algorithm(IGMBACK)for large nonsymmetric linear systems.It is based on an incomplete orthogonalization of the Krylov vectors in question,and gives an approximate or quasi-minimum backward perturbation solution over the Krylov subspace.Theoretical properties of IGMBACK including finite termination,existence and uniqueness are discussed in details,and practical implementation issues associated with the IGMBACK algorithm are considered.Numerical experiments show that,the IGMBACK method is usually more efficient than GMBACK and GMRES,and IMBACK,GMBACK often have better convergence performance than GMRES.Specially,for sensitive matrices and right-hand sides being parallel to the left singular vectors corresponding to the smallest singular values of the coefficient matrices,GMRES does not necessarily converge,and IGMBACK,GMBACK usually converge and outperform GMRES.展开更多
This paper extendes the results by E.M. Kasenally([7]) on a Generalized Minimum Backward Error Algorithm for nonsymmetric linear systems Ax = b to the problem in which pertubations are simultaneously permitted on A an...This paper extendes the results by E.M. Kasenally([7]) on a Generalized Minimum Backward Error Algorithm for nonsymmetric linear systems Ax = b to the problem in which pertubations are simultaneously permitted on A and b. The approach adopted by Kasenally has been to view the approximate solution as the exact solution to a perturbed linear system in which changes are permitted to the matrix A only. The new method introduced in this paper is a Krylov subspace iterative method which minimizes the norm of the perturbations to both the observation vector b and the data matrix A and has better performance than the Kasenally's method and the restarted GMRES method([12]). The minimization problem amounts to computing the smallest singular value and the corresponding right singular vector of a low-order upper-Hessenberg matrix. Theoratical properties of the algorithm are discussed and practical implementation issues are considered. The numerical examples are also given.展开更多
The concept of the field of value to localize the spectrum of the iteration matrices of the skew-symmetric iterative methods is further exploited. Obtained formulas are derived to relate the fields of values of the or...The concept of the field of value to localize the spectrum of the iteration matrices of the skew-symmetric iterative methods is further exploited. Obtained formulas are derived to relate the fields of values of the original matrix and the iteration matrix. This allows us to determine theoretically that indefinite nonsymmetric linear systems can be solved by this class of iterative methods.展开更多
文摘The purpose of this paper is to show the preconditioned BMinPert algorithm and analyse the practical implementation. Then a posteriori backward error for BGMRES is given. Furthermore, we discuss their applications in color image restoration. The key differences between BMinPert and other methods such as BFGMRES-S(m, p<sub>f</sub>), GsGMRES and BGMRES are illustrated with numerical experiments which expound the advantages of BMinPert in the presence of sensitive data with ill-condition problems.
文摘In this paper, a new kind of light beam called off-axial elliptical cosine-Gaussian beam (ECosGBs) is defined by using the tensor method. An analytical propagation expression for the ECosGBs passing through axially nonsymmetrical optical systems is derived by using vector integration. The intensity distributions of ECosGBs on the input plane, on the output plane with the equivalent Fresnel number being equal to 0.1 and on the focal plane are respectively illustrated for the propagation properties. The results indicate that an ECosGB is eventually transformed into an elliptical cosh- Gaussian beam. In other words, ECosGBs and cosh-Gaussian beams act in a reciprocal manner after propagation.
文摘This paper gives the truncated version of the generalized minimum backward error algorithm(GMBACK)—the incomplete generalized minimum backward perturbation algorithm(IGMBACK)for large nonsymmetric linear systems.It is based on an incomplete orthogonalization of the Krylov vectors in question,and gives an approximate or quasi-minimum backward perturbation solution over the Krylov subspace.Theoretical properties of IGMBACK including finite termination,existence and uniqueness are discussed in details,and practical implementation issues associated with the IGMBACK algorithm are considered.Numerical experiments show that,the IGMBACK method is usually more efficient than GMBACK and GMRES,and IMBACK,GMBACK often have better convergence performance than GMRES.Specially,for sensitive matrices and right-hand sides being parallel to the left singular vectors corresponding to the smallest singular values of the coefficient matrices,GMRES does not necessarily converge,and IGMBACK,GMBACK usually converge and outperform GMRES.
文摘This paper extendes the results by E.M. Kasenally([7]) on a Generalized Minimum Backward Error Algorithm for nonsymmetric linear systems Ax = b to the problem in which pertubations are simultaneously permitted on A and b. The approach adopted by Kasenally has been to view the approximate solution as the exact solution to a perturbed linear system in which changes are permitted to the matrix A only. The new method introduced in this paper is a Krylov subspace iterative method which minimizes the norm of the perturbations to both the observation vector b and the data matrix A and has better performance than the Kasenally's method and the restarted GMRES method([12]). The minimization problem amounts to computing the smallest singular value and the corresponding right singular vector of a low-order upper-Hessenberg matrix. Theoratical properties of the algorithm are discussed and practical implementation issues are considered. The numerical examples are also given.
文摘The concept of the field of value to localize the spectrum of the iteration matrices of the skew-symmetric iterative methods is further exploited. Obtained formulas are derived to relate the fields of values of the original matrix and the iteration matrix. This allows us to determine theoretically that indefinite nonsymmetric linear systems can be solved by this class of iterative methods.