Equações Diferenciais Ordinárias

Universidade Eduardo Mondlane Sansão Pedro, PhD i

Sumário

I Equações do Tipo Difusão - Equações Parabólicas 1

1.1 Equações Diferenciais Parabólicas . . . . . . . . . . . . . . . . . . . . . . . 1
1.1.1 Forma canônica das equações lineares de segunda ordem com coe-
ficientes constantes . . . . . . . . . . . . . . . . . . . . . . . . . . . 1
1.1.2 Actividades Formativas . . . . . . . . . . . . . . . . . . . . . . . . . 9
UNIDADE TEMÁTICA I

Equações do Tipo Difusão -

Equações Parabólicas

1.1 Equações Diferenciais Parabólicas

Aula 2: Equações Diferenciais Parabólicas

Objectivo geral da aula
Esta aula tem como objectivo principal demonstrar a utilidade de equações diferenciais
com derivadas parciais bem como são classificadas.

Objectivos específicos da aula

• Determinar a foma canónica das EDDPs;

1.1.1 Forma canônica das equações lineares de segunda ordem

com coeficientes constantes
Na aula passada vimos que a equação linear de segunda ordem da forma

Auxx + BUxy + Cuyy + Dux + Euy + F u + G = 0 (1.1.1)

cujos coeficientes A, B, C, D, E, F, G são funções reais definidas em uma região Ω ⊂ 2 R

e A2 + B 2 + C 2 > 0 pode ser classificada em três tipos mediante o sinal do discriminante
I = B 2 − 4AC . Nesta aula vamos procurar na base do tipo determinar uma forma mais
reduizida da equação .
Da mesma forma que o caso das cônicas, as EDDP lineares de segunda ordem tem formas
canônicas, ou seja , existe um sistema de coordenadas no qual a EDDP adopta uma forma
reduzida e mais simples, chamada forma canônica. Desta maneira nos restringiremos no

1
 Equações Diferenciais Ordinárias

estudo destas equações. Lembramos que em novas coordenadas a equação (1.1.1) tomou
a seguinte forma:

A′ uξξ + B ′ uξη + C ′ uηη + F̄ = 0 (1.1.2)

donde
A′ = Aa211 + Ba11 a12 + Ca212 ,
B ′ = 2Aa11 a21 + B(a11 a22 + a21 a12 ) + 2Ca12 a22 , (1.1.3)
C ′ = Aa221 + Ba21 a22 + Ca222 ,
e F̄ somente depende do máximo das derivadas de primeira ordem de u.

Forma canonica das equações hiperbólicas

Se considera primeiro as equações do tipo hiperbólico, ou seja, equações onde B 2 −4AC >
0. Com a finalidade de simplificá-las, se podem escolher a11 , a12 , a22 de tal forma que A′ =
C ′ e B ′ 6= 0. Suponha-se que A 6= 0. Se em (1.1.3) se força A′ = Aa211 +Ba11 a12 +Ca212 = 0
e sucede que a12 = 0 se tem a11 = 0, mas isso contradiz que a11 a22 − a21 a12 6= 0, portanto
a12 6= 0 similarmente se faz para C ′ = 0 em (1.1.3) e a22 = 0 se tem a21 = 0 o que
também contradiz que o determinante seja diferente de zero. Por tanto a21 6= 0 e a22 6= 0
. Assim A′ = 0 e C ′ = 0 em (1.1.3) se podem escrever como

Aϕ2 + Bϕ + c = 0

a11 a21
onde ϕ = a22
ou ϕ = a22
, mas em qualquer dos casos temos que
√
−B ± B 2 − 4AC
ϕ=
2A

B 2 − 4AC > 0 isto nos leva a definir novas coordenadas (ξ, η)

√
ξ = −(B + B 2 − 4AC)x + 2Ay
√ (1.1.4)
η = −(B − B 2 − 4AC)x + 2Ay

observe que a condição a11 a22 − a21 a12 6= 0 se cumpre dado que
( ) ( √ )
a11 a12 −(B + B 2 − 4AC) 2A √
= √ = −4A B 2 − 4AC 6= 0.
a21 a22 −(B − B 2 − 4AC) 2A
√
mais ainda, substituindo em B ′ da equação (1.1.3) com a11 = −B − B 2 − 4AC , a21 =
√
−B + B 2 − 4AC , a12 = 2A,a22 = 2A se tem

B ′ = −4A(b2 − 4AC) 6= 0.

Universidade Eduardo Mondlane Sansão Pedro, PhD 2

 Equações Diferenciais Ordinárias

Em novas coordenadas (ξ, η) a equação (1.1.1) adopta a forma

B ′ uξη + F = 0.

Dado que B ′ 6= 0, dividindo a esta nova equação por B ′ e definindo G = F

B′
, se obtêm
a forma canonica da equação hiperbólica

uξη = G. (1.1.5)

Donde, como foi mencionado, F depende das derivadas parciais de primeira ordem de u
e de mesma forma G. Introdizindo a nova mudança de variaveis ξ ′ = ξ+η
2
, η ′ = ξ−η
2
se
obtem uma outra forma canonica das equações hiperbólicas

uξ′ ξ′ − uη′ η′ = F ′ . (1.1.6)

Se deve fazer notar que a suposição A 6= 0 não limita a generalidade da discussão e

que se A = 0, por C 6= 0 uma discussão similar a de caso A 6= 0 conduz a mesma
forma canonica (1.1.5). Finalmente se A = C = 0 então a equação (1.1.1) toma a forma
Buxy + Dux + Euy + F u + G = 0 a qual trivialmente pode ser levada a forma canonica
já pela hipótese B 6= 0.

Exemplo 1.1.1. Mostre que a equação diferencial parcial uxx − 3uxy − 10uyy = 0 é
hiperbólica e encotre a forma canónica (1.1.5) da equação.
Solução : Em este caso A = 1, B = −3, C = −10, assim que B 2 − 4AC = 49 > 0,
portanto a equação é hiperbóloca. Por medio das equações (1.1.4) encontramos as novas
coordenadas

ξ = −4x + 2y
η = 10x + 2y

Então , as derivadas de segunda ordem da equação dada são

uxx = 16uξξ − 80uξη + 100uηη

uxy = −8uξξ + 2uξη + 20uηη
uyy = 4uξξ + 8uξη + 4uηη .

Substituindo estas equações se obtém:

−166uξη = 0 ou uξη = 0

a qual é a forma canónica da equação .

Universidade Eduardo Mondlane Sansão Pedro, PhD 3

 Equações Diferenciais Ordinárias

Forma canonica das equações parabólicas

Para as equações tipo parabólicas se cumpre I = 0. Suponhamos que A 6= 0, a equação

(1.1.1) pode simplificar-se usando (1.1.4) no qual se considera a troca de variaveis:

ξ = −Bx + 2Ay
(1.1.7)
η=x

de onde η é escolhida de maneira que

( ) ( )
a11 a12 −B 2A
= = −2A 6= 0.
a21 a22 1 0

Tomando os novos valores de a11 = −B, a12 = 2A, a21 = 1, a22 = 0 se tem

A′ = 0, B ′ = 0, C ′ = A (1.1.8)

tal que ao substituir em A′ uξξ + B ′ uξη + C ′ uηη + F̄ = 0 se obém

C ′ uηη + F̄ = 0

e dividindo por C ′ encontramos a equação na forma canónica para as equações parabólicas

′
uηη = F ′ (1.1.9)

donde F ′ contém no máximo, derivadas parciais de primeira ordem de u.

Observação 1.1.1. O caso de A = 0 implica B = 0 ja que B 2 − 4AC = 0 e, portanto,

C 6= 0 dada a condição A2 + B 2 + C 2 > 0. Desta forma, a equação (1.1.1) tem a forma

Cuyy = F̄

donde F̄ contém no máximo, derivadas parciais de primeira ordem de u. Segue-se então

que (1.1.1) pode ser facilmente levado para a forma canônica da equação parabólica (1.1.9),
se ela for dividida por C , assim a discussão do caso parabólico está concluída.

Exemplo 1.1.2. Reduza a equação 2uxx + 4uxy + 2uyy = 0 na forma canónica. Neste
caso A = C = 2, B = 4 pelo que B 2 − 4AC = 0, assim, a equação é parabólica. Ao usar
as equações (1.1.7) encontramos as novas coordenadas:

ξ = −4x + 4y
η = x.

Universidade Eduardo Mondlane Sansão Pedro, PhD 4

 Equações Diferenciais Ordinárias

Depois das substituições tem-se A′ = 0, B ′ = 0, C ′ = 2, que resulta em

2uηη = 0 ou uηη = 0

a qual é a forma canónica da equação .

Forma canonica das equações eípticas

Para as equações elípticas se cumpre que B − 4AC < 0 pelo que, A 6= 0 e C 6= 0. Notamos
que as equações (1.1.4) nos levam a equações de variavel complexa:
√
ξ = −(B + i 4AC − B 2 )x + 2Ay
√ (1.1.10)
η = −(B − i 4AC − B 2 )x + 2Ay
√
onde i = −1. Uma substituição direta de tais equações em (1.1.1) leva a uma equa-
ção em derivadas parciais com variáveis complexas. No entanto, introduzindo as novas
variáveis
ξ+η η−ξ
ξ′ = , η′ = (1.1.11)
2 2i
se obtém as variaveis reais:
ξ ′ = −Bx + 2Ay
√ (1.1.12)
η = 4AC − B 2 x,
as quais constituem uma troca viável, verificando-se que
( ) ( )
a11 a12 −B 2A √
= √ = −2A 4AC − B 2 6= 0.
a21 a22 4AC − B 2 0

Fazendo substituição com os coeficientes (1.1.12) obtém-se A′ = C ′ = A(4AC −B 2 ), B ′ =

0. De onde ao dividir a equação A′ uξξ + B ′ uξη + C ′ uηη + F̄ = 0 por A(4AC − B 2 ) se
obtem a forma canónica das equações elípticas

uξ′ ξ′ + uη′ η′ = F ′ (1.1.13)

donde F ′ contém no máximo, derivadas parciais de primeira ordem.

Exemplo 1.1.3. Classifique a equação 17uxx +uxy +uyy = 0 e reduze-a a forma canónica.
Solução : Dado que A = 17, B = C = 1 teremos B 2 − 4AC = −67, portanto, a equação
é elíptica. Por meio da transformação (1.1.12) se tem a forma canónica uξξ + uηη = 0.

Abaixo apresentamos o quadro resumo das formas canónicas

Universidade Eduardo Mondlane Sansão Pedro, PhD 5

 Equações Diferenciais Ordinárias

I = B 2 − 4AC Tipo Forma canónica

I >0 hiperbólica uξη = F ′ (uξ , uη , ξ, η) ou uξξ − uηη = F ′
I =0 parabólica uηη = G′ (uξ , uη , ξ, η)
I <0 elíptica uξη = uξξ + uηη = H ′ (uξ , uη , ξ, η)

Observação 1.1.2. Notemos que as formas canónicas das equações lineares podem ser
estudadas com maior profundidade mediante a análise das curvas características, mediante
a invariabilidade do sinal do discriminante da equação característica.

Classificação de equações lineares de segunda ordem com mais de duas variá-

veis

Uma EDDP linear de segunda ordem em n variáveis é uma equação que pode ser expressar
na forma
∑n ∑
n
L[u] = aij uxi xj + bi uxi + cu = d (1.1.14)
i,j=1 i=1

onde aij, bi, c, d são funções reais de variáveis x1 , x2 , · · · , xn . Nesta sessão se estuda a
classificação de equação (1.1.14) em caso de coeficientes constanttes. A propriedade que
será utilizada para classificar a equação (1.1.14) é aquela que se reduz à invariância do sinal
do discriminante I para a equação (1.1.1). Para motivar esta propriedade, é necessário o
conceito de uma forma quadrática real.

Definição 1.1.1. Uma forma quadrática real em n variav́eis, é uma função

∑
n
H: R n
−→ R, (x1 , x2 , · · · , xn ) 7→ aij xi xj
i,j=1

com aij ∈ R.
A toda matriz simétrica, ou seja, a toda matriz A que satisfaça A = At , onde At
é a transposta de A, corresponde uma forma quadrática e, inversamente, uma matriz
simétrica pode ser associada a toda forma quadrática. Daqui em diante vamos supor que
este será o caso. Por exemplo, a forma quadrática H(x, y) = Ax2 +Bxy+Cy 2 corresponde
à matriz simétrica
( )
A 21 B
M= 1
2
B C

observe que
( )( )
1
A 2
B x
H(x, y) = (x, y) 1
2
B C y

Universidade Eduardo Mondlane Sansão Pedro, PhD 6

 Equações Diferenciais Ordinárias

Agora com a parte principal Auxx + Buxy + Cuyy da equação (1.1.1) a forma quadrática
Ax2 + Bxy + Cy 2 será associada e a matriz associada desta forma quadrática é aquela
correspondente à parte principal da equação (1.1.1). Assim M é a matriz simétrica da
parte principal Auxx +Buxy +Cuyy da equação (1.1.1). Sejam λ1 , λ2 os autovalores de M ,
ou seja, λ1 , λ2 são raízes da equação det(λIM ) = λ2 (A + C)λ(B24AC)/4 = 0 . Se deixa
como exercício verificar que λ1 e λ2 são números reais e que λ1 λ2 = (B24AC)/4 = I/4.
De onde,

• I > 0 ⇔ λ1 & λ2 são diferentes de zero e têm sinais opostos.

• I = 0 ⇔ pelo menos um dos λi , i = 1, 2, é zero.

• I < 0 ⇔ λ1 & λ2 são diferentes de zero e têm o mesmo sinal.

Desta forma, o sinal de I pode ser caracterizado pelo sinal dos autovalores de M . Para
matrizes simétricas em R
sabe-se que elas são semelhantes a uma matriz diagonal, temos
também o seguinte teorema

Teorema 1.1.1 (Silvester). Seja M uma matriz real simétrica, então o número de en-
tradas positivas e negativas de qualquer matriz diagonal semelhante a M é independente
da escolha da matriz diagonal.

Assim este teorema fornece uma propriedade intrínseca da equação (1.1.14) que per-
mite uma classificação da referida equação.

Definição 1.1.2. Sejam λ1 , λ2 , · · · , λn os autovalores da matriz da parte principal da

equação (1.1.14)

(i) Se λ1 , λ2 , · · · , λn forem diferentes de zero e pelo menos um deles tiver sinal diferente
dos demais, a equação é chamada de hiperbólica.

(ii) Se pelo menos um dos λ1 , λ2 , · · · , λn for zero, a equação é chamada de parabólica

(iii) Se λ1 , λ2 , · · · , λn são diferentes de zero e têm o mesmo sinal, a equação é chamada

elíptico.

Segue uma classificação em termos de formas quadráticas

Definição 1.1.3. Seja H : N → R R

uma forma quadrática com matriz associada A.
Diz-se que H está definido se os autovalores correspondentes ao polinômio característico:

det(λI − A) = 0 (1.1.15)

são todos positivos ou todos negativos. H é dito indefinido se a equação (1.1.15) tem
raízes que mudam de sinal. H é dito degenerado se a equação (1.1.15) tem uma solução
de λ = 0.

Universidade Eduardo Mondlane Sansão Pedro, PhD 7

 Equações Diferenciais Ordinárias

Assim, a equação é do tipo hiperbólica, se a forma quadrática (1.1.14) for indefinida,

parabólica se (1.1.14), for degenerada, elíptico se (1.1.14) for definido.

Exemplo 1.1.4. A forma geral da equação em derivadas parciais linear de segunda ordem
em três variáveis é

Auxx + BUxy + Cuxz + Duyz + Euyy + F uzz + G(x, y, z, ux , uy , uz ) = 0

a qual tem a forma quadrática

  
1 1
A 2
B 2
C x
1 1  
H(x, y, z) = (x, y, z)  2 B E 2  
D y .
1 1
2
C 2
D F z

Em particular, a equação de Laplace ∆u = uxx + Uyy + uzz = 0 tem associada a forma

quadrática
  
1 0 0 x
  
H(x, y, z) = (x, y, z) 0 1 0 y  .
0 0 1 z

Assim, o polinomio característico correspondente é dado por (λ − 1)3 = 0 pelo qual

constata-se que o operador de Laplace é elíptico.
Para a equação Lu = uxx + Uyy − uzz = 0 se tem a forma quadrática
  
1 0 0 x
  
H(x, y, z) = (x, y, z) 0 1 0  y  .
0 0 −1 z

Cujo, polinómio característico correspondente é dado por (λ − 1)2 (λ + 1) = 0, portanto,

L é hiperbólico.

Leituras obrigatórias
A leitura dos textos obrigatórios constitui a base de consolidação da aprendizagem dos
conteúdos desta unidade temática, daí que sejam indispensáveis, sendo igualmente de
suma importância a resolução dos exercícios que complementarão as leituras.

[1 ] Stanley J. Farlow Partial Differential Equations for Scientists and Engineers, DO-
VER PUBLICATIONS, INC. New York, 1993.

[2 ] Gabriel López garza & Fco. Hugo Martínez Ortiz, Equaciones Diferenciales Par-
ciales, 2013 Gabriel López.

Universidade Eduardo Mondlane Sansão Pedro, PhD 8

 Equações Diferenciais Ordinárias

1.1.2 Actividades Formativas

As actividades que se seguem complementam o estudo desta unidade temática. Resolva-as
e verifique se estão certas.

Actividade Formativa
1. Reduza as seguintes equações a forma canónica.

a) 3uxx + 2uxy − 5uyy − uy = 0.

b) 8uxx + 2uxy − 3uyy = xy .
c) 2uxx − 3uxy + uyy − ux + uy − 3 = 0
d) 9uxx − 6uxy + uyy + 10ux − 15uy − 50u + x − 2y = 0
e) uxx + 4uxy + 13uyy + 3ux + 24uy − 9u + 9(x + y) = 0

2. Resolva o exercício acima entrando a forma quadrática e o polinómio característico

correspondente.

Universidade Eduardo Mondlane Sansão Pedro, PhD 9