삼지각 행렬

Tridiagonal matrix

선형대수학에서 삼지각행렬 대각선상에 0이 아닌 원소를 갖는 밴드 행렬로, 주 대각선 위에만 첫 번째 대각선, 그리고 주 대각선 위에 첫 번째 대각선이 있다.

예를 들어, 다음 행렬은 삼지각형이다.

삼지각 행렬의 결정 인자는 원소의 연속성에 의해 주어진다.[1]

대칭(또는 에르미트어) 행렬을 3지각형 형태로 직교 변환하는 것은 Lanczos 알고리즘을 사용하여 수행할 수 있다.

특성.

3각형 행렬은 헤센베르크 행렬의 상단과 하단에 있는 행렬이다.[2] 특히 3각 행렬은 p + q/2 = n — 3각 행렬의 치수처럼 p 1 by-1과 q 2 by-2 행렬의 직접 이다. 일반적인 3각형 행렬이 반드시 대칭이거나 에르미트어인 것은 아니지만, 선형 대수 문제를 풀 때 발생하는 많은 행렬은 이러한 속성 중 하나를 가지고 있다. 또한 실제 3각형 행렬 A가 모든 k대해k,k+1k+1,k a > 0을 만족하여 입력의 기호가 대칭인 경우, 기준 행렬의 대각선 변화에 의해 은둔형 행렬과 유사하다. 그러므로, 그것의 고유값은 진짜다. 만약 우리가 엄격한 불평등을 ≥ 0으로k,k+1k+1,k 대체한다면, 그렇다면 연속성에 의해, 고유값은 여전히 실제라고 보장되지만, 행렬은 더 이상 은둔자 행렬과 유사할 필요가 없다.[3]

모든 n × n 삼지각 행렬의 집합3n-2차원 벡터 공간을 형성한다.

많은 선형 대수 알고리즘은 대각 행렬에 적용할 때 계산 노력이 현저히 덜 필요하며, 이러한 개선은 3각 행렬에도 종종 이월된다.

결정인자

순서 n의 삼지각 행렬 A의 결정 인자는 3개월의 재발 관계에서 계산할 수 있다.[4] 쓰기1 f1 = a = a(1, f1 a로만1 구성된 1 X 1 행렬의 결정인자임)로 하고, let let

시퀀스(fi)를 연속체라고 하며, 재발 관계를 만족한다.

초기 값 f0 = 1 및 f = 0−1. 이 공식을 사용하여 3각형 행렬의 결정 인자를 계산하는 비용은 n 단위의 선형인 반면, 일반 행렬의 비용은 세제곱이다.

반전

비성삼각형 행렬역행렬 T

에 의해 주어지다

θi 재발 관계를 만족하는 경우.

초기 조건 θ0 = 1, θ1 = a 1 ϕi 만족

초기 조건 ϕn+1 = 1 및 ϕn = an.[5][6]

닫힌 폼 솔루션은 모든 대각선 및 비대각 원소[8] 동일한[7] 대칭 행렬과 같은 특수한 경우와 일반 사례에 대해서도 계산할 수 있다.[9][10]

일반적으로 삼지각 행렬의 역행렬은 반음행렬이고 그 반대도 반음행렬이다.[11]

선형계 해법

등식 Ax = b for 의 시스템은 A삼지각 행렬 알고리즘이라 불리는 삼지각일 때 가우스 제거의 효율적인 형태로 해결할 수 있어 O(n) 연산이 필요하다.[12]

아이겐값

3각형 행렬이 또한 토플리츠일 때, 그것의 고유값을 위한 단순한 폐쇄형 솔루션이 있다. 즉,[13][14] 다음과 같다.

실제 대칭 3각형 행렬은 실제 고유값을 가지며, 모든 비대칭 원소가 0이 아닐 경우 모든 고유값은 구별(단순)된다.[15] 실제 대칭 3지각 행렬의 고유값을 임의의 유한 정밀도로 수치 연산하기 위해 수많은 방법이 존재하며, 일반적으로 크기 })의 행렬에 O ) N(\displaysty n 연산이 필요하다 ) n만 필요[16]

곁눈질로, 축소되지 않은 대칭 3각 행렬은 3각형의 비대칭 원소를 포함하는 행렬로, 고유값은 구별되는 반면 고유 벡터는 스케일 팩터까지 고유하며 상호 직교한다.[17]

비대칭 3각형 행렬의 경우 유사성 변환을 사용하여 Eigendecomposition을 계산할 수 있다.

대칭 3각형 행렬과의 유사성

실제 3지각, 비대칭 행렬이 지정됨

여기서

대각선 이탈 입력의 각 제품이 엄격하게 의 b i i> {\ b_이라고 가정하고 변환 매트릭스 (를) 정의하십시오.

유사성 변환 - D은(는) 다음을 기준으로 대칭[18] 3각형 J {\ J을(를) 산출함

( J {\의 고유값이 동일하다는 점에 유의하십시오.

컴퓨터 프로그래밍

일반 매트릭스를 헤센베르크 양식으로 줄이는 변환은 에르미트 매트릭스를 3각형 형태로 감소시킬 것이다. 그러므로 많은 고유값 알고리즘이, 에르미타니아 행렬에 적용되었을 때, 입력된 에르미타니아 행렬을 첫 번째 단계로 (대칭 실제) 삼지각형 형태로 감소시킨다.

3각형 매트릭스는 또한 특별한 저장 방식을 사용하여 일반 매트릭스보다 더 효율적으로 저장할 수 있다. 예를 들어, LAPACK Fortran 패키지는 대각선 요소를 포함하는 길이 n의 하나, 그리고 대각선초대각선을 포함하는 길이 n - 1의 두 개의 n - 1의 비대칭 3차원 배열로 저장된다.

참고 항목

메모들

  1. ^ Thomas Muir (1960). A treatise on the theory of determinants. Dover Publications. pp. 516–525.
  2. ^ Horn, Roger A.; Johnson, Charles R. (1985). Matrix Analysis. Cambridge University Press. p. 28. ISBN 0521386322.
  3. ^ 혼 앤 존슨, 174페이지
  4. ^ El-Mikkawy, M. E. A. (2004). "On the inverse of a general tridiagonal matrix". Applied Mathematics and Computation. 150 (3): 669–679. doi:10.1016/S0096-3003(03)00298-4.
  5. ^ Da Fonseca, C. M. (2007). "On the eigenvalues of some tridiagonal matrices". Journal of Computational and Applied Mathematics. 200: 283–286. doi:10.1016/j.cam.2005.08.047.
  6. ^ Usmani, R. A. (1994). "Inversion of a tridiagonal jacobi matrix". Linear Algebra and its Applications. 212–213: 413–414. doi:10.1016/0024-3795(94)90414-6.
  7. ^ Hu, G. Y.; O'Connell, R. F. (1996). "Analytical inversion of symmetric tridiagonal matrices". Journal of Physics A: Mathematical and General. 29 (7): 1511. doi:10.1088/0305-4470/29/7/020.
  8. ^ Huang, Y.; McColl, W. F. (1997). "Analytical inversion of general tridiagonal matrices". Journal of Physics A: Mathematical and General. 30 (22): 7919. doi:10.1088/0305-4470/30/22/026.
  9. ^ Mallik, R. K. (2001). "The inverse of a tridiagonal matrix". Linear Algebra and its Applications. 325: 109–139. doi:10.1016/S0024-3795(00)00262-7.
  10. ^ Kılıç, E. (2008). "Explicit formula for the inverse of a tridiagonal matrix by backward continued fractions". Applied Mathematics and Computation. 197: 345–357. doi:10.1016/j.amc.2007.07.046.
  11. ^ Raf Vandebril; Marc Van Barel; Nicola Mastronardi (2008). Matrix Computations and Semiseparable Matrices. Volume I: Linear Systems. JHU Press. Theorem 1.38, p. 41. ISBN 978-0-8018-8714-7.
  12. ^ Golub, Gene H.; Van Loan, Charles F. (1996). Matrix Computations (3rd ed.). The Johns Hopkins University Press. ISBN 0-8018-5414-8.
  13. ^ Noschese, S.; Pasquini, L.; Reichel, L. (2013). "Tridiagonal Toeplitz matrices: Properties and novel applications". Numerical Linear Algebra with Applications. 20 (2): 302. doi:10.1002/nla.1811.
  14. ^ 왜냐하면 cos⁡()))− 못 말리겠고 ⁡(π −)){\displaystyle \cos())=-\cos(\pi -x)}, Kulkarni, D;슈미트, D;추이, S.K.(1999년):에 수행됩니까 이것은 또한 한+2bc. 왜냐하면 ⁡(kπ/(n+1)){\displaystyle a+2{\sqrt{bc}}\cos(k\pi/{())})}같이 쓸 수 있다."3중 대각 pseudo-Toeplitz 매트릭스의 고유"(PDF).선형 대수학과 그 응용.297:63.doi:10.1016(99)00114-7.
  15. ^ Parlett, B.N. (1980). The Symmetric Eigenvalue Problem. Prentice Hall, Inc.
  16. ^ Coakley, E.S.; Rokhlin, V. (2012). "A fast divide-and-conquer algorithm for computing the spectra of real symmetric tridiagonal matrices". Applied and Computational Harmonic Analysis. 34 (3): 379–414. doi:10.1016/j.acha.2012.06.003.
  17. ^ Dhillon, Inderjit Singh. A New O(n 2 ) Algorithm for the Symmetric Tridiagonal Eigenvalue/Eigenvector Problem (PDF). p. 8.
  18. ^ "www.math.hkbu.edu.hk math lecture" (PDF).

외부 링크