Función de Green

Función de Green

En matemáticas, una función de Green es una función matemática usada como núcleo de un operador lineal integral y usada en la resolución de ecuaciones diferenciales inhomogéneas con condiciones de contorno especificadas. La función de Green recibe ese nombre por el matemático británico George Green, que desarrolló el concepto hacia 1830.

El término también aparece en física, particularmente en teoría cuántica de campos, para referirse a varios tipos de funciones de correlación y operadores integrales para ciertas magnitudes calculables a partir del operador de campo.

Contenido

Motivación intuitiva

El término función de Green se usa para designar a un operador lineal K que tiene forma de integral, siendo el núcleo de este operador integral la función de Green propiamente dicha. Para explicar que es la función de Green consideremos un operador diferencial lineal L que actúa sobre cierto espacio de funciones definidas sobre una variedad diferenciable M, y pongamos que pretendemos resolver la ecuación diferencial:

(1) L[u(x)] = f(x) \qquad x\in\Omega\subset M

La idea del método basado en la función de Green es encontrar una función de dos variables G(x, s) continua y diferenciable en el sentido de la teoría de distribuciones que cumpla:

(2) L[G(x,s)] = \delta(x-s)\,

Donde \delta()\; es la distribución delta de Dirac. Si se puede hallar una función G que cumpla la ecuación (2) entonces la solución de la ecuación (1) sea cual sea la función f puede escribirse en la forma:

(3) u(x) = K[f(x)] := \int G(x,s)f(s)\quad ds

Puede verse informalmente que la solución así calculada es solución de la ecuación (1) ya que:

L[u(x)] = \int L[G(x,s)] f(s) ds = \int \delta(x-s)f(s) ds = f(x)

Por tanto, tenemos la siguiente relación entre el operador integral dado por la función de Green y el operador diferencial asociado a la ecuación diferencial:

L\circ K = K\circ L = Id  \qquad \Rightarrow  \qquad L^{-1} = K

Conviene añadir algunas precisiones al planteamiento informal que hemos presentado:

  1. Si el núcleo de L no es trivial, entonces la función de Green no es única, aunque en la práctica una combinación de las simetrías del problema, las condiciones de contorno y otros criterios prácticos externos nos proporcionan un única función de Green.
  2. La función de Green G usualmente no es una Función matemática ordinaria sino que puede ser una distribución o función generalizada.
  3. No cualquier operador diferencial lineal L admite función de Green. En el caso más general K es sólo un inverso por la derecha de L.

Las funciones de Green son muy útiles en teoría de la materia condensada donde permiten resolver la ecuación de difusión y también en mecánica cuántica donde la función de Green del hamiltoniano es un concepto clave, para el desarrollo de la teoría cuántica de campos.

Definición formal

Para definir la función de Green que hace de núcleo integral del operador que resuelve cierta ecuación diferencial inhomogénea es necesario introducir algunos conceptos. Empezando con un operador diferencial de Sturm-Liouville \scriptstyle L de la forma:

 L = {d \over dx}\left[ p(x) {d \over dx} \right] + q(x)

Y expresando mediante operador D las condiciones de frontera de Dirichlet:

 Du = \left\{\begin{matrix} \alpha _1 u'(0) + \beta _1 u(0) \\ \alpha _2 u'(l) + \beta _2 u(l). \end{matrix}\right.

Sea \scriptstyle  f(x) una función continua en [0,l], con la cual planteamos el siguiente problema:

 \begin{matrix}Lu = f \\ Du = 0 \end{matrix}

Este es un problema regular, lo cual significa, que para la ecuación homogénea la única solución existente es la solución trivial.

Teorema Solamente existe una solución u (x) que satisface
 \begin{matrix}Lu = f \\ Du = 0 \end{matrix}

Dicha solución viene además dada por la siguiente expresión:

 u(x) = \int_0^\ell f(s) G(x,s) \, ds

En la cual G (x, s) es la función de Green que satisface las siguientes condiciones:

  1. G (x, s) es continua x y s.
  2. Para  x \ne s ,  L G( x, s ) = 0 \,.
  3. Para  s \ne 0, l ,  D G( x, s ) = 0 \,.
  4. Salto en la derivada:  G'(s_{ + 0}, s ) - G'(s_{ - 0}, s ) = 1 / p(s). \,
  5. Simetría: G(x, s) = G(s, x).

Ejemplos

Ejemplo introductorio

Dado el problema

 \begin{cases} \cfrac{d}{dx}\left[ \cfrac{d}{dx} u(x) \right] + u(x) = f(x) \\
u(0) = 0, \quad \quad u\left(\frac{\pi}{2}\right) = 0 \end{cases}

Donde la última línea representa las condiciones de contorno o frontera. Para encontrar la función de Green del problema anterior se siguen los siguienes pasos:

  • Primer paso. La función de Green para el operador lineal es definida como la solución para_

 g'' + g=\delta(x-s).\,

Si \scriptstyle x\ \ne\ s\,\!, entonces, la distribución delta asume una valor nulo y la solución general para el problema es

 g(x,s) = A \cos x  + B \sin x.\,

Para \scriptstyle x\ <\ s, la condición de frontera en \scriptstyle x=0 significa que:

 g(0,s) = c_1 \cos 0  + c_2 \sin 0 = c_1 \cdot 1  + c_2 \cdot 0 = 0, \quad c_1 = 0.

La ecuación para \scriptstyle g(\pi/2,s) = 0 se omite pues \scriptstyle x \ne \pi/2 si \scriptstyle x < s y \scriptstyle s \ne \pi/2. Para \scriptstyle x > s la condición de frontera en \scriptstyle x=\pi/2 implica que:

 g\left(\frac{\pi}{2},s\right) = c_3 \cdot 0  + c_4 \cdot 1 = 0, \quad c_4 = 0.

La ecuación \quad g(0,s) = 0 es omitida por similares razones. Combinando ambos resultados anteriores, obtenemos, finalmente:

 g(x,s)=\left\{\begin{matrix} 
c_2 \sin x, \;\; x < s \\
c_3 \cos x, \;\; s < x \end{matrix}\right.

  • Segundo paso. A continuación, vamos a encontrar c2 y c3. Debemos asegurar la continuidad de la función de green para el intervalo escogido. Cuando x=s\,\! se tiene que:

 c_2 \sin s = c_3 \cos s.\,

También debemos asegurar la discontinuidad de la primera derivada por integración de la ecuación diferencial de \scriptstyle x=s-\epsilon a \scriptstyle x=s+\epsilon y tomando el límite cuando \scriptstyle \epsilon tiende a cero. Por lo cual, derivando la igualdad anterior y garantizando la discontinuidad de esta, tenemos:

 c_3 \cdot [ - \sin s ] - c_2 \cdot \cos s = 1\,

En la cual se iguala a 1 pues \scriptstyle p(x)= 1. Resolvemos para las constantes. \scriptstyle c_2 y \scriptstyle c_3 obteniendo:

 c_2 = - \cos s  \quad ; \quad c_3 = - \sin s

Entonces, la función de Green es:

 g(x,s) = \begin{cases}
-\cos(s) \sin(x) & x < s, \\
-\sin(s) \cos(x) & s < x \end{cases}

  • Solución final, recopilando los resultados anteriores tenemos que la solución final al problema planteado es:

u(x) = \int_0^{\pi/2} f(s)g(x,s)\ ds = 
-\int_0^x f(s)\sin s\cos x\ ds - \int_x^{\pi/2} f(s)\cos s \sin x\ ds

Dicha solución existe para cualquier función \scriptstyle f \in L^1([0,\pi/2]) integrable en el intervalo \scriptstyle [0,\pi/2].

Oscilador armónico forzado

En el caso de un oscilador armónico tenemos la siguiente ecuación diferencial

 mx''(t)+bx'(t)+kx(t)=F(t) \,

Siendo \scriptstyle F(t) la fuerza que provoca la oscilación. Supondremos que la fuerza comienza actuar en t=0 de modo que:

 F(t)=\begin{cases} 0&t<0 \\ f(t) &t \geq 0 \end{cases}

Asumiendo como condiciones de contorno x(0)=x'(0)=0 la solución de la ecuación de movimiento es:

(*)  x(t)=\int_{-\infty}^{\infty}G(t,s)F(s)ds

con

 G(t,s) = \begin{cases} \frac{1}{m\omega_1}e^{-\beta (t-s)}
\sin(\omega_1(t-s)) &t\geq s\\ 0 & t<s\end{cases} \qquad\qquad \omega_1^2=\frac{k}{m}-\beta^2
\qquad \beta=\frac{b}{2m}

Eliminando la parte nula de la función de Green resulta:

 x(t)=\int_{0}^{t}\frac{1}{m\omega_1}e^{-\beta (t-s)}
\sin(\omega_1(t-s))f(s)ds \qquad\qquad\forall t\geq 0

La integral (*) se conoce como integral de Duhamel.

Ecuaciones diferenciales con coeficientes constantes

Una ecuación diferencial lineal con coeficientes constantes de orden n-ésimo se caracteriza por ser de la forma:

 \sum_{k=0}^{n}a_k\frac{d^kx(t)}{dt^k}=F(t)

Supondremos que F(t) es de la forma:

 F(t)=\begin{cases} 0&t<0 \\ f(t) &t \geq 0 \end{cases}

La solución que cumple las condiciones de contorno:

 \frac{d^{n-1}x(t)}{dt^{n-1}}|_{t=0}=\frac{d^{n-2}x(t)}{dt^{n-2}}|_{t=0}=\cdot\cdot\cdot=\frac{dx(t)}{dt}|_{t=0}=x(0)=0

viene dada por:

 x(t)=\int_{-\infty}^{\infty}G(t,s)F(s)ds

con

 G(t,s) = \begin{cases} \frac{dx_h(u)}{du}|_{t-s} &t\geq s\\ 0 & t<s\end{cases}

siendo xh la solución particular de la ecuación homogénea que verifica:

\begin{cases}
\sum_{k=0}^{n}a_k\cfrac{d^kx_{h}(t)}{dt^k}=0\\
\cfrac{d^{n-1}x_h(t)}{dt^{n-1}}|_{t=0}=\cfrac{d^{n-2}x_h(t)}{dt^{n-2}}|_{t=0}=\dots =
\cfrac{dx_h(t)}{dt}|_{t=0}=0\\ x_h(0)=-\cfrac{1}{a_0} \end{cases}

La solución de la ecuación inhogénea viene dada por tanto:

 x(t)=\int_{0}^{t}\frac{dx_h(u)}{du}\bigg|_{t-s}f(s)ds, \qquad \forall t\geq 0

Aplicaciones

El uso principal del formalismo de la función de Green en matemáticas y física es la resolución de ecuaciones diferenciales inhomogéneas con condiciones de contorno dadas. En física las funciones de Green además son usadas como propagadores en el cálculo de diagramas de Feynman.


Bibliografía

  • Eyges, Leonard, The Classical Electromagnetic Field, Dover Publications, New York, 1972. ISBN 0-486-63947-9. (Chapter 5 contains a very readable account of using Green's functions to solve boundary value problems in electrostatics.)
  • A. D. Polyanin and V. F. Zaitsev, Handbook of Exact Solutions for Ordinary Differential Equations (2nd edition), Chapman & Hall/CRC Press, 2003. ISBN 1-58488-297-2
  • A. D. Polyanin, Handbook of Linear Partial Differential Equations for Engineers and Scientists, Chapman & Hall/CRC Press, Boca Raton, 2002. ISBN 1-58488-299-9
  • Jerry B. Marion Dinámica clásica de las partículas y sistemas

Enlaces externos


Wikimedia foundation. 2010.

Игры ⚽ Поможем написать реферат

Mira otros diccionarios:

  • Función racional — Este artículo trata sobre el concepto matemático. Para la «capacidad de razonar», véase Racionalidad. Función racional de grado 2 …   Wikipedia Español

  • Identidades de Green — En matemáticas, las Identidades de Green son un conjunto de igualdades en cálculo vectorial. Nombradas así en honor del matemático George Green, el mismo que descubrió el Teorema de Green. Contenido 1 Primera Identidad de Green 2 Segunda… …   Wikipedia Español

  • Bowling Green (línea de la Avenida Lexington) — Saltar a navegación, búsqueda Bowling Green …   Wikipedia Español

  • Historia de las ecuaciones diferenciales — Una ecuación es una igualdad condicional que se cumple sólo para las soluciones de la misma. Así, en una ecuación algebraica como x – 2 = 0, la igualdad sólo se cumple para x = 2. En forma similar, una ecuación diferencial, constituida por… …   Wikipedia Español

  • Método de elementos de frontera — El método de elementos de frontera, es también conocido como método de elementos de contorno o método de elementos de borde (es mejor conocido como BEM, del inglés boundary element method) es un método numérico para resolver ecuaciones en… …   Wikipedia Español

  • Líquido de Fermi — Un líquido de Fermi[1] es un término genérico para un líquido mecánocuántico de fermiones que surge bajo ciertas condiciones físicas cuando la temperatura es suficientemente baja. La interacción entre las partículas de un sistema de muchos… …   Wikipedia Español

  • Solución fundamental — En matematicas, una solución fundamental para un operador diferencial con derivadas parciales L es una formulación en el lenguaje de la teoría de distribuciones proveniente de la antigua idea de la función de Green. En terminos de la función… …   Wikipedia Español

  • Núcleo de Poisson — En la teoría del potencial, el núcleo de Poisson o kernel de Poisson es un núcleo integral, utilizado para resolver el problema de Dirichlet en dos dimensiones. Específicamente, sirve para hallar las soluciones a la ecuación de Laplace en dos… …   Wikipedia Español

  • Ecuación de Poisson — Saltar a navegación, búsqueda En matemática y física, la ecuación de Poisson es una ecuación en derivadas parciales con una amplia utilidad en electrostática, ingeniería mecánica y física teórica. Su nombre se lo debe al matemático, geómetra y… …   Wikipedia Español

  • Operador diferencial — En matemáticas, un operador diferencial es un operador lineal definido como una función del operador de diferenciación. Ayuda, como una cuestión de notación, considerar a la diferenciación como una operación abstracta, que acepta una función y… …   Wikipedia Español

Compartir el artículo y extractos

Link directo
Do a right-click on the link above
and select “Copy Link”