Articulo de referencia

Ecuación de diferencias racionales

Una ecuación de diferencias racional es una ecuación de diferencias no lineal de la forma [ 1 ] [ 2 ] [ 3 ] [ 4 ] incógnita norte + 1 = α + ∑ i = 0 k β i incógnita norte − i A +...

Una ecuación de diferencias racional es una ecuación de diferencias no lineal de la forma [ 1 ] [ 2 ] [ 3 ] [ 4 ]

incógnitanorte+1=α+i=0kβiincógnitanorteiA+i=0kBiincógnitanortei ,{\displaystyle x_{n+1}={\frac {\alpha +\sum _{i=0}^{k}\beta _{i}x_{ni}}{A+\sum _{i=0}^{k}B_{i}x_{ni}}}~,}

donde las condiciones inicialesincógnita0,incógnita1,,incógnitak{\displaystyle x_{0},x_{-1},\dots ,x_{-k}}son tales que el denominador nunca se anula para ningún n .

Ecuación de diferencias racionales de primer orden

Una ecuación de diferencias racional de primer orden es una ecuación de diferencias no lineal de la forma

Cuandoa,b,do,d{\displaystyle a,b,c,d}y la condición inicialw0{\displaystyle w_{0}}son números reales , esta ecuación de diferencias se llama ecuación de diferencias de Riccati . [ 3 ]

Dicha ecuación se puede resolver escribiendowt{\displaystyle w_{t}}como una transformación no lineal de otra variableincógnitat{\displaystyle x_{t}}que a su vez evoluciona linealmente. Entonces se pueden utilizar métodos estándar para resolver la ecuación de diferencias lineal enincógnitat{\displaystyle x_{t}}.

Las ecuaciones de esta forma surgen del problema de la escalera de resistencias infinita. [ 5 ] [ 6 ]

Resolver una ecuación de primer orden

Primer enfoque

WLoG, la cantidad similar a un determinanteadbdo0{\displaystyle ad-bc\neq 0}y esto se puede ver fácilmente al notar que, con una división tanto en el numerador como en el denominador en la ecuación ( 2 ) , siempre se puede establecerdo=1{\displaystyle c=1}, en cuyo caso, siadbdo=0{\displaystyle ad-bc=0}, entonces el numerador y el denominador se cancelarán, sin dejar ninguna ecuación de diferencias, sino que se reducirán por completo hasta convertirse exactamente en la constante.a{\displaystyle a}. Por lo tanto, un enfoque [ 7 ] para desarrollar la variable transformadaincógnitat{\displaystyle x_{t}}, es escribir

yt+1=αβyt{\displaystyle y_{t+1}=\alpha -{\frac {\beta }{y_{t}}}}

dóndeα=(a+d)/do{\displaystyle \alpha =(a+d)/c} yβ=(adbdo)/do2{\displaystyle \beta =(ad-bc)/c^{2}}y dóndewt=ytd/do{\displaystyle w_{t}=y_{t}-d/c}.

Más textos yt=incógnitat+1/incógnitat{\displaystyle y_{t}=x_{t+1}/x_{t}}se puede demostrar que produce

incógnitat+2αincógnitat+1+βincógnitat=0.{\displaystyle x_{t+2}-\alpha x_{t+1}+\beta x_{t}=0.}

Segundo enfoque

El enfoque anterior ya es de aplicabilidad general. El siguiente enfoque [ 8 ] proporciona una ecuación de diferencias de primer orden paraincógnitat{\displaystyle x_{t}}en lugar de uno de segundo orden. Dejemosr2=(da)2+4bdo{\displaystyle r^{2}=(da)^{2}+4bc}. Para el caso en quer20{\displaystyle r^{2}\geqslant 0}, cada término será de valor real y este método puede ser conveniente de usar. De lo contrario, el método sigue funcionando, pero aparecerán números complejos y podría ser más conveniente intentar un ansatz trigonométrico en su lugar. Sustituyendoincógnitat=1/(η+wt){\displaystyle x_{t}=1/(\eta +w_{t})}, lo cual implicawt=(1ηincógnitat)/incógnitat{\displaystyle w_{t}=(1-\eta x_{t})/x_{t}}Sustituyendo en la ecuación ( 2 ) , encontramos que siempre es posible hacerincógnitat{\displaystyle x_{t}}evolucionan según la ecuación de diferencias lineal simple no homogénea de primer orden.

incógnitat+1=(dηdoηdo+a)incógnitat+doηdo+a.{\displaystyle x_{t+1}=\left({\frac {d-\eta c}{\eta c+a}}\right)\!x_{t}+{\frac {c}{\eta c+a}}.}

al elegirη{\displaystyle \eta }de tal manera quedoη2(da)ηb=0{\displaystyle c\eta ^{2}-(da)\eta -b=0}y está claro que esto siempre se puede hacer en cualquiera de las dos opciones.η±=da±r2do{\displaystyle \eta _{\pm }={\frac {da\pm r}{2c}}}, incluso cuandordo{\displaystyle r\in \mathbb {C} }

Tercer enfoque

La ecuación

wt+1=awt+bdowt+d{\displaystyle w_{t+1}={\frac {aw_{t}+b}{cw_{t}+d}}}

También se puede resolver tratándola como un caso especial de la ecuación matricial más general.

incógnitat+1=(mi+Bincógnitat)(do+Aincógnitat)1,{\displaystyle X_{t+1}=-(E+BX_{t})(C+AX_{t})^{-1},}

donde A, B, C, E y X son matrices n × n (en este caso n = 1); la solución de esto es [ 9 ].

incógnitat=nortetDt1{\displaystyle X_{t}=N_{t}D_{t}^{-1}}

dónde

(nortetDt)=(BmiAdo)t(incógnita0I).{\displaystyle {\begin{pmatrix}N_{t}\\D_{t}\end{pmatrix}}={\begin{pmatrix}-B&-E\\A&C\end{pmatrix}}^{t}{\begin{pmatrix}X_{0}\\I\end{pmatrix}}.}

Solicitud

En [ 10 ] se demostró que una ecuación de Riccati matricial dinámica de la forma

Ht1=K+AHtAAHtdo(doHtdo)1doHtA,{\displaystyle H_{t-1}=K+A'H_{t}AA'H_{t}C(C'H_{t}C)^{-1}C'H_{t}A,}

que pueden surgir en algunos problemas de control óptimo en tiempo discreto , se pueden resolver utilizando el segundo enfoque anterior si la matriz C tiene solo una fila más que una columna.

Referencias

  1. Skellam, JG (1951). “Dispersión aleatoria en poblaciones teóricas”, Biometrika 38 196−–218, ecs. (41,42)
  2. Camouzis, Elias; Ladas, G. (16 de noviembre de 2007). Dinámica de ecuaciones en diferencias racionales de tercer orden con problemas abiertos y conjeturas . CRC Press. ISBN 9781584887669 vía Google Libros.
  3. 1 2 Kulenovic, Mustafa RS; Ladas, G. (30 de julio de 2001). Dinámica de ecuaciones en diferencias racionales de segundo orden: con problemas abiertos y conjeturas . CRC Press. ISBN 9781420035384 vía Google Libros.
  4. Newth, Gerald, "Orden mundial a partir de comienzos caóticos", Mathematical Gazette 88, marzo de 2004, 39-45 ofrece un enfoque trigonométrico.
  5. "Resistencia equivalente en circuitos en escalera" . Stack Exchange . Consultado el 21 de febrero de 2022 .
  6. "Pensando recursivamente: ¡Cómo resolver el rompecabezas de la escalera de resistencias infinitas!" . Youtube . Consultado el 21 de febrero de 2022 .
  7. Brand, Louis, "Una secuencia definida por una ecuación de diferencias", American Mathematical Monthly 62 , septiembre de 1955, 489-492 . en línea
  8. Mitchell, Douglas W., "Una solución analítica de Riccati para el control discreto en tiempo con dos objetivos", Journal of Economic Dynamics and Control 24, 2000, 615 622.
  9. ^ Martin, CF y Ammar, G., "La geometría de la ecuación matricial de Riccati y el método de valores propios asociados", en Bittani, Laub y Willems (eds.), The Riccati Equation , Springer-Verlag, 1991.
  10. Balvers, Ronald J., y Mitchell, Douglas W., "Reducción de la dimensionalidad de problemas de control lineal cuadrático", Journal of Economic Dynamics and Control 31, 2007, 141 159.

Lecturas adicionales

  • Simons, Stuart, "Una ecuación de diferencias no lineal", Mathematical Gazette 93, noviembre de 2009, 500–504.