9.2

Дополнительные задачи

[33/0%]
LaTeX
Задача 9.2.1

Пусть z1,z2,z3∈Cz_1, z_2, z_3 \in \mathbb {C}. Рассмотрим матрицу 2×22 \times 2

A(z1,z2,z3)=(z1z2z30). A(z_1,z_2,z_3) = \begin{pmatrix} z_1 & z_2 \\ z_3 & 0 \end{pmatrix}.

Вычислите exp(A(z1,z2,z3))\mathrm{exp}\left(A(z_1,z_2,z_3)\right).

?
Задача 9.2.2

Пусть α∈R\alpha \in \mathbb {R}. Вычислите вектор в R2\mathbb {R}^{2}

exp⁡(α(0100))(10) \exp \left(\alpha \begin{pmatrix} 0 & 1 \\ 0 & 0 \end{pmatrix}\right) \begin{pmatrix} 1 \\ 0 \end{pmatrix}

двумя различными способами. Сравните и обсудите.

?
Задача 9.2.3

Пусть AA, BB — матрицы n×nn \times n.

?
(i)

Покажите, что

[eA,eB]=eAeB−eBeA=∑j,k=1∞1j!k![Aj,Bk]. [e^A,e^B] = e^Ae^B - e^Be^A = \sum _{j,k=1}^{\infty } \frac{1}{j!k!}[A^j,B^k].
(ii)

Предположим, что для всех собственных значений λ\lambda выполнено ℜ(λ)<0\Re (\lambda ) < 0. Пусть BB — произвольная матрица n×nn \times n над C\mathbb {C}. Пусть

R(t):=∫0∞etA∗BetAdt. R(t) := \int _0^{\infty } e^{tA^*}Be^{tA}dt.

Покажите, что матрица RR удовлетворяет R(t)A+A∗R(t)=−BR(t)A + A^*R(t) = -B.

Задача 9.2.4

Пусть AA, BB — матрицы n×nn \times n над C\mathbb {C}. Верно ли, что tr⁡(eA⊗eB)=tr⁡(eA⊗B)\operatorname {tr}(e^A \otimes e^B) = \operatorname {tr}(e^{A \otimes B})? Докажите или опровергните.

?
Задача 9.2.5
?
(i)

Пусть A1,A2,…,ApA_1, A_2, \ldots , A_p — матрицы n×nn \times n над C\mathbb {C}. Обобщённая формула Троттера задаётся как

exp⁡(∑j=1nAj)=lim⁡n→∞fn({Aj})(1) \exp \left(\sum _{j=1}^{n} A_j\right) = \lim _{n \to \infty } f_n(\left\{ A_j\right\} ) \tag {1}

где nn-й аппроксимант fn({Aj})f_n(\left\{ A_j\right\} ) определяется как

fn({Aj}):=(exp⁡(1nA1)exp⁡(1nA2)⋯exp⁡(1nAp))n. f_n(\left\{ A_j\right\} ) := \left(\exp \left(\frac{1}{n}A_1\right)\exp \left(\frac{1}{n}A_2\right) \cdots \exp \left(\frac{1}{n}A_p\right)\right)^n.

Пусть p=2p = 2 и

A1=(1001),A2=(0110). A_1 = \begin{pmatrix} 1 & 0 \\ 0 & 1 \end{pmatrix}, \qquad A_2 = \begin{pmatrix} 0 & 1 \\ 1 & 0 \end{pmatrix}.

Вычислите левую и правую части (1).

(ii)

Пусть AA, BB — матрицы n×nn \times n над C\mathbb {C} с A2=B2=InA^2 = B^2 = I_n и [A,B]+=0n[A,B]_+ = 0_n. Пусть z∈Cz \in \mathbb {C}. Формула Ли-Троттера задаётся как

exp⁡(z(A+B))=lim⁡p→∞(ezA/pezB/p)p. \exp (z(A+B)) = \lim _{p \to \infty } \left(e^{zA/p}e^{zB/p}\right)^p.

Вычислите ez(A+B)e^{z(A+B)}, используя правую часть. Заметьте, что

ezA/p=Incosh⁡(z/p)+Asinh⁡(z/p),ezB/p=Incosh⁡(z/p)+Bsinh⁡(z/p). e^{zA/p} = I_n\cosh (z/p) + A\sinh (z/p), \qquad e^{zB/p} = I_n\cosh (z/p) + B\sinh (z/p).
Задача 9.2.6
?
(i)

Пусть AA, BB — матрицы размера n×nn \times n такие, что AB=BAAB = BA. Покажите, что

cos⁡(A+B)≡cos⁡(A)cos⁡(B)−sin⁡(A)sin⁡(B). \cos (A+B) \equiv \cos (A)\cos (B) - \sin (A)\sin (B).
(ii)

Пусть x,y∈Rx, y \in \mathbb {R}. Известно, что sin⁡(x+y)≡sin⁡(x)cos⁡(y)+cos⁡(x)sin⁡(y)\sin (x+y) \equiv \sin (x)\cos (y) + \cos (x)\sin (y). Пусть

C=(1111),D=(0110). C = \begin{pmatrix} 1 & 1 \\ 1 & 1 \end{pmatrix}, \qquad D = \begin{pmatrix} 0 & 1 \\ 1 & 0 \end{pmatrix}.

Верно ли, что sin⁡(C+D)=sin⁡(C)cos⁡(D)+cos⁡(C)sin⁡(D)\sin (C+D) = \sin (C)\cos (D) + \cos (C)\sin (D)? Докажите или опровергните.

Задача 9.2.7

Пусть AA, BB — матрицы размера n×nn \times n над C\mathbb {C}. Преобразование Лапласа L\mathcal{L} матричной экспоненты — это резольвента

L(etA)=(λIn−A)−1. \mathcal{L}(e^{tA}) = (\lambda I_n - A)^{-1}.
?
(i)

Вычислите L(etA⊗etB)\mathcal{L}(e^{tA} \otimes e^{tB}).

(ii)

Пусть L−1\mathcal{L}^{-1} — обратное преобразование Лапласа. Найдите

L−1((λIn−A)−1⊗(λIn−B)−1). \mathcal{L}^{-1}((\lambda I_n-A)^{-1} \otimes (\lambda I_n-B)^{-1}).
Задача 9.2.8

Пусть AA — матрица размера n×nn \times n. Характеристический многочлен

det⁡(λIn−A)=λn+a1λn−1+⋯+an=p(λ) \operatorname {det}\left(\lambda I_n - A\right) = \lambda ^n + a_1\lambda ^{n-1} + \cdots + a_n = p(\lambda )

тесно связан с резольвентой (λIn−A)−1(\lambda I_n - A)^{-1} формулой

(λIn−A)−1=N1λn−1+N2λn−2+⋯+Nnλn+a1λn−1+⋯+an=N(λ)p(λ) (\lambda I_n - A)^{-1} = \frac{N_1\lambda ^{n-1}+N_2\lambda ^{n-2}+\cdots +N_n}{\lambda ^n+a_1\lambda ^{n-1}+\cdots +a_n} = \frac{N(\lambda )}{p(\lambda )}

где присоединённая матрица N(λ)N(\lambda ) является многочленом от λ\lambda степени n−1n-1 с постоянными матрицами коэффициентов N1,…,NnN_1,\ldots ,N_n размера n×nn \times n. Преобразование Лапласа матричной экспоненты — это резольвента

L(etA)=(λIn−A)−1. \mathcal{L}(e^{tA}) = (\lambda I_n - A)^{-1}.

Матрицы NkN_k и коэффициенты aka_k можно вычислить рекурсивно следующим образом

N1=In,a1=−11tr⁡(AN1) N_1 = I_n, \qquad a_1 = -\frac{1}{1}\operatorname {tr}(AN_1) N2=AN1+a1In,a2=−12tr⁡(AN2) N_2 = AN_1 + a_1I_n, \qquad a_2 = -\frac{1}{2}\operatorname {tr}(AN_2) Nn=ANn−1+an−1In,an=−1ntr⁡(ANn) N_n = AN_{n-1} + a_{n-1}I_n, \qquad a_n = -\frac{1}{n}\operatorname {tr}(AN_n) 0n=ANn+anIn. 0_n = AN_n + a_nI_n.

Покажите, что

tr⁡(L(etA))=1p(λ)dp(λ)dλ. \operatorname {tr}(\mathcal{L}(e^{tA})) = \frac{1}{p(\lambda )}\frac{dp(\lambda )}{d\lambda }.
?
Задача 9.2.9

Пусть AA — произвольная матрица размера n×nn \times n. Можно ли заключить, что exp(A∗)=(exp(A))∗\mathrm{exp}\left(A^*\right) = (\mathrm{exp}\left(A\right))^*?

?
Задача 9.2.10

Пусть AA — матрица размера n×nn \times n над C\mathbb {C}. Предположим, что AA эрмитова, т.е. A∗=AA^* = A. Тогда AA имеет только вещественные собственные значения. Предположим, что A5+A3+A=3InA^5+A^3+A = 3I_n. Покажите, что A=InA = I_n.

?
Задача 9.2.11

Пусть AA, BB — эрмитовы матрицы размера n×nn \times n. Существуют унитарные матрицы UU и VV размера n×nn \times n (зависящие от AA и BB) такие, что

exp(iA)exp(iB)=exp(iUAU−1+iVBV−1). \mathrm{exp}\left(iA\right)\mathrm{exp}\left(iB\right) = \mathrm{exp}\left(iUAU^{-1}+iVBV^{-1}\right).

Рассмотрим n=2n = 2 и

A=(0110),B=12(111−1)  ⇒  [A,B]=2(0−110). A = \begin{pmatrix} 0 & 1 \\ 1 & 0 \end{pmatrix}, \qquad B = \frac{1}{\sqrt{2}} \begin{pmatrix} 1 & 1 \\ 1 & -1 \end{pmatrix} \; \Rightarrow \; [A,B] = \sqrt{2} \begin{pmatrix} 0 & -1 \\ 1 & 0 \end{pmatrix}.

Найдите UU и VV. Заметьте, что AA и BB также унитарны.

?
Задача 9.2.12

Пусть AA, BB — матрицы размера n×nn \times n над C\mathbb {C} такие, что [A,B]=A[A,B] = A. Что можно сказать о коммутаторе [eA,eB][e^A,e^B]?

?
Задача 9.2.13
?
(i)

Пусть AA — нильпотентная матрица размера n×nn \times n. Является ли матрица cos⁡(A)\cos (A) обратимой?

(ii)

Пусть NN — нильпотентная матрица размера n×nn \times n с Nj=0nN^j = 0_n при j≥1j \geq 1. Верно ли, что

ln⁡(In−N)=−(N+12N+⋯+1j−1Nj−1) \ln (I_n-N) = -\left(N+\frac{1}{2}N+\cdots +\frac{1}{j-1}N^{j-1}\right)

и exp⁡(ln⁡(In−N))=In−N\exp (\ln (I_n-N)) = I_n-N.

Задача 9.2.14

Пусть H=H1+H2+H3H = H_1+H_2+H_3, где H1H_1, H2H_2, H3H_3 — эрмитовы матрицы размера n×nn \times n. Покажите, что

e−βH=e−βH1/2e−βH2/2e−βH3e−βH2/2e−βH1/2+(−12β)3S e^{-\beta H} = e^{-\beta H_1/2}e^{-\beta H_2/2}e^{-\beta H_3}e^{-\beta H_2/2}e^{-\beta H_1/2} + \left(-\frac{1}{2}\beta \right)^3 S

где

S=16([[H2+H3,H1],H1+2H2+2H3]+[[H3,H2],H2+H3]). S = \frac{1}{6}([[H_2+H_3,H_1],H_1+2H_2+2H_3] + [[H_3,H_2],H_2+H_3]).
?
Задача 9.2.15
?
(i)

Пусть rr, ss — целые числа с r≥2r \geq 2 и s≥2s \geq 2. Пусть AA, BB — матрицы размера n×nn \times n с Ar=0nA^r = 0_n и Bs=0nB^s = 0_n. Вычислите exp(A⊗B)\mathrm{exp}\left(A \otimes B\right).

(ii)

Пусть AA, BB — базис (заданный матрицами размера n×nn \times n) некоммутативной двумерной алгебры Ли с [A,B]=A[A,B] = A. Вычислите exp⁡(A⊗In+In⊗A+B⊗B)\exp (A \otimes I_n+I_n \otimes A+B \otimes B).

(iii)

Пусть X1X_1, X2X_2, X3X_3 — базис (заданный матрицами размера n×nn \times n) алгебры Ли с коммутационными соотношениями

[X1,X2]=0n,[X1,X3]=X1,[X2,X3]=0n. [X_1,X_2] = 0_n, \qquad [X_1,X_3] = X_1, \qquad [X_2,X_3] = 0_n.

Вычислите exp⁡(X1+X2+X3)\exp (X_1+X_2+X_3) и

exp⁡(X1⊗In⊗In+In⊗X2⊗In+In⊗In⊗X3+X1⊗X2⊗X3). \exp (X_1 \otimes I_n \otimes I_n+I_n \otimes X_2 \otimes I_n+I_n \otimes I_n \otimes X_3+X_1 \otimes X_2 \otimes X_3).
Задача 9.2.16

Рассмотрим матрицу 3×33 \times 3

A=(101010101)=I3+B,B=(001000100). A = \begin{pmatrix} 1 & 0 & 1 \\ 0 & 1 & 0 \\ 1 & 0 & 1 \end{pmatrix} = I_3+B, \qquad B = \begin{pmatrix} 0 & 0 & 1 \\ 0 & 0 & 0 \\ 1 & 0 & 0 \end{pmatrix}.

Вычислите exp(A)\mathrm{exp}\left(A\right), используя eA=eI3eBe^A = e^{I_3}e^B.

?
Задача 9.2.17

Пусть z∈Cz \in \mathbb {C}. Постройте все матрицы 2×22 \times 2 AA и BB над C\mathbb {C} такие, что

exp(zA)Bexp(−zA)=e−zB. \mathrm{exp}\left(zA\right)B\mathrm{exp}\left(-zA\right) = e^{-z}B.
?
Задача 9.2.18

Пусть AA — матрица размера n×nn \times n над C\mathbb {C}. Рассмотрим ряд Тейлора

(In+A)1/2=In+12A−12⋅4A2+1⋅32⋅4⋅6A3−⋯ (I_n+A)^{1/2} = I_n+\frac{1}{2}A-\frac{1}{2 \cdot 4}A^2+\frac{1 \cdot 3}{2 \cdot 4 \cdot 6}A^3-\cdots

и

(In+A)−1/2=In−12A+1⋅32⋅4A2−1⋅3⋅52⋅4⋅6A3+⋯ (I_n+A)^{-1/2} = I_n-\frac{1}{2}A+\frac{1 \cdot 3}{2 \cdot 4}A^2-\frac{1 \cdot 3 \cdot 5}{2 \cdot 4 \cdot 6}A^3+\cdots

Каково условие (на норму) для AA, при котором эти ряды Тейлора существуют? Можно ли применить его к матрице

A=(1111)? A = \begin{pmatrix} 1 & 1 \\ 1 & 1 \end{pmatrix}?

Заметьте, что при n=1n=1 мы имеем условие −1<A≤+1-1 < A \leq +1.

?
Задача 9.2.19

Пусть UU — унитарная матрица размера n×nn \times n. Пусть H=U+U∗H = U+U^*. Вычислите exp(zH)\mathrm{exp}\left(zH\right).

?
Задача 9.2.20

Пусть AA, BB — матрицы размера n×nn \times n над C\mathbb {C} и exp(A)exp(B)=exp(C)\mathrm{exp}\left(A\right)\mathrm{exp}\left(B\right) = \mathrm{exp}\left(C\right). Тогда матрица CC может быть представлена в виде бесконечного ряда коммутаторов AA и BB. Пусть z∈Cz \in \mathbb {C}. Запишем

exp⁡(zA)exp⁡(zB)=exp⁡(C(zA,zB)) \exp (zA)\exp (zB) = \exp (C(zA,zB))

где

C(zA,zB)=∑j=1∞cj(A,B)zj. C(zA,zB) = \sum _{j=1}^{\infty } c_j(A,B)z^j.

Покажите, что разложение до четвёртого порядка задаётся выражениями

c1(A,B)=A+B,c2(A,B)=12[A,B], c_1(A,B) = A+B, \qquad c_2(A,B) = \frac{1}{2}[A,B], c3(A,B)=112[A,[A,B]]−112[B,[A,B]],c4(A,B)=−124[A,[B,[A,B]]]. c_3(A,B) = \frac{1}{12}[A,[A,B]] - \frac{1}{12}[B,[A,B]], \qquad c_4(A,B) = -\frac{1}{24}[A,[B,[A,B]]].
?
Задача 9.2.21

Пусть A1,…,AmA_1, \ldots , A_m — матрицы размера n×nn \times n. Покажите, что

eA1⋯eAm=exp⁡(12∑j<k[Aj,Ak])exp⁡(A1+⋯+Am) e^{A_1} \cdots e^{A_m} = \exp \left(\frac{1}{2}\sum _{j<k}[A_j,A_k]\right)\exp (A_1+\cdots +A_m)

если матрицы AjA_j (j=1,…,mj=1,\ldots ,m) удовлетворяют условию

[[Aj,Ak],Aℓ]=0nдля всех j,k,ℓ∈{1,…,m}. [[A_j,A_k],A_\ell ] = 0_n \quad \text{для всех } j,k,\ell \in \left\{ 1,\ldots ,m\right\} .
?
Задача 9.2.22

Пусть AA, BB, CC — положительно полуопределённые матрицы размера n×nn \times n. Определим

A∘B:=A1/2BA1/2 A \circ B := A^{1/2}BA^{1/2}

где A1/2A^{1/2} обозначает единственный положительный квадратный корень из AA. Верно ли, что A∘(B∘C)=(A∘B)∘CA \circ (B \circ C) = (A \circ B) \circ C? Докажите или опровергните.

?
Задача 9.2.23
?
(i)

Пусть σ1\sigma_1 — спиновая матрица Паули. Покажите, что

exp⁡(−12iπ(σ1−I2))=σ1. \exp \left(-\frac{1}{2}i\pi (\sigma _1-I_2)\right) = \sigma _1.
(ii)

Найдите все матрицы 2×22 \times 2 AA и c∈Cc \in \mathbb {C} такие, что exp⁡(c(A−I2))=A\exp (c(A-I_2)) = A. Заметьте, что

exp⁡(c(A−I2))=exp⁡(cA)exp⁡(−cI2)=exp⁡(cA)(e−c00e−c). \exp (c(A-I_2)) = \exp (cA)\exp (-cI_2) = \exp (cA) \begin{pmatrix} e^{-c} & 0 \\ 0 & e^{-c} \end{pmatrix}.
Задача 9.2.24

Пусть AA, BB — вещественные симметричные положительно определённые матрицы. Покажите, что

1det⁡(A)det⁡(A+B)≤exp(tr⁡(A−1B)). \frac{1}{\operatorname {det}\left(A\right)}\operatorname {det}\left(A+B\right) \leq \mathrm{exp}\left(\operatorname {tr}(A^{-1}B)\right).
?
Задача 9.2.25

Пусть ϵ∈R\epsilon \in \mathbb {R}. Пусть In−ϵAI_n - \epsilon A — положительно определённая матрица. Покажите, что

exp⁡(tr⁡(ln⁡(In−ϵA)))=det⁡(In−ϵA) \exp (\operatorname {tr}(\ln (I_n-\epsilon A))) = \operatorname {det}\left(I_n-\epsilon A\right)

используя тождество det⁡(eM)≡exp(tr⁡(M))\operatorname {det}\left(e^M\right) \equiv \mathrm{exp}\left(\operatorname {tr}(M)\right).

?
Задача 9.2.26

Пусть v→\overrightarrow {v} — нормированный (столбцовый) вектор в Cn\mathbb {C}^{n}, и пусть AA — эрмитова матрица размера n×nn \times n. Верно ли, что v→∗exp(A)v→≥exp(v→∗Av→)\overrightarrow {v}^*\mathrm{exp}\left(A\right)\overrightarrow {v} \geq \mathrm{exp}\left(\overrightarrow {v}^*A\overrightarrow {v}\right) для всех нормированных v→\overrightarrow {v}? Докажите или опровергните.

?
Задача 9.2.27
?
(i)

Пусть ϵ∈R\epsilon \in \mathbb {R}. Пусть AA — матрица размера n×nn \times n над C\mathbb {C}. Найдите

lim⁡ϵ→0sinh⁡(2ϵA)sinh⁡(ϵ). \lim _{\epsilon \to 0} \frac{\sinh (2\epsilon A)}{\sinh (\epsilon )}.
(ii)

Предположим, что A2=InA^2 = I_n. Вычислите

sinh⁡(2ϵA)sinh⁡(ϵ). \frac{\sinh (2\epsilon A)}{\sinh (\epsilon )}.
(iii)

Предположим, что A2=0nA^2 = 0_n. Вычислите

sinh⁡(2ϵA)sinh⁡(ϵ). \frac{\sinh (2\epsilon A)}{\sinh (\epsilon )}.

Используйте разложение

sinh⁡(2ϵA)=2ϵA+(2ϵ)33!A3+⋯ . \sinh (2\epsilon A) = 2\epsilon A + \frac{(2\epsilon )^3}{3!}A^3 + \cdots .
Задача 9.2.28

Пусть α>0\alpha > 0. Рассмотрим матрицу 4×44 \times 4

A(α)=(000i/α00i00i00αi000). A(\alpha ) = \begin{pmatrix} 0 & 0 & 0 & i/\alpha \\ 0 & 0 & i & 0 \\ 0 & i & 0 & 0 \\ \alpha i & 0 & 0 & 0 \end{pmatrix}.

Покажите, что exp⁡(πA(α))=−I4\exp (\pi A(\alpha )) = -I_4.

?
Задача 9.2.29

Пусть VV — матрица 2×22 \times 2 вида V=v0I2+v1σ1+v2σ2+v3σ3V = v_0I_2+v_1\sigma_1+v_2\sigma_2+v_3\sigma_3, где v0,v1,v2,v3∈Rv_0, v_1, v_2, v_3 \in \mathbb {R}. Рассмотрим уравнение

exp⁡(iϵV)=(I2−iW)(I2+iW)−1 \exp (i\epsilon V) = (I_2-iW)(I_2+iW)^{-1}

где ϵ\epsilon вещественно. Покажите, что WW как функция от VV задаётся выражением

W=−tan⁡(12ϵ(v0I2+v1σ1+v2σ2+v3σ3)). W = -\tan \left(\frac{1}{2}\epsilon (v_0I_2+v_1\sigma _1+v_2\sigma _2+v_3\sigma _3)\right).
?
Задача 9.2.30

Пусть AA, BB — матрицы размера n×nn \times n. Покажите, что

eA⊗In+In⊗AeB⊗In+In⊗B=(eA⊗eB)(eB⊗eB)=(eAeB)⊗(eAeB). e^{A \otimes I_n+I_n \otimes A}e^{B \otimes I_n+I_n \otimes B} = (e^A \otimes e^B)(e^B \otimes e^B) = (e^Ae^B) \otimes (e^Ae^B).
?
Задача 9.2.31

Пусть α∈R\alpha \in \mathbb {R} и Π\Pi — матрица проектирования размера n×nn \times n. Пусть Q=Π1/2Q = \Pi^{1/2} — квадратный корень из Π\Pi, т.е. Q2=ΠQ^2 = \Pi. Покажите, что

U(α)=exp⁡(iαΠ1/2)=In+iαΠ1/2+(cos⁡(α)−1)Π+i(sin⁡(α)−α)ΠΠ1/2. U(\alpha ) = \exp (i\alpha \Pi ^{1/2}) = I_n + i\alpha \Pi ^{1/2} + (\cos (\alpha )-1)\Pi + i(\sin (\alpha )-\alpha )\Pi \Pi ^{1/2}.
?
Задача 9.2.32

Пусть AA — матрица размера n×nn \times n над R\mathbb {R}. Тогда имеет место разложение Тейлора

sin⁡(A):=∑k=0∞(−1)k(2k+1)!A2k+1,cos⁡(A):=∑k=0∞(−1)k(2k)!A2k. \sin (A) := \sum _{k=0}^{\infty } \frac{(-1)^k}{(2k+1)!}A^{2k+1}, \qquad \cos (A) := \sum _{k=0}^{\infty } \frac{(-1)^k}{(2k)!}A^{2k}.

Вычислять sin⁡(A)\sin (A) и cos⁡(A)\cos (A) по усечённому ряду Тейлора имеет смысл только вблизи начала координат. Можно использовать многократное применение формулы двойного угла

cos⁡(2A)≡2cos⁡2(A)−In,sin⁡(2A)≡2sin⁡(A)cos⁡(A). \cos (2A) \equiv 2\cos ^2(A) - I_n, \qquad \sin (2A) \equiv 2\sin (A)\cos (A).

Найти sin⁡(A)\sin (A) и cos⁡(A)\cos (A) для матрицы AA по подходяще усечённым приближениям ряда Тейлора можно следующим образом

S0=приближение Тейлора для sin⁡(A/2k),C0=приближение Тейлора для cos⁡(A/2k) S_0 = \text{приближение Тейлора для } \sin (A/2^k), \qquad C_0 = \text{приближение Тейлора для } \cos (A/2^k)

и рекурсии

Sj=2Sj−1Cj−1,Cj=2Cj−12−In S_j = 2S_{j-1}C_{j-1}, \qquad C_j = 2C_{j-1}^2 - I_n

где j=1,2,…j = 1,2,\ldots. Здесь kk — положительное целое число, выбранное так, что, скажем, ∥A∥∞≈2k\left\| A\right\|_{\infty } \approx 2^k. Примените эту рекурсию для вычисления синуса и косинуса матрицы 2×22 \times 2

A=(2112). A = \begin{pmatrix} 2 & 1 \\ 1 & 2 \end{pmatrix}.

Используйте k=2k=2.

?
Задача 9.2.33

Пусть TT — матрица 2×22 \times 2

T=13(1111). T = \frac{1}{3} \begin{pmatrix} 1 & 1 \\ 1 & 1 \end{pmatrix}.

Вычислите ln⁡(det⁡(I2−T))\ln (\operatorname {det}\left(I_2-T\right)), используя правую часть тождества

ln⁡(det⁡(I2−T))≡−∑k=1∞1ktr⁡(Tk). \ln (\operatorname {det}\left(I_2-T\right)) \equiv -\sum _{k=1}^{\infty } \frac{1}{k}\operatorname {tr}(T^k).
?