| | | | |

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=n0,n0+1,, sendo n0 um número inteiro, a:na(n) e g:ng(n) é o termo fonte. A equação é dita ser homogênea quando g0 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),nn0, (2.3)

pode ser obtida por iterações diretas. Para nn0, temos

y(n+1)=a(n)y(n) (2.4)
=a(n)a(n1)y(n1) (2.5)
=a(n)a(n1)a×(n2)y(n2) (2.6)
(2.7)
=a(n)a(n1)×a(n0)y(n0). (2.8)

Ou seja, dado o valor inicial y(n0), temos a solução22endnote: 2A demonstração por ser feita por indução matemática.

y(n)=a(n0)a(n0+1)×a(n1)y(n0). (2.9)

A fim de termos uma notação mais prática, vamos usar a notação de produtório33endnote: 3Veja mais em Wiki: Produtório.

i=n0n1a(n)=a(n0)a(n0+1)×a(n1). (2.10)

Com esta notação, a solução de (2.3) pode ser escrita como segue

y(n)=[i=n0n1a(i)]y(n0), (2.11)

assumindo a notação de que i=n+1na(i)=1.

Exemplo 2.1.1.

Vamos calcular a solução de

y(n+1)=2y(n),n0. (2.12)
  1. a)

    Por iterações diretas.

    Comparando com (2.3), temos a(n)=2 para todo n. Calculando a solução por iterações diretas, temos

    y(n+1)=2y(n) (2.13)
    =22y(n1) (2.14)
    =22y(n1) (2.15)
    =222y(n2) (2.16)
    =23y(n2) (2.17)
    =2n+1y(0) (2.18)

    ou, equivalentemente, temos a solução

    y(n)=2ny(0). (2.19)
  2. b)

    Por (2.11).

    y(n)=[i=0n12]y(0) (2.20)
    =(222n vezes)y(0) (2.21)
    =2ny(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),nn0, (2.23)

pode ser obtida por iterações diretas.

Vejamos, para nn0 temos

y(n+1)=a(n)y(n)+g(n) (2.24)
=a(n)[a(n1)y(n1)+g(n1)]+g(n) (2.25)
=a(n)a(n1)y(n1)+a(n)g(n1)+g(n) (2.26)
=a(n)a(n1)[a(n2)y(n2)+g(n2)]
+a(n)g(n1)+g(n) (2.27)
=a(n)a(n1)a×(n2)y(n2)
+a(n)a(n1)g(n2)+a(n)g(n1)+g(n) (2.28)
=a(n)a(n1)×a(n0)y(n0)
+a(n0+1)a(n0+2)×a(n)g(n0)
+a(n0+2)a(n0+3)×a(n)g(n0+1)
++a(n)g(n1)+g(n) (2.29)

Logo, podemos inferir que a solução é dada por44endnote: 4A demonstração por ser feita por indução matemática.

y(n)=a(n0)a(n0+1)×a(n1)y(n0)
+a(n0+1)a(n0+2)×a(n1)g(n0)
+a(n0+2)a(n0+3)×a(n1)g(n0+1)
++a(n1)g(n2)+g(n1) (2.30)

Aqui, por maior praticidade, vamos empregar a notação de somatório55endnote: 5Veja mais em Wiki: Somatório

i=n0na(i)=a(n0)+a(n0+1)++a(n). (2.31)

Com isso, a solução de (2.23) pode ser escrita como segue

y(n)=[i=n0n1a(i)]y(n0)
+i=n0n1[j=i+1n1a(j)]g(i). (2.32)

No último termo, consideramos a notação j=i+1ia(i)=0.

Exemplo 2.1.2.

Vamos calcular a solução de

y(n+1)=2y(n)1,n0. (2.33)

Comparando com (2.23), temos a(n)=2 e g(n)=1 para todo n.

  1. 1.

    Cálculo por iterações diretas.

    Calculando a solução por iterações diretas, temos

    y(n+1)=2y(n)1 (2.34)
    =2[2y(n1)1]1 (2.35)
    =22y(n1)21 (2.36)
    =22[2y(n2)1]21 (2.37)
    =23y(n2)2221 (2.38)
    =2n+1y(0)i=0n2i (2.39)

    Logo, temos

    y(n)=2ny(0)i=0n12i. (2.40)

    Este último termo, é a soma dos termos da progressão geométrica66endnote: 6Veja 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=0n1qi=1qn1q. (2.41)

    Portanto, a solução (2.23) é

    y(n)=2ny(0)12n12 (2.42)
    =2ny(0)2n+1. (2.43)
  2. 2.

    Cálculo por (2.32).

    y(n)=[i=n0n1a(i)]y(n0) (2.44)
    +i=n0n1[j=i+1n1a(j)]g(i) (2.45)
    =[i=0n12]y(0) (2.46)
    +i=0n1[j=i+1n12](1) (2.47)
    =2ny(0)i=0n12n1i (2.48)
    =2ny(0)2n1i=0n12i (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=0n12i=1(12)n112 (2.50)
    =2(12n). (2.51)

    Retornando a (2.49), temos

    y(n)=2ny(0)2n12(12n) (2.52)
    =2ny(0)2n+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) =2ny(0)2n+1 (2.54)
=2n[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=1nk =n(n+1)2 (2.56)
k=1nk2 =n(n+1)(2n+1)6 (2.57)
k=1nk3 =[n(n+1)2]2 (2.58)
k=1nk4 =n(6n4+15n3+10n21)30 (2.59)
k=0n1qk =(1qn)1q,q1 (2.60)
k=1nkqk =(q1)(n+1)qn+1qn+2+q(q1)2 (2.61)

Exercícios resolvidos

ER 2.1.1.

Calcule a solução da equação à diferenças

y(n+1) =12y(n),n0, (2.62)
y(0) =1. (2.63)
Solução.

De (2.11), temos

y(n) =[i=0n112]y(0) (2.64)
=(12)n1 (2.65)
=2n. (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) =2y(n)+(12)n,n0, (2.67)
y(0) =0. (2.68)
Solução.

De (2.32), temos

y(n) =[i=0n12]y(0)
+i=0n1[j=i+1n12](12)i (2.69)
=i=0n12n1i2i (2.70)
=i=0n12n122i (2.71)
=2n1i=0n1(14)i (2.72)
=2n1[1(14)n]114 (2.73)
=2n143(114n) (2.74)
=43(2n12n14n) (2.75)
=43(2n12n122n) (2.76)
=43(2n12n1) (2.77)
=23(2n2n). (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)=3y(n),n0. (2.79)
Resposta.

y(n)=3ny(0)

E. 2.1.2.

Calcule a solução de

y(n+1) =13y(n),n0, (2.80)
y(0) =1. (2.81)
Resposta.

y(n)=13n

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.

y(n+1)=1,01y(n),y(0)=100; y(n)=1001,01n; y(36)143,08

E. 2.1.4.

Calcule a solução de

y(n+1) =3y(n)3,n0, (2.82)
y(0) =2. (2.83)
Resposta.

y(n)=12(3n+3)

E. 2.1.5.

Calcule a solução de

y(n+1) =ny(n)+n!,n0, (2.84)
y(0) =1. (2.85)
Resposta.

y(n)=n!

E. 2.1.6.

Calcule a solução de

y(n+1) =2y(n)+2n,n0, (2.86)
y(0) =2. (2.87)
Resposta.

y(n)=2n(n2+2)

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.

y(n+1)=1,01y(n)1,y(0)=100; y(n)=100;

E. 2.1.8.

Calcule a solução de

y(n+1)=ay(n)+b,n0, (2.88)

onde a e b são constantes com a1.

Resposta.

y(n)=(y(0)b1a)an+b1a


Envie seu comentário

As 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!