삼지각 행렬
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, f는1 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차원 배열로 저장된다.
참고 항목
메모들
- ^ Thomas Muir (1960). A treatise on the theory of determinants. Dover Publications. pp. 516–525.
- ^ Horn, Roger A.; Johnson, Charles R. (1985). Matrix Analysis. Cambridge University Press. p. 28. ISBN 0521386322.
- ^ 혼 앤 존슨, 174페이지
- ^ 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.
- ^ 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.
- ^ 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.
- ^ 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.
- ^ 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.
- ^ 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.
- ^ 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.
- ^ 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.
- ^ Golub, Gene H.; Van Loan, Charles F. (1996). Matrix Computations (3rd ed.). The Johns Hopkins University Press. ISBN 0-8018-5414-8.
- ^ 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.
- ^ 왜냐하면 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.
- ^ Parlett, B.N. (1980). The Symmetric Eigenvalue Problem. Prentice Hall, Inc.
- ^ 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.
- ^ Dhillon, Inderjit Singh. A New O(n 2 ) Algorithm for the Symmetric Tridiagonal Eigenvalue/Eigenvector Problem (PDF). p. 8.
- ^ "www.math.hkbu.edu.hk math lecture" (PDF).
외부 링크
- LAPACK 설명서의 트리디각 및 비디각 행렬.
- Moawwad El-Mikkawy, Abdelrahman Karawia (2006). "Inversion of general tridiagonal matrices" (PDF). Applied Mathematics Letters. 19 (8): 712–720. doi:10.1016/j.aml.2005.11.012. Archived from the original (PDF) on 2011-07-20.
- 응축형(Hessenberg, 3지각형, 비다이각형) 형태로의 축소를 위한 고성능 알고리즘
- 3각 선형 시스템 용해제(C++)
