Método matemático para ecuaciones integrodiferenciales.
El método Wiener-Hopf es una técnica matemática muy utilizada en matemáticas aplicadas . Fue desarrollado inicialmente por Norbert Wiener y Eberhard Hopf como un método para resolver sistemas de ecuaciones integrales , pero ha encontrado un uso más amplio en la resolución de ecuaciones diferenciales parciales bidimensionales con condiciones de frontera mixtas en la misma frontera. En general, el método funciona explotando las propiedades analíticas complejas de las funciones transformadas. Normalmente, se utiliza la transformada de Fourier estándar , pero existen ejemplos que utilizan otras transformaciones, como la transformada de Mellin .
En general, las ecuaciones rectoras y las condiciones de contorno se transforman y estas transformaciones se utilizan para definir un par de funciones complejas (normalmente denotadas con subíndices '+' y '-') que son respectivamente analíticas en las mitades superior e inferior del plano complejo. , y no tienen un crecimiento más rápido que los polinomios en estas regiones. Estas dos funciones también coincidirán en alguna región del plano complejo , típicamente, una franja delgada que contiene la línea real . La continuación analítica garantiza que estas dos funciones definen una única función analítica en todo el plano complejo, y el teorema de Liouville implica que esta función es un polinomio desconocido , que a menudo es cero o constante. El análisis de las condiciones en los bordes y esquinas del límite permite determinar el grado de este polinomio.
Descomposición de Wiener-Hopf
La ecuación fundamental que aparece en el método de Wiener-Hopf es de la forma
![{\displaystyle A(\alpha )\Xi _ {+}(\alpha )+B(\alpha )\Psi _ {-}(\alpha )+C(\alpha )=0,}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
donde , , son funciones holomorfas conocidas , las funciones , son desconocidas y la ecuación se cumple en una franja en el plano complejo . Encontrarlo es lo que se llama el problema de Wiener-Hopf . ![{\displaystyle A}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle B}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle C}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle \Xi _ {+}(\alpha )}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle \Psi _ {-}(\alpha )}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle \tau _{-}<{\mathfrak {Estoy}}(\alpha )<\tau _{+}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle \alpha }](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle \Xi _ {+}(\alpha )}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle \Psi _ {-}(\alpha )}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
El paso clave en muchos problemas de Wiener-Hopf es descomponer una función arbitraria en dos funciones con las propiedades deseadas descritas anteriormente. En general, esto se puede hacer escribiendo![{\displaystyle \Phi}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle \Phi _{\pm }}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle \Phi _{+}(\alpha )={\frac {1}{2\pi i}}\int _{C_{1}}\Phi (z){\frac {dz}{z- \alfa }}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
y
![{\displaystyle \Phi _{-}(\alpha )=-{\frac {1}{2\pi i}}\int _{C_{2}}\Phi (z){\frac {dz}{z -\alfa }},}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
donde las curvas de nivel y son paralelas a la recta real, pero pasan por encima y por debajo del punto , respectivamente. ![{\ Displaystyle C_ {1}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\ Displaystyle C_ {2}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle z=\alpha}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
De manera similar, las funciones escalares arbitrarias se pueden descomponer en un producto de funciones +/−, es decir , tomando primero el logaritmo y luego realizando una descomposición por suma. Las descomposiciones de productos de funciones matriciales (que ocurren en sistemas multimodales acoplados, como ondas elásticas) son considerablemente más problemáticas ya que el logaritmo no está bien definido y se podría esperar que cualquier descomposición no sea conmutativa. Khrapkov obtuvo una pequeña subclase de descomposiciones conmutativas y también se han desarrollado varios métodos aproximados. [ cita necesaria ]![{\displaystyle K(\alpha )=K_{+}(\alpha )K_{-}(\alpha )}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Ejemplo
Considere la ecuación diferencial parcial lineal.
![{\displaystyle {\boldsymbol {L}}_{xy}f(x,y)=0,}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
donde es un operador lineal que contiene derivadas con respecto a x e y , sujeto a las condiciones mixtas en y = 0, para alguna función prescrita g ( x ) ,![{\displaystyle {\boldsymbol {L}}_{xy}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle f=g(x){\text{ para }}x\leq 0,\quad f_{y}=0{\text{ cuando }}x>0}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
y decaimiento en el infinito, es decir, f → 0 como .![{\displaystyle {\boldsymbol {x}}\rightarrow \infty }](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Tomar una transformada de Fourier con respecto a x da como resultado la siguiente ecuación diferencial ordinaria
![{\displaystyle {\boldsymbol {L}}_{y}{\widehat {f}}(k,y)-P(k,y){\widehat {f}}(k,y)=0,}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
donde es un operador lineal que contiene únicamente derivadas de y , P ( k,y ) es una función conocida de y y k y ![{\displaystyle {\boldsymbol {L}}_{y}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle {\widehat {f}}(k,y)=\int _{-\infty }^{\infty }f(x,y)e^{-ikx}\,{\textrm {d}} X.}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Si una solución particular de esta ecuación diferencial ordinaria que satisface la desintegración necesaria en el infinito se denota por F ( k , y ) , una solución general se puede escribir como
![{\displaystyle {\widehat {f}}(k,y)=C(k)F(k,y),}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
donde C ( k ) es una función desconocida que se determinará mediante las condiciones de frontera en y =0.
La idea clave es dividir en dos funciones separadas, que son analíticas en las mitades inferior y superior del plano complejo, respectivamente.![{\displaystyle {\widehat {f}}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle {\widehat {f}}_{+}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle {\widehat {f}}_{-}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle {\widehat {f}}_{+}(k,y)=\int _{0}^{\infty }f(x,y)e^{-ikx}\,{\textrm {d }}X,}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle {\widehat {f}}_{-}(k,y)=\int _{-\infty }^{0}f(x,y)e^{-ikx}\,{\textrm { d}}x.}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Las condiciones de contorno dan entonces
![{\displaystyle {\widehat {g\,}}(k)+{\widehat {f}}_{+}(k,0)={\widehat {f}}_{-}(k,0)+ {\widehat {f}}_{+}(k,0)={\widehat {f}}(k,0)=C(k)F(k,0)}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
y, al tomar derivados con respecto a ,![{\displaystyle y}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle {\widehat {f}}'_{-}(k,0)={\widehat {f}}'_{-}(k,0)+{\widehat {f}}'_{+ }(k,0)={\widehat {f}}'(k,0)=C(k)F'(k,0).}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Eliminando rendimientos![{\displaystyle C(k)}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle {\widehat {g\,}}(k)+{\widehat {f}}_{+}(k,0)={\widehat {f}}'_{-}(k,0) /K(k),}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
dónde
![{\displaystyle K(k)={\frac {F'(k,0)}{F(k,0)}}.}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Ahora se puede descomponer en el producto de funciones y que son analíticas en los semiplanos superior e inferior respectivamente.![{\displaystyle K(k)}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle K^{-}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle K^{+}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Para ser precisos, donde![{\displaystyle K(k)=K^{+}(k)K^{-}(k),}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle \log K^{-}={\frac {1}{2\pi i}}\int _{-\infty }^{\infty }{\frac {\log(K(z))} {zk}}\,{\textrm {d}}z,\quad \operatorname {Estoy} k>0,}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle \log K^{+}=-{\frac {1}{2\pi i}}\int _{-\infty }^{\infty }{\frac {\log(K(z)) }{zk}}\,{\textrm {d}}z,\quad \operatorname {Estoy} k<0.}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
(Tenga en cuenta que esto a veces implica escalar para que tienda a ser como .) También descomponemos en la suma de dos funciones y que son analíticas en los semiplanos inferior y superior respectivamente, es decir,![{\displaystyle K}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle 1}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle k\rightarrow \infty }](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle K^{+}{\widehat {g\,}}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle G^{+}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle G^{-}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle K^{+}(k){\widehat {g\,}}(k)=G^{+}(k)+G^{-}(k).}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Esto se puede hacer de la misma manera que factorizamos
. En consecuencia,![{\displaystyle K(k).}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle G^{+}(k)+K_{+}(k){\widehat {f}}_{+}(k,0)={\widehat {f}}'_{-}(k ,0)/K_{-}(k)-G^{-}(k).}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Ahora bien, como el lado izquierdo de la ecuación anterior es analítico en el semiplano inferior, mientras que el lado derecho es analítico en el semiplano superior, la continuación analítica garantiza la existencia de una función completa que coincide con el semiplano izquierdo. o lados derechos en sus respectivos semiplanos. Además, dado que se puede demostrar que las funciones a ambos lados de la ecuación anterior decaen en valores grandes k , una aplicación del teorema de Liouville muestra que toda esta función es idénticamente cero, por lo tanto
![{\displaystyle {\widehat {f}}_{+}(k,0)=-{\frac {G^{+}(k)}{K^{+}(k)}},}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
y entonces
![{\displaystyle C(k)={\frac {K^{+}(k){\widehat {g\,}}(k)-G^{+}(k)}{K^{+}(k )F(k,0)}}.}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Ver también
Notas
Referencias
- "Categoría: Salchicha Salchicha - WikiWaves". wikiwaves.org . Consultado el 19 de mayo de 2020 .
- "Método Wiener-Hopf", Enciclopedia de Matemáticas , EMS Press , 2001 [1994]
- Fornberg, Bengt. Variables complejas y funciones analíticas: una introducción ilustrada. Piret, Cecile. Filadelfia. ISBN 978-1-61197-597-0. OCLC 1124781689.
- Noble, Ben (1958). Métodos basados en la técnica de Wiener-Hopf para la solución de ecuaciones diferenciales parciales . Nueva York, Nueva York: Taylor & Francis Estados Unidos. ISBN 978-0-8284-0332-0.