A new beam broadening synthesis technique for Synthetic Aperture Radar(SAR) antenna array, namely Projection Matrix Algorithm(PMA) is presented. The theory of PMA is introduced firstly, and then the iterative renewed ...A new beam broadening synthesis technique for Synthetic Aperture Radar(SAR) antenna array, namely Projection Matrix Algorithm(PMA) is presented. The theory of PMA is introduced firstly, and then the iterative renewed manner is improved to resolve the unbalance problem under amplitude and phase control. In order to validate the algorithm correct and effective, an actual engineering application example is investigated. The beam synthesis results of 1.0~4.5 times broadening under the phase only control and the amplitude and phase control using improved PMA are given. The results show that the beam directivity, the beam broadening, and the side-lobe level requirements were met. It is demonstrated that the improved PMA was effective and feasible for SAR application.展开更多
A new method for the construction of bivariate matrix valued rational interpolants (BGIRI) on a rectangular grid is presented in [6]. The rational interpolants are of Thiele-type continued fraction form with scalar de...A new method for the construction of bivariate matrix valued rational interpolants (BGIRI) on a rectangular grid is presented in [6]. The rational interpolants are of Thiele-type continued fraction form with scalar denominator. The generalized inverse introduced by [3]is gen-eralized to rectangular matrix case in this paper. An exact error formula for interpolation is ob-tained, which is an extension in matrix form of bivariate scalar and vector valued rational interpola-tion discussed by Siemaszko[l2] and by Gu Chuangqing [7] respectively. By defining row and col-umn-transformation in the sense of the partial inverted differences for matrices, two type matrix algorithms are established to construct corresponding two different BGIRI, which hold for the vec-tor case and the scalar case.展开更多
In this paper, the problem of computing the free distance of Trellis Coded Modulation (TCM) signal sequence has been discussed; a new algorithm-the matrix algorithm is proposed; and the step-number estimation problem ...In this paper, the problem of computing the free distance of Trellis Coded Modulation (TCM) signal sequence has been discussed; a new algorithm-the matrix algorithm is proposed; and the step-number estimation problem for state transmission to compute the free distance of TCM signal sequence has been theoretically solved. The matrix algorithm is derived from the Viterbi algorithm, and is an implementation of Viterbi algorithm in the form of matrix. Compared with other algorithms, the matrix algorithm gains two advantages: (1) The explicit solution, and its relatively less complexity. (2) More reflexible ability to the signal space distance variation. As examples, the results of some TCM signal sequence on AWGN channel and fading channels have been presented.展开更多
We continue to consider one of the cybernetic methods in biology related to the study of DNA chains. Exactly, we are considering the problem of reconstructing the distance matrix for DNA chains. Such a matrix is forme...We continue to consider one of the cybernetic methods in biology related to the study of DNA chains. Exactly, we are considering the problem of reconstructing the distance matrix for DNA chains. Such a matrix is formed on the basis of any of the possible algorithms for determining the distances between DNA chains, as well as any specific object of study. At the same time, for example, the practical programming results show that on an average modern computer, it takes about a day to build such a 30 × 30 matrix for mnDNAs using the Needleman-Wunsch algorithm;therefore, for such a 300 × 300 matrix, about 3 months of continuous computer operation is expected. Thus, even for a relatively small number of species, calculating the distance matrix on conventional computers is hardly feasible and the supercomputers are usually not available. Therefore, we started publishing our variants of the algorithms for calculating the distance between two DNA chains, then we publish algorithms for restoring partially filled matrices, i.e., the inverse problem of matrix processing. Previously, we used the method of branches and boundaries, but in this paper we propose to use another new algorithm for restoring the distance matrix for DNA chains. Our recent work has shown that even greater improvement in the quality of the algorithm can often be achieved without improving the auxiliary heuristics of the branches and boundaries method. Thus, we are improving the algorithms that formulate the greedy function of this method only. .展开更多
During the last decade, Model Checking has proven its efficacy and power in circuit design, network protocol analysis and bug hunting. Recent research on automatic verification has shown that no single model-checking ...During the last decade, Model Checking has proven its efficacy and power in circuit design, network protocol analysis and bug hunting. Recent research on automatic verification has shown that no single model-checking technique has the edge over all others in all application areas. So, it is very difficult to determine which technique is the most suitable for a given model. It is thus sensible to apply different techniques to the same model. However, this is a very tedious and time-consuming task, for each algorithm uses its own description language. Applying Model Checking in software design and verification has been proved very difficult. Software architectures (SA) are engineering artifacts that provide high-level and abstract descriptions of complex software systems. In this paper a Direct Model Checking (DMC) method based on Kripke Structure and Matrix Algorithm is provided. Combined and integrated with domain specific software architecture description languages (ADLs), DMC can be used for computing consistency and other critical properties.展开更多
A novel sparse matrix technique for the numerical analysis of semiconductor devicesand its algorithms are presented.Storage scheme and calculation procedure of the sparse matrixare described in detail.The sparse matri...A novel sparse matrix technique for the numerical analysis of semiconductor devicesand its algorithms are presented.Storage scheme and calculation procedure of the sparse matrixare described in detail.The sparse matrix technique in the device simulation can decrease storagegreatly with less CPU time and its implementation is very easy.Some algorithms and calculationexamples to show the time and space characteristics of the sparse matrix are given.展开更多
To reduce the computational complexity of matrix inversion, which is the majority of processing in many practical applications, two numerically efficient recursive algorithms (called algorithms I and II, respectively...To reduce the computational complexity of matrix inversion, which is the majority of processing in many practical applications, two numerically efficient recursive algorithms (called algorithms I and II, respectively) are presented. Algorithm I is used to calculate the inverse of such a matrix, whose leading principal minors are all nonzero. Algorithm II, whereby, the inverse of an arbitrary nonsingular matrix can be evaluated is derived via improving the algorithm I. The implementation, for algorithm II or I, involves matrix-vector multiplications and vector outer products. These operations are computationally fast and highly parallelizable. MATLAB simulations show that both recursive algorithms are valid.展开更多
A new algorithm of structure random response numerical characteristics, namedas matrix algebra algorithm of structure analysis is presented. Using the algorithm, structurerandom response numerical characteristics can ...A new algorithm of structure random response numerical characteristics, namedas matrix algebra algorithm of structure analysis is presented. Using the algorithm, structurerandom response numerical characteristics can easily be got by directly solving linear matrixequations rather than structure motion differential equations. Moreover, in order to solve thecorresponding linear matrix equations, the numerical integration fast algorithm is presented. Thenaccording to the results, dynamic design and life-span estimation can be done. Besides, the newalgorithm can solve non-proportion damp structure response.展开更多
Consistency test of the judgment matrix is an essential step in the aplication of Analytic Hierarchy Process (AHP). This thesis presents an algerithm concerning the adjustment of a judgment matrix when it failed to pa...Consistency test of the judgment matrix is an essential step in the aplication of Analytic Hierarchy Process (AHP). This thesis presents an algerithm concerning the adjustment of a judgment matrix when it failed to pass the consistency test.展开更多
A direct algorithm is proposed by which one can distinguish whether a matrix is an M-matrix (or H-matrix) or not quickly and effectively. Numerical examples show that it is effective and convincible to distinguish M-m...A direct algorithm is proposed by which one can distinguish whether a matrix is an M-matrix (or H-matrix) or not quickly and effectively. Numerical examples show that it is effective and convincible to distinguish M-matrix (or H-matrix) by using the algorithm.展开更多
In this paper, an improved gradient iterative (GI) algorithm for solving the Lyapunov matrix equations is studied. Convergence of the improved method for any initial value is proved with some conditions. Compared wi...In this paper, an improved gradient iterative (GI) algorithm for solving the Lyapunov matrix equations is studied. Convergence of the improved method for any initial value is proved with some conditions. Compared with the GI algorithm, the improved algorithm reduces computational cost and storage. Finally, the algorithm is tested with GI several numerical examples.展开更多
An axiomatic definition for the generalized inverse matrix Pade approximation (GMPA) is introduced. The matrix rational approximants are of the form of the matrix valued numerator and the scalar denominator. By means...An axiomatic definition for the generalized inverse matrix Pade approximation (GMPA) is introduced. The matrix rational approximants are of the form of the matrix valued numerator and the scalar denominator. By means of generalized inverse for matrices, the ε algorithm for the computation of GMPA is established. The well known Wynn identity for GMPA is proved on the basis of ε algorithm. The η algorithm is defined in a similar way. The equivalence relation between ε algorithm and η algorithm is proposed. Some common examples and a numerical example are given to illustrate the methods in this paper.展开更多
Non-negative matrix factorization (NMF) is a technique for dimensionality reduction by placing non-negativity constraints on the matrix. Based on the PARAFAC model, NMF was extended for three-dimension data decompos...Non-negative matrix factorization (NMF) is a technique for dimensionality reduction by placing non-negativity constraints on the matrix. Based on the PARAFAC model, NMF was extended for three-dimension data decomposition. The three-dimension nonnegative matrix factorization (NMF3) algorithm, which was concise and easy to implement, was given in this paper. The NMF3 algorithm implementation was based on elements but not on vectors. It could decompose a data array directly without unfolding, which was not similar to that the traditional algorithms do, It has been applied to the simulated data array decomposition and obtained reasonable results. It showed that NMF3 could be introduced for curve resolution in chemometrics.展开更多
An efficient and stable structure preserving algorithm, which is a variant of the QR like (SR) algorithm due to Bunse-Gerstner and Mehrmann, is presented for computing the eigenvalues and stable invariant subspaces of...An efficient and stable structure preserving algorithm, which is a variant of the QR like (SR) algorithm due to Bunse-Gerstner and Mehrmann, is presented for computing the eigenvalues and stable invariant subspaces of a Hamiltonian matrix. In the algorithm two strategies are employed, one of which is called dis-unstabilization technique and the other is preprocessing technique. Together with them, a so-called ratio-reduction equation and a backtrack technique are introduced to avoid the instability and breakdown in the original algorithm. It is shown that the new algorithm can overcome the instability and breakdown at low cost. Numerical results have demonstrated that the algorithm is stable and can compute the eigenvalues to very high accuracy.展开更多
K-means algorithm is one of the most widely used algorithms in the clustering analysis. To deal with the problem caused by the random selection of initial center points in the traditional al- gorithm, this paper propo...K-means algorithm is one of the most widely used algorithms in the clustering analysis. To deal with the problem caused by the random selection of initial center points in the traditional al- gorithm, this paper proposes an improved K-means algorithm based on the similarity matrix. The im- proved algorithm can effectively avoid the random selection of initial center points, therefore it can provide effective initial points for clustering process, and reduce the fluctuation of clustering results which are resulted from initial points selections, thus a better clustering quality can be obtained. The experimental results also show that the F-measure of the improved K-means algorithm has been greatly improved and the clustering results are more stable.展开更多
A recursive rational algorithm for matrix exponentials was obtained by making use of the generalized inverse of a matrix in this paper. On the basis of the n th convergence of Thiele type continued fraction expa...A recursive rational algorithm for matrix exponentials was obtained by making use of the generalized inverse of a matrix in this paper. On the basis of the n th convergence of Thiele type continued fraction expansion, a new type of the generalized inverse matrix valued Padé approximant (GMPA) for matrix exponentials was defined and its remainder formula was proved. The results of this paper were illustrated by some examples.展开更多
Differential evolution algorithm based on the covariance matrix learning can adjust the coordinate system according to the characteristics of the population, which make<span style="font-family:Verdana;"&g...Differential evolution algorithm based on the covariance matrix learning can adjust the coordinate system according to the characteristics of the population, which make<span style="font-family:Verdana;">s</span><span style="font-family:Verdana;"> the search move in a more favorable direction. In order to obtain more accurate information about the function shape, this paper propose</span><span style="font-family:Verdana;">s</span><span style="font-family:;" "=""> <span style="font-family:Verdana;">covariance</span><span style="font-family:Verdana;"> matrix learning differential evolution algorithm based on correlation (denoted as RCLDE)</span></span><span style="font-family:;" "=""> </span><span style="font-family:Verdana;">to improve the search efficiency of the algorithm. First, a hybrid mutation strategy is designed to balance the diversity and convergence of the population;secondly, the covariance learning matrix is constructed by selecting the individual with the less correlation;then, a comprehensive learning mechanism is comprehensively designed by two covariance matrix learning mechanisms based on the principle of probability. Finally,</span><span style="font-family:;" "=""> </span><span style="font-family:;" "=""><span style="font-family:Verdana;">the algorithm is tested on the CEC2005, and the experimental results are compared with other effective differential evolution algorithms. The experimental results show that the algorithm proposed in this paper is </span><span style="font-family:Verdana;">an effective algorithm</span><span style="font-family:Verdana;">.</span></span>展开更多
In this paper, an iterative algorithm is presented to solve the Sylvester and Lyapunov matrix equations. By this iterative algorithm, for any initial matrix X1, a solution X* can be obtained within finite iteration s...In this paper, an iterative algorithm is presented to solve the Sylvester and Lyapunov matrix equations. By this iterative algorithm, for any initial matrix X1, a solution X* can be obtained within finite iteration steps in the absence of roundoff errors. Some examples illustrate that this algorithm is very efficient and better than that of [ 1 ] and [2].展开更多
In this paper, a practical Werner-type continued fraction method for solving matrix valued rational interpolation problem is provided by using a generalized inverse of matrices. In order to reduce the continued fracti...In this paper, a practical Werner-type continued fraction method for solving matrix valued rational interpolation problem is provided by using a generalized inverse of matrices. In order to reduce the continued fraction form to rational function form of the interpolants, an efficient forward recurrence algorithm is obtained.展开更多
文摘A new beam broadening synthesis technique for Synthetic Aperture Radar(SAR) antenna array, namely Projection Matrix Algorithm(PMA) is presented. The theory of PMA is introduced firstly, and then the iterative renewed manner is improved to resolve the unbalance problem under amplitude and phase control. In order to validate the algorithm correct and effective, an actual engineering application example is investigated. The beam synthesis results of 1.0~4.5 times broadening under the phase only control and the amplitude and phase control using improved PMA are given. The results show that the beam directivity, the beam broadening, and the side-lobe level requirements were met. It is demonstrated that the improved PMA was effective and feasible for SAR application.
文摘A new method for the construction of bivariate matrix valued rational interpolants (BGIRI) on a rectangular grid is presented in [6]. The rational interpolants are of Thiele-type continued fraction form with scalar denominator. The generalized inverse introduced by [3]is gen-eralized to rectangular matrix case in this paper. An exact error formula for interpolation is ob-tained, which is an extension in matrix form of bivariate scalar and vector valued rational interpola-tion discussed by Siemaszko[l2] and by Gu Chuangqing [7] respectively. By defining row and col-umn-transformation in the sense of the partial inverted differences for matrices, two type matrix algorithms are established to construct corresponding two different BGIRI, which hold for the vec-tor case and the scalar case.
文摘In this paper, the problem of computing the free distance of Trellis Coded Modulation (TCM) signal sequence has been discussed; a new algorithm-the matrix algorithm is proposed; and the step-number estimation problem for state transmission to compute the free distance of TCM signal sequence has been theoretically solved. The matrix algorithm is derived from the Viterbi algorithm, and is an implementation of Viterbi algorithm in the form of matrix. Compared with other algorithms, the matrix algorithm gains two advantages: (1) The explicit solution, and its relatively less complexity. (2) More reflexible ability to the signal space distance variation. As examples, the results of some TCM signal sequence on AWGN channel and fading channels have been presented.
文摘We continue to consider one of the cybernetic methods in biology related to the study of DNA chains. Exactly, we are considering the problem of reconstructing the distance matrix for DNA chains. Such a matrix is formed on the basis of any of the possible algorithms for determining the distances between DNA chains, as well as any specific object of study. At the same time, for example, the practical programming results show that on an average modern computer, it takes about a day to build such a 30 × 30 matrix for mnDNAs using the Needleman-Wunsch algorithm;therefore, for such a 300 × 300 matrix, about 3 months of continuous computer operation is expected. Thus, even for a relatively small number of species, calculating the distance matrix on conventional computers is hardly feasible and the supercomputers are usually not available. Therefore, we started publishing our variants of the algorithms for calculating the distance between two DNA chains, then we publish algorithms for restoring partially filled matrices, i.e., the inverse problem of matrix processing. Previously, we used the method of branches and boundaries, but in this paper we propose to use another new algorithm for restoring the distance matrix for DNA chains. Our recent work has shown that even greater improvement in the quality of the algorithm can often be achieved without improving the auxiliary heuristics of the branches and boundaries method. Thus, we are improving the algorithms that formulate the greedy function of this method only. .
文摘During the last decade, Model Checking has proven its efficacy and power in circuit design, network protocol analysis and bug hunting. Recent research on automatic verification has shown that no single model-checking technique has the edge over all others in all application areas. So, it is very difficult to determine which technique is the most suitable for a given model. It is thus sensible to apply different techniques to the same model. However, this is a very tedious and time-consuming task, for each algorithm uses its own description language. Applying Model Checking in software design and verification has been proved very difficult. Software architectures (SA) are engineering artifacts that provide high-level and abstract descriptions of complex software systems. In this paper a Direct Model Checking (DMC) method based on Kripke Structure and Matrix Algorithm is provided. Combined and integrated with domain specific software architecture description languages (ADLs), DMC can be used for computing consistency and other critical properties.
文摘A novel sparse matrix technique for the numerical analysis of semiconductor devicesand its algorithms are presented.Storage scheme and calculation procedure of the sparse matrixare described in detail.The sparse matrix technique in the device simulation can decrease storagegreatly with less CPU time and its implementation is very easy.Some algorithms and calculationexamples to show the time and space characteristics of the sparse matrix are given.
文摘To reduce the computational complexity of matrix inversion, which is the majority of processing in many practical applications, two numerically efficient recursive algorithms (called algorithms I and II, respectively) are presented. Algorithm I is used to calculate the inverse of such a matrix, whose leading principal minors are all nonzero. Algorithm II, whereby, the inverse of an arbitrary nonsingular matrix can be evaluated is derived via improving the algorithm I. The implementation, for algorithm II or I, involves matrix-vector multiplications and vector outer products. These operations are computationally fast and highly parallelizable. MATLAB simulations show that both recursive algorithms are valid.
基金This project is supported by National Natural Science Foundation of China (No.59805001)
文摘A new algorithm of structure random response numerical characteristics, namedas matrix algebra algorithm of structure analysis is presented. Using the algorithm, structurerandom response numerical characteristics can easily be got by directly solving linear matrixequations rather than structure motion differential equations. Moreover, in order to solve thecorresponding linear matrix equations, the numerical integration fast algorithm is presented. Thenaccording to the results, dynamic design and life-span estimation can be done. Besides, the newalgorithm can solve non-proportion damp structure response.
文摘Consistency test of the judgment matrix is an essential step in the aplication of Analytic Hierarchy Process (AHP). This thesis presents an algerithm concerning the adjustment of a judgment matrix when it failed to pass the consistency test.
基金Foundation item: This work is supported by the Science Foundations of the Education Department of Yunnan Province (03Z169A)the Science Foundatons of Yunnan University (2003Z013B).
文摘A direct algorithm is proposed by which one can distinguish whether a matrix is an M-matrix (or H-matrix) or not quickly and effectively. Numerical examples show that it is effective and convincible to distinguish M-matrix (or H-matrix) by using the algorithm.
基金Project supported by the National Natural Science Foundation of China (Grant No.10271074), and the Special Funds for Major Specialities of Shanghai Education Commission (Grant No.J50101)
文摘In this paper, an improved gradient iterative (GI) algorithm for solving the Lyapunov matrix equations is studied. Convergence of the improved method for any initial value is proved with some conditions. Compared with the GI algorithm, the improved algorithm reduces computational cost and storage. Finally, the algorithm is tested with GI several numerical examples.
文摘An axiomatic definition for the generalized inverse matrix Pade approximation (GMPA) is introduced. The matrix rational approximants are of the form of the matrix valued numerator and the scalar denominator. By means of generalized inverse for matrices, the ε algorithm for the computation of GMPA is established. The well known Wynn identity for GMPA is proved on the basis of ε algorithm. The η algorithm is defined in a similar way. The equivalence relation between ε algorithm and η algorithm is proposed. Some common examples and a numerical example are given to illustrate the methods in this paper.
文摘Non-negative matrix factorization (NMF) is a technique for dimensionality reduction by placing non-negativity constraints on the matrix. Based on the PARAFAC model, NMF was extended for three-dimension data decomposition. The three-dimension nonnegative matrix factorization (NMF3) algorithm, which was concise and easy to implement, was given in this paper. The NMF3 algorithm implementation was based on elements but not on vectors. It could decompose a data array directly without unfolding, which was not similar to that the traditional algorithms do, It has been applied to the simulated data array decomposition and obtained reasonable results. It showed that NMF3 could be introduced for curve resolution in chemometrics.
文摘An efficient and stable structure preserving algorithm, which is a variant of the QR like (SR) algorithm due to Bunse-Gerstner and Mehrmann, is presented for computing the eigenvalues and stable invariant subspaces of a Hamiltonian matrix. In the algorithm two strategies are employed, one of which is called dis-unstabilization technique and the other is preprocessing technique. Together with them, a so-called ratio-reduction equation and a backtrack technique are introduced to avoid the instability and breakdown in the original algorithm. It is shown that the new algorithm can overcome the instability and breakdown at low cost. Numerical results have demonstrated that the algorithm is stable and can compute the eigenvalues to very high accuracy.
文摘K-means algorithm is one of the most widely used algorithms in the clustering analysis. To deal with the problem caused by the random selection of initial center points in the traditional al- gorithm, this paper proposes an improved K-means algorithm based on the similarity matrix. The im- proved algorithm can effectively avoid the random selection of initial center points, therefore it can provide effective initial points for clustering process, and reduce the fluctuation of clustering results which are resulted from initial points selections, thus a better clustering quality can be obtained. The experimental results also show that the F-measure of the improved K-means algorithm has been greatly improved and the clustering results are more stable.
文摘A recursive rational algorithm for matrix exponentials was obtained by making use of the generalized inverse of a matrix in this paper. On the basis of the n th convergence of Thiele type continued fraction expansion, a new type of the generalized inverse matrix valued Padé approximant (GMPA) for matrix exponentials was defined and its remainder formula was proved. The results of this paper were illustrated by some examples.
文摘Differential evolution algorithm based on the covariance matrix learning can adjust the coordinate system according to the characteristics of the population, which make<span style="font-family:Verdana;">s</span><span style="font-family:Verdana;"> the search move in a more favorable direction. In order to obtain more accurate information about the function shape, this paper propose</span><span style="font-family:Verdana;">s</span><span style="font-family:;" "=""> <span style="font-family:Verdana;">covariance</span><span style="font-family:Verdana;"> matrix learning differential evolution algorithm based on correlation (denoted as RCLDE)</span></span><span style="font-family:;" "=""> </span><span style="font-family:Verdana;">to improve the search efficiency of the algorithm. First, a hybrid mutation strategy is designed to balance the diversity and convergence of the population;secondly, the covariance learning matrix is constructed by selecting the individual with the less correlation;then, a comprehensive learning mechanism is comprehensively designed by two covariance matrix learning mechanisms based on the principle of probability. Finally,</span><span style="font-family:;" "=""> </span><span style="font-family:;" "=""><span style="font-family:Verdana;">the algorithm is tested on the CEC2005, and the experimental results are compared with other effective differential evolution algorithms. The experimental results show that the algorithm proposed in this paper is </span><span style="font-family:Verdana;">an effective algorithm</span><span style="font-family:Verdana;">.</span></span>
基金supported by the National Natural Science Foundation of China (No.10771073)
文摘In this paper, an iterative algorithm is presented to solve the Sylvester and Lyapunov matrix equations. By this iterative algorithm, for any initial matrix X1, a solution X* can be obtained within finite iteration steps in the absence of roundoff errors. Some examples illustrate that this algorithm is very efficient and better than that of [ 1 ] and [2].
文摘In this paper, a practical Werner-type continued fraction method for solving matrix valued rational interpolation problem is provided by using a generalized inverse of matrices. In order to reduce the continued fraction form to rational function form of the interpolants, an efficient forward recurrence algorithm is obtained.