Encontrar $\lim\limits_{t\to\infty}x(t)$ si $x'= (x-y)(1-x^2-y^2)$, $y' = (x+y)(1-x^2-y^2)$

Question

Supongamos $x_0,y_0$ son reales tal que $x_0^2+y_0^2>0.$ Supongamos $x(t)$ $y(t)$ cumplir con lo siguiente:

• $\frac{dx}{dt} = (x-y)(1-x^2-y^2),$
• $\frac{dy}{dt} = (x+y)(1-x^2-y^2),$
• $x(0) = x_0,$
• $y(0) = y_0.$

Se me pide que encuentre a $\lim_{t\to\infty}x(t).$

Dividiendo las dos ecuaciones diferenciales, tenemos $$\frac{dy}{dx} = \frac{x+y}{x-y} = \frac{1 + \frac{y}{x}}{1-\frac{y}{x}}.$$ Deje $v(x)=\frac{y}{x},$, de modo que $y = xv(x),$ y $$\frac{dy}{dx} = v + x\frac{dv}{dx}.$$ This implies that $$x\frac{dv}{dx} = \frac{1+v}{1-v} - v = \frac{1+v^2}{1-v}.$$ Rearranging and integrating, we get $$\log(x) = \arctan(v)-\frac{1}{2}\log(1+v^2)+c,$$ so $$\log(x) = \arctan\left(\frac{y}{x}\right)-\frac{1}{2}\log\left(1+\left(\frac{y}{x}\right)^2\right)+c.$$ By writing $1+\left(\frac{y}{x}\right)^2$ as $\frac{x^2+y^2}{x^2},$ we can simplify the above equation to $$\frac{1}{2}\log(x^2+y^2) = \arctan\left(\frac{y}{x}\right) + c.$$

Sin embargo, no estoy seguro de cómo proceder a partir de esta ecuación implícita. Cualquier ayuda sería muy apreciada!

Answer 1

Primero pasamos a coordenadas polares, a saber:

$x = r\cos \theta; \; y = r \sin \theta; \tag 1$

$\dot x = \dot r \cos \theta - r\dot \theta \sin \theta; \tag 2$

$\dot y = \dot r \sin \theta + r \dot \theta \cos \theta; \tag 3$

$\begin{pmatrix} \dot x \\ \dot y \end{pmatrix} = \begin{bmatrix} \cos \theta & -r \sin \theta \\ \sin \theta & r\cos \theta \end{bmatrix} \begin{pmatrix} \dot r \\ \dot \theta \end{pmatrix}; \tag 4$

$\begin{bmatrix} \cos \theta & -r \sin \theta \\ \sin \theta & r\cos \theta \end{bmatrix}^{-1} = \dfrac{1}{r} \begin{bmatrix} r\cos \theta & r\sin \theta \\ -\sin \theta & \cos \theta \end{bmatrix} ; \tag 5$

$\begin{pmatrix} \dot r \\ \dot \theta \end{pmatrix} = \dfrac{1}{r} \begin{bmatrix} r\cos \theta & r\sin \theta \\ -\sin \theta & \cos \theta \end{bmatrix}\begin{pmatrix} \dot x \\ \dot y \end{pmatrix}; \tag 6$

$\dot x = (x - y)(1 - x^2 - y^2) = r(\cos \theta - \sin \theta)(1 - r^2); \tag 7$

$\dot y = (x + y)(1 - x^2 - y^2) = r(\cos \theta + \sin \theta)(1 - r^2); \tag 8$

$\begin{pmatrix} \dot r \\ \dot \theta \end{pmatrix} = \dfrac{1}{r} \begin{bmatrix} r\cos \theta & r\sin \theta \\ -\sin \theta & \cos \theta \end{bmatrix}\begin{pmatrix} r(\cos \theta - \sin \theta)(1 - r^2) \\ r(\cos \theta + \sin \theta)(1 - r^2)\end{pmatrix}$ $= (1 - r^2) \begin{bmatrix} r\cos \theta & r\sin \theta \\ -\sin \theta & \cos \theta \end{bmatrix} \begin{pmatrix} \cos \theta - \sin \theta \\ \cos \theta + \sin \theta \end{pmatrix} = (1 - r^2) \begin{pmatrix} r \\ 1 \end{pmatrix}; \tag{10}$

$\dot r = r(1 - r^2); \tag{11}$

$\dot \theta = 1 - r^2. \tag{12}$

Desde

$x_0^2 + y_0^2 > 0, \tag{13}$

tenemos

$r_0^2 = x_0^2 + y_0^2 > 0 \Longrightarrow r_0 > 0; \tag{14}$

por (11),

$0 < r < 1 \Longrightarrow \dot r > 0; \tag{15}$

$1 < r \Longrightarrow \dot r < 0; \tag{16}$

$r = 1 \Longrightarrow \dot r = 0; \tag{17}$

es relativamente fácil ver que (14)-(17) en concierto implica que

$\displaystyle \lim_{t \to \infty} r = 1; \tag{18}$

también, (11) y (12), junto rendimiento

$r \ne 0 \Longrightarrow \dot{(\ln r)} = \dfrac{\dot r}{r} = 1 - r^2 = \dot \theta, \tag{19}$

de donde, integrando más de $t$,

$\ln r - \ln r_0 = \theta - \theta_0, \tag{20}$

o

$\theta = \ln r - \ln r_0 + \theta_0, \tag{21}$

y así, a través de (18)

$\displaystyle \lim_{t \to \infty} \theta = \lim_{t \to \infty} \ln r - \ln r_0 + \theta_0 = \theta_0 - \ln r_0; \tag{22}$

finalmente,

$\displaystyle \lim_{t \to \infty} x(t) = \lim_{t \to \infty} ( r \cos \theta) = (\lim_{t \to \infty} r)( \lim_{t \to \infty} \cos \theta) = \cos (\theta_0 - \ln r_0). \tag{23}$

Answer 2

Vamos $r^2 = x^2+y^2$, $x = rcos\theta$, y $y = rsin\theta$.

A continuación, se obtiene las siguientes ecuaciones: $$\begin{align}\frac{dr}{dt} &= r(1-r^2) \\ \frac{d\theta}{dt} &= (1-r^2) \end{align}$$

La solución de la primera ecuación es: $r^2=\frac{ke^{2t}}{1+ke^{2t}}$ donde $k=\frac{r^2_0}{1-r_0^2}$.

Entonces, la solución de la segunda ecuación es: $\theta-\theta_0=log\ \frac{r}{r_0}$,

donde$r_0^2=x_0^2+y_0^2$$tan\ \theta_0= \frac{y_0}{x_0}$.

Answer 3

Si sólo desea calcular el límite, usted ni siquiera necesita para resolver la ecuación explícita (al menos para $r$).

Seguir a @player100 sugerencia, puede transformar el sistema en $$ \begin{cases} \dot r = r(1-r^2)\\ \dot \theta = 1-r^2 \end{casos}, $$ con $r^2 = x^2+y^2$ $\theta = \arctan\frac{y}{x}.$

A continuación, podemos observar que nuestro sistema descansa en el origen, así como en el círculo unidad. Desde $\dot r > 0$ al $0 < r < 1,$ $\dot r < 0$ al $r>1$, podemos concluir que cualquier solución de partida en $(x_0,y_0)$ $x_0^2+y_0^2 > 0$ serán atraídos a la unidad de círculo, de modo que $$\lim_{t \to \infty} r(t) = 1.$$

Por lo tanto, nuestro límite se convierte en $$ \lim_{t\to \infty} x(t) = \lim_{t\to \infty} r(t)\cos\theta(t) = \lim_{t\to \infty} \cos\theta(t) = \cos \left(\lim_{t\to \infty}\theta(t)\right). $$ Ahora \begin{align} \theta(t) - \theta(t_0) &= \int_{t_0}^t \dot\theta(s) \, ds = \int_{t_0}^t 1-r(s)^2 \, ds, \end{align} y el cambio de la integración de la variable $u = r(s),$, de modo que $$du = \dot r(s) \, ds = r(s)(1-r(s)^2) \, ds = u(1-u^2) \, ds,$$ obtenemos \begin{align} \int_{t_0}^t 1-(r(s))^2 \, ds &= \int_{r(t_0)}^{r(t)} \frac{1}{u} \, du = \log(r(t))-\log(r(t_0)). \end{align} Esto le da \begin{align} \lim_{t\to \infty}\theta(t) &= \lim_{t\to \infty}\theta(t) - \theta(t_0) + \theta(t_0) = \lim_{t\to \infty} \log(r(t))-\log(r(t_0)) + \theta(t_0)\\ &= \log(1)-\log(r(t_0)) + \theta(t_0) = \arctan{\frac{y_0}{x_0}} - \log\sqrt{x_0^2+y_0^2}, \end{align} así que $$\lim_{t\to \infty} x(t) = \cos \left( \arctan{\frac{y_0}{x_0}} - \log\sqrt{x_0^2+y_0^2} \right). $$