Cómo configurar y resolver ecuaciones diferenciales para el problema de los cuatro errores

Supongo que ‘entiendes’ que las partículas se encuentran en el centro del cuadrado, y que el cuadrado en sí mismo se encoge con el tiempo, llegando finalmente al área cero cuando las partículas se encuentran.

La solución más simple (Velocidad relativa) Considera que te estás moviendo con una de las partículas. Desde su marco de referencia, una partícula se aleja de usted, mientras que otra partícula lo persigue. La partícula que se aleja de usted lo hace en una dirección perpendicular a la línea que une las partículas, por lo tanto, no aumenta la distancia entre las partículas. La partícula que se mueve hacia ti lo hace en una dirección a lo largo de la línea que une las dos partículas. El resultado neto es que ves una partícula moverse hacia ti con una velocidad de [matemáticas] v [/ matemáticas].

Como la distancia inicial entre las partículas adyacentes era [matemática] a [/ matemática], el tiempo que tardan las partículas en encontrarse es [matemática] a / v [/ matemática].

Enfoque matemático I (ecuaciones diferenciales) Considere el siguiente bosquejo. Muestra las posiciones de las partículas (en algún momento [matemáticas] t [/ matemáticas]) y las posiciones cambiadas después de algún tiempo [matemáticas] \ Delta t [/ matemáticas].

[ La imagen es solo para fines ilustrativos, por lo tanto, es aproximadamente correcta. Para una estimación más precisa de la ruta tomada por cada partícula, vea la gráfica de MATLAB a continuación. ]

Queremos descubrir cómo cambia la longitud de un lado del cuadrado con el tiempo.

Aplicando el teorema de Pitágoras en el triángulo ABC,

[matemáticas] x (t + \ Delta t) ^ 2 = (x (t) -v \ Delta t) ^ 2 + (v \ Delta t) ^ 2 [/ matemáticas]

En expansión,

[matemáticas] \ displaystyle \ implica x (t + \ Delta t) ^ 2 = x (t) ^ 2 + (v \ Delta t) ^ 2 -2v x (t) \ Delta t + (v \ Delta t) ^ 2 [/matemáticas]

Deje que [math] x (t) ^ 2 = A (t) [/ math] [math] [/ math] representa el área del cuadrado en [math] t [/ math]. Usando esto en lo anterior

[matemáticas] \ displaystyle \ implica A (t + \ Delta t) = A (t) – 2v \ sqrt {A (t)} \ Delta t + 2v ^ 2 (\ Delta t) ^ 2 [/ matemáticas]

[matemáticas] \ displaystyle \ implica A (t + \ Delta t) – A (t) = – 2v \ sqrt {A (t)} \ Delta t + 2v ^ 2 (\ Delta t) ^ 2 [/ matemática]

[matemáticas] \ displaystyle \ implica \ frac {A (t + \ Delta t) – A (t)} {\ Delta t} = – 2 \ sqrt {A (t)} v + 2v ^ 2 \ Delta t [/ math ]

Tomando el límite [matemática] \ Delta t \ rightarrow 0 [/ matemática], el LHS se convierte en una derivada y el último término en el RHS se convierte en cero.

[matemáticas] \ displaystyle \ implica \ frac {dA (t)} {dt} = -2v \ sqrt {A (t)} [/ matemáticas]

Esto produce una ecuación diferencial

[matemáticas] \ displaystyle \ implica \ frac {dA (t)} {\ sqrt {A (t)}} = -2vdt [/ matemáticas]

que puede evaluarse para el cambio en la longitud del lado del cuadrado con el tiempo.

[matemáticas] \ displaystyle \ implica \ int_ {a ^ 2} ^ {{x (t)} ^ 2} \ frac {dA (t)} {\ sqrt {A (t)}} = -2v \ int_ {0 } ^ {t} dt [/ matemáticas]

[matemáticas] \ displaystyle \ implica 2 \ sqrt {A (t)} | _ {a ^ 2} ^ {{x (t)} ^ 2} = -2vt [/ math]

[matemáticas] \ displaystyle \ implica | x (t) | – a = -vt [/ matemáticas]

Podemos encontrar el tiempo en el que las partículas se encuentran estableciendo [math] x (t = t_0) = 0 [/ math], que da [math] t_0 = a / v [/ math], lo mismo que arriba.

Enfoque matemático II (ecuaciones diferenciales acopladas)

La respuesta de Hare Krishna Singh es agradable. Este método tiene la ventaja adicional de que nos da el camino tomado por las partículas también explícitamente.

Doy el código MATLAB utilizado para resolver las cuatro ecuaciones diferenciales acopladas y el gráfico correspondiente.

función dx = fourBugs_x (t, x)

Función de espacio de estado%

dx = ceros (4,1);
dx (1) = x (2) -x (1);
dx (2) = x (3) -x (2);
dx (3) = x (4) -x (3);
dx (4) = x (1) -x (4);
final

función dy = fourBugs_y (t, y)

Función de espacio de estado%

dy = ceros (4,1);
dy (1) = y (2) -y (1);
dy (2) = y (3) -y (2);
dy (3) = y (4) -y (3);
dy (4) = y (1) -y (4);
final

% Para escribir en la ventana de comandos
% Asumo que los errores están en las esquinas de un cuadrado de 10 × 10 inicialmente
[t, x] = ode45 (@fourBugs_x, [0 100], [0 10 10 0]);
[t, y] = ode45 (@fourBugs_y, [0 100], [0 0 10 10]);
plot (x (:, 1), y (:, 1)); Espere; plot (x (:, 2), y (:, 2)); Espere; plot (x (:, 3), y (:, 3)); Espere; plot (x (:, 4), y (:, 4));

La figura generada correspondiente.

Supongamos que en el momento [math] \ displaystyle {t = 0} [/ math] los cuatro errores se ubicaron en los vértices de un cuadrado, es decir: [math] \ displaystyle {A, B, C, \ text {y} D} [/ matemáticas] respectivamente. Y se movieron a los puntos [math] \ displaystyle {A ‘, B’, C ‘} [/ math] y [math] \ displaystyle {D’} [/ math] en el momento [math] \ displaystyle {t> 0} [/matemáticas].

Debido a la simetría del problema, se puede ver que los cuatro puntos [matemática] \ displaystyle {A ‘, B’, C ‘, D’} [/ matemática] formaron un cuadrado. En particular, el vector [math] \ displaystyle {\ vec {A’B ‘}} [/ math] debe ser colineal con la velocidad del error en el punto [math] \ displaystyle {A’} [/ math]. De manera análoga, podemos concluir lo mismo con los puntos [math] \ displaystyle {B ‘, C’ \ text {and} D ‘} [/ math].

Además, nuevamente debido a la simetría, nos damos cuenta de que el punto [math] \ displaystyle {B ‘} [/ math] es solo la imagen del punto [math] \ displaystyle {A’} [/ math] al girarlo sobre el centro [math ] \ displaystyle {O} [/ math] del cuadrado [math] \ displaystyle {ABCD} [/ math] en sentido antihorario por un ángulo de [math] \ displaystyle {90 ^ o} [/ math]. Dejemos construir un sistema de coordenadas con el punto de origen en [math] \ displaystyle {O} [/ math] como se muestra en la figura a continuación.

En este sistema, suponga que [math] \ displaystyle {A ‘} [/ math] tiene coordenadas [math] \ displaystyle {A’ = (x, y)} [/ math]. Luego, a partir del argumento anterior, se pueden deducir las coordenadas del punto [matemáticas] \ displaystyle {B ‘= (- y, x)} [/ matemáticas]. Y [math] \ displaystyle {\ vec {A’B ‘} = (-xy, xy)} [/ math] debe ser tangente a la trayectoria del error en el punto [math] \ displaystyle {A’} \ qquad (1 )[/matemáticas].

Además de la velocidad del error en el punto [math] \ displaystyle {A ‘} [/ math], por definición, se da como:

[matemáticas] \ displaystyle {\ vec {v} = (x’_t, y’_t) = (x ‘, y’)} \ qquad (2) [/ math]

Combine [math] \ displaystyle {(1), (2)} [/ math], uno debe tener [math] \ displaystyle {\ vec {A’B ‘}} [/ math] es colineal con [math] \ displaystyle {\ vec {v}} [/ math] en el punto [math] \ displaystyle {A ‘= (x, y)} [/ math]. En otras palabras:

[matemáticas] \ displaystyle {\ frac {-xy} {x ‘} = \ frac {xy} {y’}} [/ matemáticas]

[matemáticas] \ displaystyle {\ Rightarrow \ frac {y ‘} {x’} = \ frac {-xy} {xy}} [/ math]

Esto simplemente se convierte en:

[matemáticas] \ displaystyle {\ frac {\ mathrm {d} y} {\ mathrm {d} x} = \ frac {yx} {x + y}} \ qquad (3) [/ math]

Resulta que la ecuación [math] \ displaystyle {(3)} [/ math] es una ODE homogénea. Al introducir [math] \ displaystyle {y = ux} [/ math], lo reescribimos de la siguiente manera:

[matemáticas] \ displaystyle {\ frac {(u + 1) \ mathrm {d} u} {u ^ 2 + 1} = – \ frac {\ mathrm {d} x} {x}} \ qquad (4) [ /matemáticas]

La ecuación de integración [matemáticas] \ displaystyle {(4)} [/ matemáticas] da:

[matemáticas] \ displaystyle {\ frac {1} {2} \ mathrm {ln} (1 + u ^ 2) + \ mathrm {arctan} u = – \ mathrm {ln} | x | + C} \ qquad (5) [/ matemáticas]

donde [math] \ displaystyle {C} [/ math] es una constante que se calculará utilizando la condición inicial.

Dado que como [matemáticas] \ displaystyle {t = 0 \ Rightarrow x = \ frac {a} {\ sqrt {2}}, y = 0} [/ math] donde [matemáticas] \ displaystyle {AB = BC = CD = DA = a} [/ matemáticas]. Luego de la ecuación (5) obtenemos:

[matemáticas] \ displaystyle {0 = – \ mathrm {ln} \ frac {a} {\ sqrt {2}} + C \ Rightarrow C = \ mathrm {ln} \ frac {a} {\ sqrt {2}}} [/matemáticas]

Inserte este resultado en la ecuación [math] \ displaystyle {(5)} [/ math] con atención a [math] \ displaystyle {u = \ frac {y} {x}} [/ math], logramos:

[matemáticas] \ displaystyle {\ frac {1} {2} \ mathrm {ln} \ left (1 + \ frac {y ^ 2} {x ^ 2} \ right) + \ mathrm {arctan} \ frac {y} {x} = – \ mathrm {ln} x + \ mathrm {ln} \ frac {a} {\ sqrt {2}}} [/ math]

Por lo tanto,

[matemáticas] \ displaystyle {y = x \ cdot \ mathrm {tan} \ left (\ mathrm {ln} \ left (\ frac {a} {\ sqrt {2 (x ^ 2 + y ^ 2)}} \ right ) \ right)} [/ math]

es la ecuación de la trayectoria del error a partir de [math] \ displaystyle {A} [/ math]. Las trayectorias de los otros errores se pueden obtener directamente al rotar esto sobre el punto [math] \ displaystyle {O} [/ math] por [math] \ displaystyle {90 ^ o, 180 ^ o \ text {y} 270 ^ o} [ / matemáticas] grado en sentido antihorario respectivamente.

También podemos expresar la trayectoria en forma de ecuaciones de parámetros de la siguiente manera:

[matemáticas] \ displaystyle {x = \ pm \ frac {a} {\ sqrt {2 (1 + u ^ 2)} e ^ {\ mathrm {arctan} u}}} [/ math]

y

[matemáticas] \ displaystyle {y = ux = \ pm \ frac {au} {\ sqrt {2 (1 + u ^ 2)} e ^ {\ mathrm {arctan} u}}} [/ math]

o en un sistema de coordenadas polares configurando [math] \ displaystyle {x = r \ mathrm {cos} \ theta, y = r \ mathrm {sin} \ theta} [/ math], obtenemos una expresión extremadamente simple de la trayectoria de la siguiente manera:

[matemáticas] \ displaystyle {r = \ frac {a} {\ sqrt {2}} e ^ {- \ theta}} [/ matemáticas]

Usando Matlab tenemos:

Y aquí hay un resultado de wolfram alpha con [math] \ displaystyle {a = 10} [/ math]:

Motor de conocimiento computacional

Ahora estamos listos para calcular el tiempo necesario para que los errores se alcancen entre sí. Lo sabemos:

[matemáticas] \ displaystyle {\ mathrm {d} s = \ sqrt {\ mathrm {d} x ^ 2 + \ mathrm {d} y ^ 2}} [/ math]

Como [math] \ displaystyle {x = r \ mathrm {cos} \ theta, y = r \ mathrm {sin} \ theta} [/ math] en el sistema polar. Luego:

[matemáticas] \ displaystyle {\ mathrm {d} s = \ sqrt {r ^ 2 + r ‘^ 2} \ mathrm {d} \ theta} [/ math]

Por lo tanto, la longitud de la trayectoria se puede obtener mediante:

[matemáticas] \ displaystyle {L = \ int_0 ^ {+ \ infty} \ sqrt {r ^ 2 + r ‘^ 2} \ mathrm {d} \ theta = \ int_0 ^ {+ \ infty} \ sqrt {\ frac { a ^ 2} {2} e ^ {- 2 \ theta} + \ frac {a ^ 2} {2} e ^ {- 2 \ theta}} \ mathrm {d} \ theta = \ int_0 ^ {+ \ infty } ae ^ {- \ theta} \ mathrm {d} \ theta = \ left. -ae ^ {- \ theta} \ right | _0 ^ {+ \ infty} = a} [/ math]

Luego:

[matemáticas] \ displaystyle {\ Delta t = \ frac {L} {v} = \ frac {a} {v}} \ qquad \ square [/ math]

Este problema puede modelarse como cuatro conjuntos de ecuaciones diferenciales acopladas.

dx1 / dt = x2-x1

dx2 / dt = x3-x2

dx3 / dt = x4-x3

dx4 / dt = x1-x4

todas las expresiones están en signo de vector. tenemos una condición inicial como x1 (t = 0) = (a, 0) y todas las demás posiciones iniciales y velocidades iniciales (como v1 (t = 0) está al lado del cuadrado). se resuelve explícitamente y todas las velocidades se normalizan a v como se indica en el problema.

Creo que esta es la respuesta que estás buscando.

4 errores persiguiéndose mutuamente la ecuación diferencial

La animación en Los vértices de un triángulo equilátero se está reduciendo entre sí.

para la versión de 3 errores del mismo problema podría ayudar en una buena visualización del problema.

Sin embargo, esto supone que desea hacerlo utilizando solo ecuaciones diferenciales.

Porque notar la simetría llevaría a una solución más pequeña y elegante como la de Los vértices de un triángulo equilátero se están reduciendo entre sí.