Fast Parallel QR Decomposition of Block-Toeplitz Matrices

A fast algorithm FBTQ is presented which computes the QR factorization a block-Toeplitz matrix A (A∈R) in O(mns3) multiplications. We prove that the QR decomposition of A and the inverse Cholesky decomposition can be computed in parallel using the sametransformation.We also prove that some kind of Toeplltz-block matrices can he transformed into the corresponding block-Toeplitz matrices.

Effective Methods of QR-Decompositions of Square Complex Matrices by Fast Discrete Signal-Induced Heap Transforms

The purpose of this work is to present an effective tool for computing different QR-decompositions of a complex nonsingular square matrix. The concept of the discrete signal-induced heap transform (DsiHT, Grigoryan 2006) is used. This transform is fast, has a unique algorithm for any length of the input vector/signal and can be used with different complex basic 2 × 2 transforms. The DsiHT is zeroing all components of the input signal while moving or heaping the energy of the signal to one component, for instance the first one. We describe three different types of QR-decompositions that use the basic transforms with the T, G, and M-type complex matrices we introduce, as well as without matrices but using analytical formulas. We also present the mixed QR-decomposition, when different type DsiHTs are used in different stages of the algorithm. The number of such decompositions is greater than 3(N-1), for an N × N complex matrix. Examples of the QR-decomposition are described in detail for the 4 × 4 and 6 × 6 complex matrices and compared with the known method of Householder transforms. The precision of the QR-decompositions of N × N matrices, when N are 6, 13, 17, 19, 21, 40, 64, 100, 128, 201, 256, and 400 is also compared. The MATLAB-based scripts of the codes for QR-decompositions by the described DsiHTs are given.

A TSE based design for MMSE and QRD of MIMO systems based on ASIP

A Taylor series expansion(TSE) based design for minimum mean-square error(MMSE) and QR decomposition(QRD) of multi-input and multi-output(MIMO) systems is proposed based on application specific instruction set processor(ASIP), which uses TSE algorithm instead of resource-consuming reciprocal and reciprocal square root(RSR) operations.The aim is to give a high performance implementation for MMSE and QRD in one programmable platform simultaneously.Furthermore, instruction set architecture(ISA) and the allocation of data paths in single instruction multiple data-very long instruction word(SIMD-VLIW) architecture are provided, offering more data parallelism and instruction parallelism for different dimension matrices and operation types.Meanwhile, multiple level numerical precision can be achieved with flexible table size and expansion order in TSE ISA.The ASIP has been implemented to a 28 nm CMOS process and frequency reaches 800 MHz.Experimental results show that the proposed design provides perfect numerical precision within the fixed bit-width of the ASIP, higher matrix processing rate better than the requirements of 5G system and more rate-area efficiency comparable with ASIC implementations.

Orthogonal-triangular decomposition ghost imaging

Ghost imaging(GI)offers great potential with respect to conventional imaging techniques.However,there are still some obstacles for reconstructing images with high quality,especially in the case that the orthogonal measurement matrix is impossible to construct.In this paper,we propose a new scheme based on the orthogonal-triangular(QR)decomposition,named QR decomposition ghost imaging(QRGI)to reconstruct a better image with good quality.In the scheme,we can change the randomly non-orthogonal measurement matrix into orthonormal matrix by performing QR decomposition in two cases.(1)When the random measurement matrix is square,it can be firstly decomposed into an orthogonal matrix Q and an upper triangular matrix R.Then let the off-diagonal values of R equal to 0.0,the diagonal elements of R equal to a constant k,where k is the average of all values of the main diagonal,so the resulting measurement matrix can be obtained.(2)When the random measurement matrix is with full rank,we firstly compute its transpose,and followed with above QR operation.Finally,the image of the object can be reconstructed by correlating the new measurement matrix and corresponding bucket values.Both experimental and simulation results verify the feasibility of the proposed QRGI scheme.Moreover,the results also show that the proposed QRGI scheme could improve the imaging quality comparing to traditional GI(TGI)and differential GI(DGI).Besides,in comparison with the singular value decomposition ghost imaging(SVDGI),the imaging quality and the reconstruction time by using QRGI are similar to those by using SVDGI,while the computing time(the time consuming on the light patterns computation)is substantially shortened.

Generalized cubature quadrature Kalman filters:derivations and extensions

A new Gaussian approximation nonlinear filter called generalized cubature quadrature Kalman filter (GCQKF) is introduced for nonlinear dynamic systems. Based on standard GCQKF, two extensions are developed, namely square root generalized cubature quadrature Kalman filter (SR-GCQKF) and iterated generalized cubature quadrature Kalman filter (I-GCQKF). In SR-GCQKF, the QR decomposition is exploited to alter the Cholesky decomposition and both predicted and filtered error covariances have been propagated in square root format to make sure the numerical stability. In I-GCQKF, the measurement update step is executed iteratively to make full use of the latest measurement and a new terminal criterion is adopted to guarantee the increase of likelihood. Detailed numerical experiments demonstrate the superior performance on both tracking stability and estimation accuracy of I-GCQKF and SR-GCQKF compared with GCQKF.

INTERFERENCE MITIGATION FOR ASYNCHRONOUS COOPERATIVE STBC-OFDM SYSTEMS

This paper addresses the problem of interference mitigation in cooperative Space Time Block Coded Orthogonal Frequency Division Multiplexing (STBC-OFDM) systems in the presence of asyn-chronism. This scheme first preprocesses the received ST codewords to convert the equivalent fading matrix into a suboptimal ordering upper triangular form based on low complexity permutation QR decomposition, and then suppresses the InterCarrier Interference (ICI) and InterSymbol Interference (ISI) by exploiting Successive Interference Cancellation (SIC) technique. Simulation results show that the performance of the proposed algorithm slightly outmatches or asymptotically approaches to that of the existing Minimum Mean Square Error (MMSE) detector depending on the magnitude of the Carrier Frequency Offsets (CFOs) but with less complexity.

On Linear Algebra for Non Mathematics Majors

To find out what knowledge in linear algebra is essential to non-mathematics students, a reverse tracking method was used. Based on practical problems likely to encountered by students in subsequent engineering courses, the minimum contents required has been determined. Rules are proposed to meet the background of most freshman students. An application oriented, easy to understand, computer based text book “Applied Popular Linear Algebra with MATLAB” [1] was published.

The structured total least squares algorithm research for passive location based on angle information

Based on the constrained total least squares （CTLS） passive location algorithm with bearing-only measurements, in this paper, the same passive location problem is transformed into the structured total least squares （STLS） problem.The solution of the STLS problem for passive location can be obtained using the inverse iteration method.It also expatiates that both the STLS algorithm and the CTLS algorithm have the same location mean squares error under certain condition.Finally, the article presents a kind of location and tracking algorithm for moving target by combining STLS location algorithm with Kalman filter （KF）.The efficiency and superiority of the proposed algorithms can be confirmed by computer simulation results.