Σελίδα 1 από 1

O e^A είναι κανονικός

Δημοσιεύτηκε: Παρ Αύγ 16, 2019 2:51 pm
από Tolaso J Kos
Να δοθεί παράδειγμα πίνακα A ο οποίος δεν είναι κανονικός αλλά ο e^A είναι κανονικός.

Δεν έχω λύση !

Re: O e^A είναι κανονικός

Δημοσιεύτηκε: Τρί Αύγ 20, 2019 3:33 pm
από ΠΑΠΑΔΟΠΟΥΛΟΣ ΣΤΑΥΡΟΣ
Tolaso J Kos έγραψε:
Παρ Αύγ 16, 2019 2:51 pm
Να δοθεί παράδειγμα πίνακα A ο οποίος δεν είναι κανονικός αλλά ο e^A είναι κανονικός.

Δεν έχω λύση !
Εστω
A=PDP^{-1}
οπου
P=\begin{pmatrix} 2 & 1 & \\ 1& 1 & \end{pmatrix}
και
D=\begin{pmatrix} 2\pi i & 0 & \\ 0& -2\pi i & \end{pmatrix}

λόγω της σχέσης

exp(PDP^{-1})=PexpDP^{-1}

εχουμε ότι expA=I κανονικός

ενώ ο A δεν είναι .(πράξεις)

Για τις ιδιότητες που χρησιμοποιήθηκαν κοίτα στο

https://en.wikipedia.org/wiki/Matrix_exponential

Re: O e^A είναι κανονικός

Δημοσιεύτηκε: Τρί Αύγ 20, 2019 7:09 pm
από Mihalis_Lambrou
Σταύρο, με πρόλαβες. Πολύ ωραίο το παράδειγμα και η τεχνική. Είχα ένα παράδειγμα με τελεστή τάξης 1, και διαφορετική (πιο πολύπλοκη) μέθοδο, που δεν ευδόκησα να γράψω. Περιληπτικά:

Έστω δύο διανύσματα p,q (που θα τα διαλέξουμε μετά) σε χώρο με εσωτερικό γινόμενο (και Hilbert αν θέλουμε). Θέτουμε Rx=<x,p>q. Στην βιβλιογραφία ο τελεστής αυτός συμβολίζεται p\otimes q, και θα ακολουθήσω τον συμβολισμό αν και μπορούμε να τον αποφύγουμε.

Εύκολα βλέπουμε επαγωγικά ότι R^n=<q,p>^{n-1}p\otimes q (πολλαπλάσιο του εαυτού του). Άρα, με ακριβώς τον ίδιο τρόπο όπως στην Άσκηση 1 εδώ έχουμε ότι

e^R= \mathbb I + \frac {1}{<q,p>} (e^{<q,p>}-1)R \, (*)

Εύκολα βρίσκουμε διανύσματα p,q με <q,p>=2\pi i οπότε, εδώ, e^R=\mathbb I, και άρα e^R κανονικός.

Όμως R^*=q\otimes p (απλό και γνωστό), οπότε RR^*=<p,p>q\otimes q και R^*R=<q,q>p\otimes p. Tα δύο τελευταία δεν είναι ίσα αν πάρουμε (που μπορούμε) τα p,q γραμμικά ανεξάρτητα. Δηλαδή R μη κανονικός, όπως θέλαμε.

Re: O e^A είναι κανονικός

Δημοσιεύτηκε: Τρί Αύγ 20, 2019 7:45 pm
από Tolaso J Kos
Για να δω αν έχω κάνει σωστά τις πράξεις... Βγάζω ότι

\displaystyle{A = P D P^{-1} = \begin{pmatrix} 
2 &1 \\ 
1 & 1 
\end{pmatrix} \begin{pmatrix} 
2 \pi i &0 \\ 
0 & -2\pi i 
\end{pmatrix} \begin{pmatrix} 
1 &-1 \\ 
-1 & 2 
\end{pmatrix} = \begin{pmatrix} 
2 \pi i & -8 \pi i \\ 
0 & - 6 \pi i 
\end{pmatrix}}

και ότι

\displaystyle{\exp A = \exp \left ( P D P^{-1} \right ) = \begin{pmatrix} 
1 & 1\\ 
1 & 1 
\end{pmatrix}}

Τα επιβεβαιώνει κανείς; Τα έκανα με το Wolfram ... :(

Re: O e^A είναι κανονικός

Δημοσιεύτηκε: Τρί Αύγ 20, 2019 10:03 pm
από Mihalis_Lambrou
Tόλη γράφω υπό πίεση γιατί έχω πολύ φόρτο εργασίας αλλά, αν δεν βιάστηκα, μου φαίνεται ότι σωστά τα λέει ο Σταύρος.
Tolaso J Kos έγραψε:
Τρί Αύγ 20, 2019 7:45 pm
Για να δω αν έχω κάνει σωστά τις πράξεις... Βγάζω ότι

\displaystyle{A = P D P^{-1} = \begin{pmatrix} 
2 &1 \\ 
1 & 1 
\end{pmatrix} \begin{pmatrix} 
2 \pi i &0 \\ 
0 & -2\pi i 
\end{pmatrix} \begin{pmatrix} 
1 &-1 \\ 
-1 & 2 
\end{pmatrix} = \begin{pmatrix} 
2 \pi i & -8 \pi i \\ 
0 & - 6 \pi i 
\end{pmatrix}}
Βγάζω το δεξί μέλος ίσο με   \begin{pmatrix} 
6 \pi i & -8 \pi i \\ 
4 \pi i & - 6 \pi i 
\end{pmatrix}}
Tolaso J Kos έγραψε:
Τρί Αύγ 20, 2019 7:45 pm
\displaystyle{\exp A = \exp \left ( P D P^{-1} \right ) = \begin{pmatrix} 
1 & 1\\ 
1 & 1 
\end{pmatrix}}
Προσοχή. Για διαγώνιο πίνακα D= \begin{pmatrix} 
p & 0 \\ 
0 & q 
\end{pmatrix}} , εύκολα βλέπουμε ότι e^D=  \begin{pmatrix} 
e^p & 0 \\ 
0 & e^q 
\end{pmatrix}}

οπότε για τον παραπάνω D, δεδομένου ότι e^{\pm 2\pi i } =1, ισχύει (όπως σωστά γράφει ο Σταύρος), e^A= \mathbb I

Re: O e^A είναι κανονικός

Δημοσιεύτηκε: Τρί Αύγ 20, 2019 10:16 pm
από Demetres
Τόλη είναι σωστό το e^A = I. Ακόμη και για το A όπως το υπολόγισες εσύ. Υποψιάζομαι ότι δεν χρησιμοποίησες τη σωστή εντολή για τον εκθετικό πίνακα. Πιο συγκεκριμένα υποψιάζομαι ότι η εντολή που χρησιμοποίησες απλά παίρνει την εκθετική συνάρτηση σε κάθε στοιχείο ξεχωριστά.

Re: O e^A είναι κανονικός

Δημοσιεύτηκε: Τετ Αύγ 21, 2019 10:15 am
από Tolaso J Kos
Tolaso J Kos έγραψε:
Παρ Αύγ 16, 2019 2:51 pm
Να δοθεί παράδειγμα πίνακα A ο οποίος δεν είναι κανονικός αλλά ο e^A είναι κανονικός.
Ωραία, ας βάλουμε τις λεπτομέρειες. Έστω P= \begin{pmatrix} 2 &1 \\ 1 & 1 \end{pmatrix} που προφανώς είναι αντιστρέψιμος με  P^{-1} = \begin{pmatrix} 1 &-1 \\ -1 & 2 \end{pmatrix}. Θεωρούμε D = \begin{pmatrix} 2 \pi i &0 \\ 0 & -2\pi i \end{pmatrix}. Τότε,

\displaystyle A = P D P^{-1} = \begin{pmatrix} 2 &1 \\ 1 & 1 \end{pmatrix} \begin{pmatrix} 2 \pi i &0 \\ 0 & -2\pi i \end{pmatrix} \begin{pmatrix} 1 &-1 \\ -1 & 2 \end{pmatrix} = \begin{pmatrix} 6 \pi i & -8 \pi i \\ 4\pi i & - 6 \pi i \end{pmatrix}
Έπεται ότι,

\displaystyle A^* = \begin{pmatrix} -6 \pi i & -4\pi i\\ 8 \pi i & 6 \pi i \end{pmatrix}
Ο A δεν είναι κανονικός διότι:

\displaystyle A A^* = \begin{pmatrix} 6 \pi i & -8 \pi i \\ 4\pi i & - 6 \pi i \end{pmatrix} \begin{pmatrix} -6 \pi i & -4\pi i\\ 8 \pi i & 6 \pi i \end{pmatrix} = \begin{pmatrix} 100 \pi^2 & 72 \pi^2 i \\ 72 \pi^2 & 52 \pi^2 i \end{pmatrix}
και

\displaystyle A^* A =\begin{pmatrix} -6 \pi i & -4\pi i\\ 8 \pi i & 6 \pi i \end{pmatrix} \begin{pmatrix} 6 \pi i & -8 \pi i \\ 4\pi i & - 6 \pi i \end{pmatrix} = \begin{pmatrix} 52 \pi^2 & -72 \pi^2 \\ -72 \pi^2 & 100 \pi^2 \end{pmatrix}
Όμως,

\displaystyle{\begin{aligned} 
 \exp A &= \exp \left ( P D P^{-1} \right ) \\  
&=P \exp (D) P^{-1} \\  
&= P \; \mathbb{I}_{2 \times 2} \; P^{-1} \\  
&= P P^{-1} \\  
&= \mathbb{I}_{2 \times 2}  
\end{aligned}}
ο οποίος είναι κανονικός.


Ευχαριστούμε Σταύρο!

Re: O e^A είναι κανονικός

Δημοσιεύτηκε: Τετ Αύγ 21, 2019 12:50 pm
από ΠΑΠΑΔΟΠΟΥΛΟΣ ΣΤΑΥΡΟΣ
Tolaso J Kos έγραψε:
Τετ Αύγ 21, 2019 10:15 am
Tolaso J Kos έγραψε:
Παρ Αύγ 16, 2019 2:51 pm
Να δοθεί παράδειγμα πίνακα A ο οποίος δεν είναι κανονικός αλλά ο e^A είναι κανονικός.
Ωραία, ας βάλουμε τις λεπτομέρειες. Έστω P= \begin{pmatrix} 2 &1 \\ 1 & 1 \end{pmatrix} που προφανώς είναι αντιστρέψιμος με  P^{-1} = \begin{pmatrix} 1 &-1 \\ -1 & 2 \end{pmatrix}. Θεωρούμε D = \begin{pmatrix} 2 \pi i &0 \\ 0 & -2\pi i \end{pmatrix}. Τότε,

\displaystyle A = P D P^{-1} = \begin{pmatrix} 2 &1 \\ 1 & 1 \end{pmatrix} \begin{pmatrix} 2 \pi i &0 \\ 0 & -2\pi i \end{pmatrix} \begin{pmatrix} 1 &-1 \\ -1 & 2 \end{pmatrix} = \begin{pmatrix} 6 \pi i & -8 \pi i \\ 4\pi i & - 6 \pi i \end{pmatrix}
Έπεται ότι,

\displaystyle A^* = \begin{pmatrix} -6 \pi i & -4\pi i\\ 8 \pi i & 6 \pi i \end{pmatrix}
Ο A δεν είναι κανονικός διότι:

\displaystyle A A^* = \begin{pmatrix} 6 \pi i & -8 \pi i \\ 4\pi i & - 6 \pi i \end{pmatrix} \begin{pmatrix} -6 \pi i & -4\pi i\\ 8 \pi i & 6 \pi i \end{pmatrix} = \begin{pmatrix} 100 \pi^2 & 72 \pi^2 i \\ 72 \pi^2 & 52 \pi^2 i \end{pmatrix}
και

\displaystyle A^* A =\begin{pmatrix} -6 \pi i & -4\pi i\\ 8 \pi i & 6 \pi i \end{pmatrix} \begin{pmatrix} 6 \pi i & -8 \pi i \\ 4\pi i & - 6 \pi i \end{pmatrix} = \begin{pmatrix} 52 \pi^2 & -72 \pi^2 \\ -72 \pi^2 & 100 \pi^2 \end{pmatrix}
Όμως,

\displaystyle{\begin{aligned} 
 \exp A &= \exp \left ( P D P^{-1} \right ) \\  
&=P \exp (D) P^{-1} \\  
&= P \; \mathbb{I}_{2 \times 2} \; P^{-1} \\  
&= P P^{-1} \\  
&= \mathbb{I}_{2 \times 2}  
\end{aligned}}
ο οποίος είναι κανονικός.


Ευχαριστούμε Σταύρο!
Να είσαι πάντα καλά Τόλη να μας διασκεδάζεις.
(εννοείται με τα προβλήματα που βάζεις)

Η ιδέα είναι απλή για την κατασκευή.

Λόγω της σχέσης

exp(PDP^{-1})=PexpDP^{-1}

αν πάρουμε τον D ετσι ώστε expD=I

τότε ο δεξιός πίνακας είναι κανονικός.

Τον D τον βρίσκουμε εύκολα.

Αρκεί να βρούμε τον P ώστε ο PDP^{-1}

να μην είναι κανονικός.

Οποίον P και να πάρουμε που οι στήλες του δεν είναι κάθετες
μας κάνει.

Μπροστά μου είχα από άλλο πρόβλημα αυτόν που πήρα και είχα υπολογίσει τον αντίστροφο του.

Πήρα λοιπόν αυτόν για να κάνω τις πράξεις ,ώστε να είμαι απόλυτα σίγουρος.

Re: O e^A είναι κανονικός

Δημοσιεύτηκε: Τετ Αύγ 21, 2019 5:04 pm
από Mihalis_Lambrou
Ωραία.
Tolaso J Kos έγραψε:
Τετ Αύγ 21, 2019 10:15 am


Ωραία, ας βάλουμε τις λεπτομέρειες. Έστω P= \begin{pmatrix} 2 &1 \\ 1 & 1 \end{pmatrix} που προφανώς είναι αντιστρέψιμος με  P^{-1} = \begin{pmatrix} 1 &-1 \\ -1 & 2 \end{pmatrix}. Θεωρούμε D = \begin{pmatrix} 2 \pi i &0 \\ 0 & -2\pi i \end{pmatrix}. Τότε,

...

Ο A δεν είναι κανονικός διότι:


\displaystyle A A^* =...  = \begin{pmatrix} 100 \pi^2 & 72 \pi^2 i \\ 72 \pi^2 & 52 \pi^2 i \end{pmatrix} και \displaystyle A^* A =... = \begin{pmatrix} 52 \pi^2 & -72 \pi^2 \\ -72 \pi^2 & 100 \pi^2 \end{pmatrix}


Όμως,


\displaystyle{ 
 \exp A=... =  \mathbb{I}_{2 \times 2}
Αν θέλεις χάριν περιέργειας να κάνεις παρόμοιες πράξεις με το παράδειγμα που έδωσα, στην ειδική περίπτωση p=\begin{pmatrix} 
2\pi i \\ 1 
 
\end{pmatrix}, \, q= \begin{pmatrix} 
1\\ 0 
 
\end{pmatrix} , θα βρεις R = \begin{pmatrix} 2 \pi i &1 \\ 0 & 0  \end{pmatrix}, οπότε  RR^*= \begin{pmatrix} -4\pi ^2  &0 \\ 0 & 0  \end{pmatrix} \ne \begin{pmatrix} -4\pi ^2  &-2\pi i \\ 0 & 0  \end{pmatrix}= R^*R . Επίσης e^R=\mathbb I