Solução via integrais sucessivas

Anterior: Introdução à equações diferenciais com atraso | Índice: Ecologia | Próximo: Estabilidade

O problema de valor inicial do sistema:

{\dot {x}}\left(t\right)=f\left(t,x_{t},u\left(t\right)\right)

para uma dada entrada ${\textstyle u\left(t\right)}$ consiste em determinar a solução contínua ${\textstyle x\left(t\right)}$ para ${\textstyle t\geq t_{0}}$ de forma que ${\textstyle x\left(t_{0}\right)=x_{0}}$ e ${\textstyle x\left(t\right)=\varphi \left(t\right)}$ para ${\textstyle t_{0}-\tau _{max}\leq t<t_{0}}$ onde ${\textstyle \tau _{max}={\text{const}}\in \left[0,\infty \right)}$ , sendo que ${\textstyle \varphi }$ é chamado muitas vezes de função inicial, ainda é comum assumir que ${\textstyle \varphi \left(t_{0}\right)=x\left(t_{0}\right)}$ . Além disso, é importante comentar que uma função de dimensão finita tem a propriedade de que o estado em qualquer tempo pode ser especificado listando um conjunto finito de valores, então de modo análogo, uma função de dimensão infinita exige um conjunto infinito de valores. Quando se fornece as condições iniciais para um sistema finito, só precisamos fornecer um pequeno conjunto de valores (basicamente os valores iniciais), porem para resolver uma equação com atraso, é preciso fornecer uma quantidade infinita, dessa equações diferenciais com atraso possuem dimensionalidade infinita.

Uma forma melhor de entender esta questão é pensar na solução da DDE (delayed differential equation) como um mapeamento de funções no intervalo ${\textstyle \left[t-\tau ,t\right]}$ em funções no intervalo ${\textstyle \left[t,t+\tau \right]}$ . Em outras palavras a solução pode ser pensada como uma sequência de funções ${\textstyle f_{0}\left(t\right),f_{1}\left(t\right),f_{2}\left(t\right)\dots .}$ definidas sobre um conjunto de intervalos contíguos de tempo com tamanho ${\textstyle \tau }$ .

Retornando à discussão da solução, a solução desejada é encontrada em intervalos sucessivos resolvendo equações diferenciais ordinários sem atraso em cada intervalo. Por exemplo, tendo um sistema:

{\begin{aligned}{\dot {x}}\left(t\right)=&f\left(t,x\left(t\right),x\left(t-\tau \right)\right)&t\geq &t_{0}\\x\left(t_{0}\right)=&x_{0}&t=&t_{0}\\x\left(t\right)=&\varphi \left(t\right)&t_{0}-\tau \leq t<&t_{0}\end{aligned}}

Então para ${\textstyle t\in \left[t_{0},t_{0}+\tau \right]}$ o sistema pode ser representado como uma equação ordinária diferencial. Para isso precisamos lembrar que em um instante qualquer ${\textstyle t\leq t_{0}+\tau }$ , o valor da função ${\textstyle ,x\left(t-\tau \right)}$ é dado no instante ${\textstyle T=t-\tau }$ , mas ${\textstyle T\leq t_{0}}$ , e nesta condição temos ${\textstyle x\left(t\right)=\varphi \left(t\right)}$ . Isto é:

{\begin{aligned}{\dot {x}}\left(t\right)=&f\left(t,x\left(t\right),\varphi \left(t-\tau \right)\right)&t_{0}\leq t\leq &t_{0}+\tau \\x\left(t_{0}\right)=&x_{0}&t=&t_{0}\end{aligned}}

Então se obtermos a solução ${\textstyle x\left(t\right)=\varphi _{1}\left(t\right)}$ no segmento ${\textstyle \left[t_{0},t_{0}+\tau \right]}$ , pode-se obter uma solução análoga para o próximo intervalo.

Exemplo:

Considerando o sistema:

{\begin{aligned}{\dot {x}}\left(t\right)=&6x\left(t-1\right)&t\geq &0\\x\left(t\right)=&t&-1\leq t<&0\end{aligned}}

Podemos observar que ${\textstyle t_{0}=0,\tau =1,\varphi \left(t\right)=t}$ e ainda ${\textstyle x\left(t_{0}\right)=\varphi \left(t_{0}\right)}$ . Então para o intervalo ${\textstyle \left[0,1\right]}$ :

{\begin{aligned}{\dot {x}}\left(t\right)=&6\varphi \left(t-1\right)&0\leq t&\leq 1\\x\left(0\right)=&0&t&=0\end{aligned}}

Resolvendo então, uma vez que ${\textstyle \varphi \left(t-1\right)=t-1}$ , então:

{\begin{aligned}{\frac {dx}{dt}}&=6\left(t-1\right)\\\int _{x\left(0\right)}^{x\left(t\right)}dx'&=6\int _{0}^{t}\left(t'-1\right)dt'\\x\left(t\right)&=3t^{2}-6t\end{aligned}}

E para o intervalo ${\textstyle \left[1,2\right]}$ , sendo ${\textstyle \varphi _{1}\left(t\right)=3t^{2}-6t}$ , agora temos o mesmo ${\textstyle \tau =1}$ , mas ${\textstyle t_{0}=1}$ , e novamente ${\textstyle x\left(t_{0}\right)=\varphi _{1}\left(t_{0}\right)}$ . Então: