To solve the symmetric positive definite linear system Ax = b on parallel and vector machines, multisplitting methods are considered. Here the s.p.d. (symmetric positive definite) matrix A need not be assumed in a spe...To solve the symmetric positive definite linear system Ax = b on parallel and vector machines, multisplitting methods are considered. Here the s.p.d. (symmetric positive definite) matrix A need not be assumed in a special form (e.g. the dissection form [11]). The main tool for deriving our methods is the diagonally compensated reduction (cf. [1]). The convergence of such methods is also discussed by using this tool. [WT5,5”HZ]展开更多
The symmetric positive definite solutions of matrix equations (AX,XB)=(C,D) and AXB=C are considered in this paper. Necessary and sufficient conditions for the matrix equations to have symmetric positive de...The symmetric positive definite solutions of matrix equations (AX,XB)=(C,D) and AXB=C are considered in this paper. Necessary and sufficient conditions for the matrix equations to have symmetric positive definite solutions are derived using the singular value and the generalized singular value decompositions. The expressions for the general symmetric positive definite solutions are given when certain conditions hold.展开更多
We report a case of remitting seronegative symmetrical synovitis with pitting edema(RS3 PE) syndrome in a 71-year-old woman. She referred to our hospital with finger stiffness, edema of both hands and feet, pain of bi...We report a case of remitting seronegative symmetrical synovitis with pitting edema(RS3 PE) syndrome in a 71-year-old woman. She referred to our hospital with finger stiffness, edema of both hands and feet, pain of bilateral shoulder, wrist, metacarpophalangeal, proximal interphalangeal, and ankle joints. Rheumatoid factor was negative, human leukocyte antigen-B7 antigen was positive. Moreover, matrix metalloproteinase 3(MMP-3) was high. She was diagnosed with RS3 PE syndrome, and treatment with prednisolone(15 mg/d) was started. One week after prednisolone treatment initiation, CRP decreased to negative, and joint pain was almost completely resolved. However, hand stiffness persisted, and MMP-3 level was still high. Thus, prednisolone dose was increased to 20 mg/d, and the stiffness resolved. Twenty days after treatment initiation, MMP-3 was normalized. MMP-3 was more indicative of RS3 PE syndrome symptoms than CRP. Thus, MMP-3 seems to be more sensitive to RS3 PE syndrome symptoms.展开更多
The following is proved: 1) The linear independence of assumed stress modes is the necessary and sufficient condition for the nonsingular flexibility matrix; 2) The equivalent assumed stress modes lead to the identica...The following is proved: 1) The linear independence of assumed stress modes is the necessary and sufficient condition for the nonsingular flexibility matrix; 2) The equivalent assumed stress modes lead to the identical hybrid element. The Hilbert stress subspace of the assumed stress modes is established. So, it is easy to derive the equivalent orthogonal normal stress modes by Schmidt's method. Because of the resulting diagonal flexibility matrix, the identical hybrid element is free from the complex matrix inversion so that the hybrid efficiency, is improved greatly. The numerical examples show that the method is effective.展开更多
Necessary and sufficient conditions are given for the existence of the general solution, the centrosymmetric solution, and the centroskewsymmetric solution to a system of linear matrix equations over an arbitrary skew...Necessary and sufficient conditions are given for the existence of the general solution, the centrosymmetric solution, and the centroskewsymmetric solution to a system of linear matrix equations over an arbitrary skew field. The representations of such the solutions of the system are also derived.展开更多
Let Ω be a finite dimensional central algebra and chart Ω≠2 .The matrix equation AXB-CXD=E over Ω is considered.Necessary and sufficient conditions for the existence of centro(skew)symmetric solutions of the matri...Let Ω be a finite dimensional central algebra and chart Ω≠2 .The matrix equation AXB-CXD=E over Ω is considered.Necessary and sufficient conditions for the existence of centro(skew)symmetric solutions of the matrix equation are given.As a particular case ,the matrix equation X-AXB=C over Ω is also considered.展开更多
This paper discusses the solutions of the linear matrix equation BT X B=Don some linear manifolds.Some necessary and sufficient conditions for the existenceof the solution and the expression of the general solution ar...This paper discusses the solutions of the linear matrix equation BT X B=Don some linear manifolds.Some necessary and sufficient conditions for the existenceof the solution and the expression of the general solution are given.And also someoptimal approximation solutions are discussed.展开更多
Clarke’s matrix has been applied as a phase-mode transformation matrix to three-phase transmission lines substituting the eigenvector matrices. Considering symmetrical untransposed three-phase lines, an actual symmet...Clarke’s matrix has been applied as a phase-mode transformation matrix to three-phase transmission lines substituting the eigenvector matrices. Considering symmetrical untransposed three-phase lines, an actual symmetrical three-phase line on untransposed conditions is associated with Clarke’s matrix for error and frequency scan analyses in this paper. Error analyses are calculated for the eigenvalue diagonal elements obtained from Clarke’s matrix. The eigenvalue off-diagonal elements from the Clarke’s matrix application are compared to the correspondent exact eigenvalues. Based on the characteristic impedance and propagation function values, the frequency scan analyses show that there are great differences between the Clarke’s matrix results and the exact ones, considering frequency values from 10 kHz to 1 MHz. A correction procedure is applied obtaining two new transformation matrices. These matrices lead to good approximated results when compared to the exact ones. With the correction procedure applied to Clarke’s matrix, the relative values of the eigenvalue matrix off-diagonal element obtained from Clarke’s matrix are decreased while the frequency scan results are improved. The steps of correction procedure application are detailed, investigating the influence of each step on the obtained two new phase-mode transformation matrices.展开更多
In this paper,the GH-congruence canonical forms of positive semidefinite and definte inite and definite(need not be self-conjugate)quaternion matrices are given,and a neccessary and sufficientcondition of GH-congruenc...In this paper,the GH-congruence canonical forms of positive semidefinite and definte inite and definite(need not be self-conjugate)quaternion matrices are given,and a neccessary and sufficientcondition of GH-congruence for two positive semidifinite(definite)quaternion matrices isgiven also.Then simultaneous GH-congruence reduced forms for two self-conjugate matri-ces and some result about the simultaneous GH-congruence diagonalization of quaternionmatrices are obtained.展开更多
We demonstrate that, when computing the LDU decomposition (a typical example of a direct solution method), it is possible to obtain the derivative of a determinant with respect to an eigenvalue of a non-symmetric matr...We demonstrate that, when computing the LDU decomposition (a typical example of a direct solution method), it is possible to obtain the derivative of a determinant with respect to an eigenvalue of a non-symmetric matrix. Our proposed method augments an LDU decomposition program with an additional routine to obtain a program for easily evaluating the derivative of a determinant with respect to an eigenvalue. The proposed method follows simply from the process of solving simultaneous linear equations and is particularly effective for band matrices, for which memory requirements are significantly reduced compared to those for dense matrices. We discuss the theory underlying our proposed method and present detailed algorithms for implementing it.展开更多
This paper establishes an improvement on the QL algorithm for a symmetric tridiagonal matrix T so that we can work out the eigenvalues of T faster. Meanwhile, the new algorithm don’t worsen the stability and precisio...This paper establishes an improvement on the QL algorithm for a symmetric tridiagonal matrix T so that we can work out the eigenvalues of T faster. Meanwhile, the new algorithm don’t worsen the stability and precision of the former algorithm.展开更多
QL(QR) method is an efficient method to find eigenvalues of a matrix. Especially we use QL(QR) method to find eigenvalues of a symmetric tridiagonal matrix. In this case it only costs O(n2) flops, to find all eigenval...QL(QR) method is an efficient method to find eigenvalues of a matrix. Especially we use QL(QR) method to find eigenvalues of a symmetric tridiagonal matrix. In this case it only costs O(n2) flops, to find all eigenvalues. So it is one of the most efficient method for symmetric tridiagonal matrices. Many experts have researched it. Even the method is mature, it still has many problems need to be researched. We put forward five problems here. They are: (1) Convergence and convergence rate; (2) The convergence of diagonal elements; (3) Shift designed to produce the eigenvalues in monotone order; (4) QL algorithm with multi-shift; (5) Error bound. We intoduce our works on these problems, some of them were published and some are new.展开更多
It is proved that the set of all symmetric real matrices of order n with eigenvalues lying in the interval(α, β), denoted by Sn(α,β), is convex in Rn×n. With this result, some known results on positive(negati...It is proved that the set of all symmetric real matrices of order n with eigenvalues lying in the interval(α, β), denoted by Sn(α,β), is convex in Rn×n. With this result, some known results on positive(negative) definiteness, and Hurwitz(Shur) stability, as well as the aperiodic property of polytopes of symmetric matrices are generalized, and a series of insightful necessary and sufficient conditions for some general set of symmetric matrices contained in Sn(α,β) are presented,which are directly available for analysis of the positive(negative) definiteness, Hurwitz(Shur) stability and the aperiodic property of a wide class of sets of symmetric matrices.展开更多
A hybrid method is presented for determining maximal eigenvalue and its eigenvector(called eigenpair)of a large,dense,symmetric matrix.Many problems require finding only a small part of the eigenpairs,and some require...A hybrid method is presented for determining maximal eigenvalue and its eigenvector(called eigenpair)of a large,dense,symmetric matrix.Many problems require finding only a small part of the eigenpairs,and some require only the maximal one.In a series of papers,efficient algorithms have been developed by Mufa Chen for computing the maximal eigenpairs of tridiagonal matrices with positive off-diagonal elements.The key idea is to explicitly construet effective initial guess of the maximal eigenpair and then to employ a self-closed iterative algorithm.In this paper we will extend Mufa Chen's algorithm to find maximal eigenpair for a large scale,dense,symmetric matrix.Our strategy is to first convert the underlying matrix into the tridiagonal form by using similarity transformations.We then handle the cases that prevent us from applying Chen's algorithm directly,e.g.,the cases with zero or negative super-or sub-diagonal elements.Serval numerical experiments are carried out to demonstrate the efficiency of the proposed hybrid method.展开更多
In this paper, we obtain a formula for the derivative of a determinant with respect to an eigenvalue in the modified Cholesky decomposition of a symmetric matrix, a characteristic example of a direct solution method i...In this paper, we obtain a formula for the derivative of a determinant with respect to an eigenvalue in the modified Cholesky decomposition of a symmetric matrix, a characteristic example of a direct solution method in computational linear algebra. We apply our proposed formula to a technique used in nonlinear finite-element methods and discuss methods for determining singular points, such as bifurcation points and limit points. In our proposed method, the increment in arc length (or other relevant quantities) may be determined automatically, allowing a reduction in the number of basic parameters. The method is particularly effective for banded matrices, which allow a significant reduction in memory requirements as compared to dense matrices. We discuss the theoretical foundations of our proposed method, present algorithms and programs that implement it, and conduct numerical experiments to investigate its effectiveness.展开更多
As the wireless sensor networks are easily deployable, the volume of sensor applications has been increased widely in various fields of military and commercial areas. In order to attain security on the data exchanged ...As the wireless sensor networks are easily deployable, the volume of sensor applications has been increased widely in various fields of military and commercial areas. In order to attain security on the data exchanged over the network, a hybrid cryptographic mechanism which includes both symmetric and asymmetric cryptographic functions is used. The public key cryptographic ECC security implementation in this paper performs a matrix mapping of data’s at the points on the elliptical curve, which are further encoded using the private symmetric cipher cryptographic algorithm. This security enhancement with the hybrid mechanism of ECC and symmetric cipher cryptographic scheme achieves efficiency in energy conservation of about 7% and 4% compared to the asymmetric and symmetric cipher security implementations in WSN.展开更多
Let DD_0(R)={A∈C^(n×#)||Rea_(ii)Rea_(jj)|≥A_iA_j,i≠j,i,j∈N}.PD_0(R)={A∈C^(n×#)||Rea_(ii)Rea_(kk)|≥A_iA_jA_k,i≠j≠k,i,j,k∈N}. In this paper,we show DD_0(R)PD_0(R),and the conditions under which the nu...Let DD_0(R)={A∈C^(n×#)||Rea_(ii)Rea_(jj)|≥A_iA_j,i≠j,i,j∈N}.PD_0(R)={A∈C^(n×#)||Rea_(ii)Rea_(kk)|≥A_iA_jA_k,i≠j≠k,i,j,k∈N}. In this paper,we show DD_0(R)PD_0(R),and the conditions under which the numbers of eigen vance of A∈PD_0(R)\DD_0(R)are equal to the numbers of a_(ii),i∈N in positive and negative real part respectively.Some couter examples are given which present the condnions can not be omitted.展开更多
The derivation of a diagonally loaded sample-matrix inversion (LSMI) algorithm on the busis of inverse matrix recursion (i.e.LSMI-IMR algorithm) is conducted by reconstructing the recursive formulation of covarian...The derivation of a diagonally loaded sample-matrix inversion (LSMI) algorithm on the busis of inverse matrix recursion (i.e.LSMI-IMR algorithm) is conducted by reconstructing the recursive formulation of covariance matrix. For the new algorithm, diagonal loading is by setting initial inverse matrix without any addition of computation. In addition, a corresponding improved recursive algorithm is presented, which is low computational complexity. This eliminates the complex multiplications of the scalar coefficient and updating matrix, resulting in significant computational savings. Simulations show that the LSMI-IMR algorithm is valid.展开更多
Generalized strictly diagonally dominant matrices play a wide and important role in computational mathematics, mathematical physics, theory of dynamical systems, etc.But it is difficult to judge a matrix is or not gen...Generalized strictly diagonally dominant matrices play a wide and important role in computational mathematics, mathematical physics, theory of dynamical systems, etc.But it is difficult to judge a matrix is or not generalized strictly diagonally dominant matrix.In this paper, by using the properties of α-chain diagonally dominant matrix, we obtain new criteria for judging generalized strictly diagonally dominant matrix, which enlarge the identification range.展开更多
文摘To solve the symmetric positive definite linear system Ax = b on parallel and vector machines, multisplitting methods are considered. Here the s.p.d. (symmetric positive definite) matrix A need not be assumed in a special form (e.g. the dissection form [11]). The main tool for deriving our methods is the diagonally compensated reduction (cf. [1]). The convergence of such methods is also discussed by using this tool. [WT5,5”HZ]
文摘The symmetric positive definite solutions of matrix equations (AX,XB)=(C,D) and AXB=C are considered in this paper. Necessary and sufficient conditions for the matrix equations to have symmetric positive definite solutions are derived using the singular value and the generalized singular value decompositions. The expressions for the general symmetric positive definite solutions are given when certain conditions hold.
文摘We report a case of remitting seronegative symmetrical synovitis with pitting edema(RS3 PE) syndrome in a 71-year-old woman. She referred to our hospital with finger stiffness, edema of both hands and feet, pain of bilateral shoulder, wrist, metacarpophalangeal, proximal interphalangeal, and ankle joints. Rheumatoid factor was negative, human leukocyte antigen-B7 antigen was positive. Moreover, matrix metalloproteinase 3(MMP-3) was high. She was diagnosed with RS3 PE syndrome, and treatment with prednisolone(15 mg/d) was started. One week after prednisolone treatment initiation, CRP decreased to negative, and joint pain was almost completely resolved. However, hand stiffness persisted, and MMP-3 level was still high. Thus, prednisolone dose was increased to 20 mg/d, and the stiffness resolved. Twenty days after treatment initiation, MMP-3 was normalized. MMP-3 was more indicative of RS3 PE syndrome symptoms than CRP. Thus, MMP-3 seems to be more sensitive to RS3 PE syndrome symptoms.
文摘The following is proved: 1) The linear independence of assumed stress modes is the necessary and sufficient condition for the nonsingular flexibility matrix; 2) The equivalent assumed stress modes lead to the identical hybrid element. The Hilbert stress subspace of the assumed stress modes is established. So, it is easy to derive the equivalent orthogonal normal stress modes by Schmidt's method. Because of the resulting diagonal flexibility matrix, the identical hybrid element is free from the complex matrix inversion so that the hybrid efficiency, is improved greatly. The numerical examples show that the method is effective.
基金Supported by the National Natural Science Foundation of China(10471085)
文摘Necessary and sufficient conditions are given for the existence of the general solution, the centrosymmetric solution, and the centroskewsymmetric solution to a system of linear matrix equations over an arbitrary skew field. The representations of such the solutions of the system are also derived.
基金Supported by the Natural Science Foundation of China(10071078)Supported by the Natural Science Foundation of Shandong Province(Q99A08)
文摘Let Ω be a finite dimensional central algebra and chart Ω≠2 .The matrix equation AXB-CXD=E over Ω is considered.Necessary and sufficient conditions for the existence of centro(skew)symmetric solutions of the matrix equation are given.As a particular case ,the matrix equation X-AXB=C over Ω is also considered.
基金This work was supposed by the National Nature Science Foundation of China
文摘This paper discusses the solutions of the linear matrix equation BT X B=Don some linear manifolds.Some necessary and sufficient conditions for the existenceof the solution and the expression of the general solution are given.And also someoptimal approximation solutions are discussed.
文摘Clarke’s matrix has been applied as a phase-mode transformation matrix to three-phase transmission lines substituting the eigenvector matrices. Considering symmetrical untransposed three-phase lines, an actual symmetrical three-phase line on untransposed conditions is associated with Clarke’s matrix for error and frequency scan analyses in this paper. Error analyses are calculated for the eigenvalue diagonal elements obtained from Clarke’s matrix. The eigenvalue off-diagonal elements from the Clarke’s matrix application are compared to the correspondent exact eigenvalues. Based on the characteristic impedance and propagation function values, the frequency scan analyses show that there are great differences between the Clarke’s matrix results and the exact ones, considering frequency values from 10 kHz to 1 MHz. A correction procedure is applied obtaining two new transformation matrices. These matrices lead to good approximated results when compared to the exact ones. With the correction procedure applied to Clarke’s matrix, the relative values of the eigenvalue matrix off-diagonal element obtained from Clarke’s matrix are decreased while the frequency scan results are improved. The steps of correction procedure application are detailed, investigating the influence of each step on the obtained two new phase-mode transformation matrices.
文摘In this paper,the GH-congruence canonical forms of positive semidefinite and definte inite and definite(need not be self-conjugate)quaternion matrices are given,and a neccessary and sufficientcondition of GH-congruence for two positive semidifinite(definite)quaternion matrices isgiven also.Then simultaneous GH-congruence reduced forms for two self-conjugate matri-ces and some result about the simultaneous GH-congruence diagonalization of quaternionmatrices are obtained.
文摘We demonstrate that, when computing the LDU decomposition (a typical example of a direct solution method), it is possible to obtain the derivative of a determinant with respect to an eigenvalue of a non-symmetric matrix. Our proposed method augments an LDU decomposition program with an additional routine to obtain a program for easily evaluating the derivative of a determinant with respect to an eigenvalue. The proposed method follows simply from the process of solving simultaneous linear equations and is particularly effective for band matrices, for which memory requirements are significantly reduced compared to those for dense matrices. We discuss the theory underlying our proposed method and present detailed algorithms for implementing it.
文摘This paper establishes an improvement on the QL algorithm for a symmetric tridiagonal matrix T so that we can work out the eigenvalues of T faster. Meanwhile, the new algorithm don’t worsen the stability and precision of the former algorithm.
文摘QL(QR) method is an efficient method to find eigenvalues of a matrix. Especially we use QL(QR) method to find eigenvalues of a symmetric tridiagonal matrix. In this case it only costs O(n2) flops, to find all eigenvalues. So it is one of the most efficient method for symmetric tridiagonal matrices. Many experts have researched it. Even the method is mature, it still has many problems need to be researched. We put forward five problems here. They are: (1) Convergence and convergence rate; (2) The convergence of diagonal elements; (3) Shift designed to produce the eigenvalues in monotone order; (4) QL algorithm with multi-shift; (5) Error bound. We intoduce our works on these problems, some of them were published and some are new.
文摘It is proved that the set of all symmetric real matrices of order n with eigenvalues lying in the interval(α, β), denoted by Sn(α,β), is convex in Rn×n. With this result, some known results on positive(negative) definiteness, and Hurwitz(Shur) stability, as well as the aperiodic property of polytopes of symmetric matrices are generalized, and a series of insightful necessary and sufficient conditions for some general set of symmetric matrices contained in Sn(α,β) are presented,which are directly available for analysis of the positive(negative) definiteness, Hurwitz(Shur) stability and the aperiodic property of a wide class of sets of symmetric matrices.
基金This work is partially supported by the Special Project on High-Performance Computing of the National Key R&D Program under No.2016YFB0200604the National Natural Science Foundation of China(NSFC)Grant No.11731006,and the NSFC/Hong Kong RRC Joint Research Scheme(NFSC/RGC 11961160718)The work of J.Yang is supported by NSFC-11871264 and Natural Science Foundation of Guangdong Province(2018A0303130123).
文摘A hybrid method is presented for determining maximal eigenvalue and its eigenvector(called eigenpair)of a large,dense,symmetric matrix.Many problems require finding only a small part of the eigenpairs,and some require only the maximal one.In a series of papers,efficient algorithms have been developed by Mufa Chen for computing the maximal eigenpairs of tridiagonal matrices with positive off-diagonal elements.The key idea is to explicitly construet effective initial guess of the maximal eigenpair and then to employ a self-closed iterative algorithm.In this paper we will extend Mufa Chen's algorithm to find maximal eigenpair for a large scale,dense,symmetric matrix.Our strategy is to first convert the underlying matrix into the tridiagonal form by using similarity transformations.We then handle the cases that prevent us from applying Chen's algorithm directly,e.g.,the cases with zero or negative super-or sub-diagonal elements.Serval numerical experiments are carried out to demonstrate the efficiency of the proposed hybrid method.
文摘In this paper, we obtain a formula for the derivative of a determinant with respect to an eigenvalue in the modified Cholesky decomposition of a symmetric matrix, a characteristic example of a direct solution method in computational linear algebra. We apply our proposed formula to a technique used in nonlinear finite-element methods and discuss methods for determining singular points, such as bifurcation points and limit points. In our proposed method, the increment in arc length (or other relevant quantities) may be determined automatically, allowing a reduction in the number of basic parameters. The method is particularly effective for banded matrices, which allow a significant reduction in memory requirements as compared to dense matrices. We discuss the theoretical foundations of our proposed method, present algorithms and programs that implement it, and conduct numerical experiments to investigate its effectiveness.
文摘As the wireless sensor networks are easily deployable, the volume of sensor applications has been increased widely in various fields of military and commercial areas. In order to attain security on the data exchanged over the network, a hybrid cryptographic mechanism which includes both symmetric and asymmetric cryptographic functions is used. The public key cryptographic ECC security implementation in this paper performs a matrix mapping of data’s at the points on the elliptical curve, which are further encoded using the private symmetric cipher cryptographic algorithm. This security enhancement with the hybrid mechanism of ECC and symmetric cipher cryptographic scheme achieves efficiency in energy conservation of about 7% and 4% compared to the asymmetric and symmetric cipher security implementations in WSN.
文摘Let DD_0(R)={A∈C^(n×#)||Rea_(ii)Rea_(jj)|≥A_iA_j,i≠j,i,j∈N}.PD_0(R)={A∈C^(n×#)||Rea_(ii)Rea_(kk)|≥A_iA_jA_k,i≠j≠k,i,j,k∈N}. In this paper,we show DD_0(R)PD_0(R),and the conditions under which the numbers of eigen vance of A∈PD_0(R)\DD_0(R)are equal to the numbers of a_(ii),i∈N in positive and negative real part respectively.Some couter examples are given which present the condnions can not be omitted.
文摘The derivation of a diagonally loaded sample-matrix inversion (LSMI) algorithm on the busis of inverse matrix recursion (i.e.LSMI-IMR algorithm) is conducted by reconstructing the recursive formulation of covariance matrix. For the new algorithm, diagonal loading is by setting initial inverse matrix without any addition of computation. In addition, a corresponding improved recursive algorithm is presented, which is low computational complexity. This eliminates the complex multiplications of the scalar coefficient and updating matrix, resulting in significant computational savings. Simulations show that the LSMI-IMR algorithm is valid.
基金Supported by the National Natural Science Foundation of China(71261010)
文摘Generalized strictly diagonally dominant matrices play a wide and important role in computational mathematics, mathematical physics, theory of dynamical systems, etc.But it is difficult to judge a matrix is or not generalized strictly diagonally dominant matrix.In this paper, by using the properties of α-chain diagonally dominant matrix, we obtain new criteria for judging generalized strictly diagonally dominant matrix, which enlarge the identification range.