Hi, I wish to solve the system where is block band Toeplitz and the unkown. In my case, has a particular form: the blocks are symmetric and quite large (. or even more) and I have around blocks on the block diagonal of , and around block lower (and upper) diagonals ( is band block), such that, for example for and : where the bar means the conjugate. The LU decomposition works quite well, but I would like to exploit this particular structure: symmetric blocks, block Toeplitz and block band matrix. However, the size of the blocks is large, and algorithms like block Levinson in the FEx are not efficient at all since they fully invert a block times. They are very efficient for small blocks. Any idea?

Solve particular block Toeplitz linear system

조회 수: 11 (최근 30일)

이전 댓글 표시

david gasperini 2021년 1월 4일

1
링크

이 질문에 대한 바로 가기 링크

https://kr.mathworks.com/matlabcentral/answers/708388-solve-particular-block-toeplitz-linear-system

댓글: Christine Tobler 2021년 1월 8일

Hi,

I wish to solve the system

where T is block band Toeplitz and X the unkown. In my case, T has a particular form: the blocks are symmetric and quite large (

or even more) and I have around n blocks on the block diagonal of T, and around m block lower (and upper) diagonals ( T is band block), such that, for example for

and

0 0

where the bar means the conjugate.

0 0

The LU decomposition works quite well, but I would like to exploit this particular structure: symmetric blocks, block Toeplitz and block band matrix.

However, the size of the blocks is large, and algorithms like block Levinson in the FEx are not efficient at all since they fully invert a block n times. They are very efficient for small blocks.

Any idea?