Tính toán số mũ ma trận
Cho trước ma trận $$M = \begin{pmatrix} 7i& -6-2i\\6-2i&-7i\end{pmatrix}$$ làm cách nào để tính toán ma trận theo cấp số nhân $e^M$?
Tôi biết tôi có thể sử dụng nó $e^A=Pe^DP^{-1}$ Ở đâu $D=P^{-1}AP$. Tôi đã tính toán đa thức đặc trưng của ma trận trên là
$$P(\lambda)=\lambda^2+89$$
Có cách nào dễ dàng hơn để làm điều này ngoài việc cố gắng tính toán ma trận đường chéo không?
Trả lời
Ma trận của bạn $M$ có thể theo đường chéo với các giá trị riêng $\pm i\sqrt{89}$. Điều này có nghĩa rằng$e^M = p(M)$ Ở đâu $p \in \Bbb{C}[x]$ là đa thức duy nhất có bậc nhỏ hơn $2$ như vậy mà $$p(i\sqrt{89}) = e^{i\sqrt{89}}, \quad p(-i\sqrt{89}) = e^{-i\sqrt{89}}.$$
Sử dụng công thức nội suy Lagrange, chúng ta thấy rằng $$p(x) = \frac{x+i\sqrt{89}}{2i\sqrt{89}}e^{i\sqrt{89}}-\frac{x-i\sqrt{89}}{2i\sqrt{89}}e^{-i\sqrt{89}} = \frac{\sin\sqrt{89}}{\sqrt{89}}x+ \cos\sqrt{89}$$ vì thế $$e^M = p(M) = \frac{\sin\sqrt{89}}{\sqrt{89}}M+ \cos\sqrt{89}I = \left( \begin{array}{cc} \frac{7 i \sin \sqrt{89}}{\sqrt{89}} +\cos\sqrt{89}& -\frac{(6+2 i) \sin \sqrt{89}}{\sqrt{89}} \\ \frac{(6-2 i) \sin \sqrt{89}}{\sqrt{89}} & -\frac{7 i \sin \sqrt{89}}{\sqrt{89}} +\cos\sqrt{89} \\ \end{array} \right).$$
Via Cayley-Hamilton, ${\rm M}^2 + 89 \, {\rm I}_2 = {\rm O}_2$. Vì thế,
$$\begin{aligned} {\rm M}^2 &= - 89 \, {\rm I}_2\\ {\rm M}^3 &= - 89 \, {\rm M}\\ {\rm M}^4 &= 89^2 {\rm I}_2\\ {\rm M}^5 &= 89^2 {\rm M}\\ &\vdots\\ {\rm M}^{2k} &= (-1)^k 89^k {\rm I}_2\\ {\rm M}^{2k+1} &= (-1)^k 89^k {\rm M} \end{aligned}$$
và
$$\exp({\rm M}) = \sum_{k=0}^{\infty} \frac{{\rm M}^k}{k!} = \cdots = \color{blue}{\cos( \sqrt{89} ) \, {\rm I}_2 +\frac{\sin( \sqrt{89} )}{\sqrt{89}} {\rm M}}$$
Bạn cũng có thể sử dụng phép chia Euclide để tính hàm mũ khi bạn có đa thức đặc trưng:
$$\forall n\in\mathbb{N}^*,\exists(P_n,R_n)\in\mathbb{R}[X], \deg(R)\leq1, X^n=P_n(X)\times(X^2+89)+R_n(X) $$
Để cho $R_n(X)=a_nX+b_n$, bạn lấy : $$M^n=a_nM+b_n1 $$
và bạn nhận được các giá trị của $a_n,b_n$ với : $$(i\sqrt{89})^n=a_ni\sqrt{89}+b_n\\ (-i\sqrt{89})^n=-a_ni\sqrt{89}+b_n$$
Vì vậy, bạn nhận được: $$b_n=\frac{(i\sqrt{89})^n+(-i\sqrt{89})^n}{2}\\ a_n=\frac{(i\sqrt{89})^n-(-i\sqrt{89})^n}{2i\sqrt{89}} $$
Và sau đó bạn sử dụng định nghĩa: $$e^M=\sum_{n=0}^{\infty}\frac{M^n}{n!}=\sum_{n=0}^{\infty}\frac{a_nM+b_n1}{n!} $$
Và sau đó bạn có thể kết luận từ đây.
Via Cayley-Hamilton,
$${\rm M}^2 = - 89 \, {\rm I}_2 = \left( i \sqrt{89} \right)^2 {\rm I}_2$$
Do đó, ma trận ${\rm A} := \frac{{\rm M}}{i \sqrt{89}}$là involutory , tức là,${\rm A}^2 = {\rm I}_2$. Sử dụng công thức của Euler ,
$$\begin{aligned} \exp({\rm M}) = \exp \left( i \sqrt{89} {\rm A} \right) &= \cos \left( \sqrt{89} {\rm A} \right) + i \sin \left( \sqrt{89} {\rm A} \right)\\ &= \cos \left( \sqrt{89} \right) {\rm I}_2 + i \sin \left( \sqrt{89} \right) {\rm A}\\ &= \color{blue}{\cos \left( \sqrt{89} \right) {\rm I}_2 + \frac{\sin \left( \sqrt{89} \right)}{\sqrt{89}} {\rm M}}\end{aligned}$$
Ở đâu $\cos \left( \sqrt{89} {\rm A} \right) = \cos \left( \sqrt{89} \right) {\rm I}_2$ và $\sin \left( \sqrt{89} {\rm A} \right)= \sin \left( \sqrt{89} \right) {\rm A}$ bởi vì ${\rm A}$là involutory .