Μια Εισαγωγή στη Γραμμική Άλγεβρα για τις Θετικές Επιστήμες

Κεφάλαιο 5Ιδιοτιμές, Ιδιοδιανύσματα

Σε αυτό το κεφάλαιο μελετούμε τις ιδιοτιμές και τα ιδιοδιανύσματα ενός τετραγωνικού πίνακα. Ως συνήθως, k είναι το σώμα ℝ ή το σώμα ℂ.

5.1 Ιδιοτιμές και Ιδιοδιανύσματα

Έστω V ένας k-διανυσματικός χώρος και ϕ:V→V μία γραμμική συνάρτηση. Σε αυτήν την ενότητα θα δούμε πως εντοπίζουμε μη μηδενικά διανύσματα v που οι εικόνες τους, ϕ⁢(v), είναι διανύσματα παράλληλα προς τα ίδια τα v. Τα διανύσματα αυτά τα ονομάζουμε ιδιοδιανύσματα της ϕ. Στο Σχήμα 5.1 με μπλε χρώμα απεικονίζεται ένα ιδιοδιάνυσμα v της ϕ, ενώ με κόκκινο χρώμα απεικονίζεται η εικόνα του.

Σχήμα 5.1: Ιδιοδιάνυσμα για την ιδιοτιμή λ
Ορισμός 5.1.1.

Έστω V ένας k-διανυσματικός χώρος και ϕ:V→V μία γραμμική συνάρτηση. To λ∈k λέγεται ιδιοτιμή (eigenvalue) της ϕ, αν υπάρχει v∈V, v≠0, έτσι ώστε

ϕ⁢(v)=λ⁢v. (5.1.1.1)

Το v της Σχέσης (5.1.1.1) λέγεται ιδιοδιάνυσμα (eigenvector) της ϕ για την ιδιοτιμή λ.

Παραδείγματα 5.1.2.

  1. 1.

    Η ταυτοτική συνάρτηση idV:V→V, idV⁡(v)=v έχει το 1 για ιδιοτιμή, ενώ κάθε v≠𝟎∈V είναι ιδιοδιάνυσμα για την ιδιοτιμή 1. Είναι φανερό ότι η γραμμική συνάρτηση idV:V→V δεν έχει άλλη ιδιοτιμή.

  2. 2.

    Κάθε μη μηδενικό στοιχείο v του V είναι ιδιοδιάνυσμα για τη μηδενική (τετριμμένη) συνάρτηση f:V→V για την ιδιοτιμή 0 αφού f⁢(v)=𝟎. Είναι φανερό ότι η μόνη ιδιοτιμή της μηδενικής συνάρτησης είναι το 0.

Έχοντας μελετήσει στο προηγούμενο κεφάλαιο τη στενή σχέση μεταξύ των γραμμικών συναρτήσεων και των πινάκων, δίνουμε στη συνέχεια τους αντίστοιχους ορισμούς για τους τετραγωνικούς πίνακες.

Ορισμός 5.1.3.

Αν Α∈ℳn⁢(k), τότε b=[b1⋯bn]T∈ℳn×1⁢(k) λέγεται ιδιοδιάνυσμα (eigenvector) του A για την ιδιοτιμή (eigenvalue) λ, αν b≠𝟎 και A⁢b=λ⁢b, δηλ.

A⁢[b1⋮bn]=λ⁢[b1⋮bn]. (5.1.3.1)

Σύμφωνα με το Θεώρημα 4.2.14, κάθε k-διανυσματικός χώρος V διάστασης n είναι ισόμορφος με τον kn . Για να απλοποιήσουμε, λοιπόν, τους συμβολισμούς μας, θα περιοριστούμε στην περίπτωση που V=kn. Σημειώνουμε την επόμενη πρόταση.

Πρόταση 5.1.4.

Έστω ϕ:kn→kn μία γραμμική συνάρτηση , D μία διατεταγμένη βάση του kn. To v∈kn είναι ιδιοδιάνυσμα της ϕ για την ιδιοτιμή λ αν και μόνο αν CD⁢(v) είναι ιδιοδιάνυσμα του πίνακα AD,Dϕ για την ιδιοτιμή λ. Ειδικότερα, αν B είναι η κανονική βάση του kn, τότε (b1,…,bn) είναι ιδιοδιάνυσμα της ϕ αν και μόνο αν b=[b1⋯bn]T είναι ιδιοδιάνυσμα του AB,Bϕ .

Απόδειξη.

Έστω ότι το v είναι ιδιοδιάνυσμα της ϕ για την ιδιοτιμή λ, δηλ. ϕ⁢(v)=λ⁢v. Τότε v≠𝟎, άρα CD⁢(v)≠𝟎. Είναι φανερό ότι CD⁢(λ⁢v)=λ⁢CD⁢(v). Από την Πρόταση 4.1.4, έχουμε ότι:

ΑD,Dϕ⁢CD⁢(v)=CD⁢(ϕ⁢(v))=CD⁢(λ⁢v)=λ⁢CD⁢(v).

Άρα το CD⁢(v) είναι ιδιοδιάνυσμα του ΑD,Dϕ. To αντίστροφο προκύπτει με τον ίδιο ακριβώς τρόπο.

Για το τελευταίο συμπέρασμα της πρότασης σημειώνουμε ότι αν B είναι η κανονική βάση του kn και v=(b1,…,bn), τότε CB⁢(v)=[b1⋯bn]T. ∎

Στα επόμενο παραδείγματα θα δούμε πως βρίσκουμε τα ιδιοδιανύσματα και τις ιδιοτιμές, χρησιμοποιώντας την παρατήρηση αυτή και τη γεωμετρική εποπτεία για το πραγματικό επίπεδο και τον πραγματικό τρισδιάστατο χώρο.

Παραδείγματα 5.1.5.

  1. 1.

    O μοναδιαίος πίνακας In έχει την ιδιοτιμή 1 με ιδιοδιάνυσμα κάθε v≠0∈kn. Σημειώνουμε ότι o In είναι o πίνακας της ταυτοτικής συνάρτησης idkn:kn→kn, id⁡(v)=v ως προς την κανονική βάση.

  2. 2.

    O μηδενικός n×n πίνακας

    𝟎=[0⋯0⋮⋱⋮0⋯0]

    έχει την ιδιοτιμή 0 με ιδιοδιάνυσμα κάθε v≠0∈kn. Σημειώνουμε ότι o 𝟎 είναι o πίνακας της μηδενικής συνάρτησης kn→kn, v↦𝟎 ως προς οποιαδήποτε βάση.

  3. 3.

    Η γραμμική συνάρτηση f:ℝ2→ℝ2 που περιστρέφει αριστερόστροφα κατά γωνία π/4 τα στοιχεία του ℝ2 (βλ. Σχήμα 5.2) δεν έχει ιδιοτιμές και ιδιοδιανύσματα. Για το συμπέρασμα αυτό αρκεί να παρατηρήσουμε τη γραφική απεικόνιση της συνάρτησης: όταν v≠0 τότε f⁢(v) και v ορίζουν δύο ευθείες που σχηματίζουν μεταξύ τους γωνία π/4. Άρα δεν υπάρχει λ έτσι ώστε f⁢(v)=λ⁢v.

    Σχήμα 5.2: Αριστερόστροφη περιστροφή κατά γωνία π/4

    Έστω

    A=AB⁢Bf=22⁢[1-111].

    Αφού η f δεν έχει ιδιοτιμές και ιδιοδιανύσματα, ούτε o A έχει ιδιοτιμές και ιδιοδιανύσματα στον ℝ.

  4. 4.

    Έστω f:ℝ2→ℝ2 η προβολή στην ευθεία y=m⁢x, όπως απεικονίζεται στο Σχήμα 4.4.

    Σχήμα 5.3: Προβολή στον ℝ2 στην ευθεία y=m⁢x

    Η ευθεία y=m⁢x περνά από την αρχή των αξόνων O⁢(0,0) και το σημείο (1,m). Είναι λοιπόν παράλληλη προς το διάνυσμα (1,m) και κάθετη προς το διάνυσμα (-m,1). Από τη γραφική απεικόνιση της συνάρτησης, συμπεραίνουμε ότι η f έχει δύο ιδιοτιμές: το 1 και το 0. Πράγματι αν v1=(1,m) τότε f⁢(v1)=v1 και άρα το v1 είναι ιδιοδιάνυσμα για το 1. Αντίστοιχα αν v2=(-m,1) τότε f⁢(v2)=0 και το v2 είναι ιδιοδιάνυσμα για την ιδιοτιμή 0. Γενικότερα, τα πολλαπλάσια του v1 είναι όλα ιδιοδιανύσματα για την τιμή 1, ενώ τα πολλαπλάσια του v2 είναι ιδιοδιανύσματα της f για την ιδιοτιμή 0.

    Στη συνέχεια θα βρούμε τον τύπο της f. Έστω D=(v1,v2). Ο πίνακας ΑD,Df είναι διαγώνιος, με τις ιδιοτιμές να εμφανίζονται ως στοιχεία της κυρίας διαγωνίου. Πράγματι f⁢(v1)=v1 και f⁢(v2)=0, επομένως

    CD⁢(f⁢(v1))=[10]⁢ και ⁢CD⁢(f⁢(v2))=[00].

    Άρα

    ΑD,Df=[1000].

    Έστω τώρα B=(e1,e2),

    S=SB←D=[1-mm1].

    Τότε

    ΑΒ,Βf=S⁢AB,Bf⁢S-1=11+m2⁢[1mmm2].

    Αφού

    ΑΒ,Βf⁢[ab]=11+m2⁢[a+m⁢bm⁢a+m2⁢b],

    έπεται ότι

    f⁢(a,b)=11+m2⁢(a+m⁢b,m⁢a+m2⁢b).
  5. 5.

    Έστω f:ℝ2→ℝ2 o αντικατοπτρισμός ως προς την ευθεία y=m⁢x, όπως απεικονίζεται στο Σχήμα 5.4.

    Σχήμα 5.4: Αντικατοπτρισμός ως προς την ευθεία y=m⁢x

    H f έχει την ιδιοτιμή 1 με ιδιοδιανύσματα τα πολλαπλάσια του v1=(1,m), και την ιδιοτιμή -1 με ιδιοδιανύσματα τα πολλαπλάσια του v2=(-m,1). Έτσι ως προς τη βάση D=(v1,v2) έχουμε ότι

    ΑD,Df=[100-1].

    Έστω Β η κανονική βάση του ℝ2 και S=SB←D. Τότε

    S=[1-mm1]

    και

    AB,Bf=S⁢AD,Df⁢S-1=11+m2⁢[1-m22⁢m2⁢mm2-1].

    Αφού

    AB,bf⁢[ab]=11+m2⁢[a⁢(1-m2)+2⁢b⁢m2⁢a⁢m+b⁢(m2-1)],

    o αναλυτικός τύπος της f είναι

    f⁢(a,b)=11+m2⁢(a⁢(1-m2)+2⁢b⁢m,2⁢a⁢m+b⁢(m2-1)).

Αν το λ είναι ιδιοτιμή της γραμμικής συνάρτησης ϕ:kn→kn , τότε o ιδιοχώρος (eigenspace) της ϕ, Vλ⁢(ϕ), είναι το σύνολο

Vλ⁢(ϕ):={v∈kn:ϕ⁢(v)=λ⁢v}.

Ο Vλ⁢(ϕ) είναι, λοιπόν, το σύνολο όλων των ιδιοδιανυσμάτων της ϕ για την ιδιοτιμή λ, μαζί με το 𝟎. Αντίστοιχα, αν A είναι ένας n×n πίνακας,

Vλ⁢(A):={b∈ℳn×1⁢(k):A⁢b=λ⁢b}.

Vλ⁢(A) Αφού σε κάθε n×n πίνακα αντιστοιχεί μία γραμμική συνάρτηση ϕ:kn→kn και αντίστροφα, είναι φανερό ότι ιδιότητες του Vλ⁢(ϕ) μεταφράζονται σε ιδιότητες του Vλ⁢(A) και τανάπαλιν. Παρατηρούμε ότι o Vλ⁢(ϕ) είναι διανυσματικός χώρος. Πράγματι:

  • •

    αν v1,v2∈Vλ⁢(ϕ), τότε ϕ⁢(v1+v2)=ϕ⁢(v1)+ϕ⁢(v2)=λ⁢v1+λ⁢v2=λ⁢(v1+v2) και επομένως v1+v2∈Vλ⁢(ϕ).

  • •

    Επίσης, αν v∈Vλ⁢(ϕ) και κ∈k, τότε ϕ⁢(κ⁢v)=κ⁢ϕ⁢(v)=κ⁢(λ⁢v)=(κ⁢λ)⁢v=λ⁢(κ⁢v)=λ⁢ϕ⁢(v) και επομένως κ⁢v∈Vλ⁢(ϕ).

Τέλος, μία ενδιαφέρουσα ιδιότητα των χώρων Vλ⁢(ϕ), είναι ότι παραμένουν αναλλοίωτοι (invariant) από τη δράση της ϕ , δηλ.

ϕ⁢(Vλ⁢(ϕ))⊂Vλ⁢(ϕ).
Παραδείγματα 5.1.6.

  1. 1.

    Δίνεται η γραμμική συνάρτηση

    ϕ:ℝ2→ℝ2,(x,y)↦(2⁢x+y,4⁢x+2⁢y).

    Έστω Β=(e1,e2) η κανονική βάση του ℝ2. Τότε

    AB,Bϕ=[2142]

    Παρατηρούμε ότι

    ϕ⁢(1,-2)=(0,0)=0⁢(1,-2),

    δηλ. το 0 είναι ιδιοτιμή και v1=(1,-2) είναι ιδιοδιάνυσμα της ϕ για την ιδιοτιμή 0 όπως βέβαια και κάθε άλλο (μη μηδενικό) στοιχείο του Ker⁡ϕ. Επίσης,

    ϕ⁢(1,2)=(4,8)=4⁢(1,2),

    δηλ. το 4 είναι ιδιοτιμή του ϕ και v2=(1,2) είναι ιδιοδιάνυσμα για την ιδιοτιμή 4. Τα v1,v2 είναι γραμμικά ανεξάρτητα και παράγουν τον ℝ2. Άρα D=(v1,v2) είναι βάση του ℝ2. Επομένως

    ΑD,Dϕ=[0004],

    και

    ΑD,Dϕ=S-1⁢AB,Bϕ⁢S⁢ όπου ⁢S=SB←D=[  11-22].

    O αναγνώστης καλείται να αποδείξει ότι V0(ϕ)=S({(v1}) και ότι V4⁢(ϕ)=S⁢({v2}).

  2. 2.

    Δίνεται η γραμμική συνάρτηση

    ϕ:ℝ2→ℝ2,(x,y)↦(2⁢x,3⁢y).

    Έστω Β=(e1,e2) η κανονική βάση του ℝ2. Παρατηρούμε ότι

    A=AB,Bϕ=[2003]

    είναι διαγώνιος πίνακας.: ϕ⁢(e1)=2⁢e1 και ϕ⁢(e2)=3⁢e2. Επομένως, e1 είναι ιδιοδιάνυσμα για την ιδιοτιμή 2, ενώ το e2 είναι ιδιοδιάνυσμα για την ιδιοτιμή 3. O αναγνώστης καλείται να αποδείξει ότι

    V2⁢(ϕ)=S⁢({e1})⁢ και ⁢V3⁢(ϕ)=S⁢({e2}).
  3. 3.

    Έστω ϕ:ℝ3→ℝ3,ϕ⁢(x,y,z)=(3⁢x,3⁢y,z) και Β=(e1,e2,e3) η κανονική βάση του ℝ3. Παρατηρούμε ότι

    AB,Bϕ=[300030001].

    O AB,Bϕ είναι διαγώνιος και τα e1, e2 είναι ιδιοδιανύσματα για την ιδιοτιμή 3, ενώ το e3 είναι ιδιοδιάνυσμα για την ιδιοτιμή 1. O αναγνώστης καλείται να αποδείξει ότι

    V3⁢(ϕ)=S⁢({e1,e2}),V1⁢(ϕ)=S⁢({e3}).
  4. 4.

    Θα γενικεύσουμε τις παρατηρήσεις μας για τους διαγώνιους πίνακες και τις ιδιοτιμές τους. Έστω ότι f:kn→kn είναι γραμμική συνάρτηση με την ιδιότητα ότι υπάρχει D=(v1,…,vn) διατεταγμένη βάση του kn έτσι ώστε το vi να είναι ιδιοδιάνυσμα της f, για i=1,…,n. Έστω ότι λi είναι η ιδιοτιμή της f που αντιστοιχεί στο vi, για i=1,…,n. Αφού f⁢(vi)=λi⁢vi, o πίνακας των συντεταγμένων του f⁢(vi) ως προς τη βάση D έχει την εξής μορφή:

    CD⁢(f⁢(vi))=[0⋯λi⋯0]T↑i-θέση

    και επομένως o ΑD,Df είναι διαγώνιος πίνακας,

    AD,Df=[λ10⋱0λn].

    Αντίστροφα, αν D=(v1,…,vn) είναι βάση του kn, f:kn→kn και

    AD,Df=[λ1O⋱Oλn],

    τότε f⁢(vi)=λi⁢vi, για i=1,…,n. Επομένως το vi είναι ιδιοδιάνυσμα της f για την ιδιοτιμή λi, για i=1,…,n.

  5. 5.

    Έστω ϕ:V→V μία γραμμική συνάρτηση του V και έστω ότι το 0 είναι ιδιοτιμή του ϕ. Υπάρχει λοιπόν 𝟎≠v∈V έτσι ώστε ϕ⁢(v)=0⁢v=𝟎. Επομένως v∈Ker⁡ϕ και Ker⁡ϕ≠{𝟎}. Αντίστροφα, αν v∈Ker⁡ϕ και v≠𝟎, τότε ϕ⁢(v)=𝟎 και το v είναι ιδιοδιάνυσμα που αντιστοιχεί στην ιδιοτιμή 0. Καταλήγουμε έτσι στο εξής συμπέρασμα:

    Το 0 είναι ιδιοτιμή της γραμμικής συνάρτησης ϕ:V→V αν και μόνο αν η ϕ δεν είναι ένα προς ένα. Ισοδύναμα, το 0 είναι ιδιοτιμή του τετραγωνικού πίνακα Α αν και μόνο αν det⁡A=0.

  6. 6.

    Έστω ότι ϕ:V→V είναι ισομορφισμός. Αν το λ είναι ιδιοτιμή της ϕ και v είναι ιδιοδιάνυσμα της ϕ για το λ , τότε (όπως θα δούμε) το λ-1 είναι ιδιοτιμή της ϕ-1 και το v είναι ιδιοδιάνυσμα της ϕ-1 για το λ-1. Πράγματι, έστω ότι το v είναι ιδιοδιάνυσμα της ϕ που αντιστοιχεί στην ιδιοτιμή λ. Αφού η γραμμική συνάρτηση ϕ:V→V είναι ισομορφισμός, έπεται ότι λ≠0. Επομένως

    ϕ-1⁢(v)=ϕ-1⁢((λ-1⁢λ)⁢v)=λ-1⁢ϕ-1⁢(ϕ⁢(v))=λ-1⁢idV⁡(v)=λ-1⁢v.

    Δείξαμε λοιπόν ότι:

    Αν ϕ:V→V είναι ισομορφισμός και το λ είναι ιδιοτιμή της ϕ με ιδιοδιάνυσμα v, τότε το λ-1 είναι ιδιοτιμή του ϕ-1:V→V με ιδιοδιάνυσμα το v.

  7. 7.

    Δίνεται η γραμμική συνάρτηση ϕ:ℝ2→ℝ2, με ϕ⁢(x,y)=(y,-2⁢x+3⁢y). O πίνακας A της ϕ ως προς την κανονική βάση B του ℝ2 είναι

    A=[0-213].

    O αναγνώστης καλείται να επιβεβαιώσει ότι το 2 είναι ιδιοτιμή του A και της ϕ και ότι κάθε μη μηδενικό διάνυσμα επί της ευθείας {t⁢(1,-1):t∈ℝ} είναι ιδιοδιάνυσμα για την ιδιοτιμή 2. H ϕ-1:ℝ2→ℝ2, δίνεται από τον τύπο ϕ-1⁢(x,y)=12⁢(3⁢x-y,2⁢x). O αναγνώστης καλείται να επιβεβαιώσει ότι το 12 είναι ιδιοτιμή της ϕ-1.

  8. 8.

    Έστω o 2×2 μιγαδικός πίνακας

    A=[i101+i].

    To e1 είναι ιδιοδιάνυσμα για την ιδιοτιμή i και το (1,1) είναι ιδιοδιάνυσμα για την ιδιοτιμή 1+ι.

  9. 9.

    Έστω o 2×2 μιγαδικός πίνακας

    A=[0i-i1].

    O αναγνώστης καλείται να επιβεβαιώσει ότι ±1 είναι ιδιοτιμές του A.

  10. 10.

    Έστω A∈ℳn⁢(ℂ) και λ ιδιοτιμή του A. Τότε το λ¯ είναι ιδιοτιμή του A¯. Πράγματι, έστω ότι b είναι ιδιοδιάνυσμα του A για την ιδιοτιμή λ. Τότε Α⁢b=λ⁢b και συνεπώς

    Α⁢b¯=λ⁢b¯⇒Α¯⁢b¯=λ¯⁢b¯,

    δηλ. το b¯ είναι ιδιοδιάνυσμα του A¯ για την ιδιοτιμή λ¯.

  11. 11.

    Έστω A∈ℳn⁢(ℂ) ένας πίνακας, τέτοιος ώστε A=A¯T. O A λέγεται αυτοπροσαρτημένος (self-adjoint). Θα αποδείξουμε ότι οι ιδιοτιμές του A είναι πραγματικοί αριθμοί.

    Απόδειξη.

    Έστω ότι το b είναι ιδιοδιάνυσμα του Α για την ιδιοτιμή λ∈ℂ. Για να δείξουμε ότι το λ είναι πραγματικός αριθμός, θα αποδείξουμε ότι λ=λ¯. Υπολογίζουμε το γινόμενο (A⁢b)T⁢b¯:

    (A⁢b)T⁢b¯=(λ⁢b)Τ⁢b¯=λ⁢(bT⁢b¯). (5.1.6.1)

    Από το προηγούμενο παράδειγμα γνωρίζουμε ότι λ¯ είναι ιδιοτιμή του A¯ και ότι b¯ είναι ιδιοδιάνυσμα για το λ¯. Από την υπόθεση, A¯T=A. Συνεπώς A¯=AT. Θα υπολογίσουμε τώρα το γινόμενο (A⁢b)T⁢b¯, χρησιμοποιώντας αυτές τις παρατηρήσεις.

    (A⁢b)T⁢b¯=(bT⁢AT)⁢b¯=bT⁢(A¯⁢b¯)=bT⁢(λ¯⁢b¯)=λ¯⁢(bT⁢b¯). (5.1.6.2)

    Από τις Εκφράσεις (5.1.6.1) και (5.1.6.2) προκύπτει ότι

    λ⁢(bT⁢b¯)=λ¯⁢(bT⁢b¯)⇒(λ-λ¯)⁢(bT⁢b¯)=0.

    Όμως, bT⁢b¯≠0 (βλ. Παράδειγμα 2.1.10.2) και άρα λ=λ¯. Αυτό σημαίνει ότι λ∈ℝ και το ζητούμενο αποδείχθηκε.

    ∎

    Αποδείξαμε λοιπόν ότι:

    Αν A∈ℳn⁢(ℂ) είναι αυτοπροσαρτημένος πίνακας τότε οι ιδιοτιμές του A είναι πραγματικοί αριθμοί. Αν A∈ℳn⁢(ℝ) είναι συμμετρικός πίνακας, τότε οι ιδιοτιμές του A είναι πραγματικοί αριθμοί.

  12. 12.

    O συμμετρικός πίνακας

    A=[i002⁢i]

    δεν είναι αυτοπροσαρτημένος, (προσοχή: A=AT, όμως A∉ℳn⁢(ℝ)). Οι ιδιοτιμές του A είναι οι μιγαδικοί αριθμοί i,2⁢i.

στη συνέχεια γενικεύουμε το συμπέρασμα του Παραδείγματος 5.1.5. 5.

Πρόταση 5.1.7.

Έστω A∈ℳn⁢(k). Το λ είναι ιδιοτιμή του Α αν και μόνο αν det⁡(A-λ⁢In)=0. Αντίστοιχα, αν ϕ:kn→kn είναι μία γραμμική συνάρτηση και A είναι o πίνακας της ϕ ως προς οποιαδήποτε διατεταγμένη βάση του kn, τότε το λ είναι ιδιοτιμή της ϕ αν και μόνο αν det⁡(A-λ⁢In)=0. Αν το λ είναι ιδιοτιμή του Α, τότε τα ιδιοδιανύσματα του Α για την ιδιοτιμή λ, είναι οι μη μηδενικές λύσεις του A-λ⁢In=0.

Απόδειξη.

Έστω ότι b είναι ιδιοδιάνυσμα για τον Α που αντιστοιχεί στην ιδιοτιμή λ. Τονίζουμε ότι b≠𝟎. Ισχύει ότι

Α⁢b=λ⁢b⇔Α⁢b-λ⁢b=𝟎⇔Α⁢b-λ⁢In=𝟎⇔(A-λ⁢In)⁢b=0.

Συνεπώς το ομογενές σύστημα (A-λ⁢In)⁢X=𝟎 έχει τη μη μηδενική λύση, b. Αυτό μπορεί να γίνει μόνο όταν A-λ⁢In δεν είναι αντιστρέψιμος και det⁡(A-λ⁢In)=0. Για την αντίστροφη κατεύθυνση, παρατηρούμε ότι όλες οι συνεπαγωγές στον συλλογισμό μας αντιστρέφονται. ∎

Παραδείγματα 5.1.8.

  1. 1.

    O πίνακας

    A=[1-11  1]∈ℳ2⁢(ℝ)

    δεν έχει ιδιοτιμές και ιδιοδιανύσματα στον ℝ. Πράγματι,

    det⁡(A-λ⁢I2)=det⁡([1-λ-11 1-λ])=(1-λ)2+1=λ2-2⁢λ+2

    δεν έχει ρίζες στον ℝ, αφού η διακρίνουσα του λ2-2⁢λ+2 είναι αρνητικός αριθμός. To πολυώνυμολ2-2⁢λ+2 έχει όμως δύο μιγαδικές ρίζες, τις 1±i. Επομένως, A∈ℳ2×2⁢(ℂ) έχει δύο ιδιοτιμές στον ℂ. O αναγνώστης καλείται να διαπιστώσει ότι (i,1)∈ℂ2 είναι ιδιοδιάνυσμα του A για την ιδιοτιμή 1+i, ενώ (-i,1)∈ℂ2 είναι ιδιοδιάνυσμα του A για την ιδιοτιμή 1-i.

  2. 2.

    Έστω ότι A,B είναι δύο όμοιοι n×n πίνακες και ότι B=P-1⁢A⁢P για κάποιον αντιστρέψιμο πίνακα P. Αν b είναι ιδιοδιάνυσμα του A για την ιδιοτιμή λ, τότε P-1⁢b είναι ιδιοδιάνυσμα του B για την ιδιοτιμή λ.

    Απόδειξη.

    Αφού b είναι ιδιοδιάνυσμα του A για την ιδιοτιμή λ, A⁢b=λ⁢b. Επομένως:

    B⁢(P-1⁢b)=(P-1⁢A⁢P)⁢(P-1⁢b)=(P-1⁢A)⁢(P⁢P-1⁢b)=(P-1⁢A)⁢(In⁢b)=P-1(Ab)=P-1(λb)=λ(P-1b).

    ∎

    Σημειώνουμε το συμπέρασμα που προκύπτει από τα προηγούμενα.

    Όμοιοι πίνακες έχουν τις ίδιες ιδιοτιμές.

Παρατηρούμε ότι σύμφωνα με την Πρόταση 5.1.7,

Vλ⁢(ϕ)=null⁢(Α-λ⁢In)

και μάλιστα, αφού o Vλ⁢(ϕ) περιέχει ένα τουλάχιστον μη μηδενικό ιδιοδιάνυσμα, ισχύει ότι

dimk⁢Vλ⁢(ϕ)≥1 .

Ασκήσεις Ενότητας 5.1

  1. 1.

    Έστω η γραμμική συνάρτηση ϕ:ℝ2→ℝ2 που περιστρέφει τα διανύσματα αριστερόστροφα κατά γωνία π. Να βρεθούν οι ιδιοτιμές και τα ιδιοδιανύσματα της ϕ.

  2. 2.

    Έστω η γραμμική συνάρτηση ϕ:ℝ3→ℝ3 που αντικατοπτρίζει τα διανύσματα ως προς το επίπεδο y=0. Να βρεθούν οι ιδιοτιμές και τα ιδιοδιανύσματα της ϕ.

  3. 3.

    Να υπολογισθούν οι ιδιοτιμές και τα ιδιοδιανύσματα των εξής γραμμικών συναρτήσεων

    1. (a)

      ϕ:ℂ2→ℂ2, (x,y)↦(x+4⁢y,2⁢x+3⁢y).

    2. (b)

      ψ:ℝ3→ℝ3, (a,b,c)↦(b,c,a).

    3. (c)

      f:R3→R3, (x,y,z)↦(2⁢x+y,y-z,2⁢y+4⁢z).

  4. 4.

    Να υπολογισθούν οι ιδιοτιμές και τα ιδιοδιανύσματα των εξής πινάκων:

    1. (a)

      [1232],

    2. (b)

      [-22-23],

    3. (c)

      [1i-i1].

  5. 5.

    Έστω Α ένας n×n πίνακας, λ μία ιδιοτιμή του Α και X είναι ένα ιδιοδιάνυσμα του Α που αντιστοιχεί στην ιδιοτιμή λ. Να αποδείξετε ότι το λ+k είναι μία ιδιοτιμή του Α+k⁢In και ότι X είναι ιδιοδιάνυσμα του Α+k⁢In. Επίσης να αποδείξετε ότι το λm είναι μία ιδιοτιμή του Am και ότι X είναι ιδιοδιάνυσμα του Am που αντιστοιχεί στο λm.

  6. 6.

    Έστω A ένας n×n πίνακας και λ είναι μία ιδιοτιμή του 3Α. Να αποδείξετε ότι το λ2-λ+2 είναι μία ιδιοτιμή του πίνακα 3⁢Α2-Α+2⁢Ιn.

5.2 Ιδιοχώροι και το Χαρακτηριστικό Πολυώνυμο

Έστω ϕ:kn→kn μία γραμμική συνάρτηση. Για να βρούμε τις ιδιοτιμές της ϕ, βασιζόμαστε στις Προτάσεις 5.1.4 και 5.1.7 και προχωράμε ως εξής:

  • •

    βρίσκουμε τον πίνακα A της ϕ , ως προς την κανονική (ή οποιαδήποτε άλλη διατεταγμένη βάση) του kn,

  • •

    υπολογίζουμε την ορίζουσα του πίνακα A-x⁢In,

  • •

    βρίσκουμε τις τιμές των λ που μηδενίζουν την ορίζουσα det⁡(A-x⁢In).

Για να βρούμε τα ιδιοδιανύσματα της ϕ όταν γνωρίζουμε τις ιδιοτιμές, μπορούμε να λύσουμε τα συστήματα A⁢X=λ⁢X για κάθε ιδιοτιμή του A. Υπάρχει μήπως κάποιος αποτελεσματικότερος τρόπος; Για να απαντήσουμε σε αυτό το ερώτημα, θα ορίσουμε το χαρακτηριστικό πολυώνυμο του A. Σημειώνουμε ότι o πίνακας A-x⁢In είναι o n×n πίνακας

A-x⁢In=[α11-xα12⋯α1⁢nα21α22-x⋯α2⁢n⋮⋮⋮αn⁢1αn⁢2⋯αn⁢n-x].

Επομένως det⁡(A-x⁢In) είναι ένα πολυώνυμο στον k⁢[x], βαθμού n. Το πολυώνυμο αυτό συμβολίζεται με PA⁢(x), δηλ.

PA⁢(x):=det⁡(A-x⁢In), (5.2.0.1)

και λέγεται χαρακτηριστικό πολυώνυμο (characteristic polynomial) του A. Γράφουμε Pϕ⁢(x) για το PA⁢(x) όταν θέλουμε να δώσουμε την έμφαση στην ϕ. Οι ιδιοτιμές του Α μηδενίζουν το PA⁢(x) και είναι οι ρίζες του PA⁢(x). PA⁢(x) Παρατηρούμε ότι το PA⁢(x) έχει την εξής μορφή:

PA⁢(x)=(-1)n⁢xn+αn-1⁢xn-1+⋯+α1⁢x+α0.

O σταθερός όρος του PA⁢(x), a0, είναι ίσος με την ορίζουσα του A. Πράγματι,

PA⁢(0)=α0=det⁡(A-0⁢In)=det⁡A.

Ενδιαφέρον επίσης έχει o συντελεστής an-1 του PA⁢(x). Από τον υπολογισμό της ορίζουσας του A-x⁢In προκύπτει ότι o an-1 είναι με προσέγγιση προσήμου το ίχνος του A:ίχνος

αn-1=(-1)n-1⁢(α11+⋯+αn⁢n)=(-1)n-1⁢Tr⁡(A).
Παραδείγματα 5.2.1.

  1. 1.

    Έστω

    A=[120i].

    Τότε

    Α-x⁢I2=[1-x20i-x]⇒PA⁢(x)=(1-x)⁢(i-x)=x2-(1+i)⁢x+i.

    Οι ιδιοτιμές του A είναι οι 1 και i. Έστω ϕ:ℂ2→ℂ2, ϕ⁢(a,b)=(a,2⁢a+i⁢b). Αν B είναι η κανονική βάση του ℂ2, τότε A=AB,Bϕ. Επομένως

    Pϕ⁢(x)=x2-(1+i)⁢x+i.
  2. 2.

    Έστω

    A=[102023023].

    Τότε

    Α-x⁢I3=[1-x0202-x3023-x]

    και det⁡(A-x⁢I3)=

    (1-x)⁢((2-x)⁢(3-x)-6)=-x⁢(x-1)⁢(x-5)=-x3+6⁢x2-6⁢x.
  3. 3.

    Έστω Α∈ℳn⁢(k). Τότε

    det⁡(A-x⁢In)=det⁡(A-x⁢In)T=det⁡(AT-x⁢In).

    Δείξαμε λοιπόν ότι

    Οι πίνακες A και AT έχουν το ίδιο χαρακτηριστικό πολυώνυμο.

  4. 4.

    Είδαμε ότι ο σταθερός όρος του πολυωνύμου PA⁢(x) είναι ίσος με την ορίζουσα του Α. Σύμφωνα με την Πρόταση 5.1.7, το 0 είναι ιδιοτιμή του Α αν και μόνο αν det⁡A=0. Επομένως

    O πίνακας A είναι αντιστρέψιμος αν και μόνο αν ο σταθερός όρος του PA⁢(x) είναι διάφορος του μηδενός.

Πόσες ιδιοτιμές έχει μία γραμμική συνάρτηση ϕ:kn→kn ; Αντίστοιχα πόσες ιδιοτιμές έχει ένας n×n πίνακας A; Σύμφωνα με τα όσα έχουμε εξετάσει προηγουμένως o αριθμός των ιδιοτιμών του ϕ και αντίστοιχα του A είναι ακριβώς o o αριθμός των ριζών του Pϕ⁢(x) και αντίστοιχα του PA⁢(x)). Πόσες λοιπόν ρίζες έχει ένα πολυώνυμο βαθμού n; H απάντηση εξαρτάται από το σώμα επάνω από το οποίο δουλεύουμε και δίνεται από την επόμενη πρόταση.

Πρόταση 5.2.2.

Αν ϕ:ℂn→ℂn είναι μία γραμμική συνάρτηση, τότε η ϕ έχει ακριβώς n μιγαδικές ιδιοτιμές. Αν ϕ:ℝn→ℝn είναι μία γραμμική συνάρτηση, τότε η ϕ έχει το πολύ n πραγματικές ιδιοτιμές. Αντίστοιχα αν Α∈ℳn⁢(ℂ) τότε ο Α έχει ακριβώς n μιγαδικές ιδιοτιμές. Αν Α∈ℳn⁢(ℝ) τότε ο Α έχει το πολύ n πραγματικές ιδιοτιμές.

Απόδειξη.

Για την απόδειξη, πρέπει να αναφερθούμε στο Θεμελιώδες Θεώρημα της Άλγεβρας. Σύμφωνα με το θεώρημα αυτό, Θεμελιώδες Θεώρημα της Άλγεβρας ένα πολυώνυμο p⁢(x)∈ℂ⁢[x] βαθμού n με συντελεστές από το ℂ έχει ακριβώς n ρίζες στο ℂ, μετρημένες σύμφωνα με την πολλαπλότητά τους. Βεβαίως, αν p⁢(x)∈ℝ⁢[x], τότε p⁢(x)∈ℂ⁢[x]. Επομένως το p⁢(x) έχει ακριβώς deg⁡p⁢(x) μιγαδικές ρίζες και κάποιες από αυτές μπορεί να είναι πραγματικοί αριθμοί, ίσως όμως όχι όλες. Συνεπώς αν p⁢(x)∈ℝ⁢[x], τότε p⁢(x) έχει το πολύ deg⁡p⁢(x) πραγματικές ρίζες. Αφού PA⁢(x) έχει βαθμό n , το συμπέρασμα της πρότασης και για τις δύο περιπτώσεις προκύπτει από τις προηγούμενες παρατηρήσεις. ∎

Στα επόμενα παραδείγματα θα κάνουμε χρήση αυτής της πληροφορίας.

Παραδείγματα 5.2.3.

  1. 1.

    Έστω ότι Α∈ℳn⁢(ℂ) και ότι οι ιδιοτιμές του A είναι οι λ1,…,λn∈ℂ. Αφού λοιπόν o αρχικός συντελεστής του PA⁢(x) είναι το (-1)n, έπεται ότι

    PA⁢(x)=(-1)n⁢(x-λ1)⁢⋯⁢(x-λn).

    Αναπτύσσοντας αυτό το γινόμενο, προκύπτει ότι o συντελεστής του όρου xn-1 είναι ίσος με το

    (-1)n-1⁢(λ1+⋯+λn),

    ενώ o σταθερός όρος του PA⁢(x), o PA⁢(0), είναι ίσος με το γινόμενο

    (-1)n⁢(0-λ1)⁢⋯⁢(0-λn)=(-1)n⁢(-1)n⁢(λ1⁢⋯⁢λn)=λ1⁢⋯⁢λn.

    Προηγουμένως είδαμε ότι o συντελεστής του όρου xn-1 είναι (-1)n-1⁢Tr⁡(A) και ότι o σταθερός όρος του PA⁢(x) είναι ίσος με την ορίζουσα του A. Επομένως αποδείξαμε την επόμενη σχέση μεταξύ της ορίζουσας του A και της ιδιοτιμής του A.

    Αν Α∈ℳn⁢(ℂ) και οι ιδιοτιμές του A είναι οι λ1,…,λn, τότε

    Tr⁡(A)=λ1+⋯+λn και det⁡A=λ1⁢⋯⁢λn.

  2. 2.

    Δίνεται ο πίνακας

    A=[3-a3⁢bccaiic+i-1]

    μαζί με την πληροφορία ότι det⁡A=2 και ότι μία ιδιοτιμή του A είναι ίση με 2. Από την Πρόταση 5.2.2 γνωρίζουμε ότι o A έχει τρεις συνολικά μιγαδικές ιδιοτιμές. Θα υπολογίσουμε τις άλλες δύο.

    Έστω ότι λ1,λ2 είναι οι άλλες δύο ιδιοτιμές του Α. Τότε

    2=det⁡A=2⁢λ1⁢λ2,  2=Tr⁡A=2+λ1+λ2.

    Επομένως

    λ1=-λ2, και ⁢λ12=-1

    και οι άλλες δύο ιδιοτιμές του Α είναι ±i.

Όπως έχουμε δει, όμοιοι πίνακες έχουν τις ίδιες ιδιοτιμές, βλ. Παράδειγμα 5.1.8.2. Ισχύει όμως κάτι πιο ισχυρό, όπως δείχνει η επόμενη πρόταση.

Πρόταση 5.2.4.

Όμοιοι πίνακες έχουν το ίδιο χαρακτηριστικό πολυώνυμο.

Απόδειξη.

Έστω ότι A,B∈ℳn⁢(k) είναι δύο όμοιοι πίνακες. Θα αποδείξουμε ότι det⁡(A-x⁢In)=det⁡(B-x⁢In). Αφού A,B είναι όμοιοι, υπάρχει αντιστρέψιμος πίνακας S τέτοιος ώστε

B=S-1⁢A⁢S.

Θα δείξουμε ότι det⁡(Β-x⁢In)=det⁡(A-x⁢In). Πράγματι

det⁡(Β-x⁢In)=det⁡(S-1⁢A⁢S-x⁢S-1⁢S)=det⁡(S-1⁢A⁢S-S-1⁢x⁢In⁢S)=
=det(S-1(A-xIn)S)=det(S-1|det(A-xIn)det(S)=det(A-xIn).

Αποδείξαμε λοιπόν ότι οι όμοιοι πίνακες A και B έχουν το ίδιο χαρακτηριστικό πολυώνυμο. ∎

Έστω ϕ:kn→kn μία γραμμική συνάρτηση. Οι πίνακες της ϕ ως προς διαφορετικές βάσεις του kn είναι όμοιοι, (Πόρισμα 4.3.4). Σύμφωνα με την προηγούμενη πρόταση, μπορούμε να επιλέξουμε οποιαδήποτε βάση του kn και θα βρούμε ακριβώς το ίδιο πολυώνυμο Pϕ⁢(x). O επόμενος αλγόριθμος συγκεντρώνει τα βήματα για τον υπολογισμό των ιδιοτιμών και των ιδιοδιανυσμάτων της ϕ.

Αλγόριθμος 5.2.1 Αλγόριθμος υπολογισμού ιδιοχώρων της γραμμικής συνάρτησης ϕ:k𝐧→k𝐧.
Είσοδος: Μια γραμμική συνάρτηση ϕ:k𝐧→k𝐧
Έξοδος: Μια βάση για τους ιδιοχώρους της ϕ
Βήμα 1 Έστω A = [ ϕ ( e 1 ) ⋯ ϕ ( e n ) ] .
Βήμα 2 Βρίσκουμε το πολυώνυμο P A ( x ) = P ϕ ( x ) υπολογίζοντας την ορίζουσα d e t ( A − x I n ) .
Βήμα 3 Βρίσκουμε τις ιδιοτιμές του Α στον k, λύνοντας την εξίσωση P A ( x ) = 0 .
Βήμα 4 Για κάθε λ που είναι ιδιοτιμή του A , βρίσκουμε μια βάση για το μηδενοχώρο του A − λ I n : V λ ( ϕ ) = n u l l ( A − λ I n ) .

Στα παραδείγματα που ακολουθούν θα δούμε τον αλγόριθμο στη πράξη.

Παραδείγματα 5.2.5.

  1. 1.

    Έστω ϕ:ℂ2→ℂ2, ϕ⁢(a,b)=(a+2⁢b,i⁢b). Υπολογίζουμε τις εικόνες ϕ⁢(e1) και ϕ⁢(e2). Αφού ϕ⁢(e1)=(1,i) και ϕ⁢(e2)=(2,i), αν B είναι η κανονική βάση του ℂ2 και Α=ΑB,Bϕ, τότε

    A=[120i].

    Όπως είδαμε στο Παράδειγμα 5.2.1.1,

    PA⁢(x)=(1-x)⁢(i-x)

    και A έχει τις ιδιοτιμές 1 και i. Έτσι

    V1⁢(A)=null⁡(A-I2)=null⁡[020i-1]=null⁡[0100]={t⁢(1,0):t∈ℂ}

    ενώ

    Vi⁢(A)=null⁡(A-i⁢I2)=null⁡[1-i200]=null⁡[11+i00]={t⁢(-1-i,1):t∈ℂ}.

    Παρατηρούμε ότι αν v1=(1,0), v2=(-1-i,1), τότε D=(v1,v2) είναι διατεταγμένη βάση του ℂ2 από ιδιοδιανύσματα της ϕ και άρα

    ΑD,Dϕ=[100i].

    Αφού ο πίνακας ΑD,Dϕ είναι διαγώνιος, η γραμμική συνάρτηση ϕ είναι διαγωνιοποιήσιμη.

  2. 2.

    Θα προσδιορίσουμε τους ιδιοχώρους της γραμμικής συνάρτησης

    ϕ:ℝ2→ℝ2,(x,y)↦(2⁢x+y,x+2⁢y).

    Ο πίνακας του ϕ ως προς τη συνήθη βάση του ℝ2 είναι ο

    A=[2112].

    Αφού A=AT, o A έχει ακριβώς δύο πραγματικές ιδιοτιμές, βλ. Παράδειγμα 5.1.5.11. Πράγματι, υπολογίζουμε και παραγοντοποιούμε το χαρακτηριστικό πολυώνυμο της ϕ:

    Pϕ⁢(x)=det⁡[2-x112-x]=(2-x)2-1=x2-4⁢x+3=(x-1)⁢(x-3).

    Επομένως οι ιδιοτιμές του A είναι το 1 και το 3. Φέρνουμε τον πίνακα A-1⁢I2 σε ελαττωμένη κλιμακωτή μορφή γραμμών:

    [1111]⟶[1100].

    O ιδιοχώρος του ϕ που αντιστοιχεί στην ιδιοτιμή 1, είναι o

    V1⁢(ϕ)={(-x2,x2):x2∈k}=S⁢({(-1,1)})=S⁢({v}),

    όπου v=(-1,1).

    Για τον χώρο V3⁢(ϕ), φέρνουμε τον πίνακα A-3⁢I2 σε ελαττωμένη κλιμακωτή μορφή γραμμών

    [-111-1]⟶[1-100],

    άρα

    V3⁢(ϕ)={(x2,x2):x2∈k}=S⁢({(1,1)})=S⁢({u}),

    όπου u=(1,1).

    Τα ιδιοδιανύσματα {v,u} είναι βάση για τον ℝ2 και o Α είναι διαγωνιοποιήσιμος και όμοιος με τον

    [1003].
  3. 3.

    Θα προσδιορίσουμε τις ιδιοτιμές και τους ιδιοχώρους του πίνακα

    A=[1110]∈ℳ2⁢(ℝ).

    Το χαρακτηριστικό πολυώνυμο του A είναι το

    |1-x11-x|=-x+x2-1=(x-1+52)⁢(x-1-52).

    Οι ιδιοτιμές του Α είναι λ1=1+52 και λ2=1-52.

    Για να βρούμε τον ιδιοχώρο Vλ1 λύνουμε το σύστημα (Α-λ1⁢Ι2)⁢X=𝐎. H ελαττωμένη κλιμακωτή μορφή γραμμών του Α-λ1⁢Ι2 είναι:

    [1-1-5200]

    και άρα

    Vλ1⁢(A)={(1+52⁢x2,x2):x2∈ℝ}=S⁢({(1+52,1)})
    =S⁢({(1+5,2)}).

    Αντίστοιχα

    Vλ2⁢(A)=S⁢({(1-52,1)})=S⁢({(1-5,2)}).

    Ο πίνακας Α είναι διαγωνιοποιήσιμος και όμοιος με τον πίνακα

    [λ100λ2].
  4. 4.

    Ο πίνακας

    A=[0-110]∈ℳ2⁢(ℝ)

    έχει χαρακτηριστικό πολυώνυμο το

    PA⁢(x)=|-x-11-x|=x2+1 .

    Το πολυώνυμο PA⁢(x) δεν έχει ρίζες στον ℝ και ο Α δεν είναι διαγωνιοποιήσιμος στον ℝ. Τα ±i∈ℂ είναι ρίζες του PA⁢(x) στον ℂ και άρα ο Α έχει ιδιοτιμές τους μιγαδικούς αριθμούς ±i και είναι διαγωνιοποιήσιμος στον ℂ. Οι αντίστοιχοι ιδιοχώροι είναι

    Vi⁢(Α)={(i⁢x2,x2):x2∈ℂ}=S⁢({(i,1)}),

    και

    V-i⁢(Α)={(i⁢x2,x2):x2∈ℂ}=S⁢({(-i,1)}).
  5. 5.

    Έστω ϕ:ℝ3→ℝ3, ϕ⁢(a,b,c)=(a+2⁢c,2⁢b+3⁢c,2⁢b+3⁢c), B η κανονική βάση του ℝ3 και Α=ΑB,Bϕ. Τότε

    A=[102023023]

    και όπως είδαμε στο Παράδειγμα 5.2.1.2, PA⁢(x)=-x⁢(x-1)⁢(x-5). Επομένως οι ιδιοτιμές του Α είναι 0, 1, 5.

    V0⁢(A)=null⁡(A-0⁢I3)=null⁡A=null⁡[102023023]=null⁡[1020132000].

    Επομένως

    V0⁢(A)={t⁢(-2,-32,1):t∈ℝ}.

    Επίσης

    V1⁢(A)=null⁡(A-1⁢I3)=null⁡[002013022]=null⁡[010001000].

    Επομένως

    V1⁢(A)={t⁢(1,0,0):t∈ℝ}.

    Τέλος

    V5⁢(A)=null⁡(A-5⁢I3)=null⁡[-4020-3302-2]=null⁡[10-1201-1000]

    και

    V5⁢(A)={t⁢(12,1,1):t∈ℝ}.

    Παρατηρούμε ότι an

    v1=(-2,-32,1),v2=(1,0,0),v3=(12,1,1)

    τότε D=(v1,v2,v3) είναι διατεταγμένη βάση του ℝ3 από ιδιοδιανύσματα της ϕ και

    ΑD,Dϕ=[000010005].

    Αφού ο πίνακας ΑD,Dϕ είναι διαγώνιος, η ϕ είναι διαγωνιοποιήσιμη.

  6. 6.

    Έστω

    Α=[401-210-201]∈M3⁢(ℝ).

    Το χαρακτηριστικό πολυώνυμο PA⁢(x) του Α είναι το

    PA⁢(x)=|4-x01-21-x0-201-x|=(1-x)⁢[(4-x)⁢(1-x)+2]⇒
    ⇒PA(x)=(1-x)(x-2)(x-3).

    Επομένως οι ιδιοτιμές του A είναι είναι οι

    λ1=1,λ2=2,λ3=3.

    Αφού ο πίνακας A έχει τρεις διακεκριμένες ιδιοτιμές είναι διαγωνιοποιήσιμος. Πράγματι, οι ιδιοχώροι του A είναι οι

    V1={(0,κ,0):κ∈ℝ}=S⁢({(0,1,0)}),V2={(-12⁢κ,κ,κ):κ∈ℝ}=S⁢({(-12,1,1)}),V3={(-κ,κ,κ):κ∈ℝ}=S⁢({(-1,1,1)}).

    Τα ιδιοδιανύσματα

    ε1=(0,1,0),ε2=(-12,1,1),ε3=(-1,1,1)

    αποτελούν βάση του ℝ3. Έστω D=(ε1,ε2,ε3). Ο A είναι όμοιος με τον πίνακα

    Ε=[100020003],

    και

    Ε=S-1⁢A⁢S

    όπου

    S=SB←D=[0-1/2-1111011].

Έστω λ ιδιοτιμή του A∈ℳn⁢(k). Τότε PA⁢(λ)=0 και PA⁢(x)=(x-λ)⁢q1⁢(x) για κάποιο q1⁢(x)∈k⁢[x]. Αν q1⁢(λ)=0, τότε q1⁢(x)=(x-λ)⁢q2⁢(x) για κάποιο q2⁢(x)∈k⁢[x] και PA⁢(x)=(x-λ)2⁢q2⁢(x). Συνεχίζοντας καταυτόν τον τρόπο, καταλήγουμε ότι υπάρχει ένας φυσικός αριθμός mλ∈ℕ τέτοιος ώστε

PA⁢(x)=(x-λ)mλ⁢q⁢(x), όπου ⁢q⁢(x)∈k⁢[x], και ⁢q⁢(λ)≠0 .

Σημειώνουμε ότι 1≤mλ≤n. O φυσικός αυτός αριθμός mλ λέγεται αλγεβρική πολλαπλότητα (algebraic multiplicity) της ιδιοτιμής λ. Αν λοιπόν λ1,…,λk είναι οι διακεκριμένες ιδιοτιμές του A στο k, τότε υπάρχει h⁢(x)∈k⁢[x] τέτοιο ώστε

deg⁡PA⁢(x)=(x-λ1)mλ1⁢⋯⁢(x-λk)mλk⁢h⁢(x)

και άρα

mλ1+⋯+mλk≤n.

Παρατηρούμε επίσης ότι αφού Vλ⁢(A)≠{𝟎}, mλnλ έπεται ότι

nλ:=dimk⁢Vλ⁢(A)≥1 .

O φυσικός αριθμός nλ λέγεται η γεωμετρική πολλαπλότητα (geometric multiplicity) της λ . Με χρήση της μεθόδου της μαθηματικής επαγωγής αποδεικνύεται ότι

Θεώρημα 5.2.6.

Αν το λ είναι ιδιοτιμή του A, τότε mλ≥nλ≥1, δηλ.

 η αλγεβρική πολλαπλότητα του λ ≥ τη γεωμετρική πολλαπλότητα του λ.

Για την απόδειξη της παραπάνω ανισότητας παραπέμπουμε στο σύγγραμμα [2, Πρόταση 5.2.5].

Παράδειγμα 5.2.7.

Έστω

A=[40000041000045000005].

Τότε PA⁢(x)=(x-2)3⁢(x-3)2 και επομένως η αλγεβρική πολλαπλότητα του 4 είναι το 3 ενώ η αλγεβρική πολλαπλότητα του 5 είναι το 2. O αναγνώστης καλείται να διαπιστώσει ότι η γεωμετρική πολλαπλότητα του 4 είναι ίση με το 2 ενώ η γεωμετρική πολλαπλότητα του 5 είναι το 2.

Ασκήσεις Ενότητας 5.2

  1. 1.

    Δίνεται ο πίνακας

    A=[ab0c1-a2⁢iic+11].

    Αν det⁡A=8 και μία ιδιοτιμή του A είναι το 2 να βρεθούν οι υπόλοιπες ιδιοτιμές.

  2. 2.

    Να βρείτε τους ιδιοχώρους του πίνακα

    Α=[21001-1024]∈ℳ3⁢(ℝ).
  3. 3.

    Να αποδείξετε ότι όμοιοι πίνακες έχουν το ίδιο ίχνος και την ίδια ορίζουσα χρησιμοποιώντας την Πρόταση 5.2.4.

  4. 4.

    Να βρείτε τις αλγεβρικές και γεωμετρικές πολλαπλότητες των ιδιοτιμών του πίνακα

    A=[1120001100001000001300031].

5.3 Διαγωνιοποιήσιμοι Πίνακες

Θυμίζουμε ότι o n×n πίνακας A είναι διαγωνιοποιήσιμος αν o A είναι όμοιος με διαγώνιο πίνακα. Σε αυτήν την ενότητα θα εξερευνήσουμε πότε μία γραμμική συνάρτηση ϕ:kn→kn είναι διαγωνιοποιήσιμη. Σύμφωνα με το Παράδειγμα 5.1.6.3, η ϕ είναι διαγωνιοποιήσιμη αν και μόνο αν o kn έχει μία διατεταγμένη βάση D από ιδιοδιανύσματα της ϕ. Τονίζουμε ότι αν το ιδιοδιάνυσμα vi κατέχει την i θέση στο διατεταγμένο σύνολο D, τότε το λi εμφανίζεται στην i γραμμή και στήλη του διαγώνιου πίνακα ΑD,Dϕ.

Παραδείγματα 5.3.1.

  1. 1.

    Έστω ότι ϕ:ℂ2→ℂ2, ϕ⁢(a,b)=(a+2⁢b,i⁢b). Όπως είδαμε στο Παράδειγμα 5.2.5.1, η ϕ είναι διαγωνιοποιήσιμη με ιδιοτιμές 1 και i και αντίστοιχα ιδιοδιανύσματα v1 και v2. Έστω τώρα D1=(v1,v2) και D2=(v2,v1). Τότε

    ΑD1,D1ϕ=[100i], ενώ ⁢ΑD2,D2ϕ=[i001].
  2. 2.

    Έστω ϕ:ℝ3→ℝ3, μία διαγωνιοποιήσιμη γραμμική συνάρτηση με ιδιοτιμές λ1, λ2, λ3 και αντίστοιχα ιδιοδιανύσματα v1, v2, v3. Αν D1=(v1,v2,v3) και D2=(v3,v1,v2), τότε

    ΑD1,D1ϕ=[λ1000λ2000λ3], ενώ ⁢ΑD2,D2ϕ=[λ3000λ1000λ2].

H επόμενη πρόταση αποδεικνύει ότι δεν υπάρχει σχέση γραμμικής εξάρτησης μεταξύ μη μηδενικών διανυσμάτων από διαφορετικούς ιδιοχώρους.

Πρόταση 5.3.2.

Έστω ϕ:kn→kn μία γραμμική συνάρτηση. Αν τα w1, w2,…,wκ είναι ιδιοδιανύσματα του V που αντιστοιχούν σε διακεκριμένες ιδιοτιμές λ1,λ2,…,λκ, τότε τα w1,w2,…,wκ είναι γραμμικά ανεξάρτητα.

Απόδειξη.

Θα δώσουμε τη βασική ιδέα της απόδειξης, αποδεικνύοντας αναλυτικά την περίπτωση δύο ιδιοδιανυσμάτων δύο διαφορετικών ιδιοτιμών και σημειώνοντας ότι η γενική περίπτωση γίνεται ανάλογα. Έστω ότι

μ1⁢w1+μ2⁢w2=𝟎, (5.3.2.1)

όπου μ1,μ2∈k. Θα αποδείξουμε ότι μ1=μ2=0. Πράγματι,

ϕ⁢(μ1⁢w1+μ2⁢w2)=ϕ⁢(𝟎)=ϕ⁢(𝟎)⇒μ1⁢ϕ⁢(w1)+μ2⁢ϕ⁢(w2)=𝟎⇒
μ1⁢λ1⁢w1+μ2⁢λ2⁢w2=0. (5.3.2.2)

Πολλαπλασιάζουμε τη Σχέση (5.3.2.1) με λ2 και την αφαιρούμε από τη Σχέση (5.3.2.2). Έτσι

μ1⁢(λ1-λ2)⁢w1=0 .

Όμως w1 είναι ιδιοδιάνυσμα και άρα w1≠𝟎. Έπεται ότι μ1⁢(λ1-λ2)= 0. Άρα μ1=0 αφού λ1≠λ2. Αντικαθιστούμε στη Σχέση (5.3.2.1) και με τον ίδιο τρόπο καταλήγουμε ότι μ2=0. ∎


Άμεση συνέπεια της προηγούμενης πρότασης είναι το εξής θεώρημα.

Θεώρημα 5.3.3.

Έστω A∈ℳn⁢(k) και λ1,…,λk∈k οι διακεκεριμένες ιδιοτιμές του A στο k. O Α είναι διαγωνιοποιήσιμος αν και μόνο αν

dimk⁢Vλ1⁢(A)+⋯+dimk⁢Vλk⁢(A)=n.

Σημειώνουμε δύο ειδικές περιπτώσεις του Θεωρήματος 5.3.3.

Πόρισμα 5.3.4.

Έστω ότι A∈ℳn⁢(k).

  1. i)

    Αν o A έχει n διαφορετικές ιδιοτιμές, τότε ο Α είναι διαγωνιοποιήσιμος.

  2. ii)

    Αν o A είναι διαγωνιοποιήσιμος, τότε η γεωμετρική πολλαπλότητα κάθε ιδιοτιμής του A είναι ίση με την αλγεβρική πολλαπλότητα της ιδιοτιμης.

Απόδειξη.

Από το Θεώρημα 5.2.6 προκύπτει ότι

k≤nλ1+⋯+mλk=mλ1+⋯+mλk≤n.

Επομένως αν o A είναι διαγωνιοποιήσιμος τότε

mλ1=nλ1,…,mλk=nλk,

ενώ αν n=k, τότε nλi=mλi=1, για i=1,…,n και o A είναι διαγωνιοποιήσιμος. ∎

Στα παραδείγματα που ακολουθούν θα εφαρμόσουμε τα παραπάνω.

Παραδείγματα 5.3.5.

  1. 1.

    Θεωρούμε τη γραμμική συνάρτηση

    ϕ:ℝ3→ℝ3,(x,y,z)↦4⁢(x,y,z)

    του ℝ-διανυσματικού χώρου ℝ3. Ο πίνακας της ϕ ως προς τη συνήθη βάση του ℝ3 είναι ο

    A=[400040004]⁢,

    και το χαρακτηριστικό πολυώνυμο είναι το PΑ⁢(x)=(4-x)3. Έτσι η μόνη ιδιοτιμή της ϕ είναι το 4 και αφού A-4⁢I3=𝟎 έπεται ότι V4⁢(ϕ)=null⁡(A-4⁢I3)=ℝ3.

  2. 2.

    Έστω

    Α=[401040004].

    Τότε PA⁢(x)=(4-x)3 και η μόνη ιδιοτιμή του Α είναι το 4. Επίσης

    Α-4⁢Ι3=[001000000]

    και V4⁢(A)={(s,t,0):s,t∈ℝ}=S⁢({e1,e2}). O πίνακας Α δεν είναι διαγωνιοποιήσιμος γιατί δεν υπάρχει βάση του ℝ3 που να αποτελείται από ιδιοδιανύσματα του Α.

  3. 3.

    Θα υπολογίσουμε τις ιδιοτιμές και τα ιδιοδιανύσματα της γραμμικής συνάρτησης f:ℝ3→ℝ3,f⁢(a,b,c)=(3⁢a+2⁢b+2⁢c,a+4⁢b+2⁢c,-a-2⁢b). Έστω Β η κανονική βάση του ℝ3 και έστω

    A=ΑΒ,Βf=[322142-1-20]⁢.

    Το χαρακτηριστικό πολυώνυμο PA⁢(x) είναι το

    |3-x  22  14-x2-1-2-x|=12-16⁢x+7⁢x2-x3=(2-x)2⁢(3-x),

    και οι ιδιοτιμές του A είναι το 2 και το 3. Φέρνουμε τον πίνακα A-2⁢I3 σε ελαττωμένη κλιμακωτή μορφή γραμμών:

    A-2⁢I3=[122122-1-40]⟶[122000000]⁢,

    επομένως το σύνολο των λύσεων του συστήματος (A-2⁢I3)⁢X=𝟎 είναι το

    V2⁢(Α)={(-2⁢x2-2⁢x3,x2,x3):x2,x3∈ℝ}=
    {t⁢(-2,1,0)+s⁢(-2,0,1):t,s∈ℝ}.

    Όμοια φέρνουμε τον πίνακα A-3⁢I3 σε ελαττωμένη κλιμακωτή μορφή γραμμών:

    A-3⁢I3=[022112-1-2-3]⟶[101011000]⁢,

    και ο ιδιοχώρος που αντιστοιχεί στην ιδιοτιμή 3 είναι ο

    V3⁢(Α)=S⁢({(-1,-1,1)}).

    Αφού

    dim⁢(V2⁢(A))+dim⁢(V3⁢(A))=3

    έπεται ότι

    V2⁢(A)+V3⁢(A)=ℝ3

    και ότι η f είναι διαγωνιοποιήσιμη. Θεωρούμε τα ιδιοδιανύσματα v1=(-2,1,0),v2=(-2,0,1) (για την ιδιοτιμή 2) και v3=(-1,-1,1) (για την ιδιοτιμή 3) και τις διατεταγμένες βάσεις D=(v1,v2,v3), E=(v2,v3,v1) του ℝ3. Τότε

    AD,Df=[200020003]⁢ ενώ ⁢AE,Ef=[200030002].

    Οι τρεις πίνακες ΑB,Bf, AD,Df, AE,Ef είναι όμοιοι ανά ζεύγη. O αναγνώστης καλείται να βρει τους πίνακες μετάβασης.

  4. 4.

    Έστω o πίνακας

    A=[20000-2100-5200002].

    Τότε PA(x)==(2-x)2(x2+1). Παρατηρούμε ότι η μοναδική ιδιοτιμή του Α στον ℝ είναι το 2, ενώ στον ℂ οι ιδιοτιμές του Α είναι το 2 και το ±i. Επομένως o A δεν είναι διαγωνιοποιήσιμος στον ℝ. Εθετάζουμε στη συνέχεια αν o A είναι διαγωνιοποιήσιμος στον ℂ.

    Για να βρούμε μία βάση για τον ιδιοχώρο V2⁢(A), φέρνουμε τον πίνακα A-2⁢I4 σε ελαττωμένη κλιμακωτή μορφή γραμμών:

    [00000-4100-5000000]⟶[0100001000000000]⁢.

    Επομένως

    V2⁢(Α)=null⁡(A-2⁢I4)={(x1,0,0,x4):x1,x4∈ℂ}=S⁢({(1,0,0,0),(0,0,0,1)}).

    Για να βρούμε τον ιδιοχώρο Vi⁢(A) παρατηρούμε ότι μετά από στοιχειώδεις πράξεις γραμμών, η ελαττωμένη κλιμακωτή μορφή γραμμών του A-i⁢I4 είναι

    [1000012-i5000010000]

    και άρα

    Vi⁢(Α)={x3⁢(0,-2+i5,1,0):t∈ℂ}={k⁢(0,-2+i,5,0):k∈ℂ}.

    Αντίστοιχα, βρίσκουμε ότι o ιδιοχώρος για την ιδιοτιμή -i είναι

    V-i⁢(Α)={t⁢(0,-2-i,5,0):t∈ℂ}.

    Ο Α είναι διαγωνιοποιήσιμος στον C, αφού τα ιδιοδιανύσματα

    (1,0,0,0),(0,0,0,1),(0,-2+i,5,0),(0,-2-i,5,0)

    είναι γραμμικά ανεξάρτητα και αποτελούν βάση για τον ℂ4. Ο πίνακας Α στον ℂ είναι όμοιος με τον πίνακα

    [2000020000i0000-i].

Παραθέτουμε το επόμενο θεώρημα χωρίς απόδειξη, τονίζουμε όμως ότι πρόκειται για ένα από τα σημαντικότερα θεωρήματα της Γραμμικής Άλγεβρας. To θεώρημα αυτό δείχνει ότι μπορούμε να συνδυάσουμε γραμμικά τις δυνάμεις του A για να προκύψει o μηδενικός πίνακας. Για την απόδειξη, παραπέμπουμε τον αναγνώστη στο σύγγραμμα [2, Ενότητα 5.3].

Θεώρημα 5.3.6.

(Cayley-Hamilton) Έστω ότι A∈Mn×n⁢(k) και ότι

PA⁢(x)=(-1)n⁢xn+αn-1⁢xn-1+⋯+α1⁢x+α0,

είναι το χαρακτηριστικό πολυώνυμο του A. Τότε

(-1)n⁢Αn+αn-1⁢Αn-1+⋯+α1⁢Α+α0⁢Ιn

είναι o μηδενικός πίνακας, δηλ. PA⁢(A)=𝟎.

Στα επόμενα παραδείγματα θα εξετάσουμε διαφορετικές εφαρμογές του Θεωρήματος των Cayley-Hamilton.

Παραδείγματα 5.3.7.

  1. 1.

    Έστω

    A=[1102].

    Τότε PA⁢(x)=x2-3⁢x+2. Παρατηρούμε ότι

    A2-3⁢A+2⁢I2=[1304]-3⁢[1102]+[2002]=[0000]=𝟎,

    διαπιστώνουμε λοιπόν την ισχύ του Θεωρήματος των Cayley-Hamilton σε αυτό το συγκεκριμένο απλό παράδειγμα. Παρατηρούμε επίσης ότι

    A2-3⁢A+2⁢I2=𝟎⇒A2-3⁢A=-2⁢I2⇒(A-3⁢I2)⁢A=2⁢I2⇒12⁢(A-3⁢I2)⁢A=I2⇒A-1=12⁢(A-3⁢I2).
  2. 2.

    Έστω

    Α=[1 020-110 10].

    Τότε PΑ⁢(x)=-x3+2⁢x-1 και σύμφωνα με το Θεώρημα των Cayley-Hamilton

    -Α3+2⁢Α-Ι3=0.

    Αφού o σταθερός όρος του PA⁢(x) είναι διάφορος του μηδενός, o A είναι αντιστρέψιμος. Θα υπολογίσουμε τον A-1 χρησιμοποιώντας την προηγούμενη σχέση μεταξύ των δυνάμεων του A.

    -Α3+2⁢Α-Ι3=𝟎⇒-Α3+2⁢Α=Ι3⇒Α⁢(-Α2+2⁢Ι3)=Ι3.

    Επομένως

    A-1=-Α2+2⁢Ι3.
  3. 3.

    Έστω ότι

    PA⁢(x)=(-1)n⁢xn+αn-1⁢xn-1+⋯+α1⁢x+α0.

    Από το Θεώρημα των Cayley-Hamilton έπεται ότι

    PA⁢(A)=𝟎⇒(-1)n⁢An+αn-1⁢An-1+⋯+α1⁢A+α0⁢In=𝟎⇒An=(-1)n+1⁢αn-1⁢An-1+⋯+(-1)n+1⁢α1⁢A+(-1)n+1⁢α0⁢In.

    Επομένως ο An εκφράζεται ως γραμμικός συνδυασμός των In,A,…,An-1. Πολλαπλασιάζοντας τις δύο πλευρές της προηγούμενης ισότητας με Α και αντικαθιστώντας το An βλέπουμε ότι

    An+1=(αn-12+αn-1⁢αn-2)⁢An-1+⋯+αn-1⁢α0⁢Ιn.

    Άρα και o An+1 εκφράζεται ως γραμμικός συνδυασμός των In,A,…,An-1. Γενικότερα, όλες οι δυνάμεις Αm όπου m≥n μπορούν να γραφούν ως γραμμικοί συνδυασμοί των δυνάμεων In,A,…,An-1.

    To ίδιο ισχύει και για το A-1 όταν o A είναι αντιστρέψιμος. Πράγματι, ο A είναι αντιστρέψιμος αν και μόνο αν α0≠0. Σε αυτήν την περίπτωση,

    (-1)n⁢An+αn-1⁢An-1+⋯+α1⁢A=-α0⁢In.

    Βγάζοντας κοινό παράγοντα το Α (από το αριστερό μέρος της ισότητας) βρίσκουμε ότι

    A⁢((-1)n⁢An-1+αn-1⁢An-2+⋯+α1⁢In)=-α0⁢In,

    δηλ.

    A-1=-1α0⁢((-1)n⁢An-1+αn-1⁢An-2+⋯+α1⁢In).
  4. 4.

    Έστω

    Α=[00ii000-10].

    Το χαρακτηριστικό πολυώνυμο του πίνακα Α είναι

    det⁡[-x0ii-x00-1-x]=(-x)⁢(-x)2-i⋅i=-x3+1.

    Σύμφωνα με το θεώρημα των Cayley-Hamilton

    -Α3+Ι3=𝟎

    και άρα Α3=I3. Αν m∈ℕ, τότε o m μπορεί να γραφεί μοναδικά ως m=3⁢κ+r, όπου r=0,1,2 και κ∈ℕ. Επομένως για m∈ℕ έχουμε τις εξής περιπτώσεις:

    Am={Α3⁢κ=(Α3)κ=I3κ=I3Α3⁢κ+1=(Α3)κ⁢A=AΑ3⁢κ+2=(Α3)κ⁢A2=A2.

    Έτσι οι διακεκριμένες δυνάμεις του πίνακα Α είναι οι Ι3,Α,A2. Αφού A3=Ι3, προκύπτει ότι A⁢A2=Ι3, άρα A-1=A2.


O πίνακας J∈ℳk⁢(ℂ) λέγεται μπλοκ του Jordan (Jordan block) για το λ αν o J έχει την εξής μορφή:

J=[λ10⋯000λ1⋯0000λ⋯00⋮⋱⋮000⋯λ1000⋯0λ].

Ένα σημαντικό αποτέλεσμα αφορά την ταξινόμηση των πινάκων σε ανάλυση μπλοκ πινάκων του Jordan. To αναφέρουμε στη συνέχεια χωρίς απόδειξη.

Θεώρημα 5.3.8.

Έστω A∈ℳn⁢(ℂ) και έστω PA⁢(x)=(x-λ1)m1⁢⋯⁢(x-λs)ms, όπου λ1,…,λs∈ℂ είναι οι διακεκριμμένες ιδιοτιμές του A. Έστω επίσης ότι ni=dimℂ⁢Vλi⁢(A), για i=1,…,s. Υπάρχει αντιστρέψιμος πίνακας P τέτοιος ώστε

P-1⁢A⁢P=[J10⋯00J2⋯0⋮⋱000⋯Jr],

όπου, για j=1,…,r, o υποπίνακας Jj του J είναι μπλοκ του Jordan για κάποια από τις ιδιοτιμές του A. Ισχύουν τα εξής για την ιδιοτιμή λi (i=1,…,s):

  1. i)

    O αριθμός των μπλοκ του Jordan για την λi είναι ίσος με ni,

  2. ii)

    η λi εμφανίζεται ακριβώς mi φορές στη διαγώνιο των μπλοκ του Jordan για τη λi.

O πίνακας P-1⁢A⁢P είναι μοναδικός και λέγεται η κανονική μορφή του Jordan (Jordan canonical form) για τον A.

Στο επόμενο παράδειγμα εφαρμόζουμε το Θεώρημα 5.3.8 για έναν 3×3 πίνακα.

Παράδειγμα 5.3.9.

Έστω ότι PA⁢(x)=(x-2)3. Τότε o A είναι όμοιος με ακριβώς έναν από τους εξής πίνακες:

J1=[200020002],J2=[200021002],J3=[210021002].

Και στις τρεις περιπτώσεις, η αλγεβρική πολλαπλότητα του 2 είναι ίση με το 3. O πίνακας A είναι διαγωνιοποιήσιμος μόνο όταν η κανονική μορφή του Jordan του A είναι o J1, όταν δηλ. η γεωμετρική πολλαπλότητα του 2 είναι 3. O J1 αποτελείται από τρία μπλοκ του Jordan, όλα 1×1 πίνακες. Όταν η γεωμετρική πολλαπλότητα του 2 είναι δύο, τότε η κανονική μορφή του Jordan του A είναι o J2. O J2 αποτελείται από δύο μπλοκ του Jordan: το πρώτο είναι ένας 1×1 πίνακας, το δεύτερο είναι ένας 2×2 πίνακας. Τέλος, όταν η γεωμετρική πολλαπλότητα του 2 είναι ένα, τότε η κανονική μορφή του Jordan του A είναι o J3. O J3 αποτελείται από ένα 3×3 μπλοκ του Jordan.

Παραπέμπουμε τον αναγνώστη στο σύγγραμμα [6] για περισσότερα επί του θέματος αυτού.

Ασκήσεις Ενότητας 5.3

  1. 1.

    Έστω η γραμμική συνάρτηση ϕ:ℝ2→ℝ2 που περιστρέφει τα διανύσματα αριστερόστροφα κατά γωνία π. Είναι η ϕ διαγωνιοποιήσιμη;

  2. 2.

    Έστω η γραμμική συνάρτηση ϕ:ℝ3→ℝ3 που αντικατοπτρίζει τα διανύσματα ως προς το επίπεδο y=0. Είναι η ϕ διαγωνιοποιήσιμη;

  3. 3.

    Να υπολογίσετε για ποιές τιμές των a,b διαγωνιοποιείται ο πίνακας

    [a10b]∈ℳ2⁢(ℝ).
  4. 4.

    Να εξετάσετε αν ο πίνακας

    Α=[21001-1024]

    είναι διαγωνιοποιήσιμος.

  5. 5.

    Να εξετάσετε αν ο πίνακας

    Α=[1i-i2]

    είναι διαγωνιοποιήσιμος.

  6. 6.

    Έστω ότι το χαρακτηριστικό πολυώνυμο του πίνακα A είναι το

    x4-2⁢x3+x+4.

    Να υπολογίσετε τον αριθμό των γραμμών του A, τη βαθμίδα rank⁡(A) και να εκφράσετε τον A-1 ως γραμμικό συνδυασμό δυνάμεων του A.

  7. 7.

    Με εφαρμογή του Θεωρήματος των Cayley-Hamilton να υπολογίσετε τον αντίστροφο του πίνακα

    Α=[1-22-1013-4-3]

    και να γράψετε τον A4 ως γραμμικό συνδυασμό των πινάκων I3,A,A2.

  8. 8.

    Να βρείτε τις κανονικές μορφές του Jordan που μπορεί να έχει o πίνακας A∈ℳ5⁢(ℂ) όταν PA⁢(x)=(x-3)2⁢(x+1)3.

5.4 Ορθογώνιοι Πίνακες και Εσωτερικό Γινόμενο στον ℝn

Ο τετραγωνικός πίνακας A∈ℳn⁢(ℝ) λέγεται ορθογώνιος (orthogonal) αν

AT⁢A=In ή ισοδύναμα αν  A-1=AT.

Είναι φανερό ότι A είναι ορθογώνιος πίνακας αν και μόνο an AT είναι ορθογώνιος πίνακας, βλ. Άσκηση 5.4.2.

Παραδείγματα 5.4.1.

  1. 1.

    O In είναι ορθογώνιος πίνακας.

  2. 2.

    Έστω

    A=[1/3-2/601/31/6-1/21/31/61/2].

    O αναγνώστης καλείται να επιβεβαιώσει ότι A⁢AT=I3 και άρα o Α είναι ορθογώνιος πίνακας και

    Α-1=ΑΤ=[1/31/31/3-2/61/61/60-1/21/2].
  3. 3.

    Έστω o Α∈ℳn⁢(ℝ) ορθογώνιος πίνακας. Αφού AT⁢A=In έπεται ότι

    det⁡AT⁢A=det⁡In=1⇒det⁡AT⁢det⁡A=1⇒det⁡A⁢det⁡A=1⇒detA=±1 .

Μία σημαντική ιδιότητα των ορθογώνιων πινάκων είναι ότι οι ιδιοτιμές τους μπορούν να πάρουν δύο μόνο τιμές, τις ±1, βλ. Άσκηση 5.3.3. Για να μπορέσουμε, όμως, να αποδείξουμε το προηγούμενο, θα χρειαστεί να αναπτύξουμε τη θεωρία του εσωτερικού γινομένου στον ℝn και να περιγράψουμε το γινόμενο AT⁢A με τη χρήση του εσωτερικού γινομένου.

Το σύνηθες εσωτερικό γινόμενο (standard inner product) στον ℝn είναι η συνάρτηση

⟨,⟩:ℝn×ℝn→ℝ, όπου 
⟨(α1,…,αn),(β1,…,βn)⟩=α1⁢β1+⋯+αn⁢βn.

Για παράδειγμα, στον ℝ2

⟨(1,2),(3,-1)⟩=1⋅3+2⋅(-1)=1,⟨(1,2),(1,2)⟩=5 .

Ομοίως στον ℝn για κάθε n∈ℕ, έχουμε ότι:

⟨(α1,…,αn),(α1,…,αn)⟩=α12+⋯+αn2.

Έστω ότι v,u,w∈ℝn, κ,λ∈ℝ. Είναι εύκολο να δούμε ότι ικανοποιούνται οι ιδιότητες του επόμενου πίνακα.

Πίνακας 5.4.1: Ιδιότητες του εσωτερικού γινομένου ⟨,⟩:ℝn×ℝn→ℝ.
  1. i)

    ⟨v,v⟩≥0, ∀v∈ℝn και ⟨v,v⟩=0 αν και μόνο αν v=𝟎.

  2. ii)

    ⟨v,u⟩=⟨u,v⟩, ∀v,w∈ℝn

  3. iii)

    ⟨κ⁢v+λ⁢w,u⟩=κ⁢⟨v,u⟩+λ⁢⟨w,u⟩, ∀v,w∈ℝn, ∀κ,λ∈ℝ,

  4. iv)

    ⟨v,κ⁢u+λ⁢w⟩=κ⁢⟨v,u⟩+λ⁢⟨v,w⟩, ∀v,w∈ℝn, ∀κ,λ∈ℝ,

  5. v)

    ⟨v,𝟎⟩=0, ∀v∈ℝn.

Ως άμεση συνέπεια της πρώτης ιδιότητας του εσωτερικού γινομένου, προκύπτει ότι

 an ⁢⟨v,w⟩=0, για κάθε ⁢w∈ℝn, τότε ⁢v=0. (5.4.1.1)

Το μήκος (norm) του v norm συμβολίζεται ∥v∥∥v∥ και είναι ίσο με ⟨v,v⟩. To διάνυσμα v λέγεται κανονικό (normal) αν ∥v∥=1. Αν w∈ℝn είναι μη μηδενικό διάνυσμα, τότε o αναγνώστης μπορεί εύκολα να διαπιστώσει ότι το διάνυσμα

1∥w∥⁢w

είναι κανονικό και παράλληλο προς το w. Η απόσταση (distance) απόσταση δύο διανυσμάτων v,u∈ℝn είναι το μήκος του διανύσματος v-u και συμβολίζεται d⁢(v,u) δηλ. d⁢(v,u)=∥u-v∥ . Μπορεί να αποδειχθεί η επόμενη ανισότητα, που είναι γνωστή ως ανισότητα των Cauchy-Schwarz (Cauchy-Schwarz inequality):

-∥u∥⁢∥v∥≤⟨u,v⟩≤∥u∥⁢∥v∥.

Έτσι ισχύει ότι

-1≤⟨u,v⟩∥u∥⁢∥v∥≤1 .

Αναφέρουμε επίσης την τριγωνική ανισότητα (triangle inequality) που προκύπτει εύκολα από την ανισότητα των Cauchy-Schwarz:

∥v+u∥≤∥v∥+∥u∥.

Αν v, u≠0 τότε γωνία (angle) των v και u ορίζεται η γωνία θ για την οποία 0≤θ≤π και

cos⁡θ=⟨u,v⟩∥u∥⁢∥v∥.
Παράδειγμα 5.4.2.

Έστω v=(1,0,1), u=(-1,1,0). Τότε ∥u∥=∥u∥=2 και d⁢(u,v)=6.

cos⁡θ=⟨v,u⟩∥v∥⁢∥u∥=-12⁢2=-12

και άρα θ=π-π3=2⁢π3.

Σχήμα 5.5: Τα διανύσματα u και v στον ℝ3

Τα διανύσματα u,v δεν είναι κανονικά. Τα διανύσματα όμως

u′=12⁢(1,0,1),v′=12⁢(-1,1,0)

είναι κανονικά και σχηματίζουν την ίδια γωνία όπως τα u και v.

Δύο διανύσματα v, u λέγονται ορθογώνια (orthogonal) όταν

⟨v,u⟩=0 .

Είναι εύκολο να αποδείξει κανείς ότι όταν τα διανύσματα v, u είναι ορθογώνια, τότε ισχύει το παρακάτω συμπέρασμα, γνωστό ως το Πυθαγόρειο Θεώρημα, βλ. Άσκηση 5.4.2.

Αν u,v∈ℝn είναι ορθογώνια διανύσματα τότε ∥u+v∥2=∥u∥2+∥v∥2.

Αν D είναι διατεταγμένη βάση του ℝn, τότε η D λέγεται ορθογώνια βάση (orthogonal basis) του ℝn, αν τα διανύσματα της D είναι ανά δύο ορθογώνια. Αν η D είναι ορθογώνια βάση του ℝn και κάθε διάνυσμα της D είναι κανονικό, τότε η D λέγεται ορθοκανονική βάση (orthonormal basis) του ℝn .

Παράδειγμα 5.4.3.

Είναι φανερό ότι η κανονική βάση B=(e1,…,en) του ℝn είναι ορθοκανονική. Πράγματι, e1,…,en είναι κανονικά αφού ∥ei∥=1, για i=1,…,n και

⟨ei,ej⟩={0,για ⁢i≠j1,για ⁢i=j

δηλ. ei,ej είναι ανά δύο ορθογώνια για 1≤i<j≤n.

Όταν o V είναι διανυσματικός υποχώρος του ℝn γενικεύουμε τους προηγούμενους ορισμούς. Έτσι, αν D είναι μία διατεταγμένη βάση του V τότε η D λέγεται ορθογώνια βάση (orthogonal basis) του V αν τα διανύσματα της D είναι ανά δύο ορθογώνια και η D λέγεται ορθοκανονική βάση (orthonormal basis) του V αν D είναι ορθογώνια βάση και κάθε διάνυσμα της D είναι κανονικό. Σημειώνουμε την επόμενη παρατήρηση.

Κάθε ορθογώνια βάση του V γίνεται ορθοκανονική αν κανονικοποιήσουμε (normalize) τα στοιχεία της, δηλ. αν διαιρέσουμε τα στοιχεία της βάσης με το μήκος τους.

Παραδείγματα 5.4.4.

  1. 1.

    Έστω v1=(1,1), v2=(1,-1). Αφού ⟨v1,v2⟩=0 και v1,v2 είναι γραμμικά ανεξάρτητα έπεται ότι D={v1,v2} είναι ορθογώνια βάση του ℝ2. Το μήκος και των δύο διανυσμάτων της βάσης είναι ∥v1∥=∥v2∥=2. Η βάση D δεν είναι ορθοκανονική. Κανονικοποιούμε, θεωρώντας τα διανύσματα w1=1/2⁢(1,1), w2=1/2⁢(1,-1). Το σύνολο {w1,w2} είναι τώρα ορθοκανονική βάση του ℝ2.

  2. 2.

    Έστω v1=1/3⁢(1,1,1), v2=1/6⁢(-2,1,1), v3=1/2⁢(0,-1,1). Είναι εύκολο να δούμε ότι τα προηγούμενα διανύσματα είναι ανά δύο ορθογώνια. Θα δείξουμε ότι v1, v2, v3 είναι γραμμικά ανεξάρτητα. Πράγματι, ο πίνακας

    A=[1/3-2/601/31/6-1/21/31/61/2]

    είναι ορθογώνιος (Παράδειγμα 5.4.1), και επομένως αντιστρέψιμος. Άρα τα v1,v2,v3 είναι γραμμικά ανεξάρτητα, βλ. Πρόταση 3.3.5. Επίσης ∥v1∥=∥v2∥=∥v3∥=1. Έπεται ότι το σύνολο {v1,v2,v3} είναι ορθοκανονική βάση του ℝ3.

  3. 3.

    Έστω o υποχώρος V του ℝ3 που ορίζεται ως το σύνολο

    V={(x,y,z)∈ℝ3:x+y-z=0}.

    Αν v1=(1,0,1), v2=(0,1,1), τότε (v1,v2) είναι διατεταγμένη βάση για τον V αλλά όχι ορθογώνια βάση για τον V. Παρατηρούμε, όμως, ότι

    w=(-1,2,1)=-v1+2⁢v2∈V

    και ότι ⟨v1,w⟩=0, άρα {v1,w} είναι ορθογώνια βάση για τον V. Επομένως

    {12⁢v1,16⁢w}

    είναι ορθοκανονική βάση για τον V.

To επόμενο θεώρημα δείχνει τη χρησιμότητα της εύρεσης ορθογώνιας βάσης.

Θεώρημα 5.4.5.

Έστω ότι V υποχώρος του ℝn και ότι τα διανύσματα v1,…,vt∈V είναι ανά δύο ορθογώνια. Τότε τα διανύσματα v1,…,vt είναι γραμμικά ανεξάρτητα και

v=κ1⁢v1+⋯+κt⁢vt⟹κi=⟨v,vi⟩∥vi∥2⁢ για ⁢i=1,…,t.
Απόδειξη.

Έστω ότι

κ1⁢v1+⋯+κt⁢vt=𝟎 (5.4.5.1)

Θα δείξουμε ότι οι συντελεστές κ1,…,κt είναι 0. Θα πάρουμε το εσωτερικό γινόμενο με το διάνυσμα v1 και των δύο πλευρών της έκφρασης (5.4.5.1). Αφού

⟨v1,κ1⁢v1+⋯+κt⁢vt⟩=⟨v1,𝟎⟩=0

έπεται ότι

⟨v1,κ1⁢v1⟩+⋯+⟨v1,κt⁢vt⟩=0

και

κ1⁢⟨v1,v1⟩+⋯+κt⁢⟨v1,vt⟩=0,

άρα

κ1⁢1+κ2⁢0+⋯+κt⁢0=0⇒κ1=0 .

Είναι φανερό ότι μπορούμε να επαναλάβουμε τα προηγούμενα παίρνοντας το εσωτερικό γινόμενο με κάθε ένα από τα διανύσματα v1,…,vt και να δείξουμε ότι κ1=⋯=κt=0 . Επομένως τα διανύσματα v1,…,vt είναι γραμμικά ανεξάρτητα.

Έστω τώρα v=κ1⁢v1+⋯+κt⁢vt. Παίρνοντας το εσωτερικό γινόμενο με το vi προκύπτει ότι:

⟨v,vi⟩=⟨κ1⁢v1+⋯+κt⁢vt,vi⟩=κi⁢⟨vi,vi⟩⟹κi=⟨v,vi⟩⟨vi,vi⟩

από όπου προκύπτει και το ζητούμενο. ∎

Έστω Β=(e1,…,en) η κανονική βάση του ℝn. Το γινόμενο των πινάκων CB⁢(u)Τ⁢CB⁢(v) είναι ένας 1×1 πίνακας με μοναδικό στοιχείο το αποτέλεσμα του εσωτερικού γινομένου των u, v, δηλ.

CB⁢(u)Τ⁢CB⁢(v)=CB⁢(v)T⁢CB⁢(u)=[⟨v,u⟩].

Εξαιτίας αυτής της παρατήρησης και για συντομογραφία γράφουμε:

⟨v,u⟩=vT⁢u.

Γενικότερα, αν A∈ℳn⁢(ℝ) και Α=[CB⁢(v1)…CB⁢(vn)], τότε

ΑΤ⁢Α=[v1T⁢v1…v1T⁢vn⋮⋮vnT⁢v1…vnT⁢vn].

Προκύπτει λοιπόν το εξής κριτήριο:

Πρόταση 5.4.6.

To σύνολο D={v1,…,vn}⊂ℝn είναι ορθογώνια βάση του ℝn αν και μόνο αν o AT⁢A είναι διαγώνιος πίνακας, όπου

Α=[CB⁢(v1)⋯CB⁢(vn)].

To σύνολο D είναι ορθοκανονική βάση του ℝn αν και μόνο αν AT⁢A=In, δηλ. αν και μόνο αν ο Α είναι ορθογώνιος πίνακας.

H προηγούμενη πρόταση γενικεύεται ως εξής για το διανυσματικό υποχώρο V του ℝn.

Πρόταση 5.4.7.

Έστω B η κανονική βάση του ℝn, V ένας διανυσματικός υποχώρος του ℝn, με dimℝ⁢V=m. Ένα υποσύνολο D={v1,…,vm}⊂V είναι ορθογώνια βάση του V αν και μόνο αν o AT⁢A είναι διαγώνιος m×m πίνακας, όπου

Α=[CB⁢(v1)…CB⁢(vm)].

To σύνολο D είναι ορθοκανονική βάση του ℝn αν και μόνο αν AT⁢A=Im, δηλ. αν και μόνο αν ο Α είναι ορθογώνιος πίνακας.

Είναι πλέον φανερό ότι αν οι στήλες του A είναι ορθοκανονική βάση του ℝn, τότε και οι γραμμές του A είναι ορθοκανονική βάση του ℝn.

Παραδείγματα 5.4.8.

  1. 1.

    Στο Παράδειγμα 5.4.1 είδαμε ότι o

    A=[1/3-2/601/31/6-1/21/31/61/2]

    είναι ορθογώνιος πίνακας. Συνεπώς

    (13⁢(1,1,1),16⁢(-2,1,1),12⁢(0,-1,1))

    και

    ((13,26,0),(13,16,-12),(13,16,12)),

    είναι δύο ορθοκανονικές βάσεις του ℝ3.

  2. 2.

    Στο Παράδειγμα 5.4.4.3 είδαμε ότι

    (12⁢(1,0,1),16⁢(-1,2,1))

    είναι ορθοκανονική βάση για το επίπεδο

    V={(x,y,z)∈ℝ3:x+y-z=0}.

    Έστω B η κανονική βάση του ℝ3. O πίνακας A της Πρότασης 5.4.7 είναι o 3×2 πίνακας

    A=[1/2-1/602/61/21/6].

    O αναγνώστης καλείται να διαπιστώσει ότι AT⁢A=I2.

Πως γενικεύονται τα προηγούμενα όταν o V είναι ένας ℝ-διανυσματικός χώρος; Αν V είναι ένας ℝ-διανυσματικό χώρος, τότε λέμε ότι o V είναι Ευκλείδειος διανυσματικός χώρος (Euclidean vector space) αν μπορεί να ορισθεί ένα εσωτερικό γινόμενο (inner product) στον V , εσωτερικό γινόμενο δηλ. μία συνάρτηση ⟨,⟩:V×V→ℝ για την οποία ισχύουν ιδιότητες, αντίστοιχες με αυτές του εσωτερικού γινομένου που μελετήσαμε στον ℝn και τους υποχώρους του, βλ. Πίνακα 5.4.1. Στον επόμενο πίνακα έχουμε συγκεντρώσει αυτές τις ιδιότητες.

Πίνακας 5.4.2: Ιδιότητες του εσωτερικού γινομένου ⟨,⟩:V×V→ℝ.
  1. i.)

    ⟨v,v⟩≥0, ∀v∈V και ⟨v,v⟩=0 αν και μόνο αν v=𝟎.

  2. ii.)

    ⟨v,u⟩=⟨u,v⟩, ∀v,w∈V

  3. iii.)

    ⟨κ⁢v+λ⁢w,u⟩=κ⁢⟨v,u⟩+λ⁢⟨w,u⟩, ∀v,w∈V, ∀κ,λ∈ℝ,

  4. iv.)

    ⟨v,κ⁢u+λ⁢w⟩=κ⁢⟨v,u⟩+λ⁢⟨v,w⟩, ∀v,w∈V, ∀κ,λ∈ℝ,

  5. v.)

    ⟨v,𝟎⟩=0, ∀v∈V.

Παραδείγματα 5.4.9.

  1. 1.

    O ℝn με το σύνηθες εσωτερικό γινόμενο.

  2. 2.

    To σύνολο V των συνεχών πραγματικών συναρτήσεων με πεδίο ορισμού το κλειστό διάστημα [a,b] είναι Ευκλείδειος διανυσματικός χώρος με εσωτερικό γινόμενο

    ⟨,⟩:V×V→ℝ,⟨f,g⟩=∫abf(t)g(t)dt.

Όλα τα συμπεράσματα αυτής της ενότητας επεκτείνονται με το φυσικότερο τρόπο σε Ευκλείδειους διανυσματικούς χώρους. Παραπέμπουμε τον αναγνώστη στα συγγράμματα της βιβλιογραφίας για περισσότερες λεπτομέρειες.

Μπορούμε, όμως, γενικότερα να ορίσουμε το σύνηθες εσωτερικό γινόμενο στον ℂn και την έννοια του εσωτερικού γινομένου για τους ℂ-διανυσματικούς χώρους. Το σύνηθες εσωτερικό γινόμενο (standard inner product) στον ℂn ορίζεται ως εξής:

⟨(α1,…,αn),(β1,…,βn)⟩=α1⁢β1¯+⋯+αn⁢βn¯.

Λέμε ότι o ℂ-διανυσματικός χώρος V είναι Ερμητιανός διανυσματικός χώρος (Hermitian vector space) αν σε αυτόν ορίζεται ένα εσωτερικό ή ερμητιανό γινόμενο (inner product), δηλ.  μία συνάρτηση ⟨,⟩:V×V→ℂ για την οποία ισχύουν οι ιδιότητες του επόμενου πίνακα.

Πίνακας 5.4.3: Ιδιότητες του Ερμητιανού εσωτερικού γινομένου ⟨,⟩:V×V→ℂ.
  1. i.)

    ⟨v,v⟩∈ℝ+, ∀v∈V και ⟨v,v⟩=0 αν και μόνο αν v=𝟎.

  2. ii.)

    ⟨v,u⟩=⟨u,v⟩¯, ∀v,w∈V.

  3. iii.)

    ⟨κ⁢v+λ⁢w,u⟩=κ⁢⟨v,u⟩+λ⁢⟨w,u⟩, ∀v,w∈V, ∀κ,λ∈ℝ.

  4. iv.)

    ⟨v,κ⁢u+λ⁢w⟩=κ¯⁢⟨v,u⟩+λ¯⁢⟨v,w⟩, ∀v,w∈V, ∀κ,λ∈ℝ.

  5. v.)

    ⟨v,𝟎⟩=0, ∀v∈V.

Ένας πίνακας A∈ℳn⁢(ℂ) λέγεται ορθογώνιος (orthogonal) αν

A⁢(AT)¯=Ιn.

Όπως σε έναν Ευκλείδειο διανυσματικό χώρο, έτσι και σε έναν Ερμητιανό διανυσματικό χώρο έχουμε τις έννοιες μίας ορθογώνιας και ορθοκανονικής βάσης. Τονίζουμε ότι αν o A∈ℳn⁢(ℝ) είναι ορθογώνιος στον ℝ, τότε είναι ορθογώνιος και στον ℂ. Έτσι, σύνολα που είναι ορθογώνιες (ορθοκανονικές) βάσεις στον ℝn είναι ορθογώνιες (ορθοκανονικές) βάσεις στον ℂn. Είναι όμως απαραίτητο για την ισχύ της θεωρίας μας, να χρησιμοποιούμε τους συζυγείς αναστρόφους πίνακες, όταν δουλεύουμε με μιγαδικούς αριθμούς.

Παράδειγμα 5.4.10.

O αναγνώστης καλείται να επιβεβαιώσει ότι o πίνακας

A=12⁢[1ii1]

είναι ορθογώνιος και ότι οι στήλες του A δίνουν μία ορθοκανονική βάση του ℂ2.

Ασκήσεις Ενότητας 5.4

  1. 1.

    Έστω v=(1,2,3,4). Να βρείτε ∥v∥, ∥3⁢v∥, ∥15⁢v∥, d⁢(u,v).

  2. 2.

    Αν u,v∈ℝn είναι ορθογώνια, τότε από τις ιδιότητες του εσωτερικού γινομένου, να αποδείξετε το Πυθαγόρειο Θεώρημα, δηλ. ότι

    ∥u+v∥2=∥u∥2+∥v∥2.
  3. 3.

    Να αποδείξετε ότι τα v1=(-3,0,4), v2=(4,0,-3), v3=(0,5,0) δίνουν μία ορθογώνια βάση του ℝ3.

  4. 4.

    Αν Α είναι ορθογώνιος πίνακας τότε να αποδείξετε ότι AT είναι ορθογώνιος πίνακας.

  5. 5.

    Για τον πίνακα

    A=[1551]

    να βρείτε ένα διαγώνιο πίνακα D και ένα πίνακα P έτσι ώστε D=PT⁢A⁢P.

  6. 6.

    Να διαγωνιοποιήσετε ορθογώνια τον πίνακα

    [100032026].
  7. 7.

    Έστω V o χώρος των συνεχών πραγματικών συναρτήσεων στο διάστημα [-π,π] με εσωτερικό γινόμενο

    ⟨,⟩:V×V→ℝ,⟨f,g⟩=∫-ππf(t)g(t)dt.

    Να αποδείξετε ότι τα στοιχεία του (άπειρου) συνόλου

    {1,cosx,sinx,cos(2x),sin(2x),…,cos(nx),sin(nx),…,}

    είναι ανά δύο ορθογώνια. Να συμπεράνετε ότι dimℝ⁢(V)=∞.

  8. 8.

    Να βρείτε ορθογώνιους πίνακες A1∈ℳ2⁢(ℝ) (αν υπάρχουν) τέτοιοι ώστε

    A1=[1/21/2bd],A2=[1/4abd].
  9. 9.

    Να βρείτε μία ορθογώνια βάση στον ℂ3 που να περιέχει το διάνυσμα (1,0,i).

  10. 10.

    Έστω ότι {u1,…,un} βάση για τον υποχώρο U του ℝm και έστω ότι w∈ℝm. Να αποδείξετε ότι ⟨w,u⟩=0, για κάθε u∈U , αν και μόνο αν ⟨w,ui⟩=0, για i=1,…,n.

  11. 11.

    Έστω ότι h:ℝn→ℝn είναι γραμμική συνάρτηση. An

    ⟨h⁢(v),w⟩=0, για ⁢v,w∈ℝn,

    τότε h=0ℝn και h⁢(v)=𝟎, για v∈ℝn. Να συμπεράνετε ότι αν g:ℝn→ℝn είναι γραμμική συνάρτηση και

    ⟨g⁢(v),w⟩=⟨v,w⟩, για ⁢v,w∈ℝn,

    τότε g=idℝn και g⁢(v)=v, για v∈ℝn.

5.5 Φασματικό Θεώρημα

Έστω ότι ϕ:ℝn→ℝn είναι μία γραμμική συνάρτηση, Β είναι η κανονική βάση του ℝn , και Α=AB,Bϕ . Ορίζουμε την προσαρτημένη γραμμική συνάρτηση (adjoint linear map) του ϕ , να είναι η γραμμική συνάρτηση ϕ*ϕ*:ℝn→ℝn , με πίνακα

AB, Bϕ*=ΑΤ.
Παράδειγμα 5.5.1.

Έστω η γραμμική συνάρτηση ϕ:ℝ2→ℝ2, ϕ⁢(x,y)=(x+2⁢y,3⁢x-y) . Αν Β είναι η κανονική βάση του ℝ2 , τότε

AB,Bϕ=[123-1].

Επομένως

AB,Bϕ*=[132-1]⁢ και⁢ϕ*:ℝ2→ℝ2,ϕ*⁢(x,y)=(x+3⁢y,2⁢x-y).

Μία σημαντική ιδιότητα του εσωτερικού γινομένου αφορά τις γραμμικές συναρτήσεις και τις προσαρτημένες συναρτήσεις.

Θεώρημα 5.5.2.

Έστω ότι ϕ:ℝn→ℝn είναι γραμμική συνάρτηση. Τότε

⟨ϕ⁢(v),w⟩=⟨v,ϕ*⁢(w)⟩⁢𝑔𝑖𝑎⁢v,w∈ℝn.
Απόδειξη.

Θα δείξουμε ότι

ϕ⁢(v)T⁢w=vT⁢ϕ*⁢(w)

και ισοδύναμα ότι

CB⁢(ϕ⁢(v))T⁢CB⁢(w)=CB⁢(v)T⁢CB⁢(ϕ*⁢(w)).

Από την Πρόταση 4.1.4, έπεται ότι

CB⁢(ϕ⁢(v))=Α⁢CB⁢(v)⁢και⁢CB⁢(ϕ*⁢(w))=AT⁢CB⁢(w),

όπου A=AB,Bϕ . Άρα

(Α⁢CB⁢(v))Τ⁢CB⁢(w)=(CB⁢(v)T⁢AT)⁢CB⁢(w)=CB⁢(v)T⁢(AT⁢CB⁢(w))

και το ζητούμενο αποδείχθηκε.∎

Στη γλώσσα των πινάκων, μπορούμε να διατυπώσουμε το Θεώρημα 5.5.2 ως εξής:

Έστω ότι A∈ℳn⁢(ℝ) και v∈ℝn. Τότε ⟨ A⁢v,v⟩=⟨v,AT⁢v⟩.

Είναι φανερό ότι (ϕ*)*=ϕ . Άμεση λοιπόν συνέπεια του Θεωρήματος 5.5.2 είναι το επόμενο πόρισμα.

Πόρισμα 5.5.3.

Έστω ότι ϕ:ℝn→ℝn είναι γραμμική συνάρτηση. Τότε

⟨ϕ*⁢(v),w⟩=⟨v,ϕ⁢(w)⟩⁢𝑔𝑖𝑎⁢v,w∈ℝn.

Ειδικότερα

⟨ϕ*⁢ϕ⁢(v),w⟩=⟨ϕ⁢(v),ϕ⁢(w)⟩⁢𝑔𝑖𝑎⁢v,w∈ℝn.

Στη συνέχεια εξετάζουμε υπό ποιες συνθήκες μία γραμμική συνάρτηση ϕ:ℝn→ℝn έχει n ορθογώνια ιδιοδιανύσματα. Όταν αυτό συμβαίνει λέμε ότι η ϕ (και ο αντίστοιχος πίνακας ως προς οποιαδήποτε βάση του ℝn ) είναι ορθογώνια διαγωνιοποιήσιμη (orthogonal diagonalizable). Τα n ορθογώνια ιδιοδιανύσματα του ℝn καθορίζουν ένα νέο σύστημα συντεταγμένων που διευκολύνει τη μελέτη της συνάρτησης. Από όσα έχουμε εξετάσει νωρίτερα και το Παράδειγμα 5.1.6.11 προκύπτει το εξής συμπέρασμα.

Πρόταση 5.5.4.

Ένας πίνακας Α∈ℳn⁢(ℝ) είναι ορθογώνια διαγωνιοποιήσιμος αν υπάρχει ορθογώνιος πίνακας P∈ℳn⁢(ℝ) έτσι ώστε

PT⁢A⁢P=[λ10⋱0λn].

Τα λ1,…,λn στην προηγούμενη έκφραση είναι οι ιδιοτιμές του Α και είναι πραγματικοί αριθμοί. Οι στήλες του P αντιστοιχούν σε ιδιοδιανύσματα του Α.

Ο Α είναι λοιπόν ορθογώνια διαγωνιοποιήσιμος όταν υπάρχει μία ορθοκανονική βάση του ℝn που να αποτελείται από ιδιοδιανύσματα του Α . To επόμενο θεώρημα δίνει μία ικανή και αναγκαία συνθήκη για να είναι ένας πίνακας ορθογώνια διαγωνιοποιήσιμος. To θεώρημα είναι γνωστόv ως το Φασματικό θεώρημα (spectral theorem). Για την αναλυτική απόδειξη του θεωρήματος παραπέμπουμε στο σύγγραμμα [2 , Θεώρημα 7.2.18]

Θεώρημα 5.5.5 (Φασματικό Θεώρημα).

Ο πίνακας Α∈ℳn⁢(ℝ) είναι ορθογώνια διαγωνιοποιήσιμος αν και μόνο αν Α=ΑΤ.

Στη συνέχεια, στο Θεώρημα 5.5.6, θα αποδείξουμε μία ειδική περίπτωση του Φασματικού Θεωρήματος.

Θεώρημα 5.5.6.

Έστω ότι Α∈ℳ2⁢(ℝ) και A=AT . Τότε o A έχει δύο ορθογώνια ιδιοδιανύσματα.

Απόδειξη.

Έστω ότι Α∈ℳ2⁢(ℝ) και A=AT . Στο Παράδειγμα 5.1.6.11 είδαμε ότι οι ιδιοτιμές του Α είναι πραγματικοί αριθμοί. Έστω ότι λ1 είναι μία από τις ιδιοτιμές του Α και ότι X1=[a1b1]T είναι ιδιοδιάνυσμα του Α για την ιδιοτιμή λ1 . Θεωρούμε τον ορθογώνιο πίνακα

P=[a1-b1b1a1].

Θα δείξουμε ότι X2=[-b111]T είναι ιδιοδιάνυσμα του A . Έστω v1=(a1,b1),v2=(-b1,a1) και B η κανονική βάση του ℝ2 . Έτσι, X1=CB⁢(v1),X2=CB⁢(v2) . Για να δείξουμε ότι X2 είναι ιδιοδιάνυσμα του A αρκεί να δείξουμε ότι A⁢X2=t⁢X2 για κάποιον πραγματικό αριθμό t . Θέτουμε Y=A⁢X2 Αφού (v1,v2) είναι βάση του ℝ2 , υπάρχουν κ,t∈ℝ έτσι ώστε Y=κ⁢X1+t⁢X2 . Ο στόχος μας, λοιπόν, είναι να δείξουμε ότι κ=0 . Για να τον πετύχουμε, θα υπολογίσουμε το γινόμενο YT⁢X1=(A⁢X2)T⁢X1 . Κάνοντας τις πράξεις βρίσκουμε ότι

(A⁢X2)T⁢X1=X2T⁢(AT⁢X1)=X2T⁢(A⁢X1)=X2T⁢(λ1⁢X1)=λ1⁢(X2T⁢X1)=0 .

Αντικαθιστώντας κ⁢X1+t⁢X2 στη θέση του A⁢X2 βρίσκουμε ότι

(A⁢X2)T⁢X1=(κ⁢X1+t⁢X2)Τ⁢X1=(κ⁢X1T+t⁢X2T)⁢X1
=κ(X1TX1)+t(X2TX1)=κ(X1TX1)+t𝟎=κ(X1TX1).

Αφού X1T⁢X1≠𝟎 και κ⁢(X1T⁢X1)=𝟎 έπεται ότι κ=0 . Άρα A⁢X2=τ⁢X2 και το X2 είναι ιδιοδιάνυσμα του Α για την ιδιοτιμή t . Επομένως ο 2×2 πίνακας Α έχει δύο ορθογώνια ιδιοδιανύσματα και είναι ορθογώνια διαγωνιοποιήσιμος.∎

Έστω ότι B είναι η κανονική βάση του ℝn και έστω ότι ϕ:ℝn→ℝn είναι μία γραμμική συνάρτηση. Παρατηρούμε ότι

AB,Bϕ=(AB,Bϕ)T⁢ αν και μόνο an⁢ϕ=ϕ*.

H γραμμική συνάρτηση ϕ:ℝn→ℝn λέγεται αυτοπροσαρτημένη (self-adjoint) an ϕ=ϕ* ενώ o πίνακας A∈ℳn⁢(ℝ) λέγεται αυτοπροσαρτημένος αν A=AT . Σύμφωνα με το Θεώρημα 5.5.5 η συνάρτηση ϕ:ℝn→ℝn είναι ορθογώνια διαγωνιοποιήσιμη αν και μόνο αν η ϕ είναι αυτοπροσαρτημένη.

Παραδείγματα 5.5.7.

  1. 1.

    Στο Παράδειγμα 5.2.5.2 δείξαμε ότι η συνάρτηση

    ϕ:ℝ2→ℝ2,(x,y)↦(2⁢x+y,x+2⁢y)

    είναι διαγωνιοποιήσιμη και ότι

    PT⁢Ab,Bϕ⁢P=[1003],

    όπου

    P=12⁢[-11  11].

    Τα ιδιοδιανύσματα για τις ιδιοτιμές 1 και 3 της ϕ είναι ορθογώνια και ϕ=ϕ*.

  2. 2.

    Θεωρούμε τον πίνακα

    Α=[1-10-12-10-11].

    Παρατηρούμε ότι Α=ΑΤ . Υπολογίζουμε το χαρακτηριστικό πολυώνυμο του Α.

    PA⁢(x)=|1-x-1  0-12-x-1  0-11-x|=(1-x)⁢[(2-x)⁢(1-x)-1]+[-(1-x)]
    ⇒PA(x)=(1-x)x(x-3).

    Οι ρίζες λ1=1,λ2=0,λ3=3 του PA⁢(x) είναι οι ιδιοτιμές του A . Υπολογίζουμε τους ιδιοχώρους V1⁢(A),V0⁢(A),V3⁢(A) του Α . Έχουμε ότι

    • –

      V1⁢(A)= {(-κ,0,κ):κ∈ℝ}=S⁢({(-1,0,1)}),

    • –

      V0⁢(A)= {(κ,κ,κ):κ∈ℝ}=S⁢({(1,1,1)}) , και

    • –

      V3(A)={(κ,-2κ,κ):κ∈ℝ}=S({(1,-2,1)}.

    Επιλέγουμε ένα διάνυσμα από κάθε ιδιοχώρο και βρίσκουμε μία βάση του ℝ3 . An ε1=(-1,0,1),ε2=(1,1,1),ε3=(1,-2,1) , τότε (ε1,ε2,ε3) είναι ορθογώνια βάση του ℝ3 . Για i=1,…,3 , διαιρούμε τα εi με το μήκος τους και βρίσκουμε μία ορθοκανονική βάση του ℝ3 από ιδιοδιανύσματα του Α:

    (12⁢(-1,0,1),13⁢(1,1,1),16⁢(1,-2,1)).

    Ο πίνακας

    S=[-1/21/31/601/3-2/61/21/31/6],

    είναι ορθογώνιος και

    ST⁢A⁢S=[100000003].
  3. 3.

    Έστω ϕ:ℝ3→ℝ3,ϕ⁢(x,y,z)=(x,2⁢y+4⁢z,4⁢y+3⁢z) . O πίνακας της ϕ ως προς την κανονική βάση B του ℝ3 είναι o

    Α=[100024043].

    Αφού Α=ΑΤ , o Α είναι ορθογώνια διαγωνιοποιήσιμος και η ϕ=ϕ* είναι αυτοπροσαρτημένη γραμμική συνάρτηση. Παρατηρούμε ότι

    PA⁢(x)=(1-x)⁢(x2-5⁢x-10)

    και ότι οι ιδιοτιμές του Α είναι οι

    λ1=1,λ2=5+652,λ3=5-652.

    H ελαττωμένη κλιμακωτή μορφή γραμμών του A-I3 είναι o πίνακας

    [010001000].

    Έπεται ότι V1⁢(A)= {t⁢(1,0,0):t∈ℝ} και ότι v1=(1,0,0) είναι ιδιοδιάνυσμα για την ιδιοτιμή 1. Στοιχειώδεις πράξεις γραμμών φέρνουν τον πίνακα Α-λ2⁢I3 στη μορφή

    [1000-1+65/24000].

    Ο πίνακας αυτός δεν είναι σε ελαττωμένη κλιμακωτή μορφή γραμμών. Παρόλα αυτά μας δίνει όλες τις πληροφορίες που χρειαζόμαστε για να υπολογίσουμε τον ιδιοχώρο Vλ2⁢(A) . Πράγματι αν (x,y,z)∈Vλ2⁢(A) τότε από τη πρώτη γραμμή του πίνακα αυτού προκύπτει ότι x=0 , ενώ από τη δεύτερη γραμμή προκύπτει ότι

    -1+652⁢y+4⁢z=0⇒z=1-658⁢y.

    Άρα

    Vλ2⁢(A)= {t⁢(0,1,1-658):t∈ℝ}={s⁢(0,8,1-65):s∈ℝ}

    και v2=(0,8,1-65) είναι ιδιοδιάνυσμα για την ιδιοτιμή λ2 . Ομοίως

    Vλ3⁢(A)= {t⁢(0,8,1+65):t∈ℝ}

    και v3=(0,8,1+65) είναι ιδιοδιάνυσμα για την ιδιοτιμή λ3 . Το σύνολο

    (v1,v2,v3)

    είναι ορθογώνια βάση του ℝ3 που αποτελείται από ιδιοδιανύσματα του A . Κανονικοποιούμε και βρίσκουμε μία ορθοκανονική βάση του ℝ3:

    ((1,0,0),1130-2⁢65⁢(0,8,1-65),1130+2⁢65⁢(0,8,1+65)).

    Ο πίνακας

    P=[10008130-2⁢658130+2⁢6501-65130-2⁢65)1+65130+2⁢65]

    είναι ορθογώνιος και

    [1000λ2000λ3]=PT⁢A⁢P.

To Θεώρημα 5.5.5 προσδιορίζει πότε ακριβώς ένας πίνακας με στοιχεία πραγματικούς αριθμούς είναι ορθογώνια διαγωνιοποιήσιμος. Θέτουμε το αντίστοιχο ερώτημα για πίνακες με στοιχεία μιγαδικούς αριθμούς. An A∈ℳn⁢(ℂ) , πότε έχειο ℂn μία βάση ορθογώνιων ιδιοδιανυσμάτων του A; O A∈ℳn⁢(ℂ) λέγεται ορθογώνια διαγωνιοποιήσιμος (orthogonal diagonalizable), αν υπάρχει ένας πίνακας P έτσι ώστε P-1=(PT)¯ και

(PT)¯⁢A⁢P=[λ1⋯00⋮λn].

Οι στήλες του P είναι ιδιοδιανύσματα του A και λ1,…,λn είναι τα ιδιοδιανύσματα του A . O A λέγεται αυτοπροσαρτημένος ή Ερμητιανός (Hermitian) an A=(AT)¯ . Στο Παράδειγμα 5.1.6.11 δείξαμε ότι an A είναι αυτοπροσαρτημένος, τότε ο Α έχει πραγματικές ιδιοτιμές. Μπορεί να αποδειχθεί ότι ισχύει το αντίστοιχο του Θεωρήματος 5.5.2 για τους μιγαδικούς n×n πίνακες και το σύνηθες εσωτερικό γινόμενο στον ℂn.

Θεώρημα 5.5.8.

Αν A∈ℳn⁢(ℂ) και v∈ℂn , τότε ⟨A⁢v, v⟩=⟨v,AT¯⁢v⟩.

Όπως στην απόδειξη του Φασματικού Θεωρήματος προκύπτει ότι αν o A είναι αυτοπροσαρτημένος τότε o A είναι ορθογώνια διαγωνιοποιήσιμος. Το αντίστροφο όμως δεν ισχύει. Υπάρχουν πίνακες στο ℳn⁢(ℂ) που είναι ορθογώνια διαγωνιοποιήσιμοι και που δεν είναι αυτοπροσαρτημένοι. Ποια είναι η κοινή ιδιότητα αυτών των πινάκων; Ένας πίνακας A∈ℳn⁢(ℂ) λέγεται κανονικός (normal) αν ισχύει η ακόλουθη συνθήκη:

A⁢A¯T=A¯T⁢A.

Είναι φανερό ότι an A=(AT)¯ τότε o A είναι κανονικός. Μπορεί, λοιπόν, κανείς να αποδείξει την εξής σημαντική πρόταση.

Θεώρημα 5.5.9.

Ένας πίνακας Α∈ℳn⁢(ℂ) είναι ορθογώνια διαγωνιοποιήσιμος αν και μόνο αν o A είναι κανονικός. Ειδκότερα, an A=(AT)¯ , τότε o A είναι ορθογώνια διαγωνιοποιήσιμος.

Για την απόδειξη της προηγούμενης πρότασης, παραπέμπουμε στο σύγγραμμα [2, Θεώρημα 7.3.11].

Παραδείγματα 5.5.10.

O πίνακας

A=[110011101]

είναι κανονικός και δεν είναι αυτοπροσαρτημένος. Στην Άσκηση 5.5.7, o αναγνώστης καλείται να διαγωνιοποιήσει ορθογώνια ton A.

Σημειώνουμε ότι οι ιδιοτιμές ενός κανονικού πίνακα δεν είναι κατ' ανάγκη πραγματικοί αριθμοί. Παρατηρούμε, όμως, ότι οι ορθογώνιοι πίνακες είναι κανονικοί. Επίσης, οι ιδιοτιμές ενός ορθογώνιου πίνακα έχουν μέτρο 1 , βλ. Άσκηση 5.5.5. Ισχύει, λοιπόν, το επόμενο συμπέρασμα.

Πόρισμα 5.5.11.

Έστω A ένας ορθογώνιος n×n πίνακας. Τότε o A είναι ορθογώνια διαγωνιοποιήσιμος. Οι ιδιοτιμές του A έχουν μέτρο 1.

Ασκήσεις Eνότητας 5.5

  1. 1.

    Έστω ϕ:ℝ3→ℝ3, ϕ⁢(x,y,z)=(x+2⁢y,x+y+z,y-x). Να βρεθεί o τύπος της ϕ*.

  2. 2.

    Έστω A∈ℳm×n⁢(ℝ). Να αποδείξετε ότι null⁡(A)=null⁡(AT⁢A). Να συμπεράνετε ότι rank⁡(AT⁢A)=rank⁡(A).

  3. 3.

    Έστω A∈ℳn⁢(ℝ) ορθογώνιος πίνακας, λ ιδιοτιμή του Α και v ιδιοδιάνυσμα του Α για την ιδιοτιμή λ. Εξετάζοντας το εσωτερικό γινόμενο ⟨A⁢v,A⁢v⟩ να αποδείξετε ότι λ=±1.

  4. 4.

    Να βρείτε έναν ορθογώνια διαγωνιοποιήσιμο πίνακα A στον ℳ3⁢(ℝ) με πρώτη γραμμή [123]. Να βρείτε τις ιδιοτιμές του A και βάσεις για τους ιδιοχώρους του.

  5. 5.

    Έστω A∈ℳn⁢(ℂ) ορθογώνιος πίνακας και λ μία ιδιοτιμή του Α. Να αποδείξετε ότι ∥λ∥=1. (Υπόδειξη: Άσκηση 3)

  6. 6.

    Να βρείτε έναν ορθογώνια διαγωνιοποιήσιμο πίνακα A στον ℳ2⁢(ℂ) με πρώτη γραμμή [11+3⁢i]. Να βρείτε τις ιδιοτιμές του A και βάσεις για τους ιδιοχώρους του.

  7. 7.

    Να διαγωνιοποιήσετε ορθογώνια τον

    A=[110011101].

5.6 Σύντομα Ιστορικά Στοιχεία

Οι ιδιοτιμές εμφανίστηκαν για πρώτη φορά στη μελέτη των τετραγωνικών μορφών και των διαφορικών εξισώσεων. H σημασία της ύπαρξης των ιδιοδιανυσμάτων έγινε αντιληπτή τον 18ο αιώνα από τους Euler και Lagrange. Στις αρχές του 19ου αιώνα o Cauchy αντιλήφθηκε ότι τα αποτελέσματά τους μπορούν να χρησιμοποιηθούν για την ταξινόμηση των κωνικών επιφανειών, ενώ o Fourier εφάρμοσε τις ιδέες τους στο βιβλίο του Αναλυτική Θεωρία της Θερμότητας. O όρος χαρακτηριστικό πολυώνυμο οφείλεται στον Cauchy, o οποίος και απέδειξε το Φασματικό Θεώρημα. O Hermite εισήγαγε τους Ερμητιανούς πίνακες το 1855 και στις αρχές του 20ου αιώνα o Hilbert επέκτεινε τη θεωρία σε άπειρους πίνακες. O όρος eigenvectors οφείλεται σε αυτόν. H κανονική μορφή του Jordαν εμφανίστηκε στη δουλειά του Jordan το 1870.

Βιβλιογραφία

  1. 1.

    H. Anton, C. Rorres, Elementary Linear Algebra, Applications Version, John Wiley and Sons, 1994.

  2. 2.

    Θ. Θεοχάρη-Αποστολίδη, Χ. Χαραλάμπους, Β. Βαβατσούλας, Εισαγωγή στη Γραμμική Άλγεβρα, Θεσσαλονίκη 2006.

  3. 3.

    V. Katz, Ιστορία των Μαθηματικών, Πανεπιστημιακές Εκδόσεις Κρήτης, 2013.

  4. 4.

    K. Nicholoson, Elementary Linear Algebra, McGraw-Hill, 2001.

  5. 5.

    Th. Shiffrin and M. R.Adams, Linear Algebra, a Geometric Approach, W. H. Freeman and Company, 2002.

  6. 6.

    S. H. Weintraub, Jordan Canonical Form, Theory and Practice, Morgan & Claypool Publishers, 2009.