The main aim of this paper is to discuss the following two problems: Problem I: Given X ∈ Hn×m (the set of all n×m quaternion matrices), A = diag(λ1,…, λm) EEEEE Hm×m, find A ∈ BSHn×n≥such th...The main aim of this paper is to discuss the following two problems: Problem I: Given X ∈ Hn×m (the set of all n×m quaternion matrices), A = diag(λ1,…, λm) EEEEE Hm×m, find A ∈ BSHn×n≥such that AX = X(?), where BSHn×n≥ denotes the set of all n×n quaternion matrices which are bi-self-conjugate and nonnegative definite. Problem Ⅱ2= Given B ∈ Hn×m, find B ∈ SE such that ||B-B||Q = minAE∈=sE ||B-A||Q, where SE is the solution set of problem I , || ·||Q is the quaternion matrix norm. The necessary and sufficient conditions for SE being nonempty are obtained. The general form of elements in SE and the expression of the unique solution B of problem Ⅱ are given.展开更多
is gained by deleting the k<sup>th</sup> row and the k<sup>th</sup> column (k=1,2,...,n) from T<sub>n</sub>.We put for-ward an inverse eigenvalue problem to be that:If we don’t k...is gained by deleting the k<sup>th</sup> row and the k<sup>th</sup> column (k=1,2,...,n) from T<sub>n</sub>.We put for-ward an inverse eigenvalue problem to be that:If we don’t know the matrix T<sub>1,n</sub>,but weknow all eigenvalues of matrix T<sub>1,k-1</sub>,all eigenvalues of matrix T<sub>k+1,k</sub>,and all eigenvaluesof matrix T<sub>1,n</sub> could we construct the matrix T<sub>1,n</sub>.Let μ<sub>1</sub>,μ<sub>2</sub>,…,μ<sub>k-1</sub>,μ<sub>k</sub>,μ<sub>k+1</sub>,…,μ<sub>n-1</sub>,展开更多
In this paper we first consider the existence and the general form of solution to the following generalized inverse eigenvalue problem(GIEP): given a set of n-dimension complex vectors {x j}m j=1 and a set of co...In this paper we first consider the existence and the general form of solution to the following generalized inverse eigenvalue problem(GIEP): given a set of n-dimension complex vectors {x j}m j=1 and a set of complex numbers {λ j}m j=1, find two n×n centrohermitian matrices A,B such that {x j}m j=1 and {λ j}m j=1 are the generalized eigenvectors and generalized eigenvalues of Ax=λBx, respectively. We then discuss the optimal approximation problem for the GIEP. More concretely, given two arbitrary matrices, , ∈C n×n, we find two matrices A and B such that the matrix (A*,B*) is closest to (,) in the Frobenius norm, where the matrix (A*,B*) is the solution to the GIEP. We show that the expression of the solution of the optimal approximation is unique and derive the expression for it.展开更多
The inverse design method of a dynamic system with linear parameters has been studied. For some specified eigenvalues and eigenvectors, the design parameter vector which is often composed of whole or part of coefficie...The inverse design method of a dynamic system with linear parameters has been studied. For some specified eigenvalues and eigenvectors, the design parameter vector which is often composed of whole or part of coefficients of spring and mass of the system can be obtained and the rigidity and mass matrices of an initially designed structure can be reconstructed through solving linear algebra equations. By using implicit function theorem, the conditions of existence and uniqueness of the solution are also deduced. The theory and method can be used for inverse vibration design of complex structure system.展开更多
An Inverse perturbation method is described for solving the general inverse eigenvalue problem. By taking the analysis of the rotor system as example based upon FEM, the new inverse perturbation method for structural ...An Inverse perturbation method is described for solving the general inverse eigenvalue problem. By taking the analysis of the rotor system as example based upon FEM, the new inverse perturbation method for structural design with specified low-order natural frequencies or frequency constraint bands is detailed as well as its complete theoretical basis. Moreover, formulations to calculate the inverse perturbation parameter ε and method to select the corresponding ε's value properly are also proposed. The proposed method is characterized in reducing frequency analysis and suitable for large and small structrual changes alike. Finally, several different numerical examples for inverse cigenvalue problem are discussed to illustrate the method, which show that this inverse perturbation method Is general and can be applied to other type of structure or dement.展开更多
1 IntroductionLet R<sup>n×n</sup> be the set of all n×n real matrices.R<sup>n</sup>=R<sup>n×1</sup>.C<sup>n×n</sup>denotes the set of all n×n co...1 IntroductionLet R<sup>n×n</sup> be the set of all n×n real matrices.R<sup>n</sup>=R<sup>n×1</sup>.C<sup>n×n</sup>denotes the set of all n×n complex matrices.We are interested in solving the following inverse eigenvalue prob-lems:Problem A (Additive inverse eigenvalue problem) Given an n×n real matrix A=(a<sub>ij</sub>),and n distinct real numbers λ<sub>1</sub>,λ<sub>2</sub>,…,λ<sub>n</sub>,find a real n×n diagonal matrix展开更多
We present a differential geometric perspective of the IEP for symmetric matrices in the framework of a fibre bundle with structure group SO(n). In particular, a Newton type algorithm is developed to construct a non...We present a differential geometric perspective of the IEP for symmetric matrices in the framework of a fibre bundle with structure group SO(n). In particular, a Newton type algorithm is developed to construct a non singular symmetric matrix for given target eigenvalues using a singular symmetric matrix as the initial matrix for the iteration. Explicit computations are performed for 2 x 2 non singular symmetric matrix to illustrate the result.展开更多
Let T1,n be an n x n unreduced symmetric tridiagonal matrix with eigenvaluesand is an (n - 1) x (n - 1) submatrix by deleting the kth row and kth column, k = 1, 2,be the eigenvalues of T1,k andbe the eigenvalues of Tk...Let T1,n be an n x n unreduced symmetric tridiagonal matrix with eigenvaluesand is an (n - 1) x (n - 1) submatrix by deleting the kth row and kth column, k = 1, 2,be the eigenvalues of T1,k andbe the eigenvalues of Tk+1,nA new inverse eigenvalues problem has put forward as follows: How do we construct anunreduced symmetric tridiagonal matrix T1,n, if we only know the spectral data: theeigenvalues of T1,n, the eigenvalues of Ti,k-1 and the eigenvalues of Tk+1,n?Namely if we only know the data: A1, A2, An,how do we find the matrix T1,n? A necessary and sufficient condition and an algorithm ofsolving such problem, are given in this paper.展开更多
In this paper, we discuss an inverse eigenvalue problem for constructing a 2n × 2n Jacobi matrix T such that its 2n eigenvalues are given distinct real values and its leading principal submatrix of order n is a g...In this paper, we discuss an inverse eigenvalue problem for constructing a 2n × 2n Jacobi matrix T such that its 2n eigenvalues are given distinct real values and its leading principal submatrix of order n is a given Jacobi matrix. A new sufficient and necessary condition for the solvability of the above problem is given in this paper. Furthermore, we present a new algorithm and give some numerical results.展开更多
In this paper,the generalized inverse eigenvalue problem for the(P,Q)-conjugate matrices and the associated approximation problem are discussed by using generalized singular value decomposition(GSVD).Moreover,the ...In this paper,the generalized inverse eigenvalue problem for the(P,Q)-conjugate matrices and the associated approximation problem are discussed by using generalized singular value decomposition(GSVD).Moreover,the least residual problem of the above generalized inverse eigenvalue problem is studied by using the canonical correlation decomposition(CCD).The solutions to these problems are derived.Some numerical examples are given to illustrate the main results.展开更多
Utilizing the properties of the smallest singular value of a matrix, we propose a new, efficient and reliable algorithm for solving nonsymmetric matrix inverse eigenvalue problems, and compare it with a known method. ...Utilizing the properties of the smallest singular value of a matrix, we propose a new, efficient and reliable algorithm for solving nonsymmetric matrix inverse eigenvalue problems, and compare it with a known method. We also present numerical experiments which illustrate our results.展开更多
In this note, we consider the backward errors for more general inverse eigenvalue problems by extending Sun's approach. The optimal backward errors are defined for diagonal-ization matrix inverse eigenvalue proble...In this note, we consider the backward errors for more general inverse eigenvalue problems by extending Sun's approach. The optimal backward errors are defined for diagonal-ization matrix inverse eigenvalue problem with respect to an approximate solution, and the upper and lower bounds are derived for the optimal backward errors. The results may be useful for testing the stability of practical algorithms.展开更多
This paper describes inverse eigenvalue problems that arise in studying qualitative dynamics in systems biology models.An algorithm based on lift-andproject iterations is proposed,where the lifting step entails solvin...This paper describes inverse eigenvalue problems that arise in studying qualitative dynamics in systems biology models.An algorithm based on lift-andproject iterations is proposed,where the lifting step entails solving a constrained matrix inverse eigenvalue problem.In particular,prior to carrying out the iterative steps,a-priori bounds on the entries of the Jacobian matrix are computed by relying on the reaction network structure as well as the form of the rate law expressions for the model under consideration.Numerical results on a number of models show that the proposed algorithm can be used to computationally explore the possible dynamical scenarios while identifying the important mechanisms via the use of sparsity-promoting regularization.展开更多
In recent years,numerical solutions of the inverse eigenvalue problems with multiple eigenvalues have attracted the attention of some researchers,and there have been a few algorithms with quadratic convergence.We prop...In recent years,numerical solutions of the inverse eigenvalue problems with multiple eigenvalues have attracted the attention of some researchers,and there have been a few algorithms with quadratic convergence.We propose here an extended two-step method for solving the inverse eigenvalue problems with multiple eigenvalues.Under appropriate assumptions,the convergence analysis of the extended method is presented and the cubic root-convergence rate is proved.Numerical experiments are provided to confirm the theoretical results and comparisons with the inexact Cayley transform method are made.Our extended method and convergence result in the present paper may enrich the results of numerical solutions of the inverse eigenvalue problems with multiple eigenvalues.展开更多
The tridiagonal coefficient matrix for the "fixed-fixed" spring-mass system was obtained by changing spring length. And then a new algorithm of the inverse problem was designed to construct the masses and the spring...The tridiagonal coefficient matrix for the "fixed-fixed" spring-mass system was obtained by changing spring length. And then a new algorithm of the inverse problem was designed to construct the masses and the spring constants from the natural frequencies of the "fixed-fixed" and "fixed-fres" spring-mass systems. An example was given to illustrate the results.展开更多
In this paper, we are concerned with the inverse transmission eigenvalue problem to recover the shape as well as the constant refractive index of a penetrable medium scatterer. The linear sampling method is employed t...In this paper, we are concerned with the inverse transmission eigenvalue problem to recover the shape as well as the constant refractive index of a penetrable medium scatterer. The linear sampling method is employed to determine the transmission eigenvalues within a certain wavenumber interval based on far-field measurements. Based on a prior information given by the linear sampling method, the neural network approach is proposed for the reconstruction of the unknown scatterer. We divide the wavenumber intervals into several subintervals, ensuring that each transmission eigenvalue is located in its corresponding subinterval. In each such subinterval, the wavenumber that yields the maximum value of the indicator functional will be included in the input set during the generation of the training data. This technique for data generation effectively ensures the consistent dimensions of model input. The refractive index and shape are taken as the output of the network. Due to the fact that transmission eigenvalues considered in our method are relatively small,certain super-resolution effects can also be generated. Numerical experiments are presented to verify the effectiveness and promising features of the proposed method in two and three dimensions.展开更多
The Bayesian method of statistical analysis has been applied to the parameter identification problem. A method is presented to identify parameters of dynamic models with the Bayes estimators of measurement frequencies...The Bayesian method of statistical analysis has been applied to the parameter identification problem. A method is presented to identify parameters of dynamic models with the Bayes estimators of measurement frequencies. This is based on the solution of an inverse generalized evaluate problem. The stochastic nature of test data is considered and a normal distribution is used for the measurement frequencies. An additional feature is that the engineer's confidence in the measurement frequencies is quantified and incorporated into the identification procedure. A numerical example demonstrates the efficiency of the method.展开更多
In this paper we first present a CG-type method for inverse eigenvalue problem of constructing real and symmetric matrices M,D and K for the quadratic pencil Q(λ)=λ^(2)M+λD+K,so that Q(λ)has a prescribed subset of...In this paper we first present a CG-type method for inverse eigenvalue problem of constructing real and symmetric matrices M,D and K for the quadratic pencil Q(λ)=λ^(2)M+λD+K,so that Q(λ)has a prescribed subset of eigenvalues and eigenvectors.This method can determine the solvability of the inverse eigenvalue problem automatically.We then consider the least squares model for updating a quadratic pencil Q(λ).More precisely,we update the model coefficient matrices M,C and K so that(i)the updated model reproduces the measured data,(ii)the symmetry of the original model is preserved,and(iii)the difference between the analytical triplet(M,D,K)and the updated triplet(M_(new),D_(new),K_(new))is minimized.In this paper a computationally efficient method is provided for such model updating and numerical examples are given to illustrate the effectiveness of the proposed method.展开更多
Least squares solution of F=PG with respect to positive semidefinite symmetric P is considered,a new necessary and sufficient condition for solvablity is given,and the expression of solution is derived in the some spe...Least squares solution of F=PG with respect to positive semidefinite symmetric P is considered,a new necessary and sufficient condition for solvablity is given,and the expression of solution is derived in the some special cases. Based on the expression, the least spuares solution of an inverse eigenvalue problem for positive semidefinite symmetric matrices is also given.展开更多
基金This work is supported by the NSF of China (10471039, 10271043) and NSF of Zhejiang Province (M103087).
文摘The main aim of this paper is to discuss the following two problems: Problem I: Given X ∈ Hn×m (the set of all n×m quaternion matrices), A = diag(λ1,…, λm) EEEEE Hm×m, find A ∈ BSHn×n≥such that AX = X(?), where BSHn×n≥ denotes the set of all n×n quaternion matrices which are bi-self-conjugate and nonnegative definite. Problem Ⅱ2= Given B ∈ Hn×m, find B ∈ SE such that ||B-B||Q = minAE∈=sE ||B-A||Q, where SE is the solution set of problem I , || ·||Q is the quaternion matrix norm. The necessary and sufficient conditions for SE being nonempty are obtained. The general form of elements in SE and the expression of the unique solution B of problem Ⅱ are given.
基金Project 19771020 supported by National Science Foundation of China
文摘is gained by deleting the k<sup>th</sup> row and the k<sup>th</sup> column (k=1,2,...,n) from T<sub>n</sub>.We put for-ward an inverse eigenvalue problem to be that:If we don’t know the matrix T<sub>1,n</sub>,but weknow all eigenvalues of matrix T<sub>1,k-1</sub>,all eigenvalues of matrix T<sub>k+1,k</sub>,and all eigenvaluesof matrix T<sub>1,n</sub> could we construct the matrix T<sub>1,n</sub>.Let μ<sub>1</sub>,μ<sub>2</sub>,…,μ<sub>k-1</sub>,μ<sub>k</sub>,μ<sub>k+1</sub>,…,μ<sub>n-1</sub>,
文摘In this paper we first consider the existence and the general form of solution to the following generalized inverse eigenvalue problem(GIEP): given a set of n-dimension complex vectors {x j}m j=1 and a set of complex numbers {λ j}m j=1, find two n×n centrohermitian matrices A,B such that {x j}m j=1 and {λ j}m j=1 are the generalized eigenvectors and generalized eigenvalues of Ax=λBx, respectively. We then discuss the optimal approximation problem for the GIEP. More concretely, given two arbitrary matrices, , ∈C n×n, we find two matrices A and B such that the matrix (A*,B*) is closest to (,) in the Frobenius norm, where the matrix (A*,B*) is the solution to the GIEP. We show that the expression of the solution of the optimal approximation is unique and derive the expression for it.
基金Science Developing Plan of Beijing Educational Committee, Beijing Natural Science Fund (No. 3022003), and NationalNatural Science Fund of China(No.50375002)
文摘The inverse design method of a dynamic system with linear parameters has been studied. For some specified eigenvalues and eigenvectors, the design parameter vector which is often composed of whole or part of coefficients of spring and mass of the system can be obtained and the rigidity and mass matrices of an initially designed structure can be reconstructed through solving linear algebra equations. By using implicit function theorem, the conditions of existence and uniqueness of the solution are also deduced. The theory and method can be used for inverse vibration design of complex structure system.
基金This research is supported by China National Natural Science Foundation (CNNSF), Research Grant No. 50128504
文摘An Inverse perturbation method is described for solving the general inverse eigenvalue problem. By taking the analysis of the rotor system as example based upon FEM, the new inverse perturbation method for structural design with specified low-order natural frequencies or frequency constraint bands is detailed as well as its complete theoretical basis. Moreover, formulations to calculate the inverse perturbation parameter ε and method to select the corresponding ε's value properly are also proposed. The proposed method is characterized in reducing frequency analysis and suitable for large and small structrual changes alike. Finally, several different numerical examples for inverse cigenvalue problem are discussed to illustrate the method, which show that this inverse perturbation method Is general and can be applied to other type of structure or dement.
文摘1 IntroductionLet R<sup>n×n</sup> be the set of all n×n real matrices.R<sup>n</sup>=R<sup>n×1</sup>.C<sup>n×n</sup>denotes the set of all n×n complex matrices.We are interested in solving the following inverse eigenvalue prob-lems:Problem A (Additive inverse eigenvalue problem) Given an n×n real matrix A=(a<sub>ij</sub>),and n distinct real numbers λ<sub>1</sub>,λ<sub>2</sub>,…,λ<sub>n</sub>,find a real n×n diagonal matrix
文摘We present a differential geometric perspective of the IEP for symmetric matrices in the framework of a fibre bundle with structure group SO(n). In particular, a Newton type algorithm is developed to construct a non singular symmetric matrix for given target eigenvalues using a singular symmetric matrix as the initial matrix for the iteration. Explicit computations are performed for 2 x 2 non singular symmetric matrix to illustrate the result.
基金Project 19771020 supported by National Science Foundation of China.
文摘Let T1,n be an n x n unreduced symmetric tridiagonal matrix with eigenvaluesand is an (n - 1) x (n - 1) submatrix by deleting the kth row and kth column, k = 1, 2,be the eigenvalues of T1,k andbe the eigenvalues of Tk+1,nA new inverse eigenvalues problem has put forward as follows: How do we construct anunreduced symmetric tridiagonal matrix T1,n, if we only know the spectral data: theeigenvalues of T1,n, the eigenvalues of Ti,k-1 and the eigenvalues of Tk+1,n?Namely if we only know the data: A1, A2, An,how do we find the matrix T1,n? A necessary and sufficient condition and an algorithm ofsolving such problem, are given in this paper.
基金This work was supported by The National Natural Science Foundation of China, under grant 10271074.
文摘In this paper, we discuss an inverse eigenvalue problem for constructing a 2n × 2n Jacobi matrix T such that its 2n eigenvalues are given distinct real values and its leading principal submatrix of order n is a given Jacobi matrix. A new sufficient and necessary condition for the solvability of the above problem is given in this paper. Furthermore, we present a new algorithm and give some numerical results.
基金Supported by the Key Discipline Construction Project of Tianshui Normal University
文摘In this paper,the generalized inverse eigenvalue problem for the(P,Q)-conjugate matrices and the associated approximation problem are discussed by using generalized singular value decomposition(GSVD).Moreover,the least residual problem of the above generalized inverse eigenvalue problem is studied by using the canonical correlation decomposition(CCD).The solutions to these problems are derived.Some numerical examples are given to illustrate the main results.
文摘Utilizing the properties of the smallest singular value of a matrix, we propose a new, efficient and reliable algorithm for solving nonsymmetric matrix inverse eigenvalue problems, and compare it with a known method. We also present numerical experiments which illustrate our results.
文摘In this note, we consider the backward errors for more general inverse eigenvalue problems by extending Sun's approach. The optimal backward errors are defined for diagonal-ization matrix inverse eigenvalue problem with respect to an approximate solution, and the upper and lower bounds are derived for the optimal backward errors. The results may be useful for testing the stability of practical algorithms.
文摘This paper describes inverse eigenvalue problems that arise in studying qualitative dynamics in systems biology models.An algorithm based on lift-andproject iterations is proposed,where the lifting step entails solving a constrained matrix inverse eigenvalue problem.In particular,prior to carrying out the iterative steps,a-priori bounds on the entries of the Jacobian matrix are computed by relying on the reaction network structure as well as the form of the rate law expressions for the model under consideration.Numerical results on a number of models show that the proposed algorithm can be used to computationally explore the possible dynamical scenarios while identifying the important mechanisms via the use of sparsity-promoting regularization.
基金supported by the National Natural Science Foundation of China(Grant 12071441).
文摘In recent years,numerical solutions of the inverse eigenvalue problems with multiple eigenvalues have attracted the attention of some researchers,and there have been a few algorithms with quadratic convergence.We propose here an extended two-step method for solving the inverse eigenvalue problems with multiple eigenvalues.Under appropriate assumptions,the convergence analysis of the extended method is presented and the cubic root-convergence rate is proved.Numerical experiments are provided to confirm the theoretical results and comparisons with the inexact Cayley transform method are made.Our extended method and convergence result in the present paper may enrich the results of numerical solutions of the inverse eigenvalue problems with multiple eigenvalues.
基金Project supported by the National Natural Science Foundation of China(Grant No.10271074)
文摘The tridiagonal coefficient matrix for the "fixed-fixed" spring-mass system was obtained by changing spring length. And then a new algorithm of the inverse problem was designed to construct the masses and the spring constants from the natural frequencies of the "fixed-fixed" and "fixed-fres" spring-mass systems. An example was given to illustrate the results.
基金supported by the Jilin Natural Science Foundation,China(No.20220101040JC)the National Natural Science Foundation of China(No.12271207)+2 种基金supported by the Hong Kong RGC General Research Funds(projects 11311122,12301420 and 11300821)the NSFC/RGC Joint Research Fund(project N-CityU 101/21)the France-Hong Kong ANR/RGC Joint Research Grant,A_CityU203/19.
文摘In this paper, we are concerned with the inverse transmission eigenvalue problem to recover the shape as well as the constant refractive index of a penetrable medium scatterer. The linear sampling method is employed to determine the transmission eigenvalues within a certain wavenumber interval based on far-field measurements. Based on a prior information given by the linear sampling method, the neural network approach is proposed for the reconstruction of the unknown scatterer. We divide the wavenumber intervals into several subintervals, ensuring that each transmission eigenvalue is located in its corresponding subinterval. In each such subinterval, the wavenumber that yields the maximum value of the indicator functional will be included in the input set during the generation of the training data. This technique for data generation effectively ensures the consistent dimensions of model input. The refractive index and shape are taken as the output of the network. Due to the fact that transmission eigenvalues considered in our method are relatively small,certain super-resolution effects can also be generated. Numerical experiments are presented to verify the effectiveness and promising features of the proposed method in two and three dimensions.
文摘The Bayesian method of statistical analysis has been applied to the parameter identification problem. A method is presented to identify parameters of dynamic models with the Bayes estimators of measurement frequencies. This is based on the solution of an inverse generalized evaluate problem. The stochastic nature of test data is considered and a normal distribution is used for the measurement frequencies. An additional feature is that the engineer's confidence in the measurement frequencies is quantified and incorporated into the identification procedure. A numerical example demonstrates the efficiency of the method.
基金Research supported by National Natural Science Foundation of China(10571047 and 10861005)Provincial Natural Science Foundation of Guangxi(0991238)。
文摘In this paper we first present a CG-type method for inverse eigenvalue problem of constructing real and symmetric matrices M,D and K for the quadratic pencil Q(λ)=λ^(2)M+λD+K,so that Q(λ)has a prescribed subset of eigenvalues and eigenvectors.This method can determine the solvability of the inverse eigenvalue problem automatically.We then consider the least squares model for updating a quadratic pencil Q(λ).More precisely,we update the model coefficient matrices M,C and K so that(i)the updated model reproduces the measured data,(ii)the symmetry of the original model is preserved,and(iii)the difference between the analytical triplet(M,D,K)and the updated triplet(M_(new),D_(new),K_(new))is minimized.In this paper a computationally efficient method is provided for such model updating and numerical examples are given to illustrate the effectiveness of the proposed method.
文摘Least squares solution of F=PG with respect to positive semidefinite symmetric P is considered,a new necessary and sufficient condition for solvablity is given,and the expression of solution is derived in the some special cases. Based on the expression, the least spuares solution of an inverse eigenvalue problem for positive semidefinite symmetric matrices is also given.