2.1 Equações lineares Ajude a manter o site livre, gratuito e sem propagandas. Colabore!
Nesta seção, discutimos sobre equações a diferenças de ordem 1 e lineares. Tais equações podem ser escritas na seguinte forma
y ( n + 1 ) = a ( n ) y ( n ) + g ( n ) ,
(2.2)
onde n = n 0 , n 0 + 1 , … , sendo n 0 um número inteiro, a : n ↦ a ( n ) e g : n ↦ g ( n ) é o termo fonte. A equação é dita ser homogênea quando g ≡ 0 e, caso contrário, é dita ser não homogênea .
2.1.1 Equação homogênea
A solução de uma equação a diferenças de ordem 1, linear e homogênea
y ( n + 1 ) = a ( n ) y ( n ) , n ≥ n 0 ,
(2.3)
pode ser obtida por iterações diretas. Para n ≥ n 0 , temos
y ( n + 1 ) = a ( n ) y ( n )
(2.4)
= a ( n ) a ( n − 1 ) y ( n − 1 )
(2.5)
= a ( n ) a ( n − 1 ) a × ( n − 2 ) y ( n − 2 )
(2.6)
⋮
(2.7)
= a ( n ) a ( n − 1 ) ⋯ × a ( n 0 ) y ( n 0 ) .
(2.8)
Ou seja, dado o valor inicial y ( n 0 ) , temos a solução2 2 endnote: 2 A demonstração por ser feita por indução matemática.
y ( n ) = a ( n 0 ) a ( n 0 + 1 ) ⋯ × a ( n − 1 ) y ( n 0 ) .
(2.9)
A fim de termos uma notação mais prática, vamos usar a notação de produtório3 3 endnote: 3 Veja mais em Wiki: Produtório .
∏ i = n 0 n − 1 a ( n ) = a ( n 0 ) a ( n 0 + 1 ) × ⋯ a ( n − 1 ) .
(2.10)
Com esta notação, a solução de (2.3 ) pode ser escrita como segue
y ( n ) = [ ∏ i = n 0 n − 1 a ( i ) ] y ( n 0 ) ,
(2.11)
assumindo a notação de que ∏ i = n + 1 n a ( i ) = 1 .
Exemplo 2.1.1 .
Vamos calcular a solução de
y ( n + 1 ) = 2 y ( n ) , n ≥ 0 .
(2.12)
a)
Comparando com (2.3 ), temos a ( n ) = 2 para todo n . Calculando a solução por iterações diretas, temos
y ( n + 1 ) = 2 y ( n )
(2.13)
= 2 ⋅ 2 y ( n − 1 )
(2.14)
= 2 2 y ( n − 1 )
(2.15)
= 2 2 ⋅ 2 y ( n − 2 )
(2.16)
= 2 3 y ( n − 2 )
(2.17)
⋮
= 2 n + 1 y ( 0 )
(2.18)
ou, equivalentemente, temos a solução
b)
Por (2.11 ).
y ( n ) = [ ∏ i = 0 n − 1 2 ] y ( 0 )
(2.20)
= ( 2 ⋅ 2 ⋅ ⋯ ⋅ 2 ⏟ n vezes ) y ( 0 )
(2.21)
= 2 n y ( 0 ) .
(2.22)
A solução vale para qualquer valor inicial y ( 0 ) .
No Python , podemos computar a solução da equação a diferenças (2.12 ) com os seguintes comandos:
1 In : from sympy import *
2 In : n = symbols ( 'n' , integer = true )
3 In : y = symbols ( 'y' , cls = Function )
4 In : ead = Eq ( y ( n +1),2* y ( n ))
5 In : rsolve ( ead , y ( n ))
6 Out : 2** n * C0
2.1.2 Equação não homogênea
A solução de uma equação a diferenças de ordem 1, linear e não homogênea
y ( n + 1 ) = a ( n ) y ( n ) + g ( n ) , n ≥ n 0 ,
(2.23)
pode ser obtida por iterações diretas.
Vejamos, para n ≥ n 0 temos
y ( n + 1 ) = a ( n ) y ( n ) + g ( n )
(2.24)
= a ( n ) [ a ( n − 1 ) y ( n − 1 ) + g ( n − 1 ) ] + g ( n )
(2.25)
= a ( n ) a ( n − 1 ) y ( n − 1 ) + a ( n ) g ( n − 1 ) + g ( n )
(2.26)
= a ( n ) a ( n − 1 ) [ a ( n − 2 ) y ( n − 2 ) + g ( n − 2 ) ]
+ a ( n ) g ( n − 1 ) + g ( n )
(2.27)
= a ( n ) a ( n − 1 ) a × ( n − 2 ) y ( n − 2 )
+ a ( n ) a ( n − 1 ) g ( n − 2 ) + a ( n ) g ( n − 1 ) + g ( n )
(2.28)
⋮
= a ( n ) a ( n − 1 ) ⋯ × a ( n 0 ) y ( n 0 )
+ a ( n 0 + 1 ) a ( n 0 + 2 ) × ⋯ a ( n ) g ( n 0 )
+ a ( n 0 + 2 ) a ( n 0 + 3 ) × ⋯ a ( n ) g ( n 0 + 1 )
+ ⋯ + a ( n ) g ( n − 1 ) + g ( n )
(2.29)
Logo, podemos inferir que a solução é dada por4 4 endnote: 4 A demonstração por ser feita por indução matemática.
y ( n ) = a ( n 0 ) a ( n 0 + 1 ) ⋯ × a ( n − 1 ) y ( n 0 )
+ a ( n 0 + 1 ) a ( n 0 + 2 ) × ⋯ a ( n − 1 ) g ( n 0 )
+ a ( n 0 + 2 ) a ( n 0 + 3 ) ⋯ × a ( n − 1 ) g ( n 0 + 1 )
+ ⋯ + a ( n − 1 ) g ( n − 2 ) + g ( n − 1 )
(2.30)
Aqui, por maior praticidade, vamos empregar a notação de somatório5 5 endnote: 5 Veja mais em Wiki: Somatório
∑ i = n 0 n a ( i ) = a ( n 0 ) + a ( n 0 + 1 ) + ⋯ + a ( n ) .
(2.31)
Com isso, a solução de (2.23 ) pode ser escrita como segue
y ( n ) = [ ∏ i = n 0 n − 1 a ( i ) ] y ( n 0 )
+ ∑ i = n 0 n − 1 [ ∏ j = i + 1 n − 1 a ( j ) ] g ( i ) .
(2.32)
No último termo, consideramos a notação ∑ j = i + 1 i a ( i ) = 0 .
Exemplo 2.1.2 .
Vamos calcular a solução de
y ( n + 1 ) = 2 y ( n ) − 1 , n ≥ 0 .
(2.33)
Comparando com (2.23 ), temos a ( n ) = 2 e g ( n ) = − 1 para todo n .
1.
Cálculo por iterações diretas.
Calculando a solução por iterações diretas, temos
y ( n + 1 ) = 2 y ( n ) − 1
(2.34)
= 2 ⋅ [ 2 y ( n − 1 ) − 1 ] − 1
(2.35)
= 2 2 y ( n − 1 ) − 2 − 1
(2.36)
= 2 2 ⋅ [ 2 y ( n − 2 ) − 1 ] − 2 − 1
(2.37)
= 2 3 y ( n − 2 ) − 2 2 − 2 − 1
(2.38)
⋯
= 2 n + 1 y ( 0 ) − ∑ i = 0 n 2 i
(2.39)
Logo, temos
y ( n ) = 2 n y ( 0 ) − ∑ i = 0 n − 1 2 i .
(2.40)
Este último termo, é a soma dos termos da progressão geométrica6 6 endnote: 6 Veja mais em Wiki:Progressão geométrica. de razão q = 2 e termo inicial 1 (veja Subseção 2.1.3 ), i.e.
∑ i = 0 n − 1 q i = 1 − q n 1 − q .
(2.41)
Portanto, a solução (2.23 ) é
y ( n ) = 2 n y ( 0 ) − 1 − 2 n 1 − 2
(2.42)
= 2 n y ( 0 ) − 2 n + 1 .
(2.43)
2.
y ( n ) = [ ∏ i = n 0 n − 1 a ( i ) ] y ( n 0 )
(2.44)
+ ∑ i = n 0 n − 1 [ ∏ j = i + 1 n − 1 a ( j ) ] g ( i )
(2.45)
= [ ∏ i = 0 n − 1 2 ] y ( 0 )
(2.46)
+ ∑ i = 0 n − 1 [ ∏ j = i + 1 n − 1 2 ] ( − 1 )
(2.47)
= 2 n y ( 0 ) − ∑ i = 0 n − 1 2 n − 1 − i
(2.48)
= 2 n y ( 0 ) − 2 n − 1 ∑ i = 0 n − 1 2 − i
(2.49)
Este último somatório é a soma dos termos da progressão geométrica de razão q = 1 / 2 e termo inicial 1 ((veja Subseção 2.1.3 ), equação (2.60 )). Logo,
∑ i = 0 n − 1 2 − i = 1 − ( 1 2 ) n 1 − 1 2
(2.50)
= 2 ( 1 − 2 − n ) .
(2.51)
Retornando a (2.49 ), temos
y ( n ) = 2 n y ( 0 ) − 2 n − 1 ⋅ 2 ⋅ ( 1 − 2 − n )
(2.52)
= 2 n y ( 0 ) − 2 n + 1 .
(2.53)
A solução vale para qualquer valor inicial y ( 0 ) .
No Python , podemos computar a solução da equação a diferenças (2.12 ) com os seguintes comandos:
1 In : from sympy import *
2 In : n = symbols ( 'n' , integer = true )
3 In : y = symbols ( 'y' , cls = Function )
4 In : ead = Eq ( y ( n +1),2* y ( n )-1)
5 In : rsolve ( ead , y ( n ))
6 Out : 2** n * C0 + 1
Observamos que esta solução é equivalente à (2.53 ), pois
y ( n )
= 2 n y ( 0 ) − 2 n + 1
(2.54)
= 2 n [ y ( 0 ) − 1 ] + 1 ,
(2.55)
onde y ( 0 ) é um valor inicial arbitrário.
2.1.3 Somas definidas
Seguem algumas somas definidas que podem ser úteis na resolução de equações a diferenças.
∑ k = 1 n k
= n ( n + 1 ) 2
(2.56)
∑ k = 1 n k 2
= n ( n + 1 ) ( 2 n + 1 ) 6
(2.57)
∑ k = 1 n k 3
= [ n ( n + 1 ) 2 ] 2
(2.58)
∑ k = 1 n k 4
= n ( 6 n 4 + 15 n 3 + 10 n 2 − 1 ) 30
(2.59)
∑ k = 0 n − 1 q k
= ( 1 − q n ) 1 − q , q ≠ 1
(2.60)
∑ k = 1 n k q k
= ( q − 1 ) ( n + 1 ) q n + 1 − q n + 2 + q ( q − 1 ) 2
(2.61)
Exercícios resolvidos
ER 2.1.1 .
Calcule a solução da equação à diferenças
y ( n + 1 )
= 1 2 y ( n ) , n ≥ 0 ,
(2.62)
y ( 0 )
= 1 .
(2.63)
Solução .
De (2.11 ), temos
y ( n )
= [ ∏ i = 0 n − 1 1 2 ] y ( 0 )
(2.64)
= ( 1 2 ) n ⋅ 1
(2.65)
= 2 − n .
(2.66)
No Python , podemos computar a solução deste exercício com os seguintes comandos:
1 In : from sympy import *
2 In : n = symbols ( 'n' , integer = true )
3 In : y = symbols ( 'y' , cls = Function )
4 In : ead = Eq ( y ( n +1), S (1)/2* y ( n ))
5 In : rsolve ( ead , y ( n ), { y (0):1})
6 Out : 0.5** n
ER 2.1.2 .
Calcule a solução de
y ( n + 1 )
= 2 y ( n ) + ( 1 2 ) n , n ≥ 0 ,
(2.67)
y ( 0 )
= 0 .
(2.68)
Solução .
De (2.32 ), temos
y ( n )
= [ ∏ i = 0 n − 1 2 ] y ( 0 )
+ ∑ i = 0 n − 1 [ ∏ j = i + 1 n − 1 2 ] ⋅ ( 1 2 ) i
(2.69)
= ∑ i = 0 n − 1 2 n − 1 − i ⋅ 2 − i
(2.70)
= ∑ i = 0 n − 1 2 n − 1 ⋅ 2 − 2 i
(2.71)
= 2 n − 1 ∑ i = 0 n − 1 ( 1 4 ) i
(2.72)
= 2 n − 1 ⋅ [ 1 − ( 1 4 ) n ] 1 − 1 4
(2.73)
= 2 n − 1 ⋅ 4 3 ⋅ ( 1 − 1 4 n )
(2.74)
= 4 3 ( 2 n − 1 − 2 n − 1 4 n )
(2.75)
= 4 3 ( 2 n − 1 − 2 n − 1 2 − 2 n )
(2.76)
= 4 3 ( 2 n − 1 − 2 − n − 1 )
(2.77)
= 2 3 ( 2 n − 2 − n ) .
(2.78)
No Python , podemos computar a solução deste exercício com os seguintes comandos:
1 In : from sympy import *
2 In : n = symbols ( 'n' , integer = true )
3 In : y = symbols ( 'y' , cls = Function )
4 In : ead = Eq ( y ( n +1),2* y ( n )+(1/2)** n )
5 In : rsolve ( ead , y ( n ), { y (0):0})
6 Out : 2*2** n /3 - 2*2**(- n )/3
Exercícios
E. 2.1.1 .
Calcule a solução de
y ( n + 1 ) = 3 y ( n ) , n ≥ 0 .
(2.79)
Resposta
E. 2.1.2 .
Calcule a solução de
y ( n + 1 )
= 1 3 y ( n ) , n ≥ 0 ,
(2.80)
y ( 0 )
= − 1 .
(2.81)
Resposta
E. 2.1.3 .
Considere um empréstimo de $ 100 a uma taxa mensal de 1 % . Considerando y ( 0 ) = 100 , qual o valor de y ( n ) no n -ésimo mês? Modele o problema como uma equação à diferenças e calcule sua solução. Então, calcule o valor da dívida no 36º mês.
Resposta
Resposta .
y ( n + 1 ) = 1,01 ⋅ y ( n ) , y ( 0 ) = 100 ; y ( n ) = 100 ⋅ 1,01 n ; y ( 36 ) ≈ 143,08
E. 2.1.4 .
Calcule a solução de
y ( n + 1 )
= 3 y ( n ) − 3 , n ≥ 0 ,
(2.82)
y ( 0 )
= 2 .
(2.83)
Resposta
E. 2.1.5 .
Calcule a solução de
y ( n + 1 )
= n y ( n ) + n ! , n ≥ 0 ,
(2.84)
y ( 0 )
= 1 .
(2.85)
Resposta
E. 2.1.6 .
Calcule a solução de
y ( n + 1 )
= 2 y ( n ) + 2 n , n ≥ 0 ,
(2.86)
y ( 0 )
= 2 .
(2.87)
Resposta
E. 2.1.7 .
Considere um empréstimo de $ 100 a uma taxa mensal de 1 % e com parcelas mensais fixas de $ 1 . Considerando y ( 0 ) = 100 , qual o valor de y ( n ) no n -ésimo mês? Modele o problema como uma equação à diferenças e calcule sua solução.
Resposta
Resposta .
y ( n + 1 ) = 1,01 ⋅ y ( n ) − 1 , y ( 0 ) = 100 ; y ( n ) = 100 ;
E. 2.1.8 .
Calcule a solução de
y ( n + 1 ) = a y ( n ) + b , n ≥ 0 ,
(2.88)
onde a e b são constantes com a ≠ 1 .
Resposta
Resposta .
y ( n ) = ( y ( 0 ) − b 1 − a ) a n + b 1 − a
Envie seu comentárioAs informações preenchidas são enviadas por e-mail para o desenvolvedor do site e tratadas de forma privada. Consulte a Política de Uso de Dados para mais informações. Aproveito para agradecer a todas/os que de forma assídua ou esporádica contribuem enviando correções, sugestões e críticas!