MAT1041 · 14. hafta · - Symmetric Matrices and Orthogonal Diagonalization

Simetrik matrisler ve ortogonal köşegenleştirme

Ortogonal matrisleri (Pᵀ = P⁻¹), simetrik matrislerin özel özelliklerini (reel özdeğerler, dik özuzaylar, her zaman köşegenleşme) ve bir simetrik matrisi PᵀAP = D olacak şekilde ortogonal bir P ile köşegenleştirmeyi, tekrarlı özdeğerde Gram-Schmidt kullanmayı öğrenirsin.

symmetric matrixorthogonal matrixorthonormal setorthogonal vectorsnormalizeorthogonally diagonalizablespectral theoremGram-Schmidt processquadratic formprincipal axes

Dönemin son konusu, simetrik matrislerin ne kadar “iyi huylu” olduğunu gösteriyor. Genel bir matris köşegenleşmeyebilir, reel özdeğeri olmayabilir, özvektörleri birbirine eğik durabilir. Simetrik matriste bunların hiçbiri olmaz: özdeğerler reel, farklı özdeğerlerin özvektörleri dik, ve matris her zaman köşegenleşir. Üstelik PP öyle seçilebilir ki sütunları birbirine dik birim vektörler olur; o zaman P−1P^{-1} hesaplamaya bile gerek kalmaz, P−1=PTP^{-1} = P^T. Bu sonuca spektral teorem denir ve istatistikten (PCA) fiziğe kadar her yerde kullanılır.

Altın kurallar:

  1. Ortogonal matris: sütunlar birim uzunlukta ve ikişer ikişer dik. Tersi transpozu.
  2. Simetrik matriste farklı özdeğerlerin özvektörleri otomatik dik.
  3. Tekrarlı özdeğerin özuzayında Gram-Schmidt uygula.
  4. En sonda her vektörü normalize et. Unutulan adım bu.

Ortogonal vektörler ve ortogonal matrisler

  1. haftadan: u⋅v=0\mathbf{u} \cdot \mathbf{v} = 0 ise u\mathbf{u} ve v\mathbf{v} diktir (orthogonal). Bir küme ikişer ikişer dik vektörlerden oluşuyorsa ortogonal küme, ayrıca hepsi birim uzunluktaysa ortonormal kümedir. Sıfır içermeyen ortogonal küme her zaman lineer bağımsızdır.

Kare QQ matrisi QTQ=IQ^TQ = I sağlıyorsa ortogonal matristir (orthogonal matrix). QTQQ^TQ‘nun (i,j)(i, j) elemanı ii. sütun ile jj. sütunun iç çarpımı olduğu için:

QTQ=I  ⟺  su¨tunlar ortonormal  ⟺  Q−1=QT.Q^TQ = I \iff \text{sütunlar ortonormal} \iff Q^{-1} = Q^T.

ÖrnekOrtogonal mi?
[cos⁡θ−sin⁡θsin⁡θcos⁡θ]\begin{bmatrix} \cos\theta & -\sin\theta \\ \sin\theta & \cos\theta \end{bmatrix} (dönme)evet, det⁡=1\det = 1
[100−1]\begin{bmatrix} 1 & 0 \\ 0 & -1 \end{bmatrix} (yansıma)evet, det⁡=−1\det = -1
[3/5−4/54/53/5]\begin{bmatrix} 3/5 & -4/5 \\ 4/5 & 3/5 \end{bmatrix}evet (Q5)
[111−1]\begin{bmatrix} 1 & 1 \\ 1 & -1 \end{bmatrix}hayır: sütunlar dik ama uzunluk 2\sqrt2

Son satır klasik tuzak: dik sütunlar yetmez, birim de olmalı. 12[111−1]\frac{1}{\sqrt2}\begin{bmatrix} 1 & 1 \\ 1 & -1 \end{bmatrix} ortogonaldir.

Ortogonal matrisler uzunluğu ve açıyı korur: ∥Qx∥=∥x∥\|Q\mathbf{x}\| = \|\mathbf{x}\|, çünkü ∥Qx∥2=xTQTQx=xTx\|Q\mathbf{x}\|^2 = \mathbf{x}^TQ^TQ\mathbf{x} = \mathbf{x}^T\mathbf{x}. Geometrik olarak sadece dönmeler ve yansımalar.

Simetrik matrislerin üç özelliği

AA reel ve simetrik (AT=AA^T = A) ise:

  1. Bütün özdeğerleri reeldir.
  2. Farklı özdeğerlere ait özvektörler diktir. İspat: Au=λuA\mathbf{u} = \lambda\mathbf{u}, Av=μvA\mathbf{v} = \mu\mathbf{v}, λ≠μ\lambda \neq \mu. λ(u⋅v)=(Au)Tv=uTATv=uTAv=μ(u⋅v)\lambda(\mathbf{u} \cdot \mathbf{v}) = (A\mathbf{u})^T\mathbf{v} = \mathbf{u}^TA^T\mathbf{v} = \mathbf{u}^TA\mathbf{v} = \mu(\mathbf{u} \cdot \mathbf{v}). (λ−μ)(u⋅v)=0(\lambda - \mu)(\mathbf{u} \cdot \mathbf{v}) = 0, yani u⋅v=0\mathbf{u} \cdot \mathbf{v} = 0. (AT=AA^T = A tam olarak burada kullanıldı.)
  3. Her özdeğer için geometrik katlılık = cebirsel katlılık. Sonuç: her zaman köşegenleşir.

Spektral teorem. AA n×nn \times n reel bir matris olsun. Şu iki ifade denktir:

  • AA simetrik.
  • PTAP=DP^TAP = D olacak şekilde ortogonal bir PP ve köşegen bir DD vardır (ortogonal köşegenleşme).

Yani “ortogonal köşegenleşir mi?” sorusunun cevabı tek bakışta: matris simetrik mi?

Ortogonal köşegenleştirme tarifi

  1. Karakteristik denklemi çöz.
  2. Her özdeğer için özuzayın bir bazını bul.
  3. Bir özuzay 2 veya daha fazla boyutluysa, bazına Gram-Schmidt uygulayıp ortogonal yap. (Tek boyutlu özuzaylar ve farklı özuzaylar arası diklik zaten bedava.)
  4. Her vektörü normalize et: v/∥v∥\mathbf{v}/\|\mathbf{v}\|.
  5. PP‘nin sütunları bu ortonormal vektörler, DD‘de aynı sırayla özdeğerler. A=PDPTA = PDP^T.

Gram-Schmidt (iki vektör). {v1,v2}\{\mathbf{v}_1, \mathbf{v}_2\} bağımsızsa

u1=v1,u2=v2−v2⋅u1u1⋅u1u1.\mathbf{u}_1 = \mathbf{v}_1, \qquad \mathbf{u}_2 = \mathbf{v}_2 - \frac{\mathbf{v}_2 \cdot \mathbf{u}_1}{\mathbf{u}_1 \cdot \mathbf{u}_1}\mathbf{u}_1.

u2\mathbf{u}_2, v2\mathbf{v}_2‘den u1\mathbf{u}_1 yönündeki bileşeni çıkarır; kalan dik parça. Üç vektörde u3\mathbf{u}_3‘ten hem u1\mathbf{u}_1 hem u2\mathbf{u}_2 bileşenleri çıkarılır.

Örnek 1 (2×2). A=[3113]A = \begin{bmatrix} 3 & 1 \\ 1 & 3 \end{bmatrix}. İz 6, det 8: λ2−6λ+8=(λ−4)(λ−2)\lambda^2 - 6\lambda + 8 = (\lambda - 4)(\lambda - 2).

λ=4:(1,1),λ=2:(1,−1),P=12[111−1],D=[4002].\lambda = 4: (1, 1), \qquad \lambda = 2: (1, -1), \qquad P = \frac{1}{\sqrt2}\begin{bmatrix} 1 & 1 \\ 1 & -1 \end{bmatrix}, \quad D = \begin{bmatrix} 4 & 0 \\ 0 & 2 \end{bmatrix}.

Özvektörler kendiliğinden dik çıktı; tek iş normalize etmekti.

Örnek 2 (3×3, tekrarlı özdeğer). A=[211121112]A = \begin{bmatrix} 2 & 1 & 1 \\ 1 & 2 & 1 \\ 1 & 1 & 2 \end{bmatrix}.

Karakteristik polinom λ3−6λ2+9λ−4=(λ−1)2(λ−4)\lambda^3 - 6\lambda^2 + 9\lambda - 4 = (\lambda - 1)^2(\lambda - 4). Kontrol: iz 6=1+1+46 = 1 + 1 + 4.

  • λ=4\lambda = 4: A−4IA - 4I‘nın satırları (−2,1,1),(1,−2,1),(1,1,−2)(-2, 1, 1), (1, -2, 1), (1, 1, -2); çözüm x=y=zx = y = z: (1,1,1)(1, 1, 1).
  • λ=1\lambda = 1: A−IA - I‘nın bütün satırları (1,1,1)(1, 1, 1): tek denklem x+y+z=0x + y + z = 0. Baz v1=(−1,1,0)\mathbf{v}_1 = (-1, 1, 0), v2=(−1,0,1)\mathbf{v}_2 = (-1, 0, 1). Bunlar dik değil (v1⋅v2=1\mathbf{v}_1 \cdot \mathbf{v}_2 = 1).

Gram-Schmidt (özuzay içinde):

u1 = (-1, 1, 0)
v2·u1 = 1 + 0 + 0 = 1,   u1·u1 = 2
u2 = (-1, 0, 1) - (1/2)(-1, 1, 0) = (-1/2, -1/2, 1)   → 2 ile çarp: (-1, -1, 2)
kontrol: u1·u2 = 1 - 1 + 0 = 0;  A u2 = (-2-1+2, -1-2+2, -1-1+4) = (-1, -1, 2) = 1·u2

(1,1,1)(1, 1, 1) ikisine de zaten dik (farklı özdeğer). Normalize: uzunluklar 2\sqrt2, 6\sqrt6, 3\sqrt3.

P=[−12−161312−161302613],D=[100010004].P = \begin{bmatrix} -\frac{1}{\sqrt2} & -\frac{1}{\sqrt6} & \frac{1}{\sqrt3} \\ \frac{1}{\sqrt2} & -\frac{1}{\sqrt6} & \frac{1}{\sqrt3} \\ 0 & \frac{2}{\sqrt6} & \frac{1}{\sqrt3} \end{bmatrix}, \qquad D = \begin{bmatrix} 1 & 0 & 0 \\ 0 & 1 & 0 \\ 0 & 0 & 4 \end{bmatrix}.

TUZAK · tekrarlı özdeğerde Gram-Schmidt’i atlamak. Aynı özuzaydan gelen iki özvektör dik olmak zorunda değil. Örnek 2’de (−1,1,0)(-1, 1, 0) ve (−1,0,1)(-1, 0, 1) ikisi de doğru özvektör ama aralarındaki açı 60∘60^\circ. Bunları direkt normalize edip PP‘ye koyarsan PTP≠IP^TP \neq I olur.

İPUCU · kontrol. PTP=IP^TP = I‘yı kontrol etmek için her sütun çiftinin iç çarpımına bak (0 olmalı) ve her sütunun uzunluğuna (1 olmalı). PTAPP^TAP‘yi tamamen çarpmaktan çok daha kısa.

Uygulama: kuadratik formlar

Q(x,y)=ax2+bxy+cy2Q(x, y) = ax^2 + bxy + cy^2 ifadesi bir kuadratik formdur ve simetrik bir matrisle yazılır:

Q(x)=xTAx,A=[ab/2b/2c].Q(\mathbf{x}) = \mathbf{x}^TA\mathbf{x}, \qquad A = \begin{bmatrix} a & b/2 \\ b/2 & c \end{bmatrix}.

A=PDPTA = PDP^T ise x=Py\mathbf{x} = P\mathbf{y} değişkeniyle Q=λ1y12+λ2y22Q = \lambda_1y_1^2 + \lambda_2y_2^2 olur: çapraz terim yok. PP‘nin sütunları asal eksenlerdir (principal axes). Q11’de 3x2+4xy+3y2=5x′2+y′23x^2 + 4xy + 3y^2 = 5x'^2 + y'^2; (1,1)(1, 1) ve (−1,1)(-1, 1) doğrultuları elipsin eksenleri. Özdeğerlerin işaretleri formun hep pozitif olup olmadığını (pozitif tanımlı) söyler.

Dönem özeti: köşegenleşme hiyerarşisi

Matris türüKöşegenleşir mi?PP
Genel kare matrisbazen (nn bağımsız özvektör varsa)tersinir
nn farklı özdeğerher zamantersinir
Simetrikher zamanortogonal (P−1=PTP^{-1} = P^T)

Sık yapılan hatalar

YanlışDoğrusu
Dik sütunlu matrisi ortogonal saymakSütunlar birim de olmalı
Normalize etmeyi unutmakHer sütunu kendi uzunluğuna böl
Tekrarlı özdeğerin özvektörlerini olduğu gibi kullanmakÖnce Gram-Schmidt
Farklı özuzaylar arasında Gram-Schmidt yapmakGerek yok, zaten dikler
A=PDP−1A = PDP^{-1} yazıp P−1P^{-1} hesaplamakP−1=PTP^{-1} = P^T
Kuadratik formda xyxy katsayısını köşegen dışına aynen yazmakYarısı: b/2b/2
Simetrik olmayan matrisi ortogonal köşegenleştirmeye çalışmakSadece simetrikler
det⁡Q\det Q her zaman 1±1\pm1

Kendini test et

Sınavdaki gibi İngilizce sorular. Önce kâğıt kalemle çöz, sonra cevabı aç.

Soru 1kolay

Which of the matrices A=[1221]A = \begin{bmatrix} 1 & 2 \\ 2 & 1 \end{bmatrix} and Q=12[11−11]Q = \frac{1}{\sqrt2}\begin{bmatrix} 1 & 1 \\ -1 & 1 \end{bmatrix} is orthogonal?

Cevabı göster
Only QQ

Ortogonal: QTQ=IQ^TQ = I, yani sütunlar birim uzunlukta ve birbirine dik.QQ: sütunlar 12(1,−1)\frac{1}{\sqrt2}(1, -1) ve 12(1,1)\frac{1}{\sqrt2}(1, 1), iç çarpım 12(1−1)=0\frac12(1 - 1) = 0, uzunluklar 1.AA: ATA=[5445]≠IA^TA = \begin{bmatrix} 5 & 4 \\ 4 & 5 \end{bmatrix} \neq I. Simetrik olmak ortogonal olmak değil.

Soru 2kolay

Find the eigenvalues of the symmetric matrix A=[5222]A = \begin{bmatrix} 5 & 2 \\ 2 & 2 \end{bmatrix}.

Cevabı göster
66 and 11

İz 7, det 10−4=610 - 4 = 6: λ2−7λ+6=(λ−6)(λ−1)\lambda^2 - 7\lambda + 6 = (\lambda - 6)(\lambda - 1).
Simetrik olduğu için reel olacakları baştan belliydi.

Soru 3orta

Orthogonally diagonalize A=[5222]A = \begin{bmatrix} 5 & 2 \\ 2 & 2 \end{bmatrix}: find an orthogonal PP and diagonal DD with PTAP=DP^TAP = D.

Cevabı göster
P=15[211−2], D=[6001]P = \frac{1}{\sqrt5}\begin{bmatrix} 2 & 1 \\ 1 & -2 \end{bmatrix},\ D = \begin{bmatrix} 6 & 0 \\ 0 & 1 \end{bmatrix}

λ=6\lambda = 6: A−6I=[−122−4]A - 6I = \begin{bmatrix} -1 & 2 \\ 2 & -4 \end{bmatrix}, x=2yx = 2y: (2,1)(2, 1).λ=1\lambda = 1: A−I=[4221]A - I = \begin{bmatrix} 4 & 2 \\ 2 & 1 \end{bmatrix}, 2x+y=02x + y = 0: (1,−2)(1, -2).
İkisi dik: 2−2=02 - 2 = 0. İkisinin de uzunluğu 5\sqrt5; normalize edip sütun yap.

Soru 4kolay

Find kk such that (1,k,2)(1, k, 2) and (3,1,−1)(3, 1, -1) are orthogonal.

Cevabı göster
k=−1k = -1

Dik   ⟺  \iff iç çarpım 0: 3+k−2=03 + k - 2 = 0, k=−1k = -1.

Soru 5orta

Find the inverse of Q=[3/5−4/54/53/5]Q = \begin{bmatrix} 3/5 & -4/5 \\ 4/5 & 3/5 \end{bmatrix}.

Cevabı göster
[3/54/5−4/53/5]\begin{bmatrix} 3/5 & 4/5 \\ -4/5 & 3/5 \end{bmatrix}

Sütunlar (3/5,4/5)(3/5, 4/5) ve (−4/5,3/5)(-4/5, 3/5): uzunluklar 9+1625=1\frac{9 + 16}{25} = 1, iç çarpım −1225+1225=0-\frac{12}{25} + \frac{12}{25} = 0. QQ ortogonal.
Ortogonal matrisin tersi transpozudur: Q−1=QTQ^{-1} = Q^T. Hesap yok.

Soru 6orta

True or false: Every real symmetric matrix is diagonalizable.

Cevabı göster
True

Spektral teorem: reel simetrik matrisin özdeğerleri reeldir ve her özdeğer için geometrik katlılık cebirsel katlılığa eşittir. Hatta ortogonal bir PP ile köşegenleşir.
Simetrik olmayan [2102]\begin{bmatrix} 2 & 1 \\ 0 & 2 \end{bmatrix} köşegenleşmiyordu; simetrikte bu olamaz.

Soru 7zor

AA is a 3×33 \times 3 symmetric matrix. (1,2,−1)(1, 2, -1) is an eigenvector for λ=3\lambda = 3 and (2,a,4)(2, a, 4) is an eigenvector for λ=−1\lambda = -1. Find aa.

Cevabı göster
a=1a = 1

Simetrik matriste farklı özdeğerlerin özvektörleri diktir.(1,2,−1)⋅(2,a,4)=2+2a−4=0⇒a=1(1, 2, -1) \cdot (2, a, 4) = 2 + 2a - 4 = 0 \Rightarrow a = 1.

Soru 8zor

Orthogonally diagonalize A=[102030201]A = \begin{bmatrix} 1 & 0 & 2 \\ 0 & 3 & 0 \\ 2 & 0 & 1 \end{bmatrix}.

Cevabı göster
P=[1/201/20101/20−1/2], D=diag⁡(3,3,−1)P = \begin{bmatrix} 1/\sqrt2 & 0 & 1/\sqrt2 \\ 0 & 1 & 0 \\ 1/\sqrt2 & 0 & -1/\sqrt2 \end{bmatrix},\ D = \operatorname{diag}(3, 3, -1)

Özdeğerler 3,3,−13, 3, -1 ((λ−3)2(λ+1)(\lambda - 3)^2(\lambda + 1)).λ=3\lambda = 3: A−3I=[−20200020−2]A - 3I = \begin{bmatrix} -2 & 0 & 2 \\ 0 & 0 & 0 \\ 2 & 0 & -2 \end{bmatrix}, x=zx = z, yy serbest: (1,0,1)(1, 0, 1) ve (0,1,0)(0, 1, 0). Bunlar zaten dik, Gram-Schmidt gerekmez.λ=−1\lambda = -1: A+I=[202040202]A + I = \begin{bmatrix} 2 & 0 & 2 \\ 0 & 4 & 0 \\ 2 & 0 & 2 \end{bmatrix}, y=0y = 0, x=−zx = -z: (1,0,−1)(1, 0, -1).
Normalize: (1,0,1)/2(1, 0, 1)/\sqrt2, (0,1,0)(0, 1, 0), (1,0,−1)/2(1, 0, -1)/\sqrt2.

Soru 9orta

What are the possible values of det⁡Q\det Q for an orthogonal matrix QQ?

Cevabı göster
11 or −1-1

QTQ=I⇒det⁡QTdet⁡Q=1⇒(det⁡Q)2=1Q^TQ = I \Rightarrow \det Q^T\det Q = 1 \Rightarrow (\det Q)^2 = 1.
Dönmeler +1+1 (Q5'teki matris), yansımalar −1-1 (diag⁡(1,−1)\operatorname{diag}(1, -1)).

Soru 10zor

Apply the Gram-Schmidt process to {(1,1,0),(1,0,1)}\{(1, 1, 0), (1, 0, 1)\} to get an orthogonal basis of their span.

Cevabı göster
{(1,1,0), (1,−1,2)}\{(1, 1, 0),\ (1, -1, 2)\}

u1=(1,1,0)\mathbf{u}_1 = (1, 1, 0). u2=(1,0,1)−(1,0,1)⋅u1u1⋅u1u1=(1,0,1)−12(1,1,0)=(12,−12,1)\mathbf{u}_2 = (1, 0, 1) - \frac{(1, 0, 1) \cdot \mathbf{u}_1}{\mathbf{u}_1 \cdot \mathbf{u}_1}\mathbf{u}_1 = (1, 0, 1) - \frac12(1, 1, 0) = \left(\tfrac12, -\tfrac12, 1\right).
2 ile çarp (diklik bozulmaz): (1,−1,2)(1, -1, 2). Kontrol: 1−1+0=01 - 1 + 0 = 0.

Soru 11zor

Write the quadratic form Q(x,y)=3x2+4xy+3y2Q(x, y) = 3x^2 + 4xy + 3y^2 as xTAx\mathbf{x}^TA\mathbf{x} with AA symmetric, and find the eigenvalues of AA.

Cevabı göster
A=[3223]A = \begin{bmatrix} 3 & 2 \\ 2 & 3 \end{bmatrix}; eigenvalues 55 and 11

Köşegen: x2x^2 ve y2y^2 katsayıları. xyxy katsayısı 4, iki köşegen dışı elemana yarı yarıya: 2 ve 2.
Özdeğerler: iz 6, det 5, λ2−6λ+5\lambda^2 - 6\lambda + 5: 55 ((1,1)(1, 1) yönü) ve 11 ((−1,1)(-1, 1) yönü).
Yeni eksenlerde Q=5x′2+y′2Q = 5x'^2 + y'^2: çapraz terim yok olur.

Doğrulama: Özdeğerler ve özuzaylar node'da BigInt kesirlerle karakteristik polinom kökleri ve tam null space ile bulundu; karekök içeren normalizasyonlardan kaçınmak için ölçeksiz sütun matrisi S ile SᵀS = diag(‖vᵢ‖²) (diklik) ve SᵀAS = diag(‖vᵢ‖²λᵢ) (köşegenleşme) tam olarak kontrol edildi; Gram-Schmidt adımları ve iç çarpımlar tam kesirlerle hesaplandı (verify/MAT1041/14.mjs).

Deftere not kartı

Simetrik: A^T = A. ORTOGONAL MATRİS: Q^T Q = I ⇔ Q⁻¹ = Q^T ⇔ sütunlar ORTONORMAL (birim + dik). det Q = ±1. Uzunluk ve açıyı korur: ||Qx|| = ||x||. SİMETRİK MATRİS (spektral teorem): (1) özdeğerler REEL (2) farklı özdeğerlerin özvektörleri DİK (3) HER ZAMAN köşegenleşir, hem de ORTOGONAL P ile: P^T A P = D. Tersi de doğru: ortogonal köşegenleşen matris simetriktir. TARİF: özdeğerler → her E_λ bazı → (tekrarlı λ'da) Gram-Schmidt → NORMALİZE → P'nin sütunları. Gram-Schmidt: u1 = v1; u2 = v2 - (v2·u1)/(u1·u1) u1. Ölçeklemek serbest. Farklı özuzaylardaki vektörler zaten dik; GS sadece AYNI özuzay içinde. Normalize: v/||v||. UNUTMA: normalize edilmemiş P ortogonal değil. A = PDP^T (P⁻¹ yerine P^T). Kuadratik form ax² + bxy + cy² → A = [a b/2; b/2 c]; özvektörler = asal eksenler.