Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                

Ejercicio 121

Descargar como pdf o txt
Descargar como pdf o txt
Está en la página 1de 2

Capı́tulo 11.

Problemas de contorno para ecuaciones diferenciales 293

Para encontrar la solución al problema de contorno basta resolver la ecuación

u′ (1) = s cos λ = 0,

despejando la variable s. Para cualquier valor de λ siempre existe la solución


trivial u = 0. Pero, si cos λ = 0, s puede ser escogido arbitrariamente. En
otras palabras, si λ = k + 12 π para algún numero entero k, la función
  
1
u(x) = sen k+ πx
2

(ó cualquier múltiplo de ella) es solución del problema de contorno homogéneo


de Neumann. ⋄

La aproximación numérica de un problema de contorno de Neumann sigue


las mismas pautas que en el caso de los problemas de Dirichlet. La única
diferencia está en cómo se aproxima la condición de contorno de Neumann.
El modo más simple consiste en aproximar la derivada de la solución mediante
una diferencia dividida progresiva en el caso del extremo lateral izquierdo o
una retrógrada, en el caso del extremo lateral derecho.
Ejercicio 122 Aproximar las soluciones del problema de contorno

−u′′ = 0, u(0) = 1, u′ (1) = 0

mediante un método de diferencias finitas que use una diferencia dividida


retrógrada para aproximar la condición de contorno en el extremo lateral
derecho. Determinar si la solución del problema aproximado es única.
Solución: Si se utiliza el esquema usual para discretizar la derivada segunda,
se obtiene el siguiente sistema lineal
ui+1 − 2ui + ui−1
− = 0, i = 1, 2, · · · , N − 1, u0 = 1, uN = uN −1 .
h2
Para analizar el problema discreto, se expresa en forma matricial Ax = b
donde
 
2 −1 0 · · · 0 0  
 −1 2 −1 · · · 1
0  0

 0 −1 2 · · ·
  0 
0  0 
..

A =  .. .. .. . . .. ,.. b= .
   
.
 . . . . .  .  
   0 
 0 0 0 ··· 2 −1 
0
0 0 0 ··· −1 1
294 Introducción al Cálculo Numérico

Si se desarrolla el cálculo del determinante de A por la última fila, se com-


prueba que está dado por la fórmula

det (A) = det (AN−2 ) − det (AN−3 ),

donde Aj representa la matriz formada por las j primeras filas y columnas


de la matriz A para j = 1, 2, · · · , N − 2. De acuerdo con el ejercicio 105, se
tiene que
det (A) = N − 1 − (N − 2) = 1 > 0,
lo que prueba que el sistema tiene una única solución. ⋄

11.6 Ejercicios
Ejercicio 123 Determinar si el siguiente problema de contorno

−u′′ + λ2 u = 0, u(0) = u(1) = 1,

tiene solución y en caso afirmativo, encontrarla.

Solución: La solución general de la ecuación diferencial es

u(x) = c1 eλx + c2 e−λx .

La solución del problema de valor inicial u(0) = u0 y u′ (0) = u′0 es

u0 + u′0 /λ λx u0 − u′0 /λ −λx


u(x) = e + e .
2 2
De este modo, las soluciones v y w están dadas por
1
(λ + s)eλx + (λ − s)e−λx ,

v(x; s) =

1 λx 
w(x) = e − e−λx .

Ahora, se busca una solución de la forma
1
(λ + 2s)eλx + (λ − 2s)e−λx

u(x) = v(x; s) + sw(x) =

Para encontrar la solución al problema de contorno basta resolver la ecuación
1
(λ + 2s)eλ + (λ − 2s)e−λ = 1

v(1; s) + sw(1) =

También podría gustarte