Efficient Inverse Cholesky Factorization for Alamouti Matrices in G-STBC and Alamouti-like Matrices in OMP
|
|
- Piers Glenn
- 5 years ago
- Views:
Transcription
1 Efficient Inverse Cholesky Factorization for Alamouti Matrices in G-STBC and Alamouti-like Matrices in OMP Hufei Zhu, Ganghua Yang Communications Technology Laboratory Huawei Technologies Co Ltd, P R China zhuhufei@huaweicom, yangyangganghua@huaweicom Wen Chen the Department of Electronic Engineering Shanghai Jiao Tong University, Shanghai, P R China wenchen@sjtueducn Abstract We apply the efficient inverse Cholesky factorization to Alamouti matrices in Groupwise Space-time Block Code (G-STBC) and Alamouti-like matrices in Orthogonal Matching Pursuit (OMP) for the sub-nyquist sampling system By utilizing some good properties of Alamouti or Alamouti-like matrices, we save about half the complexity Then we propose the whole square-root algorithms for G-STBC and OMP, respectively The proposed square-root G-STBC algorithm has an average speedup of with respect to the sub-optimal G-STBC algorithm On the other hand, when comparing the complexities of all the steps except the projection step, the complexity for the proposed square-root OMP algorithm is about 30% of that for the fast OMP algorithm by matrix inverse update I INTRODUCTION The equivalent channel matrix in Group-wise spacetime block code (G-STBC) is an Alamouti matrix ] 3], ie, a block] matrix involving 2 2 Alamouti sub-blocks c c 2 c 2 c (where c and c 2 are complex scalars) Alamouti matrices have some good properties 3], which have been utilized to develop efficient matrix manipulations ] and a fast recursive G-STBC algorithm 2] On the other hand, in the sub-nyquist sampling system 4], the Orthogonal Matching Pursuit (OMP) algorithm 5] 7] computes the inversion of Alamouti-like matrices, ie, block ] c c matrices involving 2 2 Alamouti-like sub-blocks 2 c 2 c (that is obtained by deleting the negative sign in the Alamouti sub-block) Similar to Alamouti matrices, Alamouti-like matrices have also some good properties that can be utilized to reduce the computational complexity, as shown in this paper In this paper, the efficient inverse Cholesky factorization proposed by us recently 8] is applied to Alamouti matrices in G-STBC and Alamouti-like matrices in OMP We exploit the properties of Alamouti or Alamouti-like matrices to reduce the complexity dramatically Then we develop efficient squareroot algorithms for G-STBC and OMP, respectively In the following sections, ( ) T, ( ) and ( ) H denote matrix transposition, matrix conjugate, and matrix conjugate transposition, respectively I M is the identity matrix with size M, and 0 M is the zero column vector with length M II SYSTEM MODELS In this section, we introduce the system model for G-STBC and that for OMP in the sub-nyquist sampling system 4] A System Model for G-STBC The considered G-STBC system 2] consists of 2M transmit antennas and N( M) receive antennas in a rich-scattering and flat-fading wireless channel The transmit side de-multiplexes a single data stream a,a 2,,a M,a M2 into M layers Each layer a m,a m2 (m =, 2,,M) is encoded by an independent Alamouti encoder, and then it is transmitted by its respective 2 transmit antennas in 2 time slots Let h n,i (n =, 2,,N and i =, 2,, 2M) denote the channel coefficient from transmit antenna i to receive antenna n Let x n and x n2 (n =, 2,,N) denote the signal received at antenna n in the st time slot and the 2 nd time slot, respectively Then the equivalent channel model 2] for G- STBC can be represented as x = Ha + w, () where w is the complex Gaussian noise vector with statistically independent entries, the equivalent receive signal vector x = ] x x 2 x N x T N2, (2) the equivalent transmit signal vector a = a a 2 a M a M2 ] T, (3) and the 2N 2M equivalent channel matrix h h 2 h,2m h,2m h 2 h h,2m h,2m H = h N h N2 h N,2M h N,2M (4) h N2 h N h N,2M h N,2M Assume E(aa H )=σai 2 2M and E(ww H )=σwi 2 2N The minimum mean-square error (MMSE) estimation of a is 2] where â = R H H x, (5) R = H H H + σ 2 w/σ 2 ai 2M (6) /2/$ IEEE
2 The recursive G-STBC algorithm 2] detects M Alamouti codes iteratively with the optimal ordering In each iteration, the undetected Alamouti code with the highest signal-to-noise ratio (SNR) is detected by (5), and then its effect is cancelled B System Model for OMP in Sub-Nyquist Sampling System Compressive sensing (CS) theory 9] asserts that a signal vector s C N with at most K nonzero components can be exactly recovered from a measurement vector y C M with M N The signal vector s is called K-sparse The measurement vector y = Ts, (7) where the full rank M N matrix T = ] t t 2 t N is commonly called the dictionary, and its i th column, ie t i, is commonly called an atom 6] The signal recovery problem is to find the sparse solution to the underdetermined linear system (7), ie, to represent y as a weighted sum of the least number of atoms The signal recovery algorithms include OMP 5] 7], a greedy algorithm that selects atoms from T iteratively to approximate y gradually In each iteration, OMP selects the atom best matching the current residual (ie the approximation error), renews the weights for all the already selected atoms (to obtain the least squares approximation of y), and then updates the residual accordingly In 4], the signal vector (ie X(f) in 4, left line 8 on page 387]) is conjugate symmetry, and so is the dictionary Thus the OMP algorithm is slightly modified to select a symmetric pair of atoms in each iteration 4, left line 6 on page 387] Accordingly we modify the OMP algorithm described in 6, left lines on page 2] to summarize the OMP algorithm in 4] We also follow these conventions in 6]: the set Γ n is a set containing n indices, and the matrix T Γ n is a sub-matrix of T containing only those columns of T with indices in Γ n Finally, we summarize the OMP algorithm in 4] as follows The OMP Algorithm in 4] Initialise Γ 0 =, the residual ỹ 0 = y 2 for k =, 2,,K/2 (2a) Projection: k = T H ỹ k (2b) i k = argmax N ( ) k i i= (2c) Select a symmetric pair of atoms by Γ = Γ i k ( N i k + ) (2d) Renew the weights for the selected atoms: s = ( T H Γ T ) Γ T H Γy (8) The dictionary T ] = S FD 4, equation (9)] F = F L0 F L0 4, right line 9 on page 380] is conjugate symmetry since F i = F i 4, equation (8)] S is real, since S ik = α ik 4, right line 22 on page 380] and α ik is real 4, the line below equation (4)] Moreover, D = diag ( ) d L0 d L0 4, right line 23 on page 380] has d i = d i 4, equation (6)] Thus we can verify that T = S FD = ] S F L0 d L0 S F L0 d L0 has S F i d i =( S F i d i ), ie, T = S FD is conjugate symmetry (2e) Update the residual by ỹ k = y T Γ s 3 Output: ỹ K/2, s 2 K/2 = s K III EFFICIENT INVERSE CHOLESKY FACTORIZATIONS FOR ALAMOUTI MATRICES AND ALAMOUTI-LIKE MATRICES We need to compute matrix inversions in (5) and (8) R in (5) is an Alamouti matrix 2] On the other hand, in (8) let R Γ = T H Γ T Γ (9) Obviously, R Γ =t i t i t i k t i k ]H t i t i t i k t i ]= k t H i t i t H i t i t H i t i k t H i t i k t T i t i t T i t i t T i t i k t T i t i k is an Alamouti-like matrix t H i t k i t H i t H i t k i k t H i t k i k t T i t k i t T i t T i t k i k t T i To reuse the already listed equations k for brevity, let us use (i) Γ to express equation (i) with the subscript k of each matrix replaced by Γ k and with overlines added to all vectors and scalars Then R (ie the submatrix in the first rows and columns of R) and R Γ can be represented as 2] R v v R = v H β r, (0) v H r, β and (0) Γ, respectively The inverse Cholesky factor 8] of R is F that satisfies F F H = R () We can simply apply equations () and (7) in 8] twice to compute F from F Then we have F = F u u 0 H λ f, 0 H 0 λ, (2) where we can employ equation (7) in 8] to compute λ =/ β v H F F H v, (3a) u = λ F F H v (3b) We can also employ equation (7) in 8] to compute the last column of F However, we will find some algorithms that are more efficient, to reduce the computational complexity A Inverse Cholesky Factorization for Alamouti Matrices R is an Alamouti matrix, and so are R 3, Lemma (Invariance Under Inversion)] Then the LDL T factors of R are Alamouti matrices 3, Lemma 3 (Invariance Under Triangular Factorization)], which include the unit lower triangular L and the diagonal D = diag ( ) d d 2 d Accordingly L D = L diag ( ) d d is also an Alamouti matrix, which is the unique Cholesky factor of R 0, Theorem 425 (Cholesky Factorization)], and is equivalent to F (since the permuted F is the Cholesky
3 factor of the correspondingly permuted R 8]) Thus F must be an Alamouti matrix Correspondingly in (2), f, =0, (4a) λ = λ, (4b) u = Alamouti(u ), (4c) where (4c) means that u is obtained from u directly since ] u u is an Alamouti matrix B Inverse Cholesky Factorization for Alamouti-like Matrices It is easy to verify that the sum, difference, or product of two Alamouti-like matrices is another Alamouti-like matrix Then we can prove that the inversion of an Alamouti-like matrix is another Alamouti-like matrix, as the proof for 3, Lemma (Invariance Under Inversion)] Now let us consider (2) Γ We employ (3b) Γ to obtain ] ] ū = λ f F Γ F H Γ (5), r, FΓ ū where F Γ = ] 0 H From (5) we deduce λ { ū = λ û λ ˆf, ū, (6a) f, = λ λ ˆf,, (6b) where { ˆf, = ū H + λ r,, (7a) û = F Γ F H Γ (7b) We can avoid computing (7b) and obtain û from F Γ F H Γ = ū ] / λ directly, since F Γ F H Γ = R ] Γ is an Alamouti-like matrix It can be seen from (2) Γ that the Alamoutilike sub-block in the last two rows ] and columns] of λ F Γ F H f, λ 0 Γ is 0 λ f, = λ λ2 + ] f, 2 f, λ f, λ, from which we deduce λ2 λ 2 = λ 2 + f, 2 (8) Finally we substitute (6b) into (8) to obtain λ =/ λ 2 2 ˆf, (9) Now we only need to compute (7a), (9) and (6) to obtain the last column of F IV EFFICIENT SQUARE-ROOT ALGORITHMS FOR G-STBC AND OMP IN SUB-NYQUIST SAMPLING SYSTEM A Efficient square-root algorithm for G-STBC We use F 2m (ie the equivalent Cholesky factor of R 2m ) instead of R 2m in the recursive G-STBC algorithm 2], to develop a square-root algorithm for G-STBC Similar to the square-root BLAST algorithm 8], in the initialization phase we permute the channel matrix H to be H 2M = ] h 2t h 2t h 2tM h 2tM (20) with t M,t M,,t to be the assumed successive detection order of M Alamouti codes, while in the iterative detection phase, we perform interference cancellation in and we block-upper triangularize F 2m by z 2M = H H 2M x, (2) F 2m Σ = F 2m 2 u 2m 2 u 2m 2 0 H 2m 2 f2m,2m f2m,2m (22) 0 H 2m 2 f 2m,2m f 2m,2m where Σ is unitary From (), (2) and (5) we obtain â 2m = F 2m F H 2mz 2m, ie, â 2m =(F 2m Σ)(F 2m Σ) H z 2m, into which we substitute (22) to estimate the last two entries of â 2m by ] ] â2m f = 2m,2m f2m,2m â 2m f 2m,2m f 2m,2m ] u H f 2m,2m f2m,2m (23) u H f 2m,2m f2m,2m z 2m The estimations â 2m and â 2m are quantized into a 2m and a 2m, respectively, whose effect is then cancelled in z 2m by 2] z 2m 2 = z 2] 2m v 2m 2 v 2m 2 ] a2m a 2m ] T, (24) where z 2] 2m is z 2m with the last two entries removed, and v 2m 2 and v 2m 2 are in R 2m, as shown in (0) Table I summarizes the proposed square-root G-STBC algorithm Notice that the transformation Σ in step (b) can be performed by the Givens rotations for Alamouti matrices that include inner rotations and outer rotations ] Generally, in the two minimum length rows (found in step (a)) that can be represented as rows and, the entries in the fist columns are all zeros Let Ω n i,j denote a complex Givens rotation matrix that rotates the i th entry and the j th entry in each row of F m and zeros the i th entry in the n th row 8] Then the transformation Σ in step (b) can be performed by Σ = m Ω 2j,2j m ( Ω 2j,2j+ Ω 2j,2j+2), (25) where Ω 2j,2j is an inner rotation that zeros two entries (ie f,2j and f,2j ) in an Alamouti sub-block and rotates the other two entries (ie f,2j and f,2j ) into real values, and Ω 2j,2j+ Ω 2j,2j+2 is an outer rotation that zeros the remaining two real values (ie f,2j and f,2j )ina previously inner-rotated Alamouti sub-block ] Since F 2m is still an Alamouti matrix after an inner or outer rotation, only half of the entries need to be computed ]
4 Initialization TABLE I THE PROPOSED SQUARE-ROOT ALGORITHM FOR G-STBC Permute H into H 2M by (20), and permute a 2M = a accordingly Compute R 2M and z 2M by (6) and (2), respectively Compute F 2 =(/ r, )I 2 Then compute F 2m from F 2m 2 iteratively by (3), (4) and (2) for m =2, 3,,M, to obtain F 2M Setm = M Iterative (a) Find the minimum length odd row and even row of F 2m, and permute them to the last odd row and Detection even row, respectively Correspondingly permute a 2m, z 2m, and rows and columns in R 2m (b) Block-upper triangularize F 2m by (22) to obtain F 2m 2 (c) Compute â 2m and â 2m by (23), which are quantized into a 2m and a 2m, respectively (d) Cancel the effect of a 2m and a 2m in z 2m by (24), to obtain z 2m 2 Remove a 2m and a 2m in a 2m to obtain a 2m 2 Also remove the last two rows and columns in R 2m to obtain R 2m 2 (e) If m>, go back to step (a) with F 2m 2, R 2m 2, z 2m 2 and a 2m 2 Letm = m B Efficient square-root algorithm for OMP in Sub-Nyquist Sampling System From (8), (9) and () Γ we can deduce In (26), let which satisfies y Γ = s = F Γ F H Γ TH Γy (26) y Γ = TH Γy, (27) y T Γ th γ y t H γ y ] T (28) where t γi is the i th column in T Γ Then we subsitute (27) and (2) into (26), to finally obtain s = ] s T T 0T 2 +f f ]f f ] H y Γ, (29) where f and f are the ( ) th and th columns of F Γ, respectively Moreover, from () Γ we deduce F Γ 2 = r, r,2 r, r, /(r, 2 r,2 2 ) 0 r, /(r, 2 r (30),2 2 ) To obtain the square-root algorithm for OMP, we only need to use step (2d ) instead of step (2d) in the OMP Algorithm in 4] (summarized in section IIB) Step (2d ) is as follow: if(k =): compute y Γ 2, F Γ 2 and s 2 by (27), (30) and (26) (where T H Γ y = y Γ2), respectively else: compute F Γ from F Γ by (3) Γ, (7), (9), (6) and (2) Γ (where (7b) in (7) does not need any computations, as above-mentioned) Compute y Γ and s by (28) and (29), respectively, from y Γ and s V COMPLEXITY EVALUATION In the proposed inverse Cholesky factorization for Alamouti matrices and that for Alamouti-like matrices, the complexities for the computations of the even columns (ie (4), (7a), (9), (6)) can be neglected, with respect to those for the computations of the odd columns (ie (3), (3) Γ ) Thus compared to the inverse Cholesky factorization 8] directly applied to Alamouti or Alamouti-like matrices, the proposed inverse Cholesky factorization can save about half the complexity TABLE II COMPLEXITY COMPARISON BETWEEN THE V-BLAST ALGORITHM IN 8] AND THE PROPOSED G-STBC ALGORITHM Step The V-BLAST Alg 8] The G-STBC Alg Compute R ( 2 M 2 N) (2M 2 N) Compute F ( 3 M 3 ) ( 4 3 M 3 ) Givens Rot ( 3 M 3, 9 M 3 )to(0) ( 8 9 M 3, 2 9 M 3 )to(0) Sum from ( 2 3 M M 2 N, from ( 20 9 M 3 +2M 2 N, 4 9 M M 2 4 N) 9 M 3 +2M 2 N) to ( 3 M M 2 N) to ( 4 3 M 3 +2M 2 N) Let (j, k) denote the complexity of j complex multiplications and k complex additions, and it is simplified to be (j) if j = k Table I compares the computational complexity of the proposed square-root algorithm for G-STBC and that of the square-root algorithm for Vertical Bell Laboratories Layered Space-Time architecture (V- BLAST) in 8] Notice that in the Hermitian Alamouti matrix R 2M and the upper triangular Alamouti matrix F 2m, only 4 entries need to be computed In (25), Ω 2j,2j and Ω 2j,2j+ Ω 2j,2j+2 only rotates non-zero entries in the first 2j and 2j+2 rows, respectively The complexities of a complex inner rotation and a real outer rotation are (4,2) and (2,0), ( respectively ( Thus the average complexity of ) (25) ) is M m m= k= ( 8M 3 9, 2M 3 9 m m (4, 2)(2j)+2 m (2, 0)(2j +2) ), where the probabilities for k =, 2,,m are assumed to be equal The complexity of (25) is zero when the assumed detection order is the same as the actual order 8] When the square-root V-BLAST algorithm 8] is directly applied to the equivalent channel model () that can be regarded as a 2M 2N V-BLAST system, the required complexity is from ( 6 3 M 3 +4M 2 N, 32 9 M 3 +4M 2 N) to ( 8 3 M 3 +4M 2 N) Thus when M = N, the average speedup in floating-point operations (flops) of the proposed square-root G-STBC algorithm over the square-root V-BLAST algorithm 8] directly applied to () ranges from 29 to 2 On the other hand, the complexity of the MMSE sub-optimal ordered successive interference cancelation (OSIC) detector for G-STBC using sorted QR decomposition (SQRD) is (8M 3 +4M 2 N) /2
5 FLOPS 7 x sqrt BLAST Alg (min mean value) sqrt BLAST Alg (max mean value) proposed sqrt G STBC Alg (min mean value) proposed sqrt G STBC Alg (max mean value) recursive G STBC Alg directly apply sqrt BLAST Alg to G STBC (min) directly apply sqrt BLAST Alg to G STBC (max) sub optimal OSIC SQRD G STBC Alg FLOPS Number of Receive Antennas OMP Alg by matrix inversion OMP Alg by QR decomposition 0 3 OMP Alg by matrix inverse update Apply inverse Cholesky factoriztion directly in OMP Alg proposed square root OMP Alg N: Dimension of the Original Signal Fig Comparison of computational complexity among the sub-optimal OSIC G-STBC algorithm in ], the OSIC G-STBC algorithm obtained by directly applying the square-root V-BLAST algorithm 8], the proposed square-root OSIC G-STBC algorithm, the recursive G-STBC algorithm 2], and the square-root V-BLAST algorithm 8], equation (9)] When M = N, the average speedup (in flops) of the proposed OSIC G-STBC detector over the suboptimal OSIC G-STBC detector ] ranges from 296 to 36 Assuming N = M, we carried out numerical experiments to count the average flops per time slot of the presented G- STBC or V-BLAST algorithms, for different N, the number of receive antennas The considered G-STBC system includes M Alamouti sub-streams and 2M transmitters, while the considered V-BLAST system includes M sub-streams and M transmitters All results are shown in Fig It can be seen that they are consistent with the theoretical flops calculation On the other hand, in the sub-nyquist sampling system 4], we simulate the flops of the proposed square-root OMP algorithm, the original OMP algorithm using direct matrix inversion 5], the original OMP algorithm with direct matrix inversion replaced by the inverse Cholesky factorization in 8], the fast OMP algorithm based on QR decomposition 5], 6], and the fast OMP algorithm by matrix inverse update 6], 7] Our simulations do not include the flops of the projection step (ie step (2a)) since they are the same for all the presented OMP algorithms As the simulations for Fig 5a in 2], we set K = N/2, M = K (log 2 N 2), and give flops for N = 2 i (i = 4, 5,, 0) The results are shown in Fig 2 Our simulations also show that the flops for the proposed square-root OMP algorithm are only about 30% of those for the fast OMP algorithm by matrix inverse update 6], 7] VI CONCLUSION When applying the efficient inverse Cholesky factorization 8] to Alamouti matrices in G-STBC and Alamouti-like matrices in OMP for the sub-nyquist sampling system 4], we utilize some good properties of Alamouti or Alamouti-like matrices to save about half the complexity Then we propose the whole square-root algorithms for G-STBC and OMP in Fig 2 Comparison of flops among the presented OMP algorithms 4], respectively The proposed square-root G-STBC algorithm has an average speedup of 2 29 with respect to the square-root V-BLAST algorithm 8] directly applied to the equivalent channel model for G-STBC, and has an average speedup of with respect to the sub-optimal OSIC G-STBC detector in ] On the other hand, when comparing the complexities of all the steps except the projection step, the complexity for the proposed square-root OMP algorithm is about 30% of that for the fast OMP algorithm by matrix inverse update 6], 7] REFERENCES ] J Eilert, D Wu, D Liu, D Wang, N Al-Dhahir and H Minn, Complexity reduction of matrix manipulation for multi-user STBC-MIMO decoding, IEEE Sarnoff Symposium, April 30 - May ] H Zhu, W Chen, B Li and F Gao, A Fast Recursive Algorithm for G-STBC, IEEE Trans on Communications, vol 59, no 8, 20 3] A H Sayed, W M, Younis and A Tarighat, An invariant matrix structure in multiantenna communications, Signal Processing Letters, IEEE, vol 2, Nov 2005, pp ] M Mishali and Y C Eldar, From Theory to Practice: Sub-Nyquist Sampling of Sparse Wideband Analog Signals, IEEE Journal of Selected Topics in Signal Processing, vol 4, no 2, pp 375C39, Apr 200 5] Joel Tropp and Anna Gilbert, Signal recovery from random measurements via orthogonal matching pursuit, IEEE Trans on Information Theory, vol 53, no 2, pp , December ] Thomas Blumensath and Mike E Davies, In Greedy Pursuit of New Directions: (Nearly) Orthogonal Matching Pursuit by Directional Optimisation, EUSIPCO 2007, Poznan, Poland, ] Y Fang, L Chen, J Wu and B Huang, GPU Implementation of Orthogonal Matching Pursuit for Compressive Sensing, ICPADS 20 8] H Zhu, W Chen, B Li and F Gao, An Improved Square-root Algorithm for V-BLAST Based on Efficient Inverse Cholesky Factorization, IEEE Trans on Wireless Communications, vol 0, no, pp 43-48, 20 9] D Donoho, Compressed sensing, IEEE Trans on Information Theory, vol 52, no 4, pp , April ] G H Golub and C F Van Loan, Matrix Computations, 3rd edition, Johns Hopkins University Press, 996 ] M Gomaa and A Ghrayeb, A Low Complexity MMSE Detector for Multiuser Layered Space-Time Coded MIMO Systems, CCECE 2007, Apr ] S Kunis and H Rauhut, Random sampling of sparse trigonometric polynomials II - Orthogonal matching pursuit versus basis pursuit, Preprint, 2006 Available online at
An Improved Square-Root Algorithm for V-BLAST Based on Efficient Inverse Cholesky Factorization
IEEE TRANSACTIONS ON WIRELESS COMMUNICATIONS, VOL. 0, NO., JANUARY 0 43 An Improved Square-Root Algorithm for V-BLAST Based on Efficient Inverse Cholesky Factorization Hufei Zhu, Member, IEEE, Wen Chen,
More informationLecture 5: Antenna Diversity and MIMO Capacity Theoretical Foundations of Wireless Communications 1. Overview. CommTh/EES/KTH
: Antenna Diversity and Theoretical Foundations of Wireless Communications Wednesday, May 4, 206 9:00-2:00, Conference Room SIP Textbook: D. Tse and P. Viswanath, Fundamentals of Wireless Communication
More informationNew Coherence and RIP Analysis for Weak. Orthogonal Matching Pursuit
New Coherence and RIP Analysis for Wea 1 Orthogonal Matching Pursuit Mingrui Yang, Member, IEEE, and Fran de Hoog arxiv:1405.3354v1 [cs.it] 14 May 2014 Abstract In this paper we define a new coherence
More informationAdaptive Space-Time Shift Keying Based Multiple-Input Multiple-Output Systems
ACSTSK Adaptive Space-Time Shift Keying Based Multiple-Input Multiple-Output Systems Professor Sheng Chen Electronics and Computer Science University of Southampton Southampton SO7 BJ, UK E-mail: sqc@ecs.soton.ac.uk
More informationAMS526: Numerical Analysis I (Numerical Linear Algebra for Computational and Data Sciences)
AMS526: Numerical Analysis I (Numerical Linear Algebra for Computational and Data Sciences) Lecture 1: Course Overview; Matrix Multiplication Xiangmin Jiao Stony Brook University Xiangmin Jiao Numerical
More informationSPARSE signal representations have gained popularity in recent
6958 IEEE TRANSACTIONS ON INFORMATION THEORY, VOL. 57, NO. 10, OCTOBER 2011 Blind Compressed Sensing Sivan Gleichman and Yonina C. Eldar, Senior Member, IEEE Abstract The fundamental principle underlying
More informationMaximum Achievable Diversity for MIMO-OFDM Systems with Arbitrary. Spatial Correlation
Maximum Achievable Diversity for MIMO-OFDM Systems with Arbitrary Spatial Correlation Ahmed K Sadek, Weifeng Su, and K J Ray Liu Department of Electrical and Computer Engineering, and Institute for Systems
More informationSimultaneous Sparsity
Simultaneous Sparsity Joel A. Tropp Anna C. Gilbert Martin J. Strauss {jtropp annacg martinjs}@umich.edu Department of Mathematics The University of Michigan 1 Simple Sparse Approximation Work in the d-dimensional,
More informationLecture 7 MIMO Communica2ons
Wireless Communications Lecture 7 MIMO Communica2ons Prof. Chun-Hung Liu Dept. of Electrical and Computer Engineering National Chiao Tung University Fall 2014 1 Outline MIMO Communications (Chapter 10
More informationPre-weighted Matching Pursuit Algorithms for Sparse Recovery
Journal of Information & Computational Science 11:9 (214) 2933 2939 June 1, 214 Available at http://www.joics.com Pre-weighted Matching Pursuit Algorithms for Sparse Recovery Jingfei He, Guiling Sun, Jie
More informationMatrix Factorization and Analysis
Chapter 7 Matrix Factorization and Analysis Matrix factorizations are an important part of the practice and analysis of signal processing. They are at the heart of many signal-processing algorithms. Their
More informationThe Sorted-QR Chase Detector for Multiple-Input Multiple-Output Channels
The Sorted-QR Chase Detector for Multiple-Input Multiple-Output Channels Deric W. Waters and John R. Barry School of ECE Georgia Institute of Technology Atlanta, GA 30332-0250 USA {deric, barry}@ece.gatech.edu
More informationStopping Condition for Greedy Block Sparse Signal Recovery
Stopping Condition for Greedy Block Sparse Signal Recovery Yu Luo, Ronggui Xie, Huarui Yin, and Weidong Wang Department of Electronics Engineering and Information Science, University of Science and Technology
More informationAn Efficient Optimal Algorithm for Integer-Forcing Linear MIMO Receivers Design
An Efficient Optimal Algorithm for Integer-Forcing Linear MIMO Receivers Design Jinming Wen, Lanping Li, Xiaohu Tang, Wai Ho Mow, and Chintha Tellambura Department of Electrical and Computer Engineering,
More informationof Orthogonal Matching Pursuit
A Sharp Restricted Isometry Constant Bound of Orthogonal Matching Pursuit Qun Mo arxiv:50.0708v [cs.it] 8 Jan 205 Abstract We shall show that if the restricted isometry constant (RIC) δ s+ (A) of the measurement
More informationLecture 8: MIMO Architectures (II) Theoretical Foundations of Wireless Communications 1. Overview. Ragnar Thobaben CommTh/EES/KTH
MIMO : MIMO Theoretical Foundations of Wireless Communications 1 Wednesday, May 25, 2016 09:15-12:00, SIP 1 Textbook: D. Tse and P. Viswanath, Fundamentals of Wireless Communication 1 / 20 Overview MIMO
More informationA New SLNR-based Linear Precoding for. Downlink Multi-User Multi-Stream MIMO Systems
A New SLNR-based Linear Precoding for 1 Downlin Multi-User Multi-Stream MIMO Systems arxiv:1008.0730v1 [cs.it] 4 Aug 2010 Peng Cheng, Meixia Tao and Wenjun Zhang Abstract Signal-to-leaage-and-noise ratio
More informationSparse analysis Lecture V: From Sparse Approximation to Sparse Signal Recovery
Sparse analysis Lecture V: From Sparse Approximation to Sparse Signal Recovery Anna C. Gilbert Department of Mathematics University of Michigan Connection between... Sparse Approximation and Compressed
More information5.6. PSEUDOINVERSES 101. A H w.
5.6. PSEUDOINVERSES 0 Corollary 5.6.4. If A is a matrix such that A H A is invertible, then the least-squares solution to Av = w is v = A H A ) A H w. The matrix A H A ) A H is the left inverse of A and
More informationRecovery Guarantees for Rank Aware Pursuits
BLANCHARD AND DAVIES: RECOVERY GUARANTEES FOR RANK AWARE PURSUITS 1 Recovery Guarantees for Rank Aware Pursuits Jeffrey D. Blanchard and Mike E. Davies Abstract This paper considers sufficient conditions
More information9. Numerical linear algebra background
Convex Optimization Boyd & Vandenberghe 9. Numerical linear algebra background matrix structure and algorithm complexity solving linear equations with factored matrices LU, Cholesky, LDL T factorization
More informationLecture 9: Diversity-Multiplexing Tradeoff Theoretical Foundations of Wireless Communications 1. Overview. Ragnar Thobaben CommTh/EES/KTH
: Diversity-Multiplexing Tradeoff Theoretical Foundations of Wireless Communications 1 Rayleigh Wednesday, June 1, 2016 09:15-12:00, SIP 1 Textbook: D. Tse and P. Viswanath, Fundamentals of Wireless Communication
More informationLecture 9: Diversity-Multiplexing Tradeoff Theoretical Foundations of Wireless Communications 1
: Diversity-Multiplexing Tradeoff Theoretical Foundations of Wireless Communications 1 Rayleigh Friday, May 25, 2018 09:00-11:30, Kansliet 1 Textbook: D. Tse and P. Viswanath, Fundamentals of Wireless
More informationMulti-Input Multi-Output Systems (MIMO) Channel Model for MIMO MIMO Decoding MIMO Gains Multi-User MIMO Systems
Multi-Input Multi-Output Systems (MIMO) Channel Model for MIMO MIMO Decoding MIMO Gains Multi-User MIMO Systems Multi-Input Multi-Output Systems (MIMO) Channel Model for MIMO MIMO Decoding MIMO Gains Multi-User
More informationIEEE C80216m-09/0079r1
Project IEEE 802.16 Broadband Wireless Access Working Group Title Efficient Demodulators for the DSTTD Scheme Date 2009-01-05 Submitted M. A. Khojastepour Ron Porat Source(s) NEC
More informationClosed Form Designs of Complex Orthogonal. Space-Time Block Codes of Rates (k + 1)=(2k) for
Closed Form Designs of Complex Orthogonal Space-Time Block Codes of Rates (k + 1)(k) for k 1 or k Transmit Antennas Kejie Lu, Shengli Fu, Xiang-Gen Xia Abstract In this correspondence, we present a systematic
More informationIEEE TRANSACTIONS ON INFORMATION THEORY, VOL. 49, NO. 10, OCTOBER
TRANSACTIONS ON INFORMATION THEORY, VOL 49, NO 10, OCTOBER 2003 1 Algebraic Properties of Space Time Block Codes in Intersymbol Interference Multiple-Access Channels Suhas N Diggavi, Member,, Naofal Al-Dhahir,
More informationOptimal Binary/Quaternary Adaptive Signature Design for Code-Division Multiplexing
Globecom 0 - Wireless Communications Symposium Optimal Binary/Quaternary Adaptive Signature Design for Code-Division Multiplexing ili Wei, Member, IEEE and Wen Chen, Senior Member, IEEE Department of Electronic
More informationLecture 9: Numerical Linear Algebra Primer (February 11st)
10-725/36-725: Convex Optimization Spring 2015 Lecture 9: Numerical Linear Algebra Primer (February 11st) Lecturer: Ryan Tibshirani Scribes: Avinash Siravuru, Guofan Wu, Maosheng Liu Note: LaTeX template
More informationAMS526: Numerical Analysis I (Numerical Linear Algebra)
AMS526: Numerical Analysis I (Numerical Linear Algebra) Lecture 1: Course Overview & Matrix-Vector Multiplication Xiangmin Jiao SUNY Stony Brook Xiangmin Jiao Numerical Analysis I 1 / 20 Outline 1 Course
More informationSimultaneous SDR Optimality via a Joint Matrix Decomp.
Simultaneous SDR Optimality via a Joint Matrix Decomposition Joint work with: Yuval Kochman, MIT Uri Erez, Tel Aviv Uni. May 26, 2011 Model: Source Multicasting over MIMO Channels z 1 H 1 y 1 Rx1 ŝ 1 s
More informationSparse analysis Lecture III: Dictionary geometry and greedy algorithms
Sparse analysis Lecture III: Dictionary geometry and greedy algorithms Anna C. Gilbert Department of Mathematics University of Michigan Intuition from ONB Key step in algorithm: r, ϕ j = x c i ϕ i, ϕ j
More informationMulti-Input Multi-Output Systems (MIMO) Channel Model for MIMO MIMO Decoding MIMO Gains Multi-User MIMO Systems
Multi-Input Multi-Output Systems (MIMO) Channel Model for MIMO MIMO Decoding MIMO Gains Multi-User MIMO Systems Multi-Input Multi-Output Systems (MIMO) Channel Model for MIMO MIMO Decoding MIMO Gains Multi-User
More informationAugmented Lattice Reduction for MIMO decoding
Augmented Lattice Reduction for MIMO decoding LAURA LUZZI joint work with G. Rekaya-Ben Othman and J.-C. Belfiore at Télécom-ParisTech NANYANG TECHNOLOGICAL UNIVERSITY SEPTEMBER 15, 2010 Laura Luzzi Augmented
More informationAlgebraic Multiuser Space Time Block Codes for 2 2 MIMO
Algebraic Multiuser Space Time Bloc Codes for 2 2 MIMO Yi Hong Institute of Advanced Telecom. University of Wales, Swansea, UK y.hong@swansea.ac.u Emanuele Viterbo DEIS - Università della Calabria via
More informationSGN Advanced Signal Processing Project bonus: Sparse model estimation
SGN 21006 Advanced Signal Processing Project bonus: Sparse model estimation Ioan Tabus Department of Signal Processing Tampere University of Technology Finland 1 / 12 Sparse models Initial problem: solve
More informationAMS526: Numerical Analysis I (Numerical Linear Algebra)
AMS526: Numerical Analysis I (Numerical Linear Algebra) Lecture 3: Positive-Definite Systems; Cholesky Factorization Xiangmin Jiao Stony Brook University Xiangmin Jiao Numerical Analysis I 1 / 11 Symmetric
More informationThe Effect of Ordered Detection and Antenna Selection on Diversity Gain of Decision Feedback Detector
The Effect of Ordered Detection and Antenna Selection on Diversity Gain of Decision Feedback Detector Yi Jiang Mahesh K. Varanasi Dept. Electrical & Computer Engineering, University of Colorado Boulder,
More informationExpected Error Based MMSE Detection Ordering for Iterative Detection-Decoding MIMO Systems
Expected Error Based MMSE Detection Ordering for Iterative Detection-Decoding MIMO Systems Lei Zhang, Chunhui Zhou, Shidong Zhou, Xibin Xu National Laboratory for Information Science and Technology, Tsinghua
More informationAn Adaptive Decision Feedback Equalizer for Time-Varying Frequency Selective MIMO Channels
An Adaptive Decision Feedback Equalizer for Time-Varying Frequency Selective MIMO Channels Athanasios A. Rontogiannis Institute of Space Applications and Remote Sensing National Observatory of Athens 5236,
More informationMATH 425-Spring 2010 HOMEWORK ASSIGNMENTS
MATH 425-Spring 2010 HOMEWORK ASSIGNMENTS Instructor: Shmuel Friedland Department of Mathematics, Statistics and Computer Science email: friedlan@uic.edu Last update April 18, 2010 1 HOMEWORK ASSIGNMENT
More informationSparse Solutions of Systems of Equations and Sparse Modelling of Signals and Images
Sparse Solutions of Systems of Equations and Sparse Modelling of Signals and Images Alfredo Nava-Tudela ant@umd.edu John J. Benedetto Department of Mathematics jjb@umd.edu Abstract In this project we are
More informationSpace-Time Coding for Multi-Antenna Systems
Space-Time Coding for Multi-Antenna Systems ECE 559VV Class Project Sreekanth Annapureddy vannapu2@uiuc.edu Dec 3rd 2007 MIMO: Diversity vs Multiplexing Multiplexing Diversity Pictures taken from lectures
More informationMS&E 318 (CME 338) Large-Scale Numerical Optimization
Stanford University, Management Science & Engineering (and ICME MS&E 38 (CME 338 Large-Scale Numerical Optimization Course description Instructor: Michael Saunders Spring 28 Notes : Review The course teaches
More informationGeneralized Orthogonal Matching Pursuit- A Review and Some
Generalized Orthogonal Matching Pursuit- A Review and Some New Results Department of Electronics and Electrical Communication Engineering Indian Institute of Technology, Kharagpur, INDIA Table of Contents
More informationMultipath Matching Pursuit
Multipath Matching Pursuit Submitted to IEEE trans. on Information theory Authors: S. Kwon, J. Wang, and B. Shim Presenter: Hwanchol Jang Multipath is investigated rather than a single path for a greedy
More informationRandomness-in-Structured Ensembles for Compressed Sensing of Images
Randomness-in-Structured Ensembles for Compressed Sensing of Images Abdolreza Abdolhosseini Moghadam Dep. of Electrical and Computer Engineering Michigan State University Email: abdolhos@msu.edu Hayder
More informationPhysical-Layer MIMO Relaying
Model Gaussian SISO MIMO Gauss.-BC General. Physical-Layer MIMO Relaying Anatoly Khina, Tel Aviv University Joint work with: Yuval Kochman, MIT Uri Erez, Tel Aviv University August 5, 2011 Model Gaussian
More information9. Numerical linear algebra background
Convex Optimization Boyd & Vandenberghe 9. Numerical linear algebra background matrix structure and algorithm complexity solving linear equations with factored matrices LU, Cholesky, LDL T factorization
More informationUnder sum power constraint, the capacity of MIMO channels
IEEE TRANSACTIONS ON COMMUNICATIONS, VOL 6, NO 9, SEPTEMBER 22 242 Iterative Mode-Dropping for the Sum Capacity of MIMO-MAC with Per-Antenna Power Constraint Yang Zhu and Mai Vu Abstract We propose an
More informationarxiv:cs/ v1 [cs.it] 11 Sep 2006
0 High Date-Rate Single-Symbol ML Decodable Distributed STBCs for Cooperative Networks arxiv:cs/0609054v1 [cs.it] 11 Sep 2006 Zhihang Yi and Il-Min Kim Department of Electrical and Computer Engineering
More informationA new method on deterministic construction of the measurement matrix in compressed sensing
A new method on deterministic construction of the measurement matrix in compressed sensing Qun Mo 1 arxiv:1503.01250v1 [cs.it] 4 Mar 2015 Abstract Construction on the measurement matrix A is a central
More informationOn the Role of the Properties of the Nonzero Entries on Sparse Signal Recovery
On the Role of the Properties of the Nonzero Entries on Sparse Signal Recovery Yuzhe Jin and Bhaskar D. Rao Department of Electrical and Computer Engineering, University of California at San Diego, La
More informationTHE IC-BASED DETECTION ALGORITHM IN THE UPLINK LARGE-SCALE MIMO SYSTEM. Received November 2016; revised March 2017
International Journal of Innovative Computing, Information and Control ICIC International c 017 ISSN 1349-4198 Volume 13, Number 4, August 017 pp. 1399 1406 THE IC-BASED DETECTION ALGORITHM IN THE UPLINK
More informationInteger Least Squares: Sphere Decoding and the LLL Algorithm
Integer Least Squares: Sphere Decoding and the LLL Algorithm Sanzheng Qiao Department of Computing and Software McMaster University 28 Main St. West Hamilton Ontario L8S 4L7 Canada. ABSTRACT This paper
More informationFast-Decodable MIMO HARQ Systems
1 Fast-Decodable MIMO HARQ Systems Seyyed Saleh Hosseini, Student Member, IEEE, Jamshid Abouei, Senior Member, IEEE, and Murat Uysal, Senior Member, IEEE Abstract This paper presents a comprehensive study
More informationConstrained Detection for Multiple-Input Multiple-Output Channels
Constrained Detection for Multiple-Input Multiple-Output Channels Tao Cui, Chintha Tellambura and Yue Wu Department of Electrical and Computer Engineering University of Alberta Edmonton, AB, Canada T6G
More informationDense LU factorization and its error analysis
Dense LU factorization and its error analysis Laura Grigori INRIA and LJLL, UPMC February 2016 Plan Basis of floating point arithmetic and stability analysis Notation, results, proofs taken from [N.J.Higham,
More informationInteger-Forcing Linear Receiver Design over MIMO Channels
Globecom 0 - Signal Processing for Communications Symposium Integer-Forcing inear Receiver Design over MIMO Channels ili Wei, Member, IEEE and Wen Chen, Senior Member, IEEE Department of Electronic Engineering,
More informationModel-Based Compressive Sensing for Signal Ensembles. Marco F. Duarte Volkan Cevher Richard G. Baraniuk
Model-Based Compressive Sensing for Signal Ensembles Marco F. Duarte Volkan Cevher Richard G. Baraniuk Concise Signal Structure Sparse signal: only K out of N coordinates nonzero model: union of K-dimensional
More informationBlind MIMO communication based on Subspace Estimation
Blind MIMO communication based on Subspace Estimation T. Dahl, S. Silva, N. Christophersen, D. Gesbert T. Dahl, S. Silva, and N. Christophersen are at the Department of Informatics, University of Oslo,
More informationAnalysis of Detection Methods in Massive MIMO Systems
Journal of Physics: Conference Series PAPER OPEN ACCESS Analysis of Detection Methods in Massive MIMO Systems To cite this article: Siyuan Li 208 J. Phys.: Conf. Ser. 087 042002 View the article online
More informationApplications of Lattices in Telecommunications
Applications of Lattices in Telecommunications Dept of Electrical and Computer Systems Engineering Monash University amin.sakzad@monash.edu Oct. 2013 1 Sphere Decoder Algorithm Rotated Signal Constellations
More informationNumerical Linear Algebra Primer. Ryan Tibshirani Convex Optimization /36-725
Numerical Linear Algebra Primer Ryan Tibshirani Convex Optimization 10-725/36-725 Last time: proximal gradient descent Consider the problem min g(x) + h(x) with g, h convex, g differentiable, and h simple
More information14.2 QR Factorization with Column Pivoting
page 531 Chapter 14 Special Topics Background Material Needed Vector and Matrix Norms (Section 25) Rounding Errors in Basic Floating Point Operations (Section 33 37) Forward Elimination and Back Substitution
More informationGreedy Signal Recovery and Uniform Uncertainty Principles
Greedy Signal Recovery and Uniform Uncertainty Principles SPIE - IE 2008 Deanna Needell Joint work with Roman Vershynin UC Davis, January 2008 Greedy Signal Recovery and Uniform Uncertainty Principles
More informationPolynomial matrix QR decomposition and iterative decoding of frequency selective MIMO channels
Loughborough University Institutional Repository Polynomial matrix QR decomposition and iterative decoding of frequency selective MIMO channels This item was submitted to Loughborough University's Institutional
More informationThe RF-Chain Limited MIMO System: Part II Case Study of V-BLAST and GMD
The RF-Chain Limited MIMO System: Part II Case Study of V-BLAST and GMD Yi Jiang Mahesh K. Varanasi Abstract In Part I of this paper, we have established the fundamental D-M tradeoff of a RF-chain limited
More informationNumerical Linear Algebra
Numerical Linear Algebra Decompositions, numerical aspects Gerard Sleijpen and Martin van Gijzen September 27, 2017 1 Delft University of Technology Program Lecture 2 LU-decomposition Basic algorithm Cost
More informationProgram Lecture 2. Numerical Linear Algebra. Gaussian elimination (2) Gaussian elimination. Decompositions, numerical aspects
Numerical Linear Algebra Decompositions, numerical aspects Program Lecture 2 LU-decomposition Basic algorithm Cost Stability Pivoting Cholesky decomposition Sparse matrices and reorderings Gerard Sleijpen
More informationA Computationally Efficient Block Transmission Scheme Based on Approximated Cholesky Factors
A Computationally Efficient Block Transmission Scheme Based on Approximated Cholesky Factors C. Vincent Sinn Telecommunications Laboratory University of Sydney, Australia cvsinn@ee.usyd.edu.au Daniel Bielefeld
More informationExploiting Sparsity for Wireless Communications
Exploiting Sparsity for Wireless Communications Georgios B. Giannakis Dept. of ECE, Univ. of Minnesota http://spincom.ece.umn.edu Acknowledgements: D. Angelosante, J.-A. Bazerque, H. Zhu; and NSF grants
More informationK User Interference Channel with Backhaul
1 K User Interference Channel with Backhaul Cooperation: DoF vs. Backhaul Load Trade Off Borna Kananian,, Mohammad A. Maddah-Ali,, Babak H. Khalaj, Department of Electrical Engineering, Sharif University
More informationMatrices. Chapter What is a Matrix? We review the basic matrix operations. An array of numbers a a 1n A = a m1...
Chapter Matrices We review the basic matrix operations What is a Matrix? An array of numbers a a n A = a m a mn with m rows and n columns is a m n matrix Element a ij in located in position (i, j The elements
More informationA Thesis for the Degree of Master. Efficient Detection Algorithms for MMSE-SIC in MIMO Systems
A Thesis for the Degree of Master Efficient Detection Algorithms for MMSE-SIC in MIMO Systems Junyoung Nam School of Engineering Information and Communications University 2008 Efficient Detection Algorithms
More informationSecure Degrees of Freedom of the MIMO Multiple Access Wiretap Channel
Secure Degrees of Freedom of the MIMO Multiple Access Wiretap Channel Pritam Mukherjee Sennur Ulukus Department of Electrical and Computer Engineering University of Maryland, College Park, MD 074 pritamm@umd.edu
More informationEM Channel Estimation and Data Detection for MIMO-CDMA Systems over Slow-Fading Channels
EM Channel Estimation and Data Detection for MIMO-CDMA Systems over Slow-Fading Channels Ayman Assra 1, Walaa Hamouda 1, and Amr Youssef 1 Department of Electrical and Computer Engineering Concordia Institute
More informationThis can be accomplished by left matrix multiplication as follows: I
1 Numerical Linear Algebra 11 The LU Factorization Recall from linear algebra that Gaussian elimination is a method for solving linear systems of the form Ax = b, where A R m n and bran(a) In this method
More informationBlockMatrixComputations and the Singular Value Decomposition. ATaleofTwoIdeas
BlockMatrixComputations and the Singular Value Decomposition ATaleofTwoIdeas Charles F. Van Loan Department of Computer Science Cornell University Supported in part by the NSF contract CCR-9901988. Block
More informationG1110 & 852G1 Numerical Linear Algebra
The University of Sussex Department of Mathematics G & 85G Numerical Linear Algebra Lecture Notes Autumn Term Kerstin Hesse (w aw S w a w w (w aw H(wa = (w aw + w Figure : Geometric explanation of the
More informationCHANNEL FEEDBACK QUANTIZATION METHODS FOR MISO AND MIMO SYSTEMS
CHANNEL FEEDBACK QUANTIZATION METHODS FOR MISO AND MIMO SYSTEMS June Chul Roh and Bhaskar D Rao Department of Electrical and Computer Engineering University of California, San Diego La Jolla, CA 9293 47,
More informationFeasibility Conditions for Interference Alignment
Feasibility Conditions for Interference Alignment Cenk M. Yetis Istanbul Technical University Informatics Inst. Maslak, Istanbul, TURKEY Email: cenkmyetis@yahoo.com Tiangao Gou, Syed A. Jafar University
More informationInteractive Interference Alignment
Interactive Interference Alignment Quan Geng, Sreeram annan, and Pramod Viswanath Coordinated Science Laboratory and Dept. of ECE University of Illinois, Urbana-Champaign, IL 61801 Email: {geng5, kannan1,
More information(Article begins on next page)
Chalmers Publication Library Copyright Notice 2010 IEEE. Personal use of this material is permitted. However, permission to reprint/republish this material for advertising or promotional purposes or for
More informationOn the Optimal Phase Control in MIMO Systems with Phase Quantization
On the Optimal hase Control in MIMO Systems with hase Quantization Majid Khabbazian, Kin-Kwong Leung and Mohammad A. Safari Department of Electrical and Computer Engineering University of British Columbia,
More informationInformation Theory for Wireless Communications, Part II:
Information Theory for Wireless Communications, Part II: Lecture 5: Multiuser Gaussian MIMO Multiple-Access Channel Instructor: Dr Saif K Mohammed Scribe: Johannes Lindblom In this lecture, we give the
More informationOrthogonal Transformations
Orthogonal Transformations Tom Lyche University of Oslo Norway Orthogonal Transformations p. 1/3 Applications of Qx with Q T Q = I 1. solving least squares problems (today) 2. solving linear equations
More informationLow-Complexity Encoding Algorithm for LDPC Codes
EECE 580B Modern Coding Theory Low-Complexity Encoding Algorithm for LDPC Codes Problem: Given the following matrix (imagine a larger matrix with a small number of ones) and the vector of information bits,
More informationAMS526: Numerical Analysis I (Numerical Linear Algebra)
AMS526: Numerical Analysis I (Numerical Linear Algebra) Lecture 7: More on Householder Reflectors; Least Squares Problems Xiangmin Jiao SUNY Stony Brook Xiangmin Jiao Numerical Analysis I 1 / 15 Outline
More informationDesigning Information Devices and Systems I Discussion 13B
EECS 6A Fall 7 Designing Information Devices and Systems I Discussion 3B. Orthogonal Matching Pursuit Lecture Orthogonal Matching Pursuit (OMP) algorithm: Inputs: A set of m songs, each of length n: S
More informationPartial LLL Reduction
Partial Reduction Xiaohu Xie School of Computer Science McGill University Montreal, Quebec, Canada H3A A7 Email: xiaohu.xie@mail.mcgill.ca Xiao-Wen Chang School of Computer Science McGill University Montreal,
More informationIntroduction to Compressed Sensing
Introduction to Compressed Sensing Alejandro Parada, Gonzalo Arce University of Delaware August 25, 2016 Motivation: Classical Sampling 1 Motivation: Classical Sampling Issues Some applications Radar Spectral
More informationML Detection with Blind Linear Prediction for Differential Space-Time Block Code Systems
ML Detection with Blind Prediction for Differential SpaceTime Block Code Systems Seree Wanichpakdeedecha, Kazuhiko Fukawa, Hiroshi Suzuki, Satoshi Suyama Tokyo Institute of Technology 11, Ookayama, Meguroku,
More informationMulti-Branch MMSE Decision Feedback Detection Algorithms. with Error Propagation Mitigation for MIMO Systems
Multi-Branch MMSE Decision Feedback Detection Algorithms with Error Propagation Mitigation for MIMO Systems Rodrigo C. de Lamare Communications Research Group, University of York, UK in collaboration with
More informationCME 302: NUMERICAL LINEAR ALGEBRA FALL 2005/06 LECTURE 6
CME 302: NUMERICAL LINEAR ALGEBRA FALL 2005/06 LECTURE 6 GENE H GOLUB Issues with Floating-point Arithmetic We conclude our discussion of floating-point arithmetic by highlighting two issues that frequently
More informationAn Introduction to Sparse Approximation
An Introduction to Sparse Approximation Anna C. Gilbert Department of Mathematics University of Michigan Basic image/signal/data compression: transform coding Approximate signals sparsely Compress images,
More informationSolution of Linear Equations
Solution of Linear Equations (Com S 477/577 Notes) Yan-Bin Jia Sep 7, 07 We have discussed general methods for solving arbitrary equations, and looked at the special class of polynomial equations A subclass
More informationLAPACK-Style Codes for Pivoted Cholesky and QR Updating
LAPACK-Style Codes for Pivoted Cholesky and QR Updating Sven Hammarling 1, Nicholas J. Higham 2, and Craig Lucas 3 1 NAG Ltd.,Wilkinson House, Jordan Hill Road, Oxford, OX2 8DR, England, sven@nag.co.uk,
More informationSingle-letter Characterization of Signal Estimation from Linear Measurements
Single-letter Characterization of Signal Estimation from Linear Measurements Dongning Guo Dror Baron Shlomo Shamai The work has been supported by the European Commission in the framework of the FP7 Network
More information