Edit page

Sistemas de Equações Lineares com Coeficientes Constantes

Definição

Equação do tipo:

 ⁣dy ⁣dt=Ay+B(t)\frac{\d y}{\d t} = A y + B(t)

em que

  • AA é uma matrix quadrada n×nn \times n
  • y(t)y(t) é uma função RRn\R \to \R^n
  • B(t)B(t) é uma função continua de RRn\R \to \R^n

Se B(t)0B(t) \equiv 0 a equação designa-se por homogénea.

Se tivermos, por exemplo, o sistema, que podemos facilmente resolver linha a linha:

{ ⁣dx ⁣dt=x ⁣dy ⁣dt=x+yet ⁣dz ⁣dt=x+y+z2et{x=ety=etz=et\begin{cases} \frac{\d x}{\d t} = x\\ \frac{\d y}{\d t} = x + y - e^t\\ \frac{\d z}{\d t} = x + y + z - 2e^t \end{cases} \Leftrightarrow \begin{cases} x = e^t\\ y = e^t\\ z = e^t \end{cases}

Podemos reescrever o sistema acima para obter a forma de equação linear com coeficientes constantes.

 ⁣dw ⁣dt=[xx+yetx+y+z2et]=[xx+yx+y+z]+[0et2et]=[100110111][xyz]+[0et2et]=[100110111]w+[0et2et]\begin{aligned} \frac{\d w}{\d t} &= \begin{bmatrix} x\\ x+y-e^t\\ x+y+z-2e^t \end{bmatrix}\\ & = \begin{bmatrix} x\\ x+y\\ x+y+z \end{bmatrix} + \begin{bmatrix} 0\\ -e^t\\ -2e^t \end{bmatrix}\\ & = \begin{bmatrix} 1 & 0 & 0\\ 1 & 1 & 0\\ 1 & 1 & 1 \end{bmatrix} \begin{bmatrix} x\\ y\\ z \end{bmatrix} + \begin{bmatrix} 0\\ -e^t\\ -2e^t \end{bmatrix}\\ & = \begin{bmatrix} 1 & 0 & 0\\ 1 & 1 & 0\\ 1 & 1 & 1 \end{bmatrix} w + \begin{bmatrix} 0\\ -e^t\\ -2e^t \end{bmatrix} \end{aligned}

Equações Homogéneas

Começando então pelo caso em que B(t)0B(t) \equiv 0, temos que:

 ⁣dy ⁣dt=Ay\frac{\d y}{\d t} = A y

com uma matrix AA quadrada n×nn \times n.

Para n=1n=1 sabemos que y=keAty = k e^{At}. Para n>1n>1 vamos ter exatamente a mesma coisa, mas para isso temos de definir o exponencial de uma matriz.

Exponencial de matrizes

Por definição, tal como em dimensão 1, onde

eα=n=0+αnn!e^\alpha = \sum_{n = 0}^{+\infty} \frac{\alpha^n}{n!}

passamos a ter, para uma matrix AA quadrada n×nn \times n:

eA=Id+A+12A2+13!A3+=k=0+Akk!e^A = Id + A + \frac{1}{2} A^2 + \frac{1}{3!} A^3 + \dots = \sum_{k=0}^{+\infty} \frac{A^k}{k!}

Numa explicação mais teórica, dizemos que

eA=limN+k=0N1k!Ake^A = \lim_{N \to +\infty} \sum_{k=0}^{N} \frac{1}{k!}A^k

e poderíamos concluir que este limite é absolutamente convergente (carece de demonstração).

Exemplo pela definição

Queremos calcular o valor de eAte^{At} da matriz

A=[2100]A = \begin{bmatrix} 2 & 1\\ 0 & 0 \end{bmatrix}

Vamos então tentar descobrir AkA^k, para podemos utilizar a definição.

eA=Id+A+12A2+13!A3+e^A = Id + A + \frac{1}{2} A^2 + \frac{1}{3!} A^3 + \dots
A2=[2100][2100]=[4200]A^2 = \begin{bmatrix} 2 & 1\\ 0 & 0 \end{bmatrix} \begin{bmatrix} 2 & 1\\ 0 & 0 \end{bmatrix} = \begin{bmatrix} 4 & 2\\ 0 & 0 \end{bmatrix}
A3=[4200][2100]=[8400]A^3 = \begin{bmatrix} 4 & 2\\ 0 & 0 \end{bmatrix} \begin{bmatrix} 2 & 1\\ 0 & 0 \end{bmatrix} = \begin{bmatrix} 8 & 4\\ 0 & 0 \end{bmatrix}

Podemos então concluir que AkA^k:

Ak=[2k2k100],k1A^k = \begin{bmatrix} 2^k & 2^{k-1}\\ 0 & 0 \end{bmatrix}, k \geq 1

É possível verificar este resultado por indução matemática, mas tal prova foi omitida.

Assim temos que:

eAt=Id+k=1+Aktkk!=[1001]+k=1+1k![2ktktk2k100]=[1001]+[k=1+2ktkk!k=1+2k1tkk!00]=[k=0+2ktkk!12k=1+2ktkk!01]=[e2t12(e2t1)01]\begin{aligned} e^{At} &= Id + \sum_{k=1}^{+\infty} \frac{A^k t^k}{k!}\\ &= \begin{bmatrix} 1 & 0\\ 0 & 1 \end{bmatrix} + \sum_{k=1}^{+\infty} \frac{1}{k!} \begin{bmatrix} 2^k t^k & t^k 2^{k-1}\\ 0 & 0 \end{bmatrix}\\ &= \begin{bmatrix} 1 & 0\\ 0 & 1 \end{bmatrix} + \begin{bmatrix} \sum_{k=1}^{+\infty} \frac{2^k t^k}{k!} & \sum_{k=1}^{+\infty} \frac{2^{k-1}t^k}{k!}\\ 0 & 0 \end{bmatrix}\\ &=\begin{bmatrix} \sum_{k=0}^{+\infty} \frac{2^k t^k}{k!} & \frac{1}{2}\sum_{k=1}^{+\infty} \frac{2^k t^k}{k!}\\ 0 & 1 \end{bmatrix}\\ &= \begin{bmatrix} e^{2t} & \frac{1}{2} (e^{2t} - 1)\\ 0 & 1 \end{bmatrix} \end{aligned}

Proposição

Se A=SJS1A = S J S^{-1}, em que A,S,JA, S, J são matrizes n×nn \times n e SS é invertível (detS0\det S \ne 0), então

eA=SeJS1e^A = S e^J S^{-1}
Demonstração
A2=AA=SJS1SJS1=SJJS1=SJ2S1A3=A2A=SJ2S1SJS1=SJ2JS1=SJ3S1Ak=SJkS1\begin{darray}{l} A^2 = A A = S J S^{-1} S J S^{-1} = S J J S^{-1} = S J^2 S^{-1}\\ A^3 = A^2 A = S J^2 S^{-1} S J S^{-1} = S J^2 J S^{-1} = S J^3 S^{-1}\\ A^k = S J^k S^{-1} \end{darray}

Temos assim,

eA=k=0+1k!Ak=k=0+1k!SJkS1=Sk=0+1k!JkS1=SeJS1\begin{aligned} e^A &= \sum_{k=0}^{+\infty} \frac{1}{k!}A^k\\ &= \sum_{k=0}^{+\infty} \frac{1}{k!}SJ^kS^{-1}\\ &= S \sum_{k=0}^{+\infty} \frac{1}{k!} J^k S^{-1}\\ &= S e^J S^{-1} \end{aligned}

Esta proposição irá ser útil porque se repararmos, isto significa que se tivermos uma matriz diagonal JJ, ou seja, AA for diagonalizável, podemos facilmente descobrir o exponencial de AA devido às seguinte duas propriedades:

Ak=[λ1k000λ2k00000λnk]etA=[eλ1t000eλ2t00000eλnt]\begin{darray}{cc} A^k = \begin{bmatrix} \lambda_1^k & 0 & \dots & 0\\ 0 & \lambda_2^k & \dots & 0\\ \vdots & \vdots & \ddots & 0\\ 0 & 0 & 0 & \lambda_n^k \end{bmatrix} & e^{tA} = \begin{bmatrix} e^{\lambda_1 t} & 0 & \dots & 0\\ 0 & e^{\lambda_2 t} & \dots & 0\\ \vdots & \vdots & \ddots & 0\\ 0 & 0 & 0 & e^{\lambda_n t} \end{bmatrix} \end{darray}

Se quisermos verificar que determinámos bem a exponencial de uma matrix AA, podemos usar os seguintes critérios de verificação:

  1. etAt=0=Ide^{tA} \big|_{t=0} = Id
  2.  ⁣detA ⁣dtt0=AetAt=0=A\frac{\d e^{tA}}{\d t} \big|_{t_0} = Ae^{tA} \big|_{t=0} = A
Exemplo 1

Determinar exponencial da matriz AA:

A=[3113]A = \begin{bmatrix} 3 & 1\\ 1 & 3 \end{bmatrix}

Determinar valores e vetores próprios:

det(AλI)=3λ113λ=(3λ)21\det(A - \lambda I) = \begin{vmatrix} 3 - \lambda & 1\\ 1 & 3 - \lambda \end{vmatrix} = (3-\lambda) ^2 - 1

Após calcular as raízes do polinómio, temos assim valores próprios: λ=2\lambda = 2 e λ=4\lambda = 4.
Vamos calcular agora os vetores próprios.

  • Para λ=2\lambda = 2:

    Nul(A2I)=Nul[1111]=L{(1,1)}\operatorname{Nul} (A - 2I) = \operatorname{Nul} \begin{bmatrix} 1 & 1\\ 1 & 1 \end{bmatrix} = \mathcal{L} \{(1,-1)\}
  • Para λ=4\lambda = 4:

    Nul(A4I)=Nul[1111]=L{(1,1)}\operatorname{Nul} (A - 4I) = \operatorname{Nul} \begin{bmatrix} -1 & 1\\ 1 & -1 \end{bmatrix} = \mathcal{L} \{(1,1)\}

Temos assim que A=SΛS1A = S \Lambda S^{-1} com

S=[1111]Λ=[2004]\begin{darray}{cc} S = \begin{bmatrix} 1 & 1\\ -1 & 1 \end{bmatrix} & \Lambda = \begin{bmatrix} 2 & 0\\ 0 & 4 \end{bmatrix} \end{darray}

A matriz SS corresponde aos vetores próprios, e a matrix Λ\Lambda aos valores próprios.

Então,

etA=SetΛS1=[1111][e2t00e4t]S1=[e2te4te2te4t]12[1111]=[e2t+e4t2e4te2t2e4te2t2e4t+e2t2]\begin{aligned} e^{tA} &= S e ^{t\Lambda} S^{-1}\\ &= \begin{bmatrix} 1 & 1\\ -1 & 1 \end{bmatrix} \begin{bmatrix} e^{2t} & 0\\ 0 & e^{4t} \end{bmatrix} S^{-1}\\ &= \begin{bmatrix} e^{2t} & e^{4t}\\ -e^{2t} & e^{4t} \end{bmatrix} \frac{1}{2} \begin{bmatrix} 1 & -1\\ 1 & 1 \end{bmatrix}\\ &= \begin{bmatrix} \frac{e^{2t} + e^{4t}}{2} & \frac{e^{4t} - e^{2t}}{2}\\ \frac{e^{4t} - e^{2t}}{2} & \frac{e^{4t} + e^{2t}}{2} \end{bmatrix} \end{aligned}

Vamos agora testar os critérios de verificação:

[e0+e02e0e02e0e02e0+e02]=[1001]=Id\begin{bmatrix} \frac{e^{0} + e^{0}}{2} & \frac{e^{0} - e^{0}}{2}\\ \frac{e^{0} - e^{0}}{2} & \frac{e^{0} + e^{0}}{2} \end{bmatrix} = \begin{bmatrix} 1 & 0\\ 0 & 1 \end{bmatrix} = Id
[2e2t+4e4t24e4t2e2t24e4t2e2t22e2t+4e4t2]t=0[3113]=A\begin{bmatrix} \frac{2e^{2t} + 4e^{4t}}{2} & \frac{4e^{4t} - 2e^2t}{2}\\ \frac{4e^{4t} - 2e^2t}{2} & \frac{2e^{2t} + 4e^{4t}}{2} \end{bmatrix} \underset{t=0}{\longrightarrow} \begin{bmatrix} 3 & 1\\ 1 & 3 \end{bmatrix} = A
Exemplo 2

Determinar exponencial da matriz AA:

A=[1111]A = \begin{bmatrix} 1 & -1\\ 1 & 1 \end{bmatrix}

Calculemos agora os valores e vetores próprios de AA:

det(AλI)=(1λ)2+1λ=1+iλ=1i\begin{darray}{c} \det(A - \lambda I) = (1-\lambda) ^2 + 1\\ \begin{darray}{cc} \lambda = 1 + i & \lambda = 1 - i \end{darray}\\ \end{darray}
  • Para λ=1+i\lambda = 1 + i:

    Nul(A(1+i)I)=Nul[i11i]=L{(1,i)}\operatorname{Nul} (A - (1 + i)I) = \operatorname{Nul} \begin{bmatrix} -i & -1\\ 1 & -i \end{bmatrix} = \mathcal{L} \{(1, -i)\}
  • Para λ=1i\lambda = 1 - i:

    Nul(A(1i)I)=L{(1,i)}\operatorname{Nul} (A - (1 - i)I) = \mathcal{L} \{(1, i)\}

Então, temos SS e Λ\Lambda que satisfazem SΛS1S\Lambda S^{-1}:

S=[11ii]Λ=[1+i001i]\begin{darray}{cc} S = \begin{bmatrix} 1 & 1\\ -i & i \end{bmatrix} & \Lambda = \begin{bmatrix} 1 + i & 0\\ 0 & 1 - i \end{bmatrix} \end{darray}

E finalmente:

etA=SetΛS1=[11ii][et+it00etit]12i[i1i1]=et2[11ii][eit00eit][1i11]==et2[eit+eitieitieitieit+ieiteit+eit]==et[costsintsintcost]\begin{aligned} e^{tA} &= S e^{t\Lambda} S^{-1}\\ &= \begin{bmatrix} 1 & 1\\ -i & i \end{bmatrix} \begin{bmatrix} e^{t+it} & 0\\ 0 & e^{t - it} \end{bmatrix} \frac{1}{2i} \begin{bmatrix} i & -1\\ i & 1 \end{bmatrix}\\ &= \frac{e^t}{2} \begin{bmatrix} 1 & 1\\ -i & i \end{bmatrix} \begin{bmatrix} e^{it} & 0\\ 0 & e^{-it} \end{bmatrix} \begin{bmatrix} 1 & i\\ 1 & -1 \end{bmatrix}\\ &= \dots\\ &= \frac{e^t}{2} \begin{bmatrix} e^{it} + e^{-it} & i e^{it} - i e^{-it}\\ - i e^{it} + i e^{-it} & e^{it} + e^{-it} \end{bmatrix}\\ &= \dots\\ &= e^t \begin{bmatrix} \cos t & -\sin t\\ \sin t & \cos t \end{bmatrix} \end{aligned}

Matrizes por Blocos

Por vezes as matrizes não são diagonalizável. Vamos ver mais à frente como podemos usar matrizes divididas em blocos para determinar as exponenciais destas matrizes.

Uma matriz está dividida em blocos, se tiver blocos n×nn \times n na sua diagonal, em que todos os outros valores fora da diagonal são 00. Os blocos podem ser de qualquer tamanho, inclusive 1×11 \times 1.

etA=[etB1000000000etB200000000000000000000000000000etBj]\def\marray#1{\hspace{-5pt}\begin{array}{c}#1\end{array}\hspace{-5pt}} e^{tA} = \begin{bmatrix} \boxed{e^{tB_1}} & \marray{0&0&0} & 0 & \dots &\marray{0&0}\\ \marray{0\\ 0\\ 0} & \boxed{\marray{&&\\ &e^{tB_2}&\\ &&}} & \marray{0\\0\\0} & \marray{\dots\\\dots\\\dots} & \marray{0&0\\0&0\\0&0}\\ 0 & \marray{0&0&0} & \ddots & \dots & \marray{0&0}\\ \vdots & \marray{\vdots & \vdots & \vdots} & \vdots & \ddots & \marray{0& 0}\\ \marray{0\\0} & \marray{0&0&0\\0&0&0} & \marray{0\\0} & \marray{0\\0} & \boxed{\marray{e^{tB_j} &\\&}} \end{bmatrix}

Abaixo está um exemplo de uma matriz com três blocos:

etA=[200000000031432520000000000000000110]\def\marray#1{\hspace{-5pt}\begin{array}{c}#1\end{array}\hspace{-5pt}} e^{tA} = \begin{bmatrix} \boxed{2} & \marray{0&0&0} &\marray{0&0}\\ \marray{0\\ 0\\ 0} & \boxed{\marray{0&3&1\\4&-3&2\\5&2&0}} & \marray{0&0\\0&0\\0&0}\\ \marray{0\\0} & \marray{0&0&0\\0&0&0} & \boxed{\marray{0&1\\1&0}} \end{bmatrix}
Exemplo 3

Vendo o exemplo 1 e 2 do tópico anterior, podemos fazer uma matriz que tem dois blocos, um com cada uma das matrizes dos exemplos.

A=[1100110000310013]A = \begin{bmatrix} 1 & -1 & 0 & 0\\ 1 & 1 & 0 & 0\\ 0 & 0 & 3 & 1\\ 0 & 0 & 1 & 3 \end{bmatrix}

Assim, temos que

etA=[etcostetsint00etsintetcost0000e2t+e4t2e4te2t200e4te2t2e4t+e2t2]e^{tA} = \begin{bmatrix} e^t \cos t & -e^t \sin t & 0 & 0\\ e^t \sin t & e^t \cos t & 0 & 0\\ 0 & 0 & \frac{e^{2t} + e^{4t}}{2} & \frac{e^{4t} - e^{2t}}{2}\\ 0 & 0 & \frac{e^{4t} - e^{2t}}{2} & \frac{e^{4t} + e^{2t}}{2} \end{bmatrix}

Blocos de Jordan

Dá-se o nome de bloco de Jordan de dimensão nn a uma matriz JλJ_\lambda quadrada n×nn \times n da forma:

Jλ=[λ10000λ10000λ00000λ10000λ]J_\lambda = \begin{bmatrix} \lambda & 1 & 0 & \cdots & 0 & 0\\ 0 & \lambda & 1 & \cdots & 0 & 0\\ 0 & 0 & \lambda & \cdots & 0 & 0\\ \vdots & \vdots & \vdots & \ddots & \vdots & \vdots\\ 0 & 0 & 0 & \cdots & \lambda & 1\\ 0 & 0 & 0 & \cdots & 0 & \lambda \end{bmatrix}

Portanto

Jλ=[ju,j]i,k=1,,ncomji,k={λse i=k1se i=k10nos restantes casos\begin{darray}{c} J_\lambda = [j_{u,j}]_{i,k = 1,\dots,n} & \text{com} & j_{i,k} = \begin{cases} \lambda & \text{se } i = k\\ 1 & \text{se } i=k-1\\ 0 & \text{nos restantes casos} \end{cases} \end{darray}
Exemplos

[3]\begin{bmatrix}3\end{bmatrix} é um bloco de Jordan de dimensão 11 (J3J_3)


[410041004]\begin{bmatrix}4&1&0\\0&4&1\\0&0&4\end{bmatrix} é um bloco de Jordan de dimensão 33 (J4J_4)


[0100]\begin{bmatrix}0&1\\0&0\end{bmatrix} é um bloco de Jordan de dimensão 22 (J0J_0)

É de realçar que um bloco de Jordan tem um único valor próprio e a dimensão do espaço próprio é unitária, isto é, tem apenas um vetor próprio linearmente independente.

Matrizes na forma canónica de Jordan

Uma matriz quadrada JJ é uma matriz na forma canónica de Jordan se é formada exclusivamente por blocos de Jordan sobre a diagonal:

J=[Jλ1000000000Jλ200000000000000000000000000000Jλj]\def\marray#1{\hspace{-5pt}\begin{array}{c}#1\end{array}\hspace{-5pt}} J = \begin{bmatrix} \boxed{J_{\lambda_1}} & \marray{0&0&0} & 0 & \dots &\marray{0&0}\\ \marray{0\\ 0\\ 0} & \boxed{\marray{&&\\ & J_{\lambda_2}&\\ &&}} & \marray{0\\0\\0} & \marray{\dots\\\dots\\\dots} & \marray{0&0\\0&0\\0&0}\\ 0 & \marray{0&0&0} & \ddots & \dots & \marray{0&0}\\ \vdots & \marray{\vdots & \vdots & \vdots} & \vdots & \ddots & \marray{0& 0}\\ \marray{0\\0} & \marray{0&0&0\\0&0&0} & \marray{0\\0} & \marray{0\\0} & \boxed{\marray{J_{\lambda_j} &\\&}} \end{bmatrix}

onde Jλ1,Jλ2,,JλjJ_{\lambda_1}, J_{\lambda_2}, \dots ,J_{\lambda_j} são blocos de Jordan.

Exemplos

A matriz [1000001000000002]\begin{bmatrix}1&0&0&0\\0&0&1&0\\0&0&0&0\\0&0&0&2\end{bmatrix} está na forma canónica de Jordan; é formada por 3 blocos de Jordan.


A matriz [2000020000210002]\begin{bmatrix}2&0&0&0\\0&2&0&0\\0&0&2&1\\0&0&0&2\end{bmatrix} está na forma canónica de Jordan; é formada por 3 blocos de Jordan.


A matriz [3000041000410002]\begin{bmatrix}3&0&0&0\\0&4&1&0\\0&0&4&1\\0&0&0&2\end{bmatrix} não está na forma canónica de Jordan.


A matriz [3000041000400002]\begin{bmatrix}3&0&0&0\\0&4&1&0\\0&0&4&0\\0&0&0&2\end{bmatrix} está na forma canónica de Jordan; é formada por 3 blocos de Jordan.


A matriz [3000041000410004]\begin{bmatrix}3&0&0&0\\0&4&1&0\\0&0&4&1\\0&0&0&4\end{bmatrix} está na forma canónica de Jordan; é formada por 2 blocos de Jordan.


A matriz [1201]\begin{bmatrix}1&2\\0&1\end{bmatrix} não está na forma canónica de Jordan.


A matriz [4100040000410004]\begin{bmatrix}4&1&0&0\\0&4&0&0\\0&0&4&1\\0&0&0&4\end{bmatrix} está na forma canónica de Jordan; é formada por 2 blocos de Jordan.


A matriz [1000000100000100000000001]\begin{bmatrix}1&0&0&0&0\\0&0&1&0&0\\0&0&0&1&0\\0&0&0&0&0\\0&0&0&0&1\end{bmatrix} está na forma canónica de Jordan; é formada por 3 blocos de Jordan.


A matriz [0010000000010000]\begin{bmatrix}0&0&1&0\\0&0&0&0\\0&0&0&1\\0&0&0&0\end{bmatrix} não está na forma canónica de Jordan.

Teorema

Teorema de Jordan

Seja AA uma matriz quadrada qualquer então existe SS tal que

A=SJS1A = S J S^{-1}

onde JJ está na forma canónica de Jordan formada por jj blocos de Jordan Jλ1,Jλ2,,JλjJ_{\lambda_1}, J_{\lambda_2}, \dots , J_{\lambda_j} de dimensão n1,n2,,njn_1, n_2, \dots , n_j respetivamente, onde jj é o número máximo de vetores próprios de AA linearmente independentes (λ1,λ2,,λj\lambda_1, \lambda_2, \dots, \lambda_j são valores próprios de AA).

É também importante observar que:

det(Aλ)=det(Jλ)=(λ1λ)n1(λ2λ)n2(λjλ)nj\det(A-\lambda) = \det(J - \lambda) = (\lambda_1 - \lambda)^{n_1} (\lambda_2 - \lambda)^{n_2} \dots (\lambda_j - \lambda)^{n_j}

Nota: Pode acontecer termos λa=λb\lambda_a = \lambda_b com aba \ne b.

Exponencial de Blocos de Jordan

Para um bloco de Jordan JλJ_\lambda de dimensão nn temos

etJλ=[eλtteλtt22!eλttn1(n1)!eλt0eλtteλt00eλtt22!eλtteλt000eλt]=eλt[1tt22!tn1(n1)!01t001t22!t0001]\begin{aligned} e^{tJ_\lambda} &= \begin{bmatrix} e^{\lambda t} & t e^{\lambda t} & \frac{t^2}{2!} e^{\lambda t} & \dots & \frac{t^{n-1}}{(n-1)!} e^{\lambda t}\\ 0 & e^{\lambda t} & t e^{\lambda t} & \ddots & \vdots\\ 0 & 0 & e^{\lambda t} & \ddots & \frac{t^2}{2!} e^{\lambda t}\\ \vdots & \vdots & \ddots & \ddots & t e^{\lambda t}\\ 0 & 0 & \dots & 0 & e^{\lambda t} \end{bmatrix}\\ &= e^{\lambda t} \begin{bmatrix} 1 & t & \frac{t^2}{2!} & \dots & \frac{t^{n-1}}{(n-1)!} \\ 0 & 1 & t & \ddots & \vdots\\ 0 & 0 & 1 & \ddots & \frac{t^2}{2!} \\ \vdots & \vdots & \ddots & \ddots & t \\ 0 & 0 & \dots & 0 & 1 \end{bmatrix} \end{aligned}

Funções Matriciais

Seja CC uma matrix m×nm \times n em que cada entrada é uma função escalar, e cijc_{ij} são as entradas de CC.

Podemos facilmente definir tanto a derivada como a primitiva de CC:

  • Derivada:

     ⁣dC ⁣dt=[ ⁣dcij ⁣dt]\frac{\d C}{\d t} = \begin{bmatrix}\frac{\d c_{ij}}{\d t}\end{bmatrix}

    Por exemplo:

    C=[tt2etsint] ⁣dC ⁣dt=[12tetcost]\begin{darray}{cc} C = \begin{bmatrix} t & t^2\\ e^t & \sin t \end{bmatrix} & \frac{\d C}{\d t} = \begin{bmatrix} 1 & 2t\\ e^t & \cos t \end{bmatrix} \end{darray}
  • Primitiva:

    abC ⁣dt=[abcij ⁣dt]\int_{a}^{b} C \d t = \begin{bmatrix}\int_{a}^{b} c_{ij} \d t\end{bmatrix}

Proposição

Se CC é uma matrix m×pm \times p e DD uma matrix p×np \times n, então

 ⁣d ⁣dt(CD)= ⁣dC ⁣dtD+C ⁣dD ⁣dt\frac{\d }{\d t} (CD) = \frac{\d C}{\d t}D + C \frac{\d D}{\d t}

Teorema

Seja AA uma matriz quadrada (constante). Então, a função matricial E(t)=etAE(t) = e^{tA} satisfaz as seguintes propriedades:

  •  ⁣dE(t) ⁣dt=AE(t)=E(t)A\frac{\d E(t)}{\d t} = AE(t) = E(t) A
  • E(0)=IdE(0) = Id
E(t)=Id+tA+t22A2+t33!+E(t) = Id + tA + \frac{t^2}{2}A^2 + \frac{t^3}{3!} + \dots
Demonstração
 ⁣dE(t) ⁣dt= ⁣d ⁣dt(k=0+tkk!Ak)=k=1+ktk1k!Ak=k=1+tk1(k1)!Ak=k=0+tkk!Ak+1=limn+(k=0Ntkk!Ak+1)=limn+(k=0Ntkk!Ak)A=etAA=AetA\begin{aligned} \frac{\d E(t)}{\d t} &= \frac{\d }{\d t} \left(\sum_{k=0}^{+\infty} \frac{t^k}{k!}A^k\right)\\ &=\sum_{k=1}^{+\infty} \frac{kt^{k-1}}{k!}A^k\\ &=\sum_{k=1}^{+\infty} \frac{t^{k-1}}{(k-1)!} A^k\\ &=\sum_{k=0}^{+\infty} \frac{t^k}{k!} A^{k+1}\\ &=\lim_{n \to +\infty} \left(\sum_{k=0}^{N} \frac{t^k}{k!} A^{k+1}\right)\\ &=\lim_{n \to +\infty} \left(\sum_{k=0}^{N} \frac{t^k}{k!} A^{k}\right)A\\ &=e^{tA}A = A e^{tA} \end{aligned}

Um corolário deste teorema é que

(eA)1=eA(e^A)^{-1} = e^{-A}

Resolução de Equações Homogéneas

Teorema

Seja AA matriz n×nn\times n (constante) e y0Rny_0 \in \R^n, t0Rt_0 \in \R, então o PVI (problema de valor inicial)

 ⁣dy ⁣dt=Ay com y(t0)=y0\frac{\d y}{\d t} = Ay \text{ com } y(t_0) = y_0

tem uma única solução y(t)y(t) que pode ser calculada por

y(t)=e(tt0)Ay0y(t) = e^{(t-t_0) A} y_0

Corolário: Sejam AA e BB duas matrizes n×nn\times n tais que AB=BAAB = BA, então eA+B=eAeBe^{A+B} = e^A e^B

Exemplo

Tomando agora como exemplo o seguinte sistema:

{ ⁣dx ⁣dt=2x+y ⁣dy ⁣dt=0 com y(t0)=[x0y0]\begin{cases} \frac{\d x}{\d t} = 2x + y\\ \frac{\d y}{\d t} = 0 \end{cases} \text{ com } y(t_0) = \begin{bmatrix} x_0\\ y_0 \end{bmatrix}

Podemos colocar agora na forma  ⁣dy ⁣dt=Ay\frac{\d y}{\d t} = Ay, com

y=[xy]A=[2100]\begin{darray}{c} y = \begin{bmatrix} x\\ y \end{bmatrix}& A = \begin{bmatrix} 2 & 1\\ 0 & 0 \end{bmatrix} \end{darray}

Sabemos, de um exemplo anterior, que

etA=[e2t12(e2t1)01]e^{tA}=\begin{bmatrix} e^{2t} & \frac{1}{2} (e^{2t} - 1)\\ 0 & 1 \end{bmatrix}

Então, pelo teorema acima, podemos concluir que

[x(t)y(t)]=e(tt0)A[x0y0]=[e2(tt0)12(e2(tt0)1)01][x0y0]\begin{bmatrix} x(t)\\ y(t) \end{bmatrix} = e^{(t-t_0) A} \begin{bmatrix} x_0\\ y_0 \end{bmatrix} = \begin{bmatrix} e^{2(t-t_0)} & \frac{1}{2} (e^{2(t-t_0)} - 1)\\ 0 & 1 \end{bmatrix} \begin{bmatrix} x_0\\ y_0 \end{bmatrix}
x(t)=x0e2(tt0)+y02(e2(tt0)1)y(t)=y0\begin{darray}{l} x(t) = x_0 e^{2(t-t_0)} + \frac{y_0}{2}\left(e^{2(t-t_0)} - 1 \right)\\ y(t) = y_0 \end{darray}

Teorema

Seja AA uma matrix n×nn \times n (constante). Então as soluções de

 ⁣dy ⁣dt=Ay\frac{\d y}{\d t} = Ay

constituem um espaço linear de dimensão nn.

Por outras palavras, este teorema diz que para uma equação da forma acima, com uma matrix n×nn \times n, existem nn soluções linearmente independentes.
Por exemplo, caso AA seja uma matrix 2×22 \times 2, iremos ter duas soluções linearmente independentes (isto é, como veremos abaixo, iremos encontrar dois vetores próprios).

Proposição

Um vetor vv é um vetor próprio de AA e corresponde ao valor próprio λ\lambda se e só se y(t)=etλvy(t) = e^{t \lambda} v é solução de  ⁣dy ⁣dt=Ay\frac{\d y}{\d t} = Ay.

Com esta proposição, podemos concluir que é possível resolver equações do tipo  ⁣dy ⁣dt=Ay\frac{\d y}{\d t} = Ay, se descobrirmos os vetores e valores próprios da matriz AA.

Exemplo

O exemplo seguinte corresponde ao exercício 2.1.2 a) da ficha de 2021/2022

Resolver a seguinte equação diferencial:

x=[5131]x,x(0)=[21]\begin{darray}{ccc} x' = \begin{bmatrix} 5 & -1\\ 3 & 1 \end{bmatrix} x &, & x(0) = \begin{bmatrix} 2\\ -1 \end{bmatrix} \end{darray}

Começamos então por determinar os valores e vetores próprios de AA:

det(AλI)=5λ131λ=(5λ)(1λ)+3=(λ2)(λ4)\det (A - \lambda I) = \begin{vmatrix} 5 - \lambda & -1\\ 3 & 1 - \lambda \end{vmatrix} = (5 - \lambda)(1 - \lambda) + 3 = (\lambda-2)(\lambda-4)

Temos então dois valores próprios: λ=2\lambda = 2 e λ=4\lambda = 4.

  • Para λ=2\lambda = 2:

    Temos então que

    A2I=[3131]A - 2I = \begin{bmatrix} 3 & -1\\ 3 & -1 \end{bmatrix}

    Então, os vetores próprios são os vetores que pertencem ao espaço nulo da matriz:

    Nul[3131]=L{(1,3)}\operatorname{Nul} \begin{bmatrix} 3 & -1\\ 3 & -1 \end{bmatrix} = \mathcal{L} \{(1,3)\}

    Assim, o vetor próprio do valor próprio λ=2\lambda = 2 é (1,3)(1,3).

  • Para λ=4\lambda = 4:

    Temos então que

    A4I=[1133]A - 4I = \begin{bmatrix} 1 & -1\\ 3 & -3 \end{bmatrix}

    Então, os vetores próprios são os vetores que pertencem ao espaço nulo da matriz:

    Nul[1133]=L{(1,1)}\operatorname{Nul} \begin{bmatrix} 1 & -1\\ 3 & -3 \end{bmatrix} = \mathcal{L} \{(1,1)\}

    Assim, o vetor próprio do valor próprio λ=4\lambda = 4 é (1,1)(1,1).

Assim temos que os vetores próprios de AA são (1,3)(1,3) (λ=2\lambda = 2) e (1,1)(1,1) (λ=4\lambda = 4).

De acordo com a proposição acima, temos a solução geral da nossa equação:

x(t)=c1e2t[13]+c2e4t[11]x(t) = c_1 e^{2t} \begin{bmatrix} 1\\ 3 \end{bmatrix} + c_2 e^{4t} \begin{bmatrix} 1\\ 1 \end{bmatrix}

Vamos agora descobrir os valores das constantes, de acordo com o valor inicial dado no enunciado.

x(0)=[21]c1[13]+c2[11]=[21][1131][c1c2]=[21]x(0) = \begin{bmatrix} 2\\ -1 \end{bmatrix} \Leftrightarrow c_1 \begin{bmatrix} 1\\ 3 \end{bmatrix} + c_2 \begin{bmatrix} 1\\ 1 \end{bmatrix} = \begin{bmatrix} 2\\ -1 \end{bmatrix} \Leftrightarrow \begin{bmatrix} 1 & 1\\ 3 & 1 \end{bmatrix} \begin{bmatrix} c_1\\ c_2 \end{bmatrix} = \begin{bmatrix} 2\\ -1 \end{bmatrix}
[112311][112027][10320172]\augmatrix{cc|c}{ 1 & 1 & 2\\ 3 & 1 & -1 } \to \augmatrix{cc|c}{ 1 & 1 & 2\\ 0 & -2 & -7 } \to \augmatrix{cc|c}{ 1 & 0 & -\frac{3}{2}\\ 0 & 1 & \frac{7}{2} }

Temos assim, a solução final

x(t)=32e2t[13]+72e4t[11]x(t) = -\frac{3}{2} e^{2t} \begin{bmatrix} 1\\ 3 \end{bmatrix} + \frac{7}{2} e^{4t} \begin{bmatrix} 1\\ 1 \end{bmatrix}

ou, como um sistema,

{x1(t)=32e2t+72e4tx2(t)=92e2t+72e4t\begin{cases} x_1 (t) = -\frac{3}{2} e^{2t} + \frac{7}{2} e^{4t}\\ x_2 (t) = -\frac{9}{2} e^{2t} + \frac{7}{2} e^{4t} \end{cases}

Com valores próprios complexos

Proposição

Seja AA uma matriz n×nn\times n com entradas reais.
Então y:RCny: \R \to \C^n é solução de  ⁣dy ⁣dt=Ay\frac{\d y}{\d t} = Ay se e só se

u(t)=Rey(t)ev(t)=Imy(t)\begin{darray}{ccc} u(t) = R_e y(t) & \text{e} & v(t) = I_m y(t) \end{darray}

são soluções reais da mesma equação

Esta proposição permite-nos resolver as equações quando os valores próprios da matriz são complexos.

Exemplo

Resolver a seguinte equação diferencial:

x=[1111]xx' = \begin{bmatrix} 1 & -1\\ 1 & 1 \end{bmatrix} x

Começamos então por determinar os valores e vetores próprios de AA:

det(AλI)=1λ111λ=(1λ)2+1\det (A - \lambda I) = \begin{vmatrix} 1 - \lambda & -1\\ 1 & 1 - \lambda \end{vmatrix} = (1-\lambda) ^2 + 1

Temos então dois valores próprios complexos: λ=1+i\lambda = 1 + i e λ=1i\lambda = 1 - i.

Para λ=1+i\lambda = 1 + i:

Nul[i11i]=Nul[1i00]=L{(1,i)}\operatorname{Nul} \begin{bmatrix} -i & -1\\ 1 & -i \end{bmatrix} = \operatorname{Nul} \begin{bmatrix} 1 & -i\\ 0 & 0 \end{bmatrix} = \mathcal{L} \{(1,-i)\}

Então

e(1+i)t[1i]e^{(1+i) t} \begin{bmatrix} 1\\ -i \end{bmatrix}

é uma solução do sistema de EDO.

Podemos continuar a desenvolver, de forma a separar a parte real da parte imaginária:

e(1+i)t[1i]=eteit[1i]=et(cost+isint)[1i]=et[costsint]+iet[sintcost]\begin{aligned} e^{(1+i) t} \begin{bmatrix} 1\\ -i \end{bmatrix} &= e^t e^{it} \begin{bmatrix} 1\\ -i \end{bmatrix}\\ &= e^t (\cos t + i \sin t) \begin{bmatrix} 1\\ -i \end{bmatrix}\\ &= e^t \begin{bmatrix} \cos t\\ \sin t \end{bmatrix} + i e^t \begin{bmatrix} \sin t\\ -\cos t \end{bmatrix} \end{aligned}

Logo a solução geral do sistema de EDO é

x(t)=c1et[costsint]+c2et[sintcost]x(t) = c_1 e^t \begin{bmatrix} \cos t\\ \sin t \end{bmatrix} + c_2 e^t \begin{bmatrix} \sin t\\ -\cos t \end{bmatrix}

É de notar que neste exemplo, só descobrimos os vetores próprios de um dos valores próprios complexos.
Como os números são o conjugado um do outro, os seus vetores próprios são linearmente dependentes, pelo que é inútil calcular os vetores próprios do outro.
Além disso, podemos notar que pelo teorema que diz que o as soluções constituem espaço linear um espaço linear de dimensão nn, isto é, que existem nn valores na base do espaço solução, como temos uma matrix 2×22\times 2, iremos ter apenas dois vetores nessa base, que encontramos logo com apenas um valor próprio. As restantes soluções da EDO são obtidas a partir de combinação linear das duas respostas obtidas.

Equações Não Homogéneas

Estas equações correspondem à forma

 ⁣dy ⁣dt=Ay+B(t)\frac{\d y}{\d t} = Ay + B(t)

em que AA é uma matriz constante n×nn \times n e BB é uma função IRRnI \in \R \to \R^n.

Podemos assim redefinir o teorema da variação das constantes que estava anteriormente definido.

Teorema

Variação das constantes

Seja uma equação do tipo

 ⁣dy ⁣dt=Ay+B(t)y(t0)=y0\begin{darray}{cc} \frac{\d y}{\d t} = A y + B(t) & y(t_0) = y_0 \end{darray}

com t0It_0 \in I e y0Rny_0 \in \R^n, a equação tem uma única solução:

y(t)=e(tt0)Ay0+t0te(ts)AB(s) ⁣dsy(t) = e^{(t-t_0) A} y_0 + \int_{t_0}^{t} e^{(t-s) A} B(s) \d s
Exemplo

Seja a equação

{ ⁣dx ⁣dt=2x+y ⁣dy ⁣dt=2et\begin{cases} \frac{\d x}{\d t} = 2x + y\\ \frac{\d y}{\d t} = 2 e^t \end{cases}

e x(0)=y(0)=0x(0) = y(0) = 0, determine a solução.

 ⁣d ⁣dt[xy]=[2100]A[xy]+[02et]B(t)\frac{\d }{\d t} \begin{bmatrix} x\\ y \end{bmatrix} = \underbrace{\begin{bmatrix} 2 & 1\\ 0 & 0 \end{bmatrix}}_{A} \begin{bmatrix} x\\ y \end{bmatrix} + \underbrace{\begin{bmatrix} 0\\ 2e^t \end{bmatrix}}_{B(t)}

Por um exemplo anterior, podemos calcular a exponencial da matriz:

etA=[e2t12(e2t1)01]e^{tA} = \begin{bmatrix} e^{2t} & \frac{1}{2} (e^{2t} - 1)\\ 0 & 1 \end{bmatrix}

Então, pela fórmula da variação das constantes, temos que

[xy]=etA[00]+0t[e2t12(e2t1)01][02es] ⁣ds=0t[e2tses2es] ⁣ds=[0te2tses ⁣ds0t2es ⁣ds]=[[e2tses]0t2(et1)]=[etet+e2t+12(et1)]=[(et1)22(et1)]\begin{aligned} \begin{bmatrix} x\\ y \end{bmatrix} &= e^{tA} \begin{bmatrix} 0\\ 0 \end{bmatrix} + \int_{0}^{t} \begin{bmatrix} e^{2t} & \frac{1}{2} (e^{2t} - 1)\\ 0 & 1 \end{bmatrix} \begin{bmatrix} 0\\ 2e^s \end{bmatrix} \d s\\ & = \int_{0}^{t} \begin{bmatrix} e^{2t-s}-e^s\\ 2e^s \end{bmatrix} \d s\\ &= \begin{bmatrix} \int_{0}^{t} e^{2ts} - e^s \d s\\ \int_{0}^{t} 2e^s \d s \end{bmatrix}\\ &= \begin{bmatrix} \left[ -e^{2t-s} - e^s\right]^t_0\\ 2(e^t-1) \end{bmatrix}\\ &= \begin{bmatrix} -e^t -e^t + e^{2t} + 1\\ 2(e^t-1) \end{bmatrix}\\ &= \begin{bmatrix} (e^t - 1)^2\\ 2(e^t - 1) \end{bmatrix} \end{aligned}