= I + ji lx,, Con ello se da por terminado el estudio del esquema (26) (el estudio del esquema (26) se ha ilustrado, de hecho, con Jos tres últimos ejemplos). Todo lo expuesto más arriba sirve de ejemplo ti pico de cómo se realiu el estudio de los esquemas de diferencies.
Cap. IV. Mito®• el• dl/erenc/aa de la rHohicl6n
t72
§ 2. Esquemas de diferencias homogéneos tri puntuales t. Problema de partida. Consideremos el primer problema de contorno para una ecuación düerencial ordinaria de segundo orden: Lu=
!
(k(z): )-q(.r)u "" -/(z),
k(z);;;;;i:q>O,
q(z);;;;i.0, u(0)=µ 1,
0
(1)
u(1)=112·
Una ecuación de este tipo describe la distribución estacionar ia de temperatura, es decir, una distribución que no varía en tiempo (ecuación estacionaria de conductibilidad térmica), o bien la distribución de concentración (ecuación de difusión). Si u = u (z) es la temperatura, entonces W (z) = = - k (z) : es un flujo térm ico (k (z) es el coeficiente de conduct ibilidad térmica). El problema (1) tiene una solución única, si k (z), q (.:r). f (z) son funciones continuas a trozos. Si k (z) tiene una discontinuidad de primera especie en el punto x = ~. de modo que [k] = k (~ O) - k (~ - O) -:p O, en dicho punto deben ser continuos tanto la temperatura u como el flujo térmico -(ku'): [ul = O, lku'] = O para x = ~.
+
Son posibles tambián otras condiciones de contorno para z =O, z = i: ku' = a 1 u - 111 para x = O, -ku' = a 1 u 111 para x = 1. Si a 1 >O, entonces la citada condición es de tercera especie; cuando o 1 = O, tenemos la condición de segunda especie (ku' = - 111 para x = O). Son posibles combinaciones de diferentes condiciones para x = O y x = 1. 2. &quemas de dJfenmclu trfpuntuale& Introduzcamos en el segmento O< t una red uniforme IA>h = {x1 = = lh, t = O, 1, .. ., N} de paso h = t/N y elíjamos un molde tripuntual (z,,, x 1, z 1 + 1 ). en el cual escribiremos el esquema de diferencias que aproxima el problema (1.). Cualquier ecuación en düerencias en este molde tendrá por expresión (2) b1¡¡1 +1 - c1¡¡1 + a1y,_. = - h1q¡1 ,
x<
ª"
donde b,, e, son loa coeficientes dependientes de k (z), q (z) y h. Estos coeficientes son , por ahora, deseonoeidoa. Escribamos (2) de otra forma: t
T
(b 1--,. ,,,.. --,,,--a, -,,,_,,_, -,.-- )
.J
-
u11/1 =
(3)
-cp,,
b1)/h1 •
d1 = (c1 - a 1 -
Diremos que un esquema de diferencias es homlJglneo, si sus coeficientes en todos los nodos de la red para cualesquier a coeficientes de la ecuación diferencial se calculan según unas JllÍsmas fórmulas. Asi, por ejemplo, si introducimos las funcionales A lk (1)), B [k (1)), D (k (s)J, F U-(1)), definid as para cualesquier a funciones continuas a trozos sobre el segmento -1~ s~ 1, y calculamos los coeficientes del esquema (3) por las fórmulas a1 = A lk (z1 sh)I, b 1 = B lk (z1 sh)], d1
=
Dlk (:r1
+ + &h)l,
cp1
=
F [/ (:r1
+ + sh)J, k (s) = k (z 1 + +
sh),
entonces t al esquema será homogéneo. He aqui las funcionales más simples
A (k (s)J - k(-0,5), P (i (s)) = f (0),
a,= k1 -112 -
k (:r1 -0,5h),
cp1 = / 1 = / (z1), et.e.
Si un esquema es homogéneo, resulta más cómodo servirse de l sistema de designaciones sin indices: 1
A, =-¡¡-(byx-aYi) - dy - -cp,
:rECJ>,,,
y (0) ... "" donde a = a (z), b = b (z), y,, = (y (z + h) - y (z))lh,
y
=
!/i
y (z),
= (y (:r)
y (1) = ""'
(4)
z = th E w11. - y (:r - h))/h.
Para que el problema.(4) sea resoluble, es suficiente que sea O, d;;;i. O, y en este caso In solución puede ser determinada por el método do factoriznción (véase el cap. I, § 3).
a> O, b >
Cap. IV. Milodo1 d• dl/erencltu de la rc1oluct6n
t74
3. Condiciones de aproximación. Calculemos el error de aproximación del esquema (4):
'!> = (J\v + cp)-(Lv +
/)
= (J\v- Lv) +(cp- /) =
= [} (bv,.-avi)-(kv')']-(d-q) v+(cp-/), donde v (z) es una función arbitraria suficientemente suave; k, q, j cuentan con un número de derivadas necesarias en el t ranscurso de la exposición. Hegamos uso de la fórmula de Taylor:
v (z
±
h) = v (z)
±
hv' (z)
,.. v• (z) ± ,.. v" (z) +O (h')
+2
6
y bailemos
v.,=v'+{-v·+ ": v"+O(h').
vi=v'-{-v•+ ~·
v" + O(h3).
Sustituyamos estas expresiones para para 'iJ:
'l>=(+(b-a)-k')v'+(
bta -
11.,
y v;¡ en la fórmula
k)v• +
+ h (b;-aJ v· -(d-q) v + (q¡- /)+o (h2). De aquí se ve que el esquema tiene el segundo orden de aproximación si quedan cumplidas las condiciones
b~a = k'(z)+O(h1),
bta
=k(x)+O(h1),
d=q(x)+O(h1),
q> = /(x)+O(h1).
En este caso '!> = O (h1 ) . El esquema (4) con los coeficientes b, = k1+1¡2.
ª' = k1 - 1¡2.
d, = q¡,
b1 = 1c, + 21c1+112+1c,.. , a, = k 1- 112. 4
q>¡ =
d¡- q,
¡,,
t
q>¡= ,.
(5)
f t.
f7~
E•9acm11 de dl/ncntlo• lrlpanluo/cr
satisface las condiciones (5) del segundo orden de aproximación, mientras que el esquema con los coeficientes b, = k1+1•
no satisface ni siquiera la condición del primer orden de nproximación, puest.o quo t T (b1 -a1)-kl=O(t).
§ 3. Esquemas de diferencias conservativos t. &quemas conaervativoe homogéneos. En el § 4, cap . 1 fue establecido que Ja condición necesaria y sufícíento para que un operador de diferencial! Ay sea eutoconjugedo (la matriz sea simétrica) consiste en 11.ue b1 = a 1 +i· En este caso el problema (2) del § 2 adquiere la forma t [ ,,., - ,, 111 - 1/1-1] d 1v1 = -cp, , Av -=-¡-
a,•. --,.---a,--,.-- í=t, 2, ... , N-1,
y 0 = ¡.i1 ,
y,, - µ.,.
(1)
La ecuació n
hd,v, = -
ª1+1
hcp,
(2)
es un análogo reticular de la ecuación de balance del calor sobre el intervalo (z_1,,, z1+ 1,,): "I+ 1/2
W1+ 112 -W1-112 -
J
"l+l/2
qudz= -
"t - 1/2
J /(z)dz,
W=ku',
"l - 1/ 2
(que se obtiene integrando la ecuación (1) del § 2 a lo largo z 1 + 1•1 y lleva el nombré do esdel segmento· z1-i.t z quema conservativo, es decir, esquema para el cual so cumplen los análogos de diferencies de las leyes flsicas do co nservación. El requ isi to b1 = a 1 para un esquema homogéneo significa que 8 lk (z +&Ji) = A lk (z + (s + l )/h) , o bien, B lk (s)I = A lk (s + 1)1 para cualesquiera funcio nes conlin uas a trozos k (s) en el segmento (-1, 1 1. Esto es posi-
< < +\
t76
Cc p. IV. 11110"'1• de 41/•nn
ble eólo cuando la funcional .A [k (•)I no depende de los valores de k (1) para O~ s~ 1, y B lk (1)) no depende de loa valores de k (s) para - 1 ~ s=s;;;; O, de modo que a (z) = = Alk(z + sh)I para - 1 ~8~ O. El coeficiente a (z) del esquema coll!ervativo depende eólo de los valores de k (z) en el segmento (z - h, z). Las condiciones del segundo orden de aproximación (5) del § 2 toman, para el esquema conservativo (2), la forma siguiente •(s+1'~ - •(s) - k ' (z)+O (h1). a(s+1'l+•(s) = k (z)+O(hª),
d (z) = q (z) +o (h•),
cp (z) .... , (z) +o (laJ).
(3)
(4)
De aquí, en particular, proviene que
a (z) = k (z)- 1/,/Jk' (z)+O(h 1 ) =k(z- 1 /Jl) + O(h1 ). Escribamos el esquema conservativo (2) ulllir.ando las designaciones sin indices: (aJl;r).- d (z) JI= -cp (z),
z=thECt.)11 ,
v(0)=)11 ,
y(t) = 1&1·
Exigiremos que se cumplan también las condiciones a;>c1 >0, d;;;i.O .
(5) (6)
En Ja práctica se deben emplear las fórmulas sencillas para a, d y cp, por ejemplo, a1 = k1-1, 1 • d1 = q,. epi = f,. Si la d iscontinuidad de la función k (z) se halla dentro del nodo z = z 1 de la red, calculemos los coeficientes del esquema homogéneo: a, - k1 -111 o bien a 1 "" 1/ 1 (k (z,_1 +O)+ k (z 1 - 0)), d1 - •/1 (q(z1 -0) + q(z1 +0)), cp1 =- '/1 (/(z, -0) + + f (z1 +0)). En este caso las condiciones (3) se cumplen en todo punto. mientras que las condiciones (4) se sustituyen por las condiciones d1 (91-0 + 91-+e) (h 1),
•f,
=o
Demoa a eonoeet loa ejemplos de un eequema cuyos coeficientes ae calculan por integración en los lnt.ervaloa de la
red;
a,-
t {T
"/ " 1-1
+J
) _,
;
'Iª f(z)dz= )
f(z1 +1h)dl,
-112
"l-1(2
d, - {-
d.
i(•1+11\)
- 1
"H-112
q>, ...
ºJ
ds )-' { k(z) =
.l
I
'[ª
"1+1/1
J
g(z}dz=
"1-t/2
g(z1 + 1h)dl.
- 1/1
Es evidente, pues, que las condiciones (3), (4) se cumplen. 2. Enor de aproñllll~l6n. Veamos un esquema conservativo de segundo orden de Ja aproximación. Sea u. = u.(z) la snluclón exacta del problema
Lu
=
(ku.')' - q (z) u
=- f
(z).
u (O)
=
O< z µ1 ,
<
1,
u (t)
=
µ 1,
(6)
y sea y 1 = y (z1) la 110lución del problema do contorno en diferencias (5). Analicemos el error del esquema, os decir, una función reticular
s (z) = y (z) - u (z), Al sustituir lf (z) = s (:t) = u (z) en la ecuación (5) y al suponer que u (z) ea la función dada , obtendremos pllra eJ error s (z) un problema ,'\% = (a~).r -
z (0)
= O,
ch
=-
'!> (z),
a (1) =O,
z E "'11>
a;;;. c1 > O,
d> O,
(6')
donde ;> (z) ... Au + q> (z) = (ou;;). - du + q> es el residuo del esquema (5) en la solución u = u (z) del problema diferencial de partida.
Cop. /Y. M'i<>do• cü dl/n•nclo• de lo r uoh1cl6n
178
Teniendo presente que Lu 1j> = (J\u
+ q>) =
+f
=
O, escri bamos
(Lu + /) = (J\u - fa) + (q> - /) = [(au;;)z - (ku')' 1 - (d - q) u (q> - /).
+
Por hipótesis, el esquema (5) satisface las condiciones del segundo orden de la aproximación. Esto significa que 1j> = = o (h9 ), si k E CI'>, q, I E cm, u E e<•>, y, por lo tanto, ll'l>llc=O(h2 ). Con estas suposiciones el esquema tiene el segundo orden de exactit ud. No obst ante, el mismo orden de exactitud tiene lugar también para las exigencias más débiles tespecto a la suavidad: k (.:t), q (x), f (x) E c1•>, u E CP>. (7) LB.llA.
St se cumplen las condiciones (7). queda lfcita la
/6rmu.la (ku')
-{ku')
1+113"
i - 112
= (ku')í+O(hª),
(8)
dorule u = u (.:t) es la soluci6n de la ecuaci6n (6). DEKOSTRACION. Hagamos uso de la fórmula do Taylor: t
" ' *''2 =111 ± T
lw'+ ,., 1
·+ 48 ,.. 111... (.:t1 :±: 0h),
-¡¡-111
0~ 8~ i ,
+
(111+1¡2-V1 - 1¡2) = 11í+O(hi).
Sustituyendo aquí 11 = ku' y teniendo en cuenta que (ku')·= = (qu - /)', (ku' )" = (qu - /)' , obtenemos la fórmula (8). En virt ud del lema, el error de aproximación 1j> puede ·ser representado en la forma
"11=1'),..1 + 11>:,
'11=(au;;)1-(ku')1 - 112.
1j>j = O(h2 )
bajo las condiciones (7). Ahora, teniendo presente que
a, = k1-11i +O (h1) U.X, 1 =
para
ui-;."1-1 = (u') 1_ 112 -¡- 0 (h1 )
k (x} E 0 1>, para
uEC1»,
u 1 ... ut-t/a +
'lt =O (h1 ). Efectivamente, 1 + 'IJ!uí-112 + /1""-it-t12 +O (hª).
obtenemos
"1-1 = "1-111-{ ltui-111 +{-h•ui-111 +O (h1),
ª'"i., =
"i. f = "Í- 111 +o (hl). (kt-111 +o (h')) M-111 +o (,..))=(ku')•- 112+0 (hl), '11 =0(h1 ) .
Mb abajo se obtendd la estimación aprioristica 11 s lle directam.e nte en términos de '1 y 111•. 3. &timadonee aprtorfaUca Pasemos a la estimación del error s en términos de lj). Recordemos, ante todo, la eatimación obtenida en el § 5 del cap. 1 con ayuda del método de fact.orisaci6o:
de donde se infiere t
11 s lle~ 2t, 11 ;>lle· Mostremos que para la solución del problema (as;;),, -
az =
- ""' z E "'"· 1 (O) a ;;;;i.c 1 > 0, d>O
=
i
(t)
= O,
es válida la estimación
Uslle~.!. ., (t, 11111.
(9)
N
donde se designa (!/, 11) = ~ 11111,h. f -1
Representemos 1 en forma de una suma % = w donde w y v son las soluciones de los problemas (aw:;)x = - µ,,, % E Wh• W (O) = W (1) = O; Av
=
(aui)x - du
=-
dw,
z E w,. ,
u (0)
+ u,
= u (1) = O
t80
d~
Cop. IV. Jlltodo1 .U dl/trenc/41
la rc1olucl
La función ID se hallaré. en la fonna explicita, para esti.m ar v se haré. uso del principio del máximo. De la ecuación
=
(~+µ)" = O, (aw_c)1+1 = J.11 +1 (~, + µ, se deduce que aw;; + µ - const = c0 • Realicemos las transformaciones evidente!: ID1
= ID1-1 +
..
1
(co- 111) /1
= e, ~
0 = 1DN = c0
~
•-1
N
.l!.. h + w,,
·-·
N
~ ~- ~ ~h ·
LJ
LJ
al
N •-1 N ·-·
Co =
..
1
h
~-
..
ºA
'
...
~_ah¡~~.
·-·
•-t
encontramos
De aqui proviene 1
N
... ,
•-1+1
~ ~+a1 ~ ~1 ~
¡w,i ... l-(1-a1) t
.,.- (1 - a) ~~+a ""'>
1
Jt.J
k-t
••
1
N
N
~ ~ ~ ~ ~.
Jt.J •-•+t
••
Ahora nos queda tomar en consideración que y obtenemos
¿,¡ ª" •-s
4 11 ;.. c1
>
_0
(10)
§ 3.E1qaam1JJ1 de dl/arencla• conun1atlrJ01
t8t
Con el fin de estimar v hagamos uso del teorema 4 del § 5 del cap. 1: (U) 11 V lle~ 11 W 11 e· Al reunir las desigualdades (10 y (11), tenemos 11zlle = 11 w + vlfe~211 wlle~2.(1., 11-'fl, '1 es decir, queda demostrada la estimación (9). Volvamos ahora al problema (6'), donde 1j) = T)" + i¡i•. Representemos en la forma J-1
1jl
= µ",
donde µ 1 = TJi + L] h\jl:.
--·
(12)
y hagamos uso de la estimación (9). Entonces, para la solución del problema (6') obtendremos las siguientes estimacio-
nes a prioríst icas:
llzlle~ ~
N
l -1
A-1
A-1
{c1. ITJIJ + ~hl~
h11l%[}, (13)
11 zlle~+-{(f, ITJll+(1, 111l*IJ). 1 Queda probar que tiene lugar la fórmula (12). En efecto,
·-·
al designar p1 = L] hw%. vemosquep1+ 1 - P1 = h1j>f,esdecir, A- 1
"1t = p", 1
Y IP = TJ,.+ p,.= µ,., donde µ 1 = "' + p1• 4. Convergencia y exactitud del esquema de diferencias. Pasemos a estimar la exactitud de un esquema de diferenc ias. Suponiendo que k (.i:), q (.i:), f (z) E C<'>, u (.i:) E C<'>, obtenemos 1'J (z) = O (h1 ), i¡i• = O (h'). Ahora resta por utiliu r la esti mación apriorist ica (t3), la que pod ría ser susti tu ida por una estimación más aproximada 2 11=11e~c.
lle)·
De aquí se desprende que el esquema (5) converge uniformemente con el segund o orden , es decir, 11 z lle= 11 y u ll e~ Mh1 , si se cumplen las condiciones (7).
182
C•p. IV. Mltou1 114 tll/orcncl1u tic le rc1olacl6n
Resulta mb dificil demoat.rar la convergencia del esquema en la clase de coeficientes discontinuos k (%), q (%), f (%). Para simplificar, analicemos un caso en que k (%) tiene la discontinuidad de primera especie en un punto, mientras que q (%) y f (%) son continuas y pertenecen ambas a la clase C<'>. Denotemos con QCA> [a, bl un conjunto de funciones continuas a trozos que esUn definidas en el segmento [a, b] y tienen en (a, b) /e derivadas con\inuas a t.rozos. As{ pues, sea k (%) E Q<1>, q (%), f (%) E e<'> y k (%) tiene discontinuidad de primera especie en el punto ~ del segmento lz. , %n+11, de modo que t = z. 6h, O ~ e~ i. Para z = t se cumplen las condiciones de conjugación u- = u+, (ku')- = (ku') + = w0 , donde
+
V+
=
V (~
+ 0),
V-
=
V
(t - O).
Entonces "' = O (h1 ) para t =F n + t , 11>T = O (h1 ) para todo 1 .,. 1. 2, ... , N - t, 'ln+ 1 = 4,,+ 1u,.,,. - (ku')n+lli· Sustituyendo aquí u,.+. = u (t - 0) (h'). u. = u 8hu~ + o (h'),
hu: + o
m+ m-
u,.," ... (u,.+¡-U,.)/h - eu~+(t-0) u;+o (h) =
ª (k~}- +
(1-0) (~)+ +O(h)=
=w,(:_+ (ku')n+ ll• = (ku')_ (ku')• +l/t = (ku')+
+ O (h) = w0 + O (h) + O(h) = w0 +O (h)
obtenemos ,.,,..1
= w, [ "••• (
1 9 ;.
) + O(h),
para 0 > para 0 <
1
/ 1, 1 / 1,
~ + •;.9 )- 1 +o (Ja),
J
ea decir, 'ln+t = O (t) para cualquier esquema y sólo para un eequema con coeficiente
~. - [+·~ k~) r •1-1
f 4. B19N,.... 1ob,. lcl NtUI M ••lfo,..u
188
i.eeemo11 'ln+i =O (A). En efecto, t
t
. - = -¡-
""+1
es decir,
!
•,.
t
b
•,.u
ü
8
t-8
k(s)+T ( k(s) =-¡:-+-.-+O(A),
~
•
;,,+2 {:_ + t;:e )- t +o (la), y, por conal¡uiente,
,.,,.. 1 =O {h). En el aegundo miembro de la desl¡ualdad (t3) figura la ma¡nltud N
(1,
11111=- ._,, ~ hf11tl +hl11.. +1I· 1.... +1
Con esto queda demostrado el teorema siguiente. TP.ORBXA. En la cla11 tk coeftclenlu dúcontll'UIOI k (z) E E Ql, g (z), / (z) E C<'> cualquter esquema tk diferencias lwmoglneo (5) de 1egundo orden de la aprozimact6n tlenil el primer orden de '3:actltud, mlentru gue el esquema con coeft-
clente a 1 = ; , tiene el iegundo orden de ezactltud.
§ 4. Esquemas homogéneos sobre las redes no uniformes t.F.aquema co119enaU't'o en una red no uniforme. Elij~ mos en el segmento O r ~ t una red arbitraria no uniforme
<
l=-0, t, ... , N, z.-0, zN- t}.
Para o,bteoer un 89q119m& eo1189r't'atl'f'o \ripuntual en la red no uniforme, 85Cribamos una ecnacl6o de balance en el + 111h segmento ~1-t11•
z,
•1+111 ID1+111 -
ID1- 111-
l
«c-111
•.1+111
qudz - -
J
f(z)dz, w-ku'·
•t-1/1
Dicha ecuación ee anota igual tanto para una red uniforme como para una no uniforme. Noe queda apt0ximar lu lnte-
18"
Cap. IV. Mllodo1 tk dl/trencla1 de la re1olucl6n
grales y derivadas que intervienen en la ecuación de balance: h¡ = I¡ '"I+ 1/2
J / (z) dx -
It-1'
'"1+ 1/2
J
cp11i,.
qu dz - d 1 u11i 1 ,
Zj-1/2
rlonde d 1 y cp1 son unas funciones reticulares. Como resultado, se obtiene un esquema de diferencias 1 [ l/1+ 1-111 111-111-1] d1!11 = -cp,. T, ª1+1~ - a,---- 11 1
i = 1, 2, .. ., N - 1,
Yo=l11•
llN=i-Li·
(1)
Para determinar d1 y cp1 usaremos las fórmulas m'8 sencillas cp, = / 1 , d 1 = q 1, i = i, 2, ... , N - 1.. El coeficiente se determina por los valores k(z) en el intervalo (z1_ 1, z 1), a consecuencia de lo cual puede tomarse igual al que figura sobre la red uniforme, de modo que a 1 = k 1_ 113 + O(hl) para k (z) E C<2>. 2. Error de aproximación. Introduzcamos las designaciones
a,.
JI; = 111-111-1 Z, 1 /i¡ '
l/:r, ¡ = 111.,-111' /il+I
11 · :r. 1
=
lli.1-111 /o¡
y escribamos el esquema de diferencias en la forma
(a11-) • -
. ".
d11
Al suponer
J
= -
= 11 -
cp, :r
=
Z¡
E
m
h•
l/o
=
111• 11 H
=
111· (t)
u, obtendremos para s la ecuación
(as¡¡); - ds = -;>,
x E~'"
'• = 'N =O,
(2)
donde (3)
es el residuo para el esquema (1) en la solución u= u (z).
§ l . E1q1,.m41 1ol>re lu reiü1 "" an l/ormu
185
r>.
LEMA t. Sl IJU E e~>, f E e enloncu para el error de aprozimaci6n \j> es válida la f6rmula
"' = '1~ + "'.
(4)
t
donde 'lr = (au.s}1 - (ku')i-111 - ht (qu - nf/8, i!>t =0 (h1) para cp, = f 1, d 1 = q,. Hagamos uso de la indeo tidad del p. t. escribiéndola en la la forma "1+1/2
O= w. _ J_ 11,
"· 1
r
J
(qu-f)dr , w1 =(ku')1-11z·
'"l-1/Z
Sustraye1nos esta identidad de la igualdad (3): '"1+1/ 1
1j>=((aui) 1 - (ku')1 - 112J - (du) 1 +cp1
+
!, J
(du-f)dr.
"•-111
(5)
La integral que figura en el segundo miembro se representará en forma de una suma de dos integrales: de r 1_, 1 , a z 1 y de z 1 a x 1+ 1¡,; al desarrollar después la función su bintegrsl / = qu - fea el entorno del nodo z = z., hallaremos
+,
"1+1/2
J
~' {
f(x) dx =
"1
J
1/,+(z-z,)f;Jdz+O(h~H-
"1- 1/2
"1-112
"1+112
+
J
(/,+(z-z1)fíJdx+O(h:.,>} =
=T, + s!, (h;., -hn tí +o (lll). h: + hl+ 1 < (21i1)' . La sustitución hl+ ,¡¡ =
puesto que =hl+1fl+1 + O(hl +1} nos da "1+112 I
A;
f
J
•1- 112
f(z)dz = Tr+(hZf) .. +0(1il).
"· 1
C•p. IV. Jlltodo1 .U tlJ/
18&
7
Sustituyendo esta expresión con = qu - f en (5), llegamos a la fórmula (4). Para estimar 11 1 aegún el orden veamos la diferencia (au;¡)1 - (ku'),_111 a condición de que k E e<•>, u E C<'>. Empleando la suposición a1 = k 1_11 1 O (ht) y las fórmulas u, = h1ui..vJ2 + htui...11 J8 O (h~). u1_1 = 111 = "1-111 - h1u"t-11J2 h1ut-11J8 O (h~). u;;. 1 = (u, "1-1)/h, = "1-111 + O (h1), obtenemos (au;¡)1 - (ku')1-111 =
+
u,_ +
= (k1_1, 1
+
+
+ O (hl))(ui..111 + O (h})) -
+
(ku'),_1, 1
= O (ht).
De este modo, es válida la estimación 111 = O (hl) para (le (z), q (z), J (%) E c< 1 >, u (z)
E cea>.
Se suponía que d 1 y qi1 ae determinan según las fórmulas más sencillu: d 1 = q,. cp 1 = (1• Si, en cambio, so emplean las fórmulas más complejas, por ejemplo OBSl!RVACION.
"l+l/J
epi
= ""•- 112+"1+1!1+112 2i1 •
t q>, ~ "ij"
Jr ¡ (z) clz,
"t-t/J
;ir
+
entonces la función reticular = O (h1) - (d1 - q 1) u 1 (cp 1 - f 1) puede ser representada en la forma ;i~ = = p¡, iW, donde ;it• = O (l.t). p 1 = O (111) y 11 .. en la fórmula (4) se sustituye por la suma 111 + p 1: ,¡> = ÍP1 + '11); ,¡>•. ' (4') P1 (ht) '11 (ht), ,¡>f9 (llT} para le, q, /E C'''. u EC'". 3. &umad6n de la nloddad de conYer¡encla. Para el problema (2)- (4) es válida la aiguiente estimación aprioristica
+
+
=o
=o
=o
+
ll•llo~-?-l ((t, l '11 J+(1, 1;i•1) }, donde
(v.
N
v)I ~ ~
t-t
v1v 1h 1•
(7) del§ 3, entonees
(6)
SI ee cumplen tu oondíclonea '11 -0(.\1), 'lj>f ... O(~. 1jif en (6), nos cooveoeemoa de que es
Al sustituir 11 1 y varidico el aiguieñte teorema.
e.
T.ORBllA. En la clase rk c~ficicnU• ·~ le, q, f Cfl> eO
exactitud en una sucesi6n arbitrarla. rk la.8 rerks no uniformes. Al tomar en consideración la observación del p. 2, podemos representar '1>1 en la forma 'l>t = p ;,,, + 'l>t•, donde p 1 = O (h1), 'l>iº = O (Jlt). Entonces, en lugar de (6) queda vUida la flstimación
JlsJlc<;-?{(t, 11l+PIJ+(1, J,p••IJ); l
el teorema sobre el segundo orden de exactitud sobre una red no uniforme queda en vigor. Si el coeficiente k (:r) tiene discontinuidad03 de primera especie en un número finito de puntos, siempre podemos elegir tal red no uniforme ~h (k) que los puntos de discontinuidad sean los nodos de dfoha red. En tal caso cualquier esquema tendrá el segundo orden de exactitud. Asi pues, cualquier esquema homogéneo de segundo orden de aproximación (11> = O (h')) sobre una red no uniforme y en la clase de coeficientes suaves tiene el segundo orden de exactitud con la elección especial de las redes no uniformes (le) en la clase de coeficientes discontinuos. 4. &quema e.neto. Para el problema (1) del § 2 podemos construir un esquema tri puntual exacto cuya solución en los nodos de una red arbitraria coincide con la solución exacta u = u (x) del problema de contorno para una ecuación di ferencial. Ilustremos la posibilidad de construir el esquema exacto co n un caso particular del problema para q (%) • O: (ku')' = - f (%), O < % < 1, u (O) = O, u (1) = O. (7) Al haber integrado la ecuación desde % 1 hasta :r, obtendremos una ecuación
w,.
(ku')-(ku'),
+ l" f a>~=O. "1
Dividámosla por le (r) y integremos respecto de x desde x 1 hasta r 1+ 1:
x,., •. t >+ ) . ."'(~.> ) tm~=O, .., ..,
lt'+t
u«+1-u,-(ku'), )
...
(8)
188
Cap. IV. Jlllo4'1• .U dlfertncla• dt l
(9) Jotroduttamos un.a designación "1
af = [
r
t
T," J
t1z ] - ' k(s) .
"1 -1
Multipliquemos (8) por al+ 1 /h 1 +1 , (9) por ollh,, y sustrayamos del primer resultado el segundo. Obtendremos una ecueción t ( 0 u1..-u1 l .. ,-u,_,]+ O
¡¡;-
ª'••~-a--,.-,-
ci>1= '
o bien (tO)
donde
+
8h1 para z 1_ 1 ~ z' ~ z,, Y z' = Si ponemos z' = z 1 = z 1 + 8h1+ 1 para z 1 ~ z' ~ z¡+1 , entonces dicha fórmula puede reescribirse de la manera siguiente: 1'¡af o o cp, .. -¡;"(•1~íl1> I (z, + M,) dA +
J
-1
J +~ lo 1cc,,!',¡, )lo f(z,+M1+1)d>.. •
•
•
1
•
+l
De este modo, el esquema (10) es exacta sobre una red arbitraria no uniforme y para cualesquiera funciones continuas a trozos k (z) y f (z). Por supuesto, el empleo práctico de este esquema está obstaculizado por el hecho de que los coaficientesdedicho.eaquema se expresan a través de las in-
§ 4. E•qaettuu •obre la• red.. no anl/orMu
t89
legrales de k (z) y/ (z), razón por la cual au c!lculo requiere In aplicación de las fórmulas de integración numárica. 5. Aumento del orden de ex.actitud. De lo dicho ae hace claro que para aumentar la exactitud de la solución aproximada se debe o bien disminuir el paso de la red h o bien aumenlar el orden de exact itud del esquema. No obstante, es eonveniente construir esquemas eon orden de exactitud aumentado sólo para las ecuaciones de coeficientes eonstant.es, puesto que la anotación de tales esquemas para las ecuaciones de eoeficientes va.r iables está relacionada con grandes dificultades tl!cnicas y conduce, a menudo, a los algoritmos engarrosos. Ya hemos aducido un ejemplo del esq uema O (h') para la ecuación u• = - f (z). Examinemos ahora una ecuación u: - qu
=-
q = const >O.
/(z),
Escribamos un esquema de diferencias sobre la red uniforme: /\.y
= Y;,. -
dy
=-
cp (z)
y elijamos d y cp de un modo tal que tenga la aproximación O (h'). El error de la aproximación es ,¡i - /\u + cp = (/\u -u") - (d - q) u+ cp - / = h' =u uIV - (d-q) u+
Al sustitu ir aquí uIV obtendremos
q/ -
r.
q> = - (
d-q-
~; qi) u+
por consiguiente, 'I> = O (h 4), si se pone d = q+ ~; q2 , cp =
= / +12 "' (q/ + /") .
El orden de exactitud queda intacto, en la fórmula para