Articulo de referencia

Transformación de dispersión inversa

El algoritmo de 3 pasos: transformar la solución inicial en datos de dispersión iniciales, evolucionar los datos de dispersión iniciales, transformar los datos de dispersión evo...

El algoritmo de 3 pasos: transformar la solución inicial en datos de dispersión iniciales, evolucionar los datos de dispersión iniciales, transformar los datos de dispersión evolucionados en una solución evolucionada.

En matemáticas , la transformada de dispersión inversa (o transformada de Fourier no lineal ) es un método que resuelve el problema de valor inicial para una ecuación diferencial parcial no lineal utilizando métodos matemáticos relacionados con la dispersión de ondas . [ 1 ] : 4960 La transformada de dispersión directa describe cómo una función dispersa ondas o genera estados ligados . [ 2 ] : 39–43 La transformada de dispersión inversa utiliza datos de dispersión de ondas para construir la función responsable de la dispersión de ondas. [ 2 ] : 66–67 Las transformadas de dispersión directa e inversa son análogas a las transformadas de Fourier directa e inversa que se utilizan para resolver ecuaciones diferenciales parciales lineales . [ 2 ] : 66–67

Mediante un par de operadores diferenciales , un algoritmo de tres pasos puede resolver ecuaciones diferenciales no lineales : la solución inicial se transforma en datos de dispersión (transformación de dispersión directa), los datos de dispersión evolucionan en el tiempo (evolución temporal) y, finalmente, los datos de dispersión reconstruyen la solución en el tiempo (transformación de dispersión inversa). [ 2 ] : 66–67

Este algoritmo simplifica la resolución de una ecuación diferencial parcial no lineal a la resolución de dos ecuaciones diferenciales ordinarias lineales y una ecuación integral ordinaria , un método que, en última instancia, conduce a soluciones analíticas para muchas ecuaciones diferenciales parciales no lineales que, de otro modo, serían difíciles de resolver. [ 2 ] : 72

El problema de dispersión inversa es equivalente a un problema de factorización de Riemann-Hilbert , al menos en el caso de ecuaciones de una dimensión espacial. [ 3 ] Esta formulación puede generalizarse a operadores diferenciales de orden mayor que dos y también a problemas periódicos. [ 4 ] En dimensiones espaciales superiores se tiene, en cambio, un problema de factorización de Riemann-Hilbert "no local" (con convolución en lugar de multiplicación) o un problema d-barra.

Historia

La transformada de dispersión inversa surgió del estudio de ondas solitarias. JS Russell describió una "onda de traslación" u "onda solitaria" que se produce en aguas poco profundas. [ 5 ] Primero JV Boussinesq y más tarde D. Korteweg y G. deVries descubrieron la ecuación de Korteweg-deVries (KdV) , una ecuación diferencial parcial no lineal que describe estas ondas. [ 5 ] Posteriormente, N. Zabusky y M. Kruskal, utilizando métodos numéricos para investigar el problema de Fermi-Pasta-Ulam-Tsingou , encontraron que las ondas solitarias tenían las propiedades elásticas de partículas en colisión; las amplitudes y velocidades iniciales y finales de las ondas permanecieron inalteradas después de las colisiones de ondas. [ 5 ] Estas ondas de tipo partícula se llaman solitones y surgen en ecuaciones no lineales debido a un equilibrio débil entre los efectos dispersivos y no lineales. [ 5 ]

Gardner, Greene, Kruskal y Miura introdujeron la transformada de dispersión inversa para resolver la ecuación de Korteweg-de Vries . [ 6 ] Lax, Ablowitz, Kaup, Newell y Segur generalizaron este enfoque, lo que permitió resolver otras ecuaciones no lineales, incluyendo la ecuación de Schrödinger no lineal , la ecuación de seno-Gordon , la ecuación de Korteweg-de Vries modificada , la ecuación de Kadomtsev-Petviashvili , la ecuación de Ishimori , la ecuación de la red de Toda y la ecuación de Dym . [ 5 ] [ 7 ] [ 8 ] Este enfoque también se ha aplicado a diferentes tipos de ecuaciones no lineales, incluyendo ecuaciones diferenciales-diferenciales, de diferencias parciales, multidimensionales y sistemas no lineales fraccionarios integrables. [ 5 ]

Descripción

ecuación diferencial parcial no lineal

Las variables independientes son una variable espacial.incógnita{\displaystyle x}y una variable de tiempot{\displaystyle t}. Subíndices u operadores diferenciales (incógnita,t{\textstyle \partial _ {x}, \partial _ {t}}) indican diferenciación. La función(incógnita,t){\displaystyle u(x,t)}es una solución de una ecuación diferencial parcial no lineal,t+norte()=0{\textstyle u_{t}+N(u)=0}, con condición inicial (valor)(incógnita,0){\textstyle u(x,0)}. [ 2 ] : 72

Requisitos

La solución de la ecuación diferencial cumple las condiciones de integrabilidad y Fadeev: [ 2 ] : 40

Condición de integrabilidad: |(incógnita)| dincógnita <{\displaystyle \int _{-\infty }^{\infty }\ |u(x)|\ dx\ <\infty }
Condición de Fadeev: (1+|incógnita|)|(incógnita)| dincógnita <{\displaystyle \int _{-\infty }^{\infty }\ (1+|x|)|u(x)|\ dx\ <\infty }

Par de operadores diferenciales

Los operadores diferenciales de Lax ,L{\textstyle L}yMETRO{\textstyle M}, son operadores diferenciales ordinarios lineales con coeficientes que pueden contener la función(incógnita,t){\textstyle u(x,t)}o sus derivados. El operador autoadjuntoL{\textstyle L}tiene derivada temporalLt{\textstyle L_{t}}y genera una ecuación de valores propios (espectral) con funciones propias.ψ{\textstyle \psi }y valores propios constantes en el tiempo ( parámetros espectrales )λ{\textstyle \lambda }. [ 1 ] : 4963 [ 2 ] : 98

L(ψ)=λψ, {\displaystyle L(\psi )=\lambda \psi ,\ }y Lt(ψ)=dmiF(L(ψ))tL(ψt){\textstyle \ L_{t}(\psi ){\overset {def}{=}}(L(\psi ))_{t}-L(\psi _{t})}

El operadorMETRO{\textstyle M}describe cómo evolucionan las funciones propias a lo largo del tiempo y genera una nueva función propia.ψ~{\textstyle {\widetilde {\psi }}}del operadorL{\textstyle L}a partir de la autofunciónψ{\textstyle \psi }deL{\textstyle L}. [ 1 ] : 4963

ψ~=ψtMETRO(ψ) {\displaystyle {\widetilde {\psi }}=\psi _{t}-M(\psi )\ }

Los operadores de Lax se combinan para formar un operador multiplicativo, no un operador diferencial, de las autofunciones.ψ{\textstyle \psi }. [ 1 ] : 4963

(Lt+LMETROMETROL)ψ=0{\displaystyle (L_{t}+LM-ML)\psi =0}

Los operadores de Lax se eligen para que el operador multiplicativo sea igual a la ecuación diferencial no lineal. [ 1 ] : 4963

Lt+LMETROMETROL=t+norte()=0{\displaystyle L_{t}+LM-ML=u_{t}+N(u)=0}

Los operadores diferenciales AKNS , desarrollados por Ablowitz, Kaup, Newell y Segur, son una alternativa a los operadores diferenciales de Lax y logran un resultado similar. [ 1 ] : 4964 [ 9 ] [ 10 ]

Transformación de dispersión directa

La transformada de dispersión directa genera datos de dispersión iniciales; estos pueden incluir los coeficientes de reflexión, el coeficiente de transmisión, los datos de valores propios y las constantes de normalización de las soluciones de autofunción para esta ecuación diferencial. [ 2 ] : 39–48

L(ψ)=λψ{\displaystyle L(\psi )=\lambda \psi }

Evolución temporal de los datos de dispersión

Las ecuaciones que describen cómo evolucionan los datos de dispersión a lo largo del tiempo aparecen como soluciones de una ecuación diferencial ordinaria lineal de primer orden con respecto al tiempo. Utilizando diferentes enfoques, esta ecuación diferencial lineal de primer orden puede surgir de los operadores diferenciales lineales (par de Lax, par de AKNS), una combinación de los operadores diferenciales lineales y la ecuación diferencial no lineal, o a través de operaciones adicionales de sustitución, integración o diferenciación. Ecuaciones espacialmente asintóticas (incógnita±{\textstyle x\to \pm \infty }) simplificar la resolución de estas ecuaciones diferenciales. [ 1 ] : 4967–4968 [ 2 ] : 68–72 [ 6 ]

Transformación de dispersión inversa

La ecuación de Marchenko combina los datos de dispersión en una ecuación integral lineal de Fredholm . La solución de esta ecuación integral conduce a la solución, u(x,t), de la ecuación diferencial no lineal. [ 2 ] : 48–57

Ejemplo: ecuación de Korteweg-De Vries

La ecuación diferencial no lineal de Korteweg-De Vries es [ 11 ] : 4

t6incógnita+incógnitaincógnitaincógnita=0{\displaystyle u_{t}-6uu_{x}+u_{xxx}=0}

Operadores de Lax

Los operadores de Lax son: [ 2 ] : 97–102

L=incógnita2+(incógnita,t) {\displaystyle L=-\partial _{x}^{2}+u(x,t)\ } y METRO=4incógnita3+6incógnita+3incógnita{\textstyle \ M=-4\partial _{x}^{3}+6u\partial _{x}+3u_{x}}

El operador multiplicativo es:

Lt+LMETROMETROL=t6incógnita+incógnitaincógnitaincógnita=0{\displaystyle L_{t}+LM-ML=u_{t}-6uu_{x}+u_{xxx}=0}

Transformación de dispersión directa

Las soluciones a esta ecuación diferencial

L(ψ)=ψincógnitaincógnita+(incógnita,0)ψ=λψ{\textstyle L(\psi )=-\psi _{xx}+u(x,0)\psi =\lambda \psi }

Puede incluir soluciones de dispersión con un rango continuo de valores propios ( espectro continuo ) y soluciones de estado ligado con valores propios discretos ( espectro discreto ). Los datos de dispersión incluyen coeficientes de transmisión.T(k,0){\textstyle T(k,0)}, coeficiente de reflexión izquierdaRL(k,0){\textstyle R_{L}(k,0)}, coeficiente de reflexión derechaRR(k,0){\textstyle R_{R}(k,0)}, valores propios discretosκ12,,κnorte2{\textstyle -\kappa _ {1}^{2},\ldots,-\kappa _ {N}^{2}}y constantes de normalización (normanización) de estados ligados izquierdo y derecho . [ 1 ] : 4960

do(0)Lj=( ψL2(ikj,incógnita,0) dincógnita)1/2 j=1,,norte{\displaystyle c(0)_{Lj}=\left(\int _{-\infty }^{\infty }\ \psi _{L}^{2}(ik_{j},x,0)\ dx\right)^{-1/2}\ j=1,\dots ,N}
do(0)Rj=( ψR2(ikj,incógnita,0) dincógnita)1/2 j=1,,norte{\displaystyle c(0)_{Rj}=\left(\int _{-\infty }^{\infty }\ \psi _{R}^{2}(ik_{j},x,0)\ dx\right)^{-1/2}\ j=1,\dots ,N}

Evolución temporal de los datos de dispersión

La izquierda espacialmente asintóticaψL(k,incógnita,t){\textstyle \psi _{L}(k,x,t)}y derechaψR(k,incógnita,t){\textstyle \psi _{R}(k,x,t)}Las funciones de Jost simplifican este paso. [ 1 ] : 4965–4966

ψL(incógnita,k,t)=miikincógnita+o(1), incógnita+ψL(incógnita,k,t)=miikincógnitaT(k,t)+RL(k,t)miikincógnitaT(k,t)+o(1), incógnitaψR(incógnita,k,t)=miikincógnitaT(k,t)+RR(k,t)miikincógnitaT(k,t)+o(1), incógnita+ψR(incógnita,k,t)=miikincógnita+o(1), incógnita{\displaystyle {\begin{aligned}\psi _{L}(x,k,t)&=e^{ikx}+o(1),\ x\to +\infty \\\psi _{L}(x,k,t)&={\frac {e^{ikx}}{T(k,t)}}+{\frac {R_{L}(k,t)e^{-ikx}}{T(k,t)}}+o(1),\ x\to -\infty \\\psi _{R}(x,k,t)&={\frac {e^{-ikx}}{T(k,t)}}+{\frac {R_{R}(k,t)e^{ikx}}{T(k,t)}}+o(1),\ x\to +\infty \\\psi _{R}(x,k,t)&=e^{-ikx}+o(1),\ x\to -\infty \\\end{aligned}}}

Las constantes de dependenciaγj(t){\textstyle \gamma _{j}(t)}relacionar las funciones de Jost derecha e izquierda y las constantes de normalización derecha e izquierda. [ 1 ] : 4965–4966

γj(t)=ψL(incógnita,iκj,t)ψR(incógnita,iκj,t)=(1)nortejdoRj(t)doLj(t){\displaystyle \gamma _{j}(t)={\frac {\psi _{L}(x,i\kappa _{j},t)}{\psi _{R}(x,i\kappa _{j},t)}}=(-1)^{N-j}{\frac {c_{Rj}(t)}{c_{Lj}(t)}}}

El LaxMETRO{\textstyle M}El operador diferencial genera una autofunción que puede expresarse como una combinación lineal dependiente del tiempo de otras autofunciones. [ 1 ] : 4967

tψL(k,incógnita,t)METROψL(incógnita,k,t)=aL(k,t)ψL(incógnita,k,t)+bL(k,t)ψR(incógnita,k,t){\displaystyle \partial _{t}\psi _{L}(k,x,t)-M\psi _{L}(x,k,t)=a_{L}(k,t)\psi _{L}(x,k,t)+b_{L}(k,t)\psi _{R}(x,k,t)}
tψR(k,incógnita,t)METROψR(incógnita,k,t)=aR(k,t)ψL(incógnita,k,t)+bR(k,t)ψR(incógnita,k,t){\displaystyle \partial _{t}\psi _{R}(k,x,t)-M\psi _{R}(x,k,t)=a_{R}(k,t)\psi _{L}(x,k,t)+b_{R}(k,t)\psi _{R}(x,k,t)}

Las soluciones a estas ecuaciones diferenciales, determinadas mediante funciones de Jost espacialmente asintóticas de dispersión y de estado ligado, indican un coeficiente de transmisión constante en el tiempo.T(k,t){\textstyle T(k,t)}, pero coeficientes de reflexión dependientes del tiempo y coeficientes de normalización. [ 1 ] : 4967–4968

RL(k,t)=RL(k,0)mii8k3tRR(k,t)=RR(k,0)mi+i8k3tdoLj(t)=doLj(0)mi+4κj3t, j=1,,nortedoRj(t)=doRj(0)mi4κj3t, j=1,,norte{\displaystyle {\begin{aligned}R_{L}(k,t)&=R_{L}(k,0)e^{-i8k^{3}t}\\R_{R}(k,t)&=R_{R}(k,0)e^{+i8k^{3}t}\\c_{Lj}(t)&=c_{Lj}(0)e^{+4\kappa _{j}^{3}t},\ j=1,\ldots ,N\\c_{Rj}(t)&=c_{Rj}(0)e^{-4\kappa _{j}^{3}t},\ j=1,\ldots ,N\end{aligned}}}

Transformación de dispersión inversa

El núcleo de Marchenko esF(incógnita,t){\textstyle F(x,t)}. [ 1 ] : 4968–4969

F(incógnita,t)=dmiF12πRR(k,t)miikincógnita dk+j=1nortedo(t)Lj2miκjincógnita{\displaystyle F(x,t){\overset {def}{=}}{\frac {1}{2\pi }}\int _{-\infty }^{\infty }R_{R}(k,t)e^{ikx}\ dk+\sum _{j=1}^{N}c(t)_{Lj}^{2}e^{-\kappa _{j}x}}

La ecuación integral de Marchenko es una ecuación integral lineal que se resuelve paraK(incógnita,y,t){\textstyle K(x,y,t)}. [ 1 ] : 4968–4969

K(incógnita,z,t)+F(incógnita+z,t)+incógnitaK(incógnita,y,t)F(y+z,t) dy=0{\displaystyle K(x,z,t)+F(x+z,t)+\int _{x}^{\infty }K(x,y,t)F(y+z,t)\ dy=0}

La solución a la ecuación de Marchenko,K(incógnita,y,t){\textstyle K(x,y,t)}genera la solución(incógnita,t){\textstyle u(x,t)}a la ecuación diferencial parcial no lineal. [ 1 ] : 4969

(incógnita,t)=2K(incógnita,y,t)incógnita{\displaystyle u(x,t)=-2{\frac {\partial K(x,y,t)}{\partial x}}}

Ejemplos de ecuaciones integrables

Véase también

Citas

Referencias

  • Ablowitz, MJ; Kaup, DJ; Newell, AC; Segur, H. (1973). "Método para resolver la ecuación de Sine-Gordon" . Physical Review Letters . 30 (25): 1262– 1264. Bibcode : 1973PhRvL..30.1262A . doi : 10.1103/PhysRevLett.30.1262 .
  • Ablowitz, MJ; Kaup, DJ; Newell, AC; Segur, H. (1974). "La transformada inversa de dispersión: análisis de Fourier para problemas no lineales" . Estudios en matemáticas aplicadas . 53 (4): 249– 315. doi : 10.1002/sapm1974534249 .
  • Ablowitz, Mark J.; Segur, Harvey (1981). Solitones y la transformada de dispersión inversa . SIAM. ISBN 978-0-89871-477-7.
  • Ablowitz, Mark J.; Fokas, AS (2003). Variables complejas: Introducción y aplicaciones . Cambridge University Press. pp. 604–620 . ISBN  978-0-521-53429-1.
  • Ablowitz, Mark J. (2023). "Ondas no lineales y la transformada de dispersión inversa" . Optik . 278 170710. Bibcode : 2023Optik.27870710A . doi : 10.1016/j.ijleo.2023.170710 .
  • Aktosun, Tuncay (2009). «Transformación de dispersión inversa y la teoría de los solitones». Enciclopedia de la complejidad y la ciencia de los sistemas . Springer. pp. 4960–4971 . arXiv : 0905.4746 . doi : 10.1007/978-0-387-30440-3_295 . ISBN  978-0-387-30440-3.
  • Drazin, PG; Johnson, RS (1989). Solitones: Una introducción . Cambridge University Press. ISBN 978-0-521-33655-0.
  • Gardner, Clifford S.; Greene, John M.; Kruskal, Martin D.; Miura, Robert M. (1967). "Método para resolver la ecuación de Korteweg-de Vries" . Physical Review Letters . 19 (19): 1095– 1097. Bibcode : 1967PhRvL..19.1095G . doi : 10.1103/PhysRevLett.19.1095 .
  • Konopelchenko, BG; Dubrowsky, VG (1991). «Solitones localizados para la ecuación de Ishimori». En Sattinger, David H.; Tracy, CA; Venakides, Stephanos (eds.). Dispersión inversa y aplicaciones . American Mathematical Soc. pp. 77–90 . ISBN  978-0-8218-5129-6.
  • Oono, H. (1996). «Solución de N solitones de la ecuación de Harry Dym mediante el método de dispersión inversa». En Alfinito, E.; Boiti, M.; Martina, L. (eds.). Física no lineal: teoría y experimento . World Scientific Publishing Company Pte Limited. pp. 241–248 . ISBN  978-981-02-2559-9.
  • Osborne, AR (1995). "Física de solitones y la transformada de dispersión inversa periódica" . Physica D: Nonlinear Phenomena . 86 (1): 81– 89. Bibcode : 1995PhyD...86...81O . doi : 10.1016/0167-2789(95)00089-M . ISSN 0167-2789 . 

Lecturas adicionales

  • Ablowitz, Mark J.; Clarkson, PA (12 de diciembre de 1991). Solitones, ecuaciones de evolución no lineales y dispersión inversa . Cambridge University Press. ISBN 978-0-521-38730-9.
  • Bullough, RK; Caudrey, PJ (11 de noviembre de 2013). Solitones . Springer Science & Business Media. ISBN 978-3-642-81448-8.
  • Gardner, Clifford S.; Greene, John M.; Kruskal, Martin D.; Miura, Robert M. (1974), "Ecuación de Korteweg-de Vries y generalización. VI. Métodos para la solución exacta.", Comm. Pure Appl. Math. , 27 (1): 97– 133, Bibcode : 1974CPAM...27...97G , doi : 10.1002/cpa.3160270108 , MR 0336122 
  • Gelʹfand, Izrailʹ Moiseevich (1955). Sobre la determinación de una ecuación diferencial a partir de su función espectral . American Mathematical Society. págs.  253-304.
  • Marchenko, Vladimir A. (1986). Operadores de Sturm-Liouville y aplicaciones . Teoría de operadores: avances y aplicaciones. Vol.  22. Basilea: Birkhäuser. doi : 10.1007/978-3-0348-5485-6 . ISBN 978-3-0348-5486-3.
  • Shaw, JK (1 de mayo de 2004). Principios matemáticos de la comunicación por fibra óptica . SIAM. ISBN 978-0-89871-556-9.
  • "Artículo matemático introductorio sobre IST" (PDF) . (300 KiB ) 
  • Transformada de dispersión inversa y la teoría de los solitones