Metodos Numericos Mir PDF
Metodos Numericos Mir PDF
Metodos Numericos Mir PDF
ASamarski ~
Introduccin
a los mtodos numricos
Introduccin
a los mtodos numricos
Editorial Mir
Mosc6
Traduoldo del ruso poc el ioptuo K. p . Medito
lmJlteeO eu la URSS
fk)'ll
' ' Jl3.A&Tea1>0T10
:J1!I isWTepuypu, i982
.raaaau pe,p11~R t a1110-wue1.1uaec-
t~ t.nduoc16o al Mpallol, edltorlaHMlr. 1986
Contenido
Prlogo 1
Jotroduc:ci6D 9
Captulo J
Ecuaciones en diferencias
t f. Fuoclooe. retlcuJa - . . . . . . . . . . . . . . 28
t 2. Ecuolooes en dllereoclu . . . . . . . . . St
t 3. Reaolucl611 de la11 problelllU de c:ootomo en dUennciu
Pn Ju ecuaciones de segu ndo orden . . . . . . . 40
4. Ecuaciones ea diferencias como ecuaciones operaclonalea e
5. Principio del mWmo para las ecuaclooe1 en dlfereoelu 65
Capitulo 11
lnterpolari6n e integuci6n numrica
t f. loterpolac16o y prodmaclo de 185 luncionl!I 72
2. 1otegraclo o um6rlca . . . . . . . . . . 82
Capitulo I lI
Resolucin numrica de los sistemas
de ecuaciones algebraicas lineales
t t. Si~temas de ecuaciones algebrelcas lineales .
. fOO
t 2. Metodos directos . . . . . . . . . . . . . . toe
J 3. M6t.od os l\erUvos . . . . . . . . . . . .. . . . tt3
4. Esque1na ituativo de dos capu con puimet.ros
de Ch6bl
.~.., . . . . . . . . . . . . . . . . . . . . . . . 129
S. M6todo Jtemado triangula r . . . . . 1"3
6. M6todos i\era\hos de t ipo variaciooaJ . 147
7. R~ lucin de las ecuaciones no linea les 150
6 Pr6lo10
Capitulo IV
Mtodos de diferencias de la resolucin
de los problemas de contorno para ecuaciones
diferenciales ordinarias
t l. Coru:eptos fundamenlalee de la teoria de esquemas de di-
fere11eias . . . . . . . . . . . . . . . . . . 158
t 2. Eaqllemaa de diJerenclas homogneos tripuntuales . . . t72
t 3. Eaqllemaa de dlfere.nclas conaervativos . . . . . . . 175
f 4. Eaquemu homogneos 110bre lu redes no un.ifotmes . . . t83
t 5. M6t4dos de construccin de los esquemas de diferencia t 91
Capitulo V
Problema de Cauchy para las ecuaciones
diferenciales ordinarias
t t. Mtodos de Runge-Kut t.a . . . . . . . . . . . . . . 200
t 2. Esquemas de varios paaoe. Mtodos de Adamll . . . . . 212
f 3. Aproximacin del problema de Cancb.y pra un 1C.tema de
ecnaciones dilerencial linealee ordlnar1u de primer orden 224
f 4. Estabilidad del esquema de dos u.pu . . . . . . 230
Captulo VI
Mtodos de dHerencias para las ecuaciones
elipticas
f t. Eeq11emas de dUnenclu para la ecuacin de Polaeon 241
f 2. Resolucin dll las ecuaciones en diferencias . . . 252
Captulo VJI
Mtodos de diferencias para resolver
la ecuacin de conductibilidad trmica
f t. Ecnaci6n de conductibilidad thmlca coa co.fiefeates
cona1.antea . . . . . . . . . . . . . . . . . . . . 364
t 2. Problemas mullidimenalonales de conduclibilldad trmica 277
1. 3. Esquemas econ111 lcos . . . . . . . . . . . . . . . . 285
Anexo 295
Blbllografla 302
Lta de daignclonl'S 304
lodlu alfabtico 306
Prlogo
91 = 10"''
q, ... 1()1'12, q 1 ==103"''
En este caso obtendremos en el primer paso
s, = </191 = f Q-6P/4 < M'
es decir, S 1 es un cero de mquina; todos los productos
sucesivos S 3 , S,, S& son tambin nulos; de este modo aqu
ocurre una prdida total de exactitud.
3or Mltrooo. Mezclemos estos nmeros suponiendo q1 =
tQ -="I', q 1 = 101'/I, IJs = 1Q3P/, IJc = 10-P/S, </ oc 10"1'.
Entonces hallaremos sucesivamente:
S t = </191 = iQ -P/', s, = S,q, = 10P/1,
S , = S 3 q, = 10, s.= S ,q. = 10"1',
1- ~ / (z) lb.
1
Suaihu)'&lllos f por 7y veamos Y- Jo f (z) dz y la di.leren-
1
""
r, =f(t, u), t>O, u(O) = U, (8)
Ecuaciones en diferencias
1. Funciones reticulares
t. Funclonee retlcularee y operaciones sobre ellas. Y a
se ha mencionado que en los mtodos aproximados las
funciones de un argumento continuo se sustituyen habitual-
mente por las de argumento discreto, esto es, por las funcio-
nes reticulares. La fun,ct,611 rettcular puede, pues, considerarse
como una funcin cuyo argumento es un nmero entero:
!/ (i) = "" t = o, 1, 2, ...
Podemos introducir para (i) las operaciones que represen-
tan un anlogo discreto (de diferencias) de las operaciones
de difereoclacl6n e integnci6n.
El aoilogo de la primera derivada lo constituyen las
diferencias de prtmer onkn:
!1111 = !11 +1 - ,, la dtfer.ncta ckrecha;
Vy, = V1 - llt-t la dtfermcto tsquterda;
-
N-1
H-1
N
~ y1ll111 + ~ 111Vl11=
-
-= 1) . (111111)-1)
e- o
N-1
- o
N
11u-1VY1+1
N
N
+ c-i
1) 111VV1 -
v,tw,- - ~
N- 1
,,_,
N N -1
2. Calclese SN =- ~ 1(1- 1) = ~
,_, ,_, t(I + 1). Ponemos
11 1 = t, V111 = l + 1. En eate caso 111+1 = 111 + (1 + 1) =
= 111 + (2 + 3 + . . . + (1 + 1)) = (11 1 - 1) + (1 + 1) X
X (1 + 2)/2, 111 = 111 - 1 + i (t + 1)/2. Elijamos 111 de la
condicin 11N =O, es decir, 111 = 1 - N (N + 1)/2. Apli-
eando la frmula (3) y teniendo en cuenta que !le = O,
llN = O, Vy 1 = 1, encontramos
N-1 N-t N
N-1
-o 1-1
N-1
2. Ecuaciones en diferencias
f. Ecuaciones en dllerencl.u. Una ecuacin lineal res-
pecto de la funcin reticular y 1 = y (t) (f =O, 1, 2, ...)
a 0 (1) y (1) + a (t) 11 (l + 1) + . . . + "' (1) 11 (1 + m) =
- / (1), (1)
donde a,.
(t) (k = O, 1, ... , m), f (1) son las funciones
reticulares prefijadas, a 0 ( -t) .p O, "' (1) .,:.O, lleva el
nombre de ecuact6n en dlferenct(J.$ lineal de m.-simo orden.
Dicha ecuacin contiene m + 1 valores de la funcin y (i).
Haciendo uso de las frmulas para las diferencias &,.
&,, ... , 6"'-11, podemos expresar los va lores 111+1,
!/1+1 , Ym+1 en lrminos de y 1 y las diferencias citadas:
111+1 = !11 + 6,, !11+1 = 6 1111 + 2111+1 - !11 = 6'111 +
32 Cap. J. BcU4clone1 en dl/ercnc1a1
91
e+ y'Ci=44b
2b 9
,2
= e-~
2b . (14)
Segn sean los valores del discriminante D = e - 4ab,
son posibles tres casos:
1) D = e - 4ab >O. Las ralees 91 y 9 1 son reales
y distintas. Les corresponden las soluciones particulares
l/~1) = q~ //~Z) = la;
estas soluciones son linealmente independientes, puesto que
es distinto de cero el determinante:
9' 9'+1 1
t1.i. .1t+1= ~ :~+I = ~~ (92-91),PQ
1
Ha de notarse que 91 =I= O y q 1 =I= O, pues en el caso contra-
rio a = O y la ecuacin (12) no sera ecuacin en diferencias
de segundo orden. La solucin general de la ecuacin (12)
tiene por expresin
(15)
ca,. l. Bcaac>MI In dl/renc/.41
'11 =
+l v'fDi
2b i 'la =
-1 v'fDi
2b
'
puesto que bt.- a= b ---a= D u= O. Como 6A. A+ =
Y11 = e, + e"'.
2. Hllese la solucin de la ecuacin
Y11+1 - Y1a+1 = 2y,. = O.
El discr iminante es igual a D = 1 + 8 = 9, las ralees
sern 9 1 1 = (1 3)/2, 91 = 2, 91 = -1. La solucin
general es de la forma
Y11 = C12" + C1 (-1)".
3. Hllese la solucin general de la ecuacin
Y11+1 - Y11 - 6y,._1 = 2a.1 (11:!)
La solucin general de una eruacin no homog11oa os la
~umu y,. u:
= y11 + de la solucin general i'A do la 11cuucin
homognea y de la solucin particular v:
de la ecuacin
no homognea. Hallemos primero la solucin general de le
ecuacin homognea. El discriminante es D = 1 + 24 =
= 25 > O, y las races de Ja ecuacin cuadrtica q2 - q -
- 6 = O son <JJ = 3, q1 :== -2, de suerte que y',." = 34 ,
yj." = (-2). La solucin particular u: buscaremos en la
forma y; = c2.4 , donde e = const. Sustituyendo y = c:2."'
en (18), obtendremos e (2.4 1 - 2" - 6 24 - 1) =e .zH
(-4) ~
= 24 1. e = - 1.
La solucin general de la ecuacin (18) tiene por expre-
sin
Y11. =c1 3"' +es (-2)"'-2".
6. Ecuacin en diferencias de eegundo orden con coefi-
cientes variables. Problema de Caueby y problema de contorno.
Examinemos ahora una ecuacin en diferencias con coefi-
cientes variables
b1Yt+l - C1Y1 + 1Yt-1 = /,,
a 1 *0, b1 .,PO, '=0,1 , 2, ... (19)
Dado que b 1
recurrente:
* O, de (19) obtenemos la siguiente relacin
1111-01111-1 + !1
V1+,= b1 , b1* (20)
Expresemos Y1+i e y 1_1 en t rminos de y 1 y las dilerencias
de primero y segundo rdenes. La ecuacin (19) se anotar
en este caso en la forma
!J.Vy1 +(b1 - a 1) 6y 1 - (c1 - a 1 - b1) y 1 = / 1 ,
4 * b, *
La solucin de una ecuacin en diferencias de primer
ocden depende de una constante arbitra.r ia y se determina
univoca.mente, si est prefijada uni. condicin complemen-
taria, por ejemplo, Yo = c0 La aoluci6n de Is ecuacin
de segundo orden se determina por dos constantes a.rbitrarias
y se puede hallarla, siempre que vienen dadas dos condicio-
nes complementarias. Si ambas condiciones estn dadas en
dos puntos vecinos, se trata de un problema de Cauchy.
SI J11s dos condiciones ostn dadas en dos puntos diferentes
(pero 110 vecinos), obtenemos un problema de contorno.
De mayor inters para nosotros sern precisamente . los
1 J. lfcuc"'IU' .,. ,,,,,,.11ellu
problemas de contorno. l ntroduscamos las deslgnaoiooee
Ly, = b1Y1+1 - c,y, + a1Y1-1
y enunciemos los problemas mencionados m'8 de\allada-
mente.
PROBLEMA DE CAtlCBY: hllese la solucin de la ecuacin
Ly1 = fi. t = 1, 2, (21)
con las condiciones complementarias
Yo = 1, Y1 = ,. (22)
La segunda condicin de (22) puede notarse de otro modo:
.Yo = Yt - Yo = -
1 = ji',_; podemos, entonoes, decir
que eo el caso de probleme de Cauchy vieoeo dadas en
un punto t = O las magnitudes
Yo = i, Ayo = 1 (22')
PROBLEMA DE coNTORNo:h&llese la solucin de la ecuacin
Ly1 = /,, t = 1, 2, ... , N - 1
para las condiciones complementarias
Yo = 1' y N = ,, N ;;;;. 2. (23)
En los nodos de frontera t = O e 1 = N pueden defi-
nirse no slo los valores de las funciones, sino tambin sus
diferencias y combinaciones, es decir, las expresiones
et, Ay0 + ~ 1 Yo para i = O, y et 1 Vy N + ~sY N para i = N.
Tales condiciones se pueden anotar en la forma
Yo = X1Y1 + t> Y N = x.y N-1 + , . (24)
Si x 1 = x 1 = O, obtenemos de aqu las condiciones de primer
gnero; cuando x 1 = 1, x 1 = 1, tenemos condiciones de
segundo gnero
Ayo = - 1 , VYN = ,. (25)
Si ><1 , 1 =I= O; 1, (24) llevan el nombre de condiciones de tercer
gnero:
->< 1 l:J.y 0 + (t - >< 1) !lo = 1,
x,VyN + (1 - x,) YN = ,. (26)
Adoms, pueden encontrarse problemas de contorno con
ciertas combinaciones de las citadas condiciones de con-
torno: las condiciones de un tipo para t = O, y condiciones
de otro tipo, para l = N.
La solucin del problema de Cauchy so halla directa-
mente de la ecuacin (21) seg6n la frmula recurrente (20),
tomando en consideracin los datos iniciales y 0 = .,
111 = 14,. Para la resolucin de los problemas de contorno
se emplea un m6todo mh complejo (m6todo de eliminaciones)
el cual se expondr en adelante.
Para nna ecuacin con coeficientes conatantes la solu-
cin del problema de contorno puede expresarse en la forma
explicita.
IUlDD'LO. Hllese la solucin del problema de contorno
1 = t, 2, .... N - 1, (1)
rep?919D\a el sistema de eeuaclones alpbraJeas lineal oon
matrfa trldiagonal de dimenal6n (N +t) X (N t): +
1 - x, o ... o o o ... o o o
a1 -c1 b1 O O O O O O
t-) a1P1+/1
~l+ . -
-
-- -
e1-a1a1 '
l=i, 2, ... , N - t, ~. = 111; (9)
<-)
111 = 1+1111+1 + P~ t= N -1, N-2, . .. , 2, 1, O,
11.+x.P" (tO)
VN t-aN"s
o bien 1258
1 + lt+l>f'll+l
11m e1-b1~t+l lli-t 1-1>1(1+1
Al cotejar con la frmula y 1 = ~ 1 11 1 _ 1 + ~ 1 obLendremos
(13) y (14). El vaJor de Yo hallamos de la coudicin llo =
= x 1 y1 +
A1 y de la frmula Yo = ~1 y 1 'li De la desi- +
g11Aldad 1C1 - b1i1+ 1 1;> 1c1 1- 1b, 1 1i1+1 1;;> l 111 1 +
+ 1b, 1(t - 1i1+1 1), 11 - i1X1 1>1 - 1i1 1l "11116 ve
que las condiciones (11) garantiaan que les frmulas de
far.terlzacin izquierda sean aplicables y su lculo sea
estable, puesto que 1i1 1~1 (t = t, 2, ... , N).
La combinacin de las factoriuciones izquierda y dere-
cha da el mtodo <U facto rlsacloM1 opru.W. EmplndOM
este m6todo, en la regin O ~ f ~ t 0 + 1 se ealculan, aeg6n
las frmulas (8), (9), loa coeficientes de factorlucln '"
p,, y en la regin t1 ~ t ~ N se hallan, por laa f6rmulu
(13), (14), i1 y 111 Cuando t = t,, 118 realiaa la conjugacin
de soluciones en la forma (10) y (15).
De lu frmulas 111, - cz1,+11/1,+1 + ~1.+11/1,+1- ~1.+1Pt, +
+ 111,+1 hallamos
111,+1 +a1,+1'11,
111 = t -ac.+1l1,+1
La cl\ada frmula tiene senUdo, puesto que por lo menos
una de las magnitudes 1~ 1 , + 1 1 6 1cz1,+1 1 es, en virtud
de (U), inferior a la unidad y, por lo tanto, 1 - cz1,+1E1,+1>
>O. Al conocer y 1 , podemos hallar, airvi6ndonos ae la
frmula (10), todos los valores de y1 para t < t 0 , y, por la
frmula (15), todos los valores de y1 para t > f 0 Cuando
t > 10 e t < t 0 , los c6.lculos son autnomos (se llevan a cabo
paralelamente). El m6todo de factorizaciones opuesta.s es
particularmente cmodo, si, por ejemplo, se pide hallar
y 1 slo en un nodo t = t 0
4) si x *
3) ()...x, y) = A. (z, y) para cualquier nmero complejo A.;
O, entonces (z, z) > O.
Hemos de observar que el producto escalar introducido
(x, y) engendra en H la norma
11 xn=V ex. x). (2)
Resulta vlida aqui la desigualdad deCauchy-Buniakov-
ski
1(z, 11) 1' ~ (x, .x) (y, y), (3)
la cual puede escribirse, tomand o en consideracin (2),
en la forma
1 (.x, y) 1 ~ 11 x 11-11 Y 11
2. Operadores lineales en un espacio de dimeualn finita.
Sea H un espacio lineal de dimensin finita provisto de
producto escalar (z, y). Designemos con D cierto subespacio
de H. Si a todo vector x E D se le ha puesto en corresponden-
cia, de acuerdo con una regla determinada, el vector 11 = Ax
de fJ , suele decirse que en H est dado el operador A. El
conjunto D e: U se llama dominio de definicin del opera-
Cap. l. BcaoclonH en dl/erenc/JJ1
Z = ,._,
~ ~~~ (14)
y escribamos
N
0=Az-f = l: (A.-c.- - /1t)t.-.
~-1
Al multiplicar esta igualdad escalarmente por t" y teniendo
en cuenta que (t"' ;".) = Bu. hallemos 0 = A."c" - /,.,
es decir, e" = hl'J..".
6) La norma de un operador autocoojugado A es igual
al mdulo de su valor propio mximo:
11A 11 = mx IA.-1 = IA.NI (15)
1..;A.,;N
diferencias A:
Ay, = l>1Y1+1 - c,y, +
1Y1-1 =
= I>, 6. (VY1) -(I>, - a1) (Vy1) - (c 1 - a1 - 1>1) y,,
1 = 1, 2, ... , N - 1,
es decir , Ay1 E QN-i si y, E QN+i o bien, en la notacin
reducida, A: nN+ l-+ QN-1'
Analicemos un problema de contorno en diferencias
Ay1 = - /1 , i = 1, 2, .. . , N - 1,
Yo = 111 YN = 111 (28)
y escribmosla en la forma matricial:
AY = cI>, (29)
donde; cI> = U1 +
1111 f ,, , f N- f N-1 + I> N-1111) es el
vector conocido e Y = (y1 , y 1 , . ., y 1!~ y N-i> es un
vector desconocido, ambos de dimensin N - 1; A es una
matriz t rid iagonal cuadrada de dimensin (N - 1) X
X (N - 1):
(y, V) = ~ ~ flV
1-1
y una norma
llYll = V(y, y).
Si se estudian el segundo (x1 = x. = 1) o el tercero
(x, =fo O, x 1 =fo O) de los problema.a de contorno (vase (1)
del 3), la matriz A ser cuadrada de dimensin (N +1) X
x (N + 1) y el operador A se definir de la manera siguien-
te:
Ay, = - Ay, = -(b1Y1+1 - Cy + 1Y1-1l.
t = 1, 2, .. . , N - 1,
Ay0 = - (x1y1 - Yo), AyN = - (YN - X1YN-1)
En este caso el operador A aplica el espacio de funciones
reticulares H = QN +i en s mismo: A: H - !l .
En adelante se analizar el primer problema de contorno
para una ecuacin en dcrencias de segundo orden; en
este caso, segn lo mostrado ms arriba, H = Q,,_ 1
6. Frmulas de Green de diferencias. Examinemos un ope-
rador ele diferencias D:
Ly, = b,y, + 1 - C!f + 1!11-1 l = 1, .. . N - 1. (30)
58 Cop. /. BcuocloM tlt/rncla1
1-t
~ '' (oy;),..1 h= - 1-1
~ ''i1Yi1 h+
+ (0Yis),., - 1o (oy) 1. (38')
Al sustraer (38') de (38) obtenemos la segundo /6rmula de
Green de diferencias
l't - 1 N- 1
~ y 1 (az;)r, 1 h = ~ z 1 (ayz),., 1 h+aN(Yzi-sY;;)N-
1-1 1-1
- (Yo (ozi) 1- 10 (ay).). (39)
Si estn cumplidas las condiciones
Yo= Zo =O, (40)
es decir, si y=;,
z = ; EON+I entonces en el segundo
miembro de la igualdad (39) dos ltimos sumandos se anulan
y
N- 1
~ . .
y1 (oz;),., 1 h =
,, . .
N-1
L.J
,_, -
kl z 1 (oyi),., 1 h.
suma
N-1
~ d1
1- t
y,; h,
1 ohtenemos la segundo /6rmulo de Green para
eo
N-1 N-1
~ Y1~1h = ~ ;,Ay,h (42)
1-1 1-1
crihir:
(~ Ly, ~) - (y, :. L~) =
= (Ay. ~) - (y,/\~)+(:, LiY. ~) - (y, -~Li~) =
N
=h :, (b1-a1+ 1>CV1+1V1 - y,u,. 1) h.
1-1
De aqu se vo que (Ly, ~) = (y, Lu), es decir , L = L*
s lo a la condicin de que
N- 1
~ (b1-a1+1)(Y1+1U1-ll1U1+1)h = O. (49)
1-1
Por sor arhitrari11s y 1 y u,, podemos toma r y 1 = <'l1.1,+a
v, = 61.1,. donde i 0 es un nodo fijo cualquiera (i 0 = 1, 2, . ..
. . .. N - 1) . mientras que <'l1 1, es el slmbolo de Kronecker.
62 Ct1>p. J. Bcat1>clone en dl/eronc'41
ll-1
(N- t )A t
11 sen n171%11 lll = - -2- - - 21 "'
L.J h cos 2nmz11 = 2
NA
= 2,
11-1
. , / 4 nh nh
x1mo es 1gua1 a "'N-t =V cos2 ~, ya que :ion - 2-
n nh ) nh
(N-t)=sen ( 2--2- =COS - 2-
S S. PrlnclplO del m4.zlmo 65
~}
2
Escribiendo A1 en la forma A1 = rt2 ( , ; == nh/2-
= nh/2 ~ l\/4, y teni(mdo presente que sen ti~ decrece y
tiene mnimo para ; = n/4, obtenemos A.1 ;;i: 8 para h ~ 1/2.
Para Af'l-l tenemos una estimacin A.N-i < 4/h1, y, por
consiguiente,
k = 1, 2, . .. , N - 1.
tenemos
LY1 0= bi. (Y1 0 +1 - Y10 ) - a1 0 (Y1 0 - Y1.-1) -
O,
- (c1, - a1 0 - b1,) Y1. <
lo que contradice la suposicin Ly1 ;;;is O para cualquiera de
los t = 1, 2, . .. , N - 1, inclu.so para i = i 0 La primera
afirmacin del teorema queda demostrada. La segunda afir-
macin se demuestra anlogamente (basta sustituir y 1 por
-y1 y aprovechar la afirmacin que acabamos de demostrar).
co110L.UU0 1. Si se cumplen las condlclones (2), es decir,
si Lv1 ~O (t = t, 2, ... , N - t), Vo;,;.. O, y,,;;., O, enton-
ces 1 ;;;is O (i = O, 1, ... , N).
SI L111 ;;;is O, Yo ~ O, 11 ~ O, entonces y 1 ~ O
(1 =O, i, ... , N).
DBMOBTllACION. Supongamos que Ly1 ~ O e 1 < O por
lo menos en uno de los nodos interiores i = t. ; entonces !h
alcanza el valor negativo mnimo en el nodo interior, lo
que es imposible en virtud del principio del mximo.
COROLARIO 2. SI cp >o, 111 >o, 111 >o, entonces la
aoluctn ckl probkma (1)-(2) ea no negativa: y 1 >O (t =
=o, 1, ..... N).
COROLARIO a. St quedan cumplidas las condtctones (2),
el problmul
Ly 1 = O, 1=1 , 2, ... , N-1, Yo = O, YN = O (3)
ttene '6lo una aoluci6n lrlvlal v el problema (1), (2) es resoluble
unlvocam4nte, cualesquiera que nan cp,, 111 f,
DllM08"Jl.Ac10N. Suponiendo que la so ucin y1 del pro-
blema (3) es diferente de cero por lo menos en un solo punto
l . = ' llegamos a una contradiccin con el principio del
mximo: sl Jfi. > O (111 <O), entonces 111 alcansa el valor
mximo poeitivo (m{nimo negativo) en cierto punto interior
t = 10 , lo que es imposible. Por consiguiente, 1 - O.
TBORBKA 2. (teorema de ~omparact6n). Supongamos que
111 es la soluct611 ckl problema (1), (2) e y,,
la soluct6n chl
problema
Ly1 = -iP,, 1 = 1, 2, ... , N - 1, Yo = 111
-l/ N = 111
f .S. Prl11el1 1D del mlmo 67
'
d1,,, ~ - b,, c'ii,,+. - ,.>+ cii1. -,,_,) + d1,v1. =
=d1,lw1, I
Si d1, >O, entonces ii,. < 1 w1 , 1 y obtenemos en seguida
la estimacin (12), puesto que 1 u1 ~1 Si d 1, = 0, la ecua-
.(v
cin (t 1) tomar la forma-b 1 1, + 1 -v
1, ) + ' (v1, -i.- s) =
= O, y de esta ltima se deduce que V.,+1 = ..- i = 1, . Porv
70
de modo que
f(:z:)=Yo+ ~ (:z:-:z:o)(.:t-:t1)
-l . . . (x - :z:..._,) 11 (zo, :z:h ., :z:11.)
(poltnomio de Newton). Calculadas las diferencias divididas,
el polinomio de Newton se calcular con toda la comodidad
segn el esquema de Horner
+ +
J (z) = y (z0) (z - x 0 ) ly (z0 , z1) (x - :r) X
X (y (z0 , z 1, z,) + . .. 11.
El clculo de / (z) para cada x requiere n multiplicaciones
y 2n opetaciones de sumaci6n y sus\racci6n.
Existen taJJlbin otras frmulas de interpolacin. E ntre
ella.a resulta ms aplicable la tnterpolacwn herrnitiana. Aqu
el problema se plantee del modo siguiente. Estn prefijados n
nodos {zi}, n valores de la funcin {} y n valores de la
der ivada {yi}; se plde ballar tal polinomio de grado m h mo
2n - t que se verifique
P (z1) = y,, P' (z,) = Y l = 1, 2, ... , n.
Si t.odos l os :r:1 son di11tintoa, exia\e la pit.a aoluoin que
se halla por un mtodo an6logo al de Lagraoge.
Se debe tener en cuenta que l a aplicacin de un polinomio
de alto gr ado puede conducir a Jos problemas diffoUes rela-
cionados con los errores de redondeo.
4. Spllne-ln1A!rpolaci6o. Estudiemos uo caso especial de
la inter polacin polinomial a tro~os cuando entre cuales-
q11iera nodos vecinos de la red la f110cin viene interpola da
por un polinomio cbico (spline-Interpolacin cbica). Los
coeficien\.es de dicho polinoinio se determinan en cada inter-
valo partiendo de las condiciones de conjugacin en los
nodos:
,, = y,
f' (z1 - 0) = /' (z1 + 0),
t (z1 - 0) = !" (z1 + 0), l = t , 2, ... , n - 1.
Adems, en la frontera se ponen las condiciones para z = z 0
y%= z,,
(5)
Buscaremos el polinomio cbico en la forma
f (z) = a, + b1 (z - z,_1) +e, (z - %1-1) + d, (z - z 1_J,
1
z 1_ l ~ z ~ z 1 (6)
De la condicin / 1 = v1 tenemos
I (z,_,) = ' = y,_,,
f (z1} = a. 1 + b1h1 + c,h1 + d,h\ = 111o (7)
h = Z - Zl-1' t = 1, 2, ... , n - 1.
CaJcuJemos las derivadas:
/' (z) = b1 + 2c1 (z - z 1_ 1) + 3d1 (:r: - :r:1_ 1),
r (z) = 2c, + 6d, (z - z,_,), z,_. ~ % ~ %1
78 Cop. 11. Jnterpolocl6n e lnle1rocl611 numlrlca
-
" c11q>,. (z),
11 (z) = ~
donde {q>,. (z)} es una familia de funcionO.!I ortooormalizadaa
en {o, b)
1,k = m,
6"' = { O. k + m,
(q>,., q>,,.) = 6 111
y e,.son unos coeficientes arbitrarios. Entonces, el problema
de encontrar la mejor aproximacin media cuadr6.tica se
reduce a la b<isqueda del mlnimo de la funcin de n + 1
variables c0 , c1 , . . . . en:
:t c11 (/,
(/, //) = 11-0 cp,.)- :t cJ,.,
11-0
/,. -(!, cp,.),
111111- :tel.
-o
$ 1. lnlerpolcl<ln 11 11.pto%1Macl6n de /01 /unclonu 8t
obtendremos
En este C4.90
(1.i)
Z 11 = J/(z)d.z =lt/11 1
(12)
-o
Al compar&r (H) con (12), encontramos
2. Integracin numrica
t. Planteamiento del problema. Frmulas de Integracin
numrica (de cuadratura) ms simple& El objetivo de la
lntegraei6n numirtca es hallar el valor aproximado de la
integral
b
JlfJ= J/(r)dz,
Q
(1)
L llJ - I
1
(J c (3)
a
J
-
/(z)dz-x J(1)ds->cLl/J,
t
O
"-~ -a
r-
1- 1
J 111 = ~
1-1
Jz_..
"11-1 o
etc.
Asi pues, el probluma se reduce a la construccin de la
frmu la de cuadratura para Ja integral (3) que se calcula
por ua segmento unidad. Escojamos en el segmento O< s <
< 1 los nodos O< s0 < s1 < ... < s,,.<1 (molde de Ja
frmula de cuadratura) y a la int.egral (3) le pondremos ea
correspondencia la frmula
(6)
ecuaciones
Po+ Pt + + Pm = 1,
P'9+ p,s, + ... + p,,..s,,. = t/2,
o
Las frmulas del rectngulo y de trapecio son exactas
para una funcin lineal, es decir, para un polinomio de
primer grado, de lo que es fcil convencerse inmediat amente .
. En el caao general, a titulo de P,,. (.s) puede elegirse
el polinomio de inte.r polacin de Lagrange
9 ...u o ~-~
se ve que la frmula (10) es exacta para un polinomio de
grado m, si los factores ponderales PA se determinan segn
la frmula
t
PA = ~ z1."'l (r) ch. (t2)
r
J -ni>"
1 ,(n+I) <t> dt =
o
= - <- ')"1
(n+t)I
,(',..,,<,)!'o+ Jr <- '>"1 t<"'"'>(t)dt =
(n+t)I
o
- ~/<"*ll(O)+
- (n+ t)I Jf <-
1
>n+1 t<i>(t)dt '
(n+ t)I
(14)
o
la cual demuestra la frmula (13) precisamente para n + 1.
( ntroduciendo la funcin
~/ni para ~~O,
K.m= { o para ~<O.
(\5)
(18)
De lo expuesto anter iormente (vase el p . 2) est. claro
que para un polinomio Pn (s) de grado n la frmula (10) es
exacta en dos easo11: para n ~ m + 1 = n 0 , si m es par y la
frmula es si mtrica; para n ~ m = n 0 en todos los dems
casos. Su pondremos por ahora que
/\ IPnl = L IPnl, es decir, n ~ n 0 (19)
Volvamos ahora a la diferencia (/) y sustituy8Jllos
f = P" + Rn+i en (17). Tomando en consideracin (16) y
(19), obtendremos
(/) = /\ lfl - Lf/l =
= (J\ fPnT- L (P,.1) + (A (R,.+ 1 1- L (R,.+11) -
"' 1
= A IRn il-LfR,.+al - ~ P K,. (1.- t)f<"'ll(t)dt -J
-u o
l 1
- J JKn (1-t)f("+ll (t}dtds =
o o
1 "' l
= ~ [~ p.K.(s.-t)- JK,.(s -
t)ds]ft"ll(t}dt.
o ~-o o
HacienJo uso do le expresin (15) para Kn (s - 1) , hollamos
1 l
Jr Kn (s-t)ds= Jr -,.
<-'>"
-1- di= <(n+
t-1> "''
l )l
o 1
para 1/1" 1> (l) J < M,.+i, donde M,.+1 >O es una cons-
tante, y para
1
Cn+1 = Jo JF,.+1 (t)J clt.
Si F 0 +1 (r) no cambia de signo en el segmento O~ s ~t.
entonces, en virtud del teorema del valor medio, tenemos
1
Adfl= J F ( t )/(t)dt,
1
91
'lE(O, ll
92 Cap. 11. /1tl1T/11el61t ' l1tle11'llcl611 numhlco
Jf
o
(z)dz= J
o
ip(z)dz+c,
t t
c=-z (b' -a~)+ T~(b'-a').
o
t
F, (t) - -- (K, (O-t) + K1 (1 - t)) +
+ 64 K (
1
t
T-
t) --u--
(t-1)
De aqul encontramos
P,(t) = -ff.(2t'- 3tt), t< ! ;
t t
P.(t) =n (2(1. - t)3-3(t-t)'), t > 2" ,
F&(t) >O para todos los t E{O, 1),
y, por consiguiente,
1
Jo F&(t)dt = ~
rle 11uerte que es exacta la frmula
- 1 - 1v
Aa (/) = 2Eiiif I (T)), TI E(O, 1).
114 Cap. //. Jnterpol4ct6 n tnt1rocl6n num~rtco
donde N = 2i 0 , h = 1/N.
Si es que / (z) E c1n> (n ;:;;i: 6), entonces podemos obtener
un desarrollo de la forma
DH(f)=a,h'+a..h, a 9 = 0 (1),
1
a, = t!o J/ v (z)dz= 1:0 (/"'(1)-/9(0)).
1
o
6. Aumento del orden de enctitud. Mtodo de Runge.
Para laa frmulas de cuadratura (por la aoalogia con lo
a.n teror) se puede obtooer un desa.r rollo asintt ico de la
forma
D N (/) = J N (/) - J (/) = a.Ji' + r..h + a.J + .. .,
al / (.r.} us una funcin suficientemente suave. En este caso
a.,.
1'1-A ~ . 1 es considerablemente inferior a 1 1 (k = 2, 4),
r<.:i.on por la cual el aumento del orden de exactitud de la
frmula de cuadratura resulta muy importante.
Realicemos los clculos sobre dos redes uniformes con
los pasos h y h 1 , respectivamente, y hallemos las expre-
siones
l~'[fl=Jto 1 (1J y 111'(/J=IH, (fl, h,N,="2Na=b-a.
Exijamos que el error para su combinacin lineal
Y. U> = U 1 (f) + (1- o) u <t>
sea una magnitud de orden superior en comparacin con
D'" y DA . Si para D" = D" tiene lugar la frmula del
tipo
D" = 11 (fl - J (/) = a.pliP + aqli9 + . .. , q > P
f l. /nle(,.cl4n namlrtea
y bailemos
t
p~ Jo A/ln -p
Sabiondo el valor aproximado do p, se puedo, empleando el
mtodo de Rungo expuesto ms arriba, aumentar el ord en
do exactit ud. Con esto fin formemos una co mbinacin J" =
= al" + (t - a) l" y elijamos a do una manera tal que
se verifique ahf + (t - a) h~ = (a + (t - a) pP) hP = O,
es decir, a = pJ'/ (pP - t ) = t / (1 - A). Entonces, para el
error D" = j. - J obteoomos
J N (.r"'f =~
-- o
CAZA = Jo x"' dz = :":;-11 1: =
t
= m+ t , m=O, 1, . .. , 2N-1 ,
obtendremos 2N + 2 ecuaciones para los nodos y pesos
c0 +c, + . . . +cN = i
c.zo + c,.r, + ... + CNXN "" 1/ 2,
de donde encontramos
u<'>= /<'>Ja,,, l= i. 2, .. . , N.
t02 Cop. 111. Ruoloel611 de ceaoelo11e1 ol1ebroleo1
A= [~:: ~
a,.1 a,.2 ~:J
entonces el sistema de ecuacionei1 tiene por expresin
a11u<> =/<O,
a,,u<ll + auu<I) = /12'.
u<"-11 = - -- - (f<n+I)_
On+I . n+I
~ ~ ""'<-'>},
A-1
n = 2, 3, ... , N - t.
Pera hallar el vector u = (u<1>, .. , u<N>) 11e necesitan
1 + 3 + 5 + ... + 2N - 1 = NI operaciones aritmticas.
Si A es Ja matrts triangular iuperlor, es decir , si a 1 = O
para j < t, a 11 =FO (1, j = 1, 2, ... , N),
A [11 :: ... ~: ]
= . . . ~ ... ~;_.: ~-~~ ~"~.:,.a~,, '
entonces el sistema (1) tiene por expresin
a11u<> + a 11u<~ + ... + a 111u(N) = /<O,
auu<O + ... + a,,,u<Nl = fli>,
n=N-f, N-2,
-+
... , 2, f, (2)
lo que tambin requiere N' operaciones.
La eleccin del mtodo y la eficacia de 6ste dependen
de la forma de la matriz A y, adems, del tipo del ordenador.
Para varios problemas A es una matrts 1mrarectda y la
mayora de sus elementos son ceros. Las matrices de este
tipo aparecen con frecue.n cia al realizarse la aproximacin
de d.iferencias de las ecuaciones diferenciales. Los elementos
de tal matriz se calculan, eorrientemente, por las frmulas
dadas, ratn por la cual no han de ser conservadas en la
memoria de acceso rpido del ordenador. Esta cuestin es
de mucha importancia, pues el orden de las matrices de
este tipo puede alcanzar '\lasta decenas de miles e incluso
centenares de miles.
De caso particular de una matriz enrarecida sirve la
matriz de cinta; todos sus elementos no nulos se encuentran
cerca de la diagonal principal, es decir, au = O, si 1 t - I 1>
> l, donde Z< N. Los elementos distintos de cero se dispo-
nen en 2Z + t diagonales, incluida la diagonal principal.
Como ejemplo Interviene una matriz tridlagonal
A= [~;.~;.-~ . . ~....]
O O aN-cN
El sistema correspondiente (1) tiene por expresin
-c 1u<O +b1u< -f<O,
a 1u1HI - c1u<1>+ b1u<1+11 = t<1>,
a,,u<N- 1>-c,,u<N>
t = 2, 3, ... , N - t.
t04 Cap. 111. Ruoluc16n de ecuac/onu al1ebralcu
2. Mtodos directos
t. M6todo de GaWlll. Hay varias variantes de clculo
del mtodo de GanM basado en la idea de eliminacin suce-
siva. El proe68o de resolucin del sistema de ecuaciones al-
gebraicas lineales Az = /, o
/1(
Al escoger
1-1
d 11 = slgn(a11 - ~ 1,.,d.,). (24)
l-1
hallemos
(25)
l
0 1 Cu C2N
C= 0 0 1 ... CJN
.. ........
1
. . . . ...
'
oo o 1
es decir. b, A = O para k > l, c1A = O para k < i, c11 = l.
De (27) se infiere que
4 =
__z
N
,bACA
f l. 111101 dlrtcto 113
c11 = +11
[a.,-~
,._,
1-t
3. Mtodos iterativos
t. Mtodo de iteraciones para reeoller el allltema de
ecuaciones algebraicas llneala Una atencin especial se
prestar en este capitulo a l os mtStodos iterativos, puesto
que dichos mto.dos son de amplo uso en la resolucin de
las ecuaciones en diferencias de la fiaica ma\em\ica cuyos
operadores estn en correspondencia con las matrices de
cinta A de orden superior.
t 14 CoJp. ///. RC1oluc/611 Je ecuoclone. al1cbralca'
..
lm 11 z,. llo=O,
,._ donde 11 z llo V (Dz, z), D = D* >0,
D: H-H.
En el caso general se prefija cierto error (relativo) e> O
con el que se debe hallar la solucin aproximada y 11 , los
c lcu los se clan por terminados cuando queda cu mplida la
condicin
11 Yn - u llo,.;; e 11 Yo - u llo (6)
Si n = 11 (e) es el minimo de los nmeros, para los cuales se
verifica (6). entonces el nmero total de operaciones arit m-
ti cas que han de realizarse para hallar la solucin aprox ima-
da de la ecuacin (1) es igual a On (E) = n (e) q0 , donde
q0 es el nmero de operaciones que se realizan para hallar
una iteracin, es decir, para resolver la ecuacin (4). El
problema consiste en minimizar Q,. (e) eligiendo de modo
adecuado By los parmet ros {i:,. ). Comencemos por analizar
los mtodos iterativos ms simples.
2. Mtodo de la Iteracin sim ple. Para la resolucin del
siste ma ele ecuaciones (1) puede emplearse el mtodo de la
H6 Cap. IIJ. Rolucl6n d "4,lont1 olgebralca1
iteracin 1tmple
J-1
N 1- 1
11 = - '- (111>-
111A+ 1 0 11 ~
L..J a,1111n-"'
A L..J a, 1g1n
A+ 1 ) t
;-1+1 J-1
L= 2, ... , N.
Hacie nd o uso de (10) encontramos sucesivamente para
1 = N, N - 1, ... , 1
N- 1
y\N)
Ji + t
= _NN
I _ (<N> _ "'
L.J
QN ,AJ))
JllJi '
1-1
1-1 N
t,>
A+I
=- -'-(J<ll_
0
"'
L..J 11!11i>
A
_ "'
L..J a,,yjJ>
A+ I ) '
11
J I J-1+1
l = N-1, . .. , 1.
Escribamos este mhodo en la forma matricial (operacio-
nal). Con este fin representemos le matriz A como una su ma
A=A-+ D+A ..
donde D = (a 116 11) es una matriz diagonal de dim ensin
N X N, A -= (ajj) es la matriz triangu lar (subcliagonal)
inferior co n lo diagonal principal llenada ele cer os, ii = O
ia ra j ;;> I, ii = a 11 para j < 1, A = (aj1) es la matriz
tr iangular (sobrediagonal) superior con la diagonnl principal
tt8 Cap. 111. R..ola dllll .U 1ca11c1on., a/1ebrolco1
Ay<I> =
;-1+1
"
" 1111111, (A+ D) 11<1> = } a1y<i>.
}
;-1
Por est-0, la ecuacin (10) puede anotarse en la forma
((A+D)y., 1)<1>+ (A-y.)< 1>=/<1>, i = i , 2, ... , N,
o bien, en la forma vectorial,
(A + D) l/u 1 + A -y. = / .
Rea lizadas unas transformaciones evidentes
(A + D) Yu1 + A-y, = (A+ D) (Yu1 - 11) +
+ (A- + (A + D)) Y = (A + D) (Yu 1 - 11) + Ay,
escribamos el mtodo de Seidel (10) en la forma cannica:
(D + A+) (Yu 1 - Y) + AyA = f, k = 0, 1. 2.... (11)
Comparando con (2) vemos que el mtodo do Seidel (10) co-
rresponde a
B .:... D + A +,
ea decir, el esquema (tt) es implicito. Sin embargo, por
cuanto 8 = D + A <Js una matriz tria ngular, entonces
la iteracin Y+ se determina por las frmula s explcitas.
Anlogamente se escribe la otra variante dol mtodo de
Seidel:
(D +A-) (Yu1 - y,) + Ay = /, k = O, 1, . .. , (12)
cuando 8 = D + A - ea una matriz triangular inferior.
A continuacin (en ol p. 5) se mostrar que el mtodo de
Seidel converge, si A es una matriz simtrica definida
positiva.
4. Mtodo de ret.jac16n auperior. Con el objeto de ace-
lerar un proceso iterativo, se puede reduci r el mtodo de
Seid11l al de relajacin superior , introduciendo el parmetro
3. Jlllotlol ltCNlllllOI U9
De aqu encontramos
f- 1 N
A+ 1 =
y<fl A +~
!ll 0 11 [1111 - ~
LJ ~ 111/\A11 .]
111/\'A+1t - LJ t
J-1 ;-1
i = t, 2, ... , N.
Cuando 111 = 1, ob tenem os la frmula del mtodo de Seidel.
La velocidad de convergencia del mtodo de relajacin
su perior depende del parmetro (1). En el p. 5 se mostrar que
pnra la convergencia del mtodo se ha do exigir que O < (1) <
<2.
5. Convergencia de los mtodos iterativos estacionarlos.
El mtodo de Seidel y el de relajacin su perior sir ven de
ejem plo de los esquemas implfcilos de la forma
B>TA. (16)
~11
T
Zu1- Z ll'+IJ tl+1 ll~~ll s. JI~ ()
Dado que (11 ' llA} ~ convergente, de aqul se infiere que
ex iste
llm 11 ' + 1 -
~oo
' 11 = O. ( 19)
Ai = _ _!_B(z
~
,-z.). ' = - ..!..A-B
T
(z 1 -z.),
11z11~ ~~
l
11 A-11111B112 11 '1-z111 ( . .)
122 Cap. /JI. R1olucl6n de ecuocl<mo algebratca1
para
11 11
!
- ~ >O. El mtodo de la iteracin simple
converge para todos los valores de i:. que satisfacen lo desi -
guaJdad T < 2111 A 11.
M2TODO 02 SEWEL, B = D + A - . i: =
1. En este caso
t _t _. D t _
B-2A = D+A - 2 (A +A+D) = 2 + 2 (A -A),
((8- + A) y, y} =+(Dy, y)+T((A -A-)y, y)=
t .
=T (Dy, y)> O,
siedlpre que D > O.
oe sERVACION. La desigualdad D > O proviene de Ja con-
dicin A > O. En efecto , supongamos que A > O y ~ =
= (;1. O, . .. , O); entonces (A;, ;) = (.D~. ;) = a11 (t1 ) 2 >
>O, es decir , a 11 >O. De un modo anlogo nos convence-
mos de que a11 >O, y, por consiguiente. D >O. As pues,
\\l mtodo de Seidel es siempre convergente, si A es un
operador aut()C{)njugado positivo.
Para estimar la velocidad de convergencia se deben esti-
puJar las suposiciones mss fuertes. Citemos el siguiente
'l'BORKMA 2. El mtodo de Seidel converge con la velocidad
de una progresi6n geomitrica de ra.z6n q < 1 , si A = (alJ) =
= A> O, y
l+ N
1:
.-1
la1Jl ~9 a,,, l=i, 2, .... N, q<L (20)
123
J<I J>I
Supongamos que el mh lz~1 1 118 alcanza para cierto t =10 ,
de modo que
~ fa1,1I 11 "1 lle+
11 z_., lle = lz~1+11 l. la1 1,l ll z , lle~ J<lo
~ la1,1I 11 'lle
+ J>lo
llz 1lle ~ I ~ fa1,1l/(fa1.1,I- ~ la1,1IHllzlle
J>lo J<lo
F:n vir tud de la condicin (20) tenemos
~ l a1 0 1l~qla1.1.I-~ la1 01l<q(la1,1I- ~ la1,1ll.
J>lo J<fo J<lo
y, por consiguiente,
11 '1lle~qll ' lle~</' llzolfc,
1
= (t --i ) D + ~ (A- - A)
y calculemos
A= --F .
[
1 -2 1 .. .
0
o
... i - 2
... o -2
]
.~ .
Ahora tenemos
El nmero de iteraciones:
n 0 (e) :::::: 32 + 60 para N = 10,
n 0 (e) ::::: 320 + 620 para N = 100,
ea decir, considerablemente menos que para la iteracin
s[mple.
El esquema implcito de tres capas tien.e por expreain
Brut = (1 +c:i) (B--roA) 1111 -c:iBy,._, + (t +c:i) -rJ,
""""1, 2, . . ,
Bv1 = By0 - ToAVo + -rol para cualesquiera Yo E H.
Si B = B > O y se cumplen las desigualdades (23). (24)
mientras que c:t, T 0 se calculan segn las frmulas (35). enton
cea la estimacin (36) para el nn:1ero de iteraciones queda
justa en este caso tambin.
4. B111111m11 U.Nlluo de do1 captu 129
...
,.
8 ..... - +Az.=0, k=O, 1, 2, ... , :r0 = y 0 -u,
por la correccin w
= B-1 (AY - f) (k = O, t, ... ) con
la condicin inicial Wo = B- 1 (Ayo - n.
La condicin para
que terminen las iteraciones tiene por ex presin
11 z,. llD ~ e 11 Zo llo. o bien 11 w,. llD~ e. 11 wo llo (4)
De (3) se ve que
(5)
donde Sti es el operador de paso de la capa k a la capa
k +
1. Eliminando z., '- .. . , z1 , encontramos para
k = n - 1:
z,. = T,.z 0 , T,. = S,.S,._ 1 S 1S1 ,
donde T,. es el operador de resoluci6n del esquema (3). De
aqui proviene que
'
11 z,. llo ~ q,. 11 llo. q,. = 11 T n llo (6)
La condicin para que t.erminen las iteraciones queda cum-
plida, si
(7)
Para estimar el nmero de iteraciones n = n (e) se Jebe
obtener la desigualdad (7).
Estudiemos el eaquema expllcito (2)
11u-n
'f... + Av = f . k = O 1' 2 (8)
parmetros
k = {, 2. ... , n,
t +Ptl'A '
2
'fe= Yi+Ya ' (9)
= In t +
t- ~
,. i > 2 'Jl1, sustituyamos (11) por el requisito
ms fuerte:
t 2
n (a)> "o (a) = .n In-, (12)
t,. i
De aqu se ve que
nCfJ ~ 34, n<J> ~ 200 para N = 10 (h = 1/ 10).
n<f> ~ 340, n<J> ~ 20 000 para N = 100 (h = 1./100).
2. Argumentacin de la eleccin . optima) de los parme-
troe. Demostremos la estimacin (10) en el caso de par-
metros iterativos (9). Nos hace falta hallar el mln 11 P,, (A) 11.
El polinomio
P,. (A) = rj
A- 1
(E - ;AA) =
- mx 1 T,.(1) ,.,. _ _ t_
-1..;;1..;;1 ir;;;) 1r,. (te)I '
11,1 Vt&-t=...!... .. /
P V
-A.-
Pi
t =..!_11 Vt-p&J=
Pt
= _!_
P
(t 2
1-
Yf) ... _!_ c1 Vlt(t-~>-
P1
= lt V[)t(t-~) = c1 V~vct :::V~>.
da modo que t,l+Vtl-t - _!_, (t1 1-Vti-t-P1 Y
PI
t(tPi'-+P~ )=~'""t:
IT,. lloJl =2 ~ t
t1 = COS 2n,"
2.t-f
<X.s)!'i. = -cos~n.
TABL~ 1
Surtido a 1
1 Surtld01 1
53 O,t2 t 39,6
55 27 2 2,6 tos
57 t,9-104 8,2.10
59
60
3,7107
2,610'
"7
9
S,3-tOU
t,210"
6t 2,510" it t,9t0U
62 s,s.1011 t2 Parem
63 51011
64 Parem
us,..,11 = "'"+IA-Ell = {
137
de suerte que 11s ,11 < 1 para todos los I' >O, y llS1i+ill > t
para 1'11 < - (1 -p0)/(2p0). Por cuanto
n (2n-1) n
-cos z;- ~ 11 ~ -
2n-- n = cos z;-,
cos - -
k = 1, 2, ... , n,
entonces, para la mayor cantidad de nmeros k tenemos
11 S 11 11 > 1, y si se emplean muchos parmetros 't 11 seguidos,
pera los cuales 11 S 11 11 > 1, tienen lugar la acumulacin
del error de redondeo y el crecimiento de las aproximaciones
iterativas, lo que conduce a la inestabilidad com putacional.
Con el fin de debil i tar el efecto mencionado, es natural
tratar de alternar los paramelros 't 11 , para los cualos 11 S 11 11 >
> 1 con aquellos. para los cuales 11 S 11 11<1. En este pro-
cedimient o se realiza precisamente la construccin de una
sucesin de parmetros ('t 11 }. para la cual la convergencia
de las iteraciones tiene un carcter montono y la inestabili-
dad com putacional est ausente. Existe una regla para tal
ordenacin de los ceros t 1 = - cos ~ cosn del polinomio
21 1
O,t
32
0,06
50
~ ..
s -a- + Cz =- O, k=-0, 1, ... , z. = Bllw,,
donde .i: = s112w C = B-lflAB-1/1 es un operador posi-
(25)
= ((D + wR 1) y, D - (D + wR 1 ) v) =
1
R 1 ca (riJ), riJ=
{ r 11 /2,
r,,,
,... ,,
J<t.
R 2 = (rlJ), rl1=
o,
{ r 11 /2,
r,,,
,... ,,
J> t;
J>l,
O, J<t.
Si R es una matriz sim6trica, r11 = r 11 , entonces R 1 y R 1
son reclprocament.e eonjugedos, R, = R~ .
A titulo de D = (d11 ) tomemos una matriz diago nal.
Entoneea, D + wR 1 es una matriz triangular inferior y
D + wR, , una matriz triangular superior. De este modo,
el proceso de la Iteracin se reduce a la inversin alternada
de las matrices triangulares inferior y superior (de aqu
proviene la denominacin del m6todo). Efectivamente, para
cada iteracin se debe resolver una ecuacin
B11u1 = (D + wR1} D - (D + wR,) l/Hi = F,..
1
(3)
Al denotar D - (D + wR,.} llHt = i11.+i obtenemos
1
(
~f 1+ o>6 + e>~} R ,
1
y,-6 11 (
613 +-,-
Clll 6A )-'
.
El nmero de iteraciones necesarias para Ja resolucin
de la ecuacin Ry = f depende de la rain
1=1, 2, .. ., N- 1,
Ui, = 0,
Sea B = O un espacio de fu nciones reticulares definidas en
los nodos interiores t =- 1, 2, ... , N - t de la red <i> 11;
introduzcamos un producto escalar
li- 1
(y, 11) ... ~ y 1v 1h .
-
f S. llltotlo eU.:rudo lrten1olor 145
en definitiva tenemos
y (I) = "'11 (t + t) + A' v+i (1)
HI IO+A'
i = N-1, N - 2, ... , 2, t.
es decir, A1 = A:.
Calculemos la constante ~:
V-V
'f.t.u.
+AY = /. k = O, 1., . , Para todo Yo EH. (1)
"f+iAr 11 ~ 11 rh - "foAr 11 2 ~
2
11'.+1112 = 11 ' -
~11 E - "'oA11 2 11r,.11 2
Por otra parl.O se conoce que
11 1!.' - T0A 11 = Po pura i:, 2/(y, 1 Yt).
Do aqu precisamente se deduce que 11 ' H i 11 ~ Po Ur,. lf.
Do esto modo, el mtodo do los residuos mnimos con-
verge con la misma velocidad quo el mtodo de iteracin
simple (siempre que en este ltimo se empleen lo!! valores
exactos do y 1 Y y,).
En el caso dol mtodo de residuos impl!cito o del mtodo
de c.orrecciones en vez d e (1) obtenemos una ecuacin para
la correccin:
B "' - "' + Aro. =0, k = O, 1, ... ,
1' w,. = B-1r11t,, w0 = B- 1 (Ay0 - /), (5)
donde T t i se determine por la frmula
TH1 =
(Au>1t . ... ,.
(S-lA"' AIOA) , k=
o, 1, (6)
En lugar de (4) obtenem os la estimacin
11 A11,. -lllB"~P: 11 Ayo - /11.s-
~ ms r pido. El mtod o ex plcito
2. M todo del
dll dlscenso ms rtfptdo se d iferencia d e l mtodo de los resi
duos mlnimos slo e n la fr mul a pnrn TH 1:
k -= O, 1, ... (7)
6. IUlodo1 llcrelho1 d4 llpo 11ar'4clo11cl t411
t-l"t
Pi= t+yt , ~=..h,.,, (iO)
Znt = +q>mZ~.
Do aqul y de (20) se despren de
(21)
Denotando 1111 = qz,,I. obtenemos 11,.+ 1 ~11! ~ .,!~ 1 ~ .. . ~
~ ... ~11f'~11f ... y, por eqnsiguiente,
t .....
s,.+d~'"f(qf s,I) . (22)
t56 Cop. /11. Re1olacl6n ecaoclon" ol1abra1<:01
que coincide con (21) (la diferencia consiste slo en g). Los
raaonamienio. ulteriores no. llevan a (22), (23), y (24).
iS7
L' 11=""..T
-v,.
--
+'-11n-f
- "'
J. Teorla do "I'""'' dt dl/trenela1 t63
&raNPLO 2. Primer problema de contorno:
Lu =u= -1 (z), O< z < 1, u (O) = 1 ,
u (t) = 1'1 (6)
Hagamos uso del operador de diferencias t ripuntual (2):
Lll1 = lliz 1 = (111+ 1 - 2111 + 111_1)/h1 y obtendremos un
problema de contorno en diferencias sobre la red e;,. =
= {z1 = ih, o~ t~ N, % N = t}:
L,.y, = 11;;,.,, = -/, i = 1, 2, ... , N -1 ,
Yo = l-'1 llN = 1-'1 (6')
3. Estabilidad. Nos resulta ms conveniente pasar a la
notacin del esquema de diferencias (5) en la forma operacio-
nal. Con este fin escribamos al principio la ecuacin (5) en
le forma matricial
AY h = I>11.
expresin
N-1 1 N
11 Yh - Y11 11~ t?
188 Cop. IV. Mllodo1 de dlfcr.,1cu de lo ruolucl611
U81 lle.~ e.
De aqui determinamos el paso admisible h;;;;. h 0 , es decir, h 0
Expliquemos esto con el problema concreto (t7). Para
dicho problema tenemos
.. 4 nA 4
11BA'11 = 11A1111 = u =-rcos1 2 ~ -;;r.
T
t (b 1--,.
,,,.. --,,,--a, -,,,_,,_,
-,.-- ) -
.J
u11/1 = -cp,, (3)
d1 = (c1 - a 1 - b1)/h1
Diremos que un esquema de diferencias es homlJglneo, si sus
coeficientes en todos los nodos de la red para cualesquier a
coeficientes de la ecuacin diferencial se calculan segn unas
Jllsmas frmulas. Asi, por ejemplo, si introducimos las fun-
cionales 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 frmulas
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 homogneo. He aqui las funcio-
nales ms simples
A (k (s)J - k(-0,5), a,= k1 -112 - k (:r1 -0,5h),
P (i (s)) = f (0), cp1 = / 1 = / (z1), et.e.
Si un esquema es homogneo, resulta ms cmodo servirse
de l sistema de designaciones sin indices:
1
A, =--(byx-aYi) - dy - -cp, :rECJ>,,,
J
a,- {T
t "/
.l
" 1-1
ds )-'
k(z) =
{
- 1
d.
i(1+11\)
) _,
;
"l-1(2
f(z)dz= )
-112
'I
f(z1 +1h)dl,
'[
d, - {- I
"1+1/1
g(z}dz= J g(z1 + 1h)dl.
"1-t/2 - 1/1
de donde se infiere
t
11 s lle~ 2t, 11 ;>lle
Mostremos que para la solucin del problema
(as;;),, - az = - ""' z E "'" 1 (O) = i (t) = O,
a ;;;;i.c 1 > 0, d>O
es vlida la estimacin
Uslle~.!.
., (t, 11111. (9)
N
donde se designa (!/, 11) = ~ 11111,h.
f -1
Representemos 1 en forma de una suma % = w + u,
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 (1) = O
t80 Cop. IV. Jlltodo1 .U dl/trenc/41 d~ la rc1olucl<ln
ID1 = ID1-1 + ..
(co- 111) /1 = e, ~
N
1
-1
~-
~
N
h
-
1
..
.l!.. h + w,,
0 = 1DN = c0 ~ ~- ~ ~h
LJ al LJ A '
N -1 N -
Co =
-
..
~_ah~~.
-t
...
encontramos
De aqui proviene
1 N
w,i ... l-(1-a1) ~ ~+a1 ~ ~1 ~
... , -1+1
t N N
.,.- (1 - a) ~~+a ~ ~ ~ ~ ~.
""'> 1 Jt.J
k-t
1 Jt.J
-+t ,
-s "
Ahora nos queda tomar en consideracin que 4 11 ;.. c1 > _0
y obtenemos
(10)
3.E1qaam1JJ1 de dl/arencla conun1atlrJ01 t8t
-
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 ), ii = O (h'). Ahora resta por uti-
liu r la esti macin apriorist ica (t3), la que pod ra ser sus-
ti tu ida por una estimacin ms aproximada
2
11=11e~c.<ll TJlle+lf 1j)* lle)
De aqu se desprende que el esquema (5) converge unifor-
memente con el segund o orden , es decir, 11 z lle= 11 y -
u ll e~ Mh1 , si se cumplen las condiciones (7).
182 Cp. IV. Mltou1 114 tll/orcncl1u tic le rc1olacl6n
obtenemos
,.,,..1 = w, [ " ( ~ + ;.9 )- 1 +o (Ja), J
ea decir, 'ln+t = O (t) para cualquier esquema y slo para
un eequema con coeficiente
~. - [+~ k~) r
1-1
f 4. B19N,.... 1ob,. lcl NtUI M lfo,..u 188
h = I - It-1'
(a11-) -
. ".
d11 = - cp, :r = Z E mh l/o = 111 11 H = 111 (t)
Al suponer J = 11 - u, obtendremos para s la ecuacin
(as); - ds = -;>, x E~'" ' = 'N =O, (2)
donde
(3)
7
Sustituyendo esta expresin con = qu - f en (5), llegamos
a la frmula (4).
Para estimar 11 1 aegn el orden veamos la diferencia
(au;)1 - (ku'),_111 a condicin de que k E e<>, u E C<'>.
Empleando la suposicin a1 = k 1_11 1 +
O (ht) y las frmu-
las u, = 111 u,_ +
h1ui..vJ2 + htui...11 J8 O (h~). u1_1 =+
= "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 =
= (k1_1, 1 + O (hl))(ui..111 + O (h})) - (ku'),_1, 1 = O (ht).
De este modo, es vlida la estimacin
111 = O (hl) para (le (z), q (z), J (%) E c< 1 >, u (z) E cea>.
Se supona que d 1 y qi1 ae determinan segn
OBSl!RVACION.
las frmulas ms sencillu: d 1 = q,. cp 1 = (1 Si, en cambio,
so emplean las frmulas ms complejas, por ejemplo
"l+l/J
epi = ""- 112+"1+1!1+112
2i1
t
q>, ~ "ij" Jr (z) clz,
"t-t/J
entonces la funcin reticular ;ir
= 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
frmula (4) se sustituye por la suma 111 + p 1:
,> = P1 + '11); ,>. ' + (4')
P1=o (ht) '11 =o
(ht), ,>f9 (llT} =o
para le, q, /E C'''. u EC'".
3. &umad6n de la nloddad de conYerencla. Para el
problema (2)- (4) es vlida la aiguiente estimacin aprioristica
llllo~-?-l ((t, l '11 J+(1, 1;i1) }, (6)
N
donde (v. v)I ~ ~ v1v 1h 1SI ee cumplen tu oondclonea
t-t
(7) del 3, entonees '11 -0(.\1), 'lj>f ... O(~.
Al sustituir 11 1 y 1jif en (6), nos cooveoeemoa de que es
varidico el aiguiete teorema.
T.ORBllA. En la clase rk c~ficicnU ~ le, q, f Cfl>e.
eO<Jq eStuema de la forma (1) mantiene el ugundo orden de
exactitud en una sucesi6n arbitrarla. rk la.8 rerks no uniformes.
Al tomar en consideracin la observacin del p. 2, pode-
mos 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 flstimacin
JlsJlc<;-?{(t, 11l+PIJ+(1, J,pIJ);
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 nmero finito de puntos, siempre podemos ele-
gir tal red no uniforme ~h (k) que los puntos de disconti-
nuidad sean los nodos de dfoha red. En tal caso cualquier
esquema tendr el segundo orden de exactitud.
Asi pues, cualquier esquema homogneo de segundo or-
den de aproximacin (11> = O (h')) sobre una red no uniforme
y en la clase de coeficientes suaves tiene el segundo orden de
exactitud con la eleccin especial de las redes no uniformes
w,. (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 solucin en los
nodos de una red arbitraria coincide con la solucin exacta
u = u (x) del problema de contorno para una ecuacin 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 ecuacin desde % 1 hasta :r, obtendremos
una ecuacin
(9)
af = [
t
T," J
r t1z ] - '
k(s) .
"1 -1
=u h'
uIV - (d-q) u+ <p-f +O (h').
Al sustitu ir aqu uIV = qu - r = q (qu - f) - r = qu-
q/ - r. obtendremos
~i
"' (q/ + /") .
= / +12 El orden de exactitud queda intacto,
en la frmula para <p sustituimos la derivada por su r
aprox imacin de diferencias /iz, puesto que h'f" = h'fr:c +
+o (h 4
).
El aumento de exactilud del esquema disminuyendoh vi&-
ne limitada tambin por el requisito de la cconomla del
t iempo indispensable paru la obt encin de la solucin con
t90 Cap. IV. Mltodo1 tk dl/er~ncl121 de /4 r.,olucl6n
"l+ l fl
'4>1 = ! J / (z) dz.
"1- 1/1
v,=x.v,-J.l.,
- --2h_
x1~
o,+.a, J.l1 = 1 1 +lo- 1
Con el fin de aumeuta.r la exactitud del esquema al calcular
las interales se debe utilizar la interpolacin de orden ms
elevado.
2. Mtodo de aproximacin de una funcional cuadrAtica.
Un problema de contorno
Lu = (ku')' - qu = - f (.z:), O< .:z: < 1, u (O) =O,
u (1) =o
es equivalente al problema de buscar el elemento minimi-
udor de la funcional cuadrtica
1 1
Ju - ) (k(u')+qudz-2 Jfudz.
u u
lntroduzcnmos en el segmento O~ z~ t una red (;) 11 =
= {z1 = ih, i = O, 1, ... , N) y aproximemos la funcional.
S 5. llltodo1 de con11rutcl6n d<1 101 q11cma1 I~
J (u) = ,_,
~ J,(uJ, J 1 (u)"" 1
.,_,
(k (u')1 + qu-2/u) tlr,
1lnll(IU~s de lo cnnl nproximomo!'I J,. por ejemplo, 11si:
") le (u')' dz ~ a 1 (U, 1) h,
es decir, (a;y)., - qg = - /.
4. Mtodos de Rlll y de Bubnov-Galerkl n (mtodos
varlaclonalee de dUerenclu ). El problema sobre el mlnimo
S S. Mltodo de co11alruccld11 de loa e1qucma1 197
de una funcional
/ (u) = (Au, u) - 2 (u, f),
donde A es un operador lineal autoconjugado y do!inido
positivo en el espacio de Hilbert H con el producto escalar
(x, y), es equivalente el problema sobre la resolucin de una
ecuacin
Au =f.
So introduce una sucesin de espacios do dimensin finita Vn
con base (epi"'}, i = t, 2, .. ., n.
El mtodo de Ritz consiste on que se busca un elemento
u. E v. que minimiza In funcional 1 (u) en v. La solucin
uproximada Un se busca en forma rle le suma
Un =
" YJ'Pb
~ (5)
J-1
l lunl = ~
l.J-1
a.11Y1Y1-2 ~1Y1t
1-1
~ 11Y1-~1=0,
1-1
t=t, 2, .... n,
para determinar y1 , y 1 , ., y,..
Ilustremos el mtodo de Ritz con un ejemplo del proble-
ma (lt). A titulo tlo In funcin ip1 (x) toml\mos
s < -- i<s<O,
t. > 1, 8
cp 1 (x)='l ( .r-;;-z: ) = ri,(x), ri(s)
{ 1+s,
.1 -s, 0<s<1.
t98 c.,. IV. Mllodo1 d dl/crn<lu de lo ruoluc16n
a: 11 = (Acp,, qi 1) - .~ ( k .;!_1 ~ +
4 4
q')1'1i) dz,
o
1
(5)
donde u'' ' (0) = ~~ (0) = /(O, u0), y los valores de las deri-
vadas u<- >(0) (k~2) se hallan mediante le diferenciacin
.~llC~!li Yn cJu 111 N'llBl'in ( 1)
1z:,.+ 1 1~ I to I + i;
...o
't l 'IJ, I = :i1
....o
't l ..P. \, (H)
,,,.t:/2
- 1/n -- /(t n y)
n ~
't
-/
-
(t n + ...!..2 -Yn )
t
(15)
n = O, 1, .. ., Yo = u.. (17)
c.,. V. Problema de Cuch11
(19)
J/n+1-J/n
T
+ AYn =O, Yn +I = (1 - TA} Yn. (23)
T lln +A.Yn+i = O
lln+t - 1
Yn+t = 1 +Tk Yn
(26)
ftn+tf~ft,.f + i:flj>,.f~tT
a-o
flj>hf, (33)
~"'
Yn =o;t L.J t
(b,.Tf,._. - a.y,. _,.) = ~ F(y,._,, Yn-t y,._,.).
-
Los clculos empiezan con n = m. Para hallar y,., se deben
pre fijar m valores iniciales y0 , y1 , y,._ 1 ; stos pueden
determinarse, por ejemplo, mediante el mtodo do Runge-
Kutl a en el que se utiliia solamente el va lor inicial de
!lo = "o
Si b0 . ; : O, el esquema (2) se denomina implicito (de in-
lerpolarin): al hall ar y,., se debe resol ver , para cada n,
una ecuaci n no lineal
4o!/n - bof (1,. , !In) = F (!ln-1 Yn - t ., Yn - m). (3)
Dicha ecuacin no lineal puede resolverse, por ejemplo, mc-
diente el mtodo de Newton.
El error do aproximacin del esq uema (2) en la !olucin
u = u (t) de la ecuacin (1) o el residu o se determina por la
rmule
11>,. =
m
2;
-o
b.f(i,._. , u,,_.) -+-oh "-
m
(4)
(6)
+'f '"-1
f (t, " (t)) dt,
i>,. =
~
~ bf (t,._., u,._,.)-+ '" f (t, u (t)) dt. (t2)
-o '-
3. Esquemas explcitos e Implcitos. Si bo = O, el es-
q nema (H) ser explcito y
( 13)
-1
Oo ejemplo mn!I simple del esquema explicito de Adams sirve
~I do Euler
y,. - y,._, = Tfn-i pnrn m = 1, b0 = O, b1 = 1. (t4)
2i6 Cap. V. Problema do Cauch
si ponemos
u1 = u1 + 1,t:/ (O,
1
u 0 ), u1 = u0 + t:"1
Si b_, = O, el esquema (24) ser explicito, puesto que en
el segundo miembro figuran slo valores conocidos de y,.,
!/- . . . , y,._,,.. Si b_1 = O, el esquema (24) es implcito
y para determinar y,.+ 1 se debe resolver la ecuacin
!/+I - b_if (l,.+l !/n+1) = F (y,. , !/n-1 !/n-m t.).
Para obtener el esquema de diferencias (24) calculemos una
integral
'
} u11 dt =
... u11 dt + ..,,J u11 dt =
'-t t._, ' ...
- (u'11-iw')l:: _1 +(u'11-uv')/:;_1 + u11dt, (25)
111-1
s a. B1qDMOI de uorl<ll , ..... llllodo1 "' Adam1 2t9
.,." -= ~
..,
~ n-'- u,.__) -
b/(t u,..,-2... + u,._,
,. '
---
la cual puede ser escrita , en vi r tud de la identidad (28), en
la forme
11>. =
"'
h bAfCt. __ , u. __
,. _ _ ,
'"1
)-+ /
~ (t, u (t)) 11 (t) dt . (29)
1
n..-1
+'"r'J1n-l
F (t) 11 (t) dt = 1F (t,. + n)
1
- t
(1) d1,
_ {t +1, 1<0,
F = f (t, u (1)), u(s) - t - , i >O.
De (29) se ve que el primer sumando es una fr m11l11 de
cuadratura para la integral de la funci n F (1) = f (1 , u (t))
220 CI' V. Probltma tl Caueltr
---=
""' f 1 (t, u.), t > O, o 1
u 1 ()=u,, !=.,
" ... , N .
..,
(31)
Sean u, 11 dos solucione~ del problema (30) c.on los datos
iniciales u (O) = u0 , 11 (0) 110 Para= su diferencia i =
=u - 11 obtendremos un sistema de eeuacionea lineales
ria' ~
.
di= ; a.,, (t) ' '
S- 1
donde a.,,
es el valor de la derivada ar11Jui en cierto punto
u
medio (t, 1), 1 = (111, ~ . .. , vl-1 u' +e. , 1, u'+ ....
. . . , uN) (O ::;;; e1 ::;;; 1, J = 1. 2, ... N). Por eso el modelo
lineRl del sistema de eeuaciones no lineaies (30) ser' repre-
~ntado por un sistema lineal
N
t= le
mlnAA
~ t.
A
Son rgidas, en particular, las ecuaciones que se obtienen
ni reducir las ecuaciones con derivadas parciales 11 los siste-
Crap. V. Prol>lunra do CrauchJI
,_ t, 2, . , N . (t)
Al designar por A = (1111) una matril cuadrada de dimensin
N X N con elementos 110 que no dependen de t, por u (t) =
= <u.' <1>. u. <t>, ... u." (e)) 1 1<e> = (f1 <t). r <t ...
. . ., f" (t)) loa vectores buscado y prefijado de dimensin
N, respectivamente, eseriba.mos el sistema en la forma
di>
Tt+Au=f(t), t >O, u(O) = u 0 (2)
y, por lo tanto, ": + A" =0, " (t) = ex.. (O) . -..1, de
suerte que
lf N
llu (t)ll' "" :Z cxl (O)e-U' ~ ,-a>-1 ' ::Z cxl(O)=
A- 1 A-S
==,-u, 'llu (0)111
2. Fauem de dllerend& Introduicamos una red de
paso 't segn la variable t: w~ = {t,. .. n"I", n = O, 1, 2, ... }
y designemos con 11,. = v (t,.) una funcin reticular del argu-
mento 1,. = n"I" (o bien n) con valorea en HN. Escrlbamoa wi
esquema explicito
ru1 - rn
~
+ A.11,. = /,., n ... o, t, 2, .. . ' v.=11o. (5)
(t4)
(BS, Sv) = ~
N - al,
-1 -
= 11s1111 (811, 11>.
donde
USllL = mx ( t - i:A,,)1 , (28)
l<;U;N
La desigualdad (25) es equivalente a la eondici6o
TA~ 2, le = 1, 2, .... N, (29)
la cual, a su ve~. es equivalente a la desigualdad (14). pues-
to que
5. p-establlldad.
TBOREMA a. St A = A > O, B = B > O, ernonus la
condtci6n necesaria y suficiente de la p-establlldad del eSIJue-
ma ( 2) con p > O cualquiera:
llYn+11lo~PllYnllo1 D=A,8, (30)
ser representada por las desigualdades op11raclonales
t -"t p B &
"";;;:::
A&
-..;;::::
t+p
't
B. (31)
"'
+11(1,, 1, - t)-211(~11.+11(1,, 1.+t) - / {i, fs), (3)
111 au 111
/Pu
+Thf (z,, za) 6 ar, (z,, zz) +
Cap. VI. Bcaaclonc clfpllco1
encontramos
au au ) hl 8'u - hl au -
IP= ( Bzf +Bif- /(z) + 24 Bz~ (x,, ~) + 24 8.z1 (z,,z:)
De aqu y de (1) proviene (9).
Asi pues, el esquema (4) es de segundo orden de aproxi-
macin.
4. Esquema de orden aumentado de exactitud. Haciendo
uso de un mol de nonapunt ual (z,, x,), (Xi h 1 , z 2 ), (z1),
z 1 h 1 ), (z1 h, z 1 h 2 ), podemos construir un esquema
que tenga el cuarto orden de aproximacin (y de exactitud),
si suponemos que la solucin del problema (1) - (2) u =
= u (x) E e<> (G). Dicho esquema t iene po r expresin
q> = /+ ~
hf para l
%1= 1 - h1 o<:cs<-sr
1
obtendremos
donde
0 4 nA 1 4 nll 1
1 =Tr een1 Ta", 1 = T? cos1 T,"
Al sumar estas desigualdades 5egn i 1 ""' 1, 2, ... , N 1 - 1,
obtendremos 01 (//, 11) <; (A 1 y, 11) ~ 6 1 (//, y).
Del modo an logo encontramos 6 1 (g, y)~ (A.y, g):;;;;
~ 6 1 (y, y), donde
..
u2
4
.
=To sen2
De aqu se dMprende
nh 1
-21 '
'
6 11 y11' ~(Ay, y)~ 6 11 !1 11'. (16)
donde
6 = 6, + 8i = ...!.. sen ~ + ...!.. een ~
"' u, "I "' ' (i7)
~
6 = 6 1 + 6 1 "" Tt cos1 2l. +-- cos1 21.,
4 n1a, '
"A = Tf
' ... 4 sen1 -U.
""" ,
v11,.,, (z 1) - / 2
T, sen -nlr,z,
-,- 1 , k, = 1, 2 , . . . . N 1- 1,
'... -
"" 4 sen 2
1a: u;-
""" ,
,., ( )
w11, :ri = V/ T n/r z 1
;-sen - 11- - , lci=1, 2, ... , N1 - 1,
1
(26)
b1., 1, - 0.
El operador z [y) puede anotarse de otra forma :
Z (y) == d1.1,V1,1, + b1,- 1, 1, (111,, 1, -Yi.- 1. 11 ) +
+ b1, +1.1, (111,, 1,-Y11 +1. 1,) + b1,. 1, -1(111 . 1, - 111. 1, - 1) +
+ b1,1 ,+1(ll1,1, -ll1,, 1, + 1). (30)
donde d1,1 , = a1,1,-b1,-1, 1, - b1,+1.1, -b1,. 1, - 1- b1,. 1,+ 1
2sn Cap. VI. BcD4clonu clfpllc1u
(35)
w1 = L..J
"" cr; seo --
lml , L= t, 2, ... , N-1. . (8)
-
Se ha elaborado uo algoritmo especial de transformacin
rpida de Fourier para calcular sumas, el cual permite obte-
ner la suma (8) en el transcurso de 5N log, N operaciones
aritmticas (cuando N = 2", siendo n un nmero entero) en
lugar de O (N') si se usa el modo de sumar habitual. Este al-
goritmo permite hallar la solucin del problema de partida
(2) en el transcurso de O (N1 N 1 log1 N,) operaciones. El
mtodo de separacin de vuiables puede combinarse con el
de reduccin o descomposicin que representa una modifi-
cacin del mtodo de Gauss. De resultas obtenemos un algo-
ritmo con el nmero do operaciones Q:::::: 5N1 N 1 log, N,
que es dos veces menor que el nmero correspondiente para
&l algorit mo de separacin aducido ms arriba.
2. Mtodos lteratlvoe. Los mtodos directos son ms
econmicos cuando se resuelve el problema de Dricblet en
diferencias para la ecuacin de Poissoo eo un rectngulo.
Actualmente existen programas estndar en el lenguaje aJgo-
ritmico FORTRAN y ALGOL para resolver las ecuaciones
de Poisson en un rectngulo con las condiciones de contorno
de tres tipos y tambin con las condiciones de contorno mix-
tas. No obstante, cuando el dominio no es un rectngulo o
cuando se analizan l11s ecuaciones de coef1fientes variables,
f 2. Re1olutl6n de /,u et11Attont1 en dl/erentlar "55
(t7}
" = if t ><a = if ,
(I 1 > ["';..., c1,+t, 1,i+x.i..., c11. 11 + t>+v.\" c11 , 11 1]
l/,l + t 1. I .. (f +)(1+X.J
(i9)
Para determinar yf,'' (11 , 11 ) elijamos un norlo 11 = t ,
t, = 1 en la eaquina izqu ierda del reeU.ngulo; entonces, los
dos nodos reatant.es (ti - t, 11 ) y (1 1 , 11 - t) del molde
{(11 , 11 ). (1 1 - 1, 11 ) ('1, 1,-,)} se disponen en la fronter11
y, por lo ta11to, y"'
('1 - 1, 11 ) = Y<>
('1, t, - 1) = O
son conocidos. Sabiendo vl"
para 11 = 1, 1, = 1, hallamos
sucesivamente 111" para '1 = 2, 3, .. ., N1 - 1 y 11 = 1
(en lo primera f~l a). Luego suponemos 1, ... 2 y encontramos
sucesivamente 111" en Ja segunda fila paro 11 = 1, 2, ...
. . ., N - 1. Con el fin de ballar VHi
realizamos los clcu-
los en el molde {(ti. 11), (11 + 1, 11), (11> 11 + 1)} segn
las columnas de arriba aba jo: fijamos 11 = N 1 - 1, N, -
- 2, . .. , 2, 1, 'J para cada i 1 cambiamos i 1 = N, - 1,
N, - 2, .. ., 2, 1. Empezamos la cuenta de Y+
con el
nodo (11 = N 1 - 1, t, = N, - 1) en Ja esquina derecha
superior. Se debe observar que la cuenta de Y+
puede rea-
lisaree tambi'n segn las mas de derecha a izquierda: lija-
mos 11 = N 1 - 1, N 1 - 2, .. ., 2, 1 y para cada 11 cam-
!
biamos 1 . . N, - 1, N1 - 2, .. ., 2, 1. Es mb, el cl-
culo de 111" podemos llevarlo a cabo no por las filas, sino por
las columnas de abajo arriba. Esto lo muestran las mismas
fnnulas.
Los clculos se realiian por las frmulas recurrentes (19);
la cuento es, evidentemente, estable. El algoritmo dol t ipo
semejante so llama (como ya se ha indicado) algoritmo tk
c6mpulo m6vll.
Calculemos el nmero de operaciones aritmticas que co-
rresponden a un nodo de lo red: el clculo do F,. requiere
2. R..olucl6n de ''" ec11aclone1 en dlferencl111 2.59
= ( 11 +11 ) t.-1
~ ~ Hv... > +<v...)1 l1,1,ll1h,~
1.-1
1
~
\ tl) (Ay, y),
(rt+7'f
puesto que (vluse ol ~ 1, cap. V)
t/tO 200 32 9
1/50 5000 t60 21
1/100 20000 320 29
B = (E + ioR1) (E + ioR,), R1 + R, = R , R~ ==
= R par a D = E.
1
11-1132 ,. _ 1/111
..!!.
D -11
1
D-d(x).S D-11 D.-d(l<).S
2 23 20 ~ 30
8 46 23 90 47
32 92 2S t80 53
t28 t84 26 360 57
512 387 28 720 ~9
Capitulo VJI
(6)
266 Cop. Y//. Bcuael611 ti. co11daclt6UUfad tlrmlca
donde A,. y X11. (z) son los valores propios y las funciones pro-
pias ortonormalizadaa del problema
X-+ AX= O, O< z < 1, X (O) = X (t ) =O,
iguales a
)." = k'n1 , X,. (z) = V2 son knz, (7)
con la particularidad de que
t
(X,., X,,.)= JX,. (z)X.,(z)dz =O"""
o
{ t, k=m,
0
,.,,.= O, k.pm.
Efectivamente, todas las soluciones porticulares (arm-
nicas) u11. (z, t) = c11.e->."'X,. (z) satisfacen la ecuacin y
las condiciones de contorno (5). De la condicin inicial
-
u (z, O)= u0 (z) = ~ c,.X,. (z)
.\-1
(8)
se balla11 los coeficientes e,.
= (u0 , X11) .
De (6) y (8) se infiere
11 u (1) 11' = (u (z, 1), u (z, 1)) =
puesto quo
lluoll = ~el,
11.-1
( 11u)
11.z 1
_
"1+i-2u1+"1- 1
/1
A
"xz, ( = U1,
1
Pera determina r y{+ de (13) obtenemos el problema de con-
torno
2.. yJ+ 1 - ( 1+
2
T ) yi+ 1 + h_ yi+ 1= - F~ o< i < N.
h' 1-1 "' ' 1-1 ..
de rontorno:
-'t-yl+l
2"' 1- 1
- (1 +~)
h'
yl+l
.J.' -:!,.__ yl+t =
2/t 1-1
-P{ 0< i < N,
g~+I = Ui (tt1), y~+I =~(t+ 1 ), (15)
Pl-(1-
1- - 't} Y1+u.1
J J + Y1+1
-.!._( Y1+1 J >+ cp,.
J
obienomos
Y1 = 6y< l + cp, (z1t) E wh,, y (z, O) - IL 0 (z),
Yo - 1 (t), Y N = u, (t) (t = t, = /-c. i - o. 1 . . . . ).
(17)
Sea 1t = u (z,, t 1) lo solucin exacta del problomo de parti-
Jn (3) y sea y la solucin del problema d i' ti iferencias (17).
270 Cap. VII. E cG4cl611 de conductibilidad llrmlca
= V ( %, tJ + { T) 1
- 1 av
.,;;
v = v- -z't 8I + --81'-
i: 'ti av + O ('t'),
Q>
48
,.. ou
/\u = u;;z = Lu + 12 u IV + O (h') , Lu = or ,
obtonomos
YJ+
1
1
-
( 1 - ~)
11 , YJ1 + 'f ( J
;T Y, _ , + Yl+
J ) + 't'l'
1
_, O< l<N,
(22)
y~+ - O, y~+t = 0, y~ - u0 (z,), O ~ l ~ N.
(25)
o bien
~yJ+-(t
,.. 1-1 1-~)yJ++~yi+=
,.. 1 ,.. 1+1 -P 1
1 O<i<N,
vt = 11~+ 1
=0.
Ahora, hagamos uso del teorema 3 del 5, cap. 1: para la
sol ucin del problema
A1!11-1 - C1111 + A1+11/1+1 - - F,,
C, = A1 + A1+ 1 + D,, O< i < N, lit = !IN = 0
es justa la e!timaein
11 //lle~ll ~lle
En nueetro CBllo A 1 = A 1+ 1 = -clh', D, = 1,
ll llM lle ~ 11 Fl lle ~ 11 Yl lle + T 11 Cl'l llr. (26)
De aqu, sumando segn k = O, t, . .. , J - t, obtenemos
In estimacin (25). De este modo, un esquema Implcito
puro es incondicionalmente esLeble, es decir, es estable para
cualesquiera T y h . Siendo a arbitrario, la e~ uaein en di-
f / . Bcaocl6n con co/lclcnlu corulon1.. 273
Aqu1 p = (1 - 1/2T6)/(1 +
1/2T6) = e-"'
(1 +O (T 1 )).
SI 111 condicin T ~To est perturbado, t111 decir, si
T > T 0 , entonces para 1 gran des predomina no 111 p rimera,
si no la ltima armnica:
11{ ~ c,p; sen n (N - 1) z 1 ~ c,pi ( - t)' sen nz1 ,
t/2tt. - t < e- ).1, 1o que, por su pues t o, no t 1ene
112, 6 + 1
d on d e p =
nadn en comn con la sol ucin de una ecuac in diferencial.
2. Problmu mullltllmen1lonale1 de eontlucllbllldad 211
con la frontera
Yh = {.z, = (i1h.,
l 1h 1 ): 11 =O, N 1 , O< i~ < N,:
t, = O, N ,, O< i 1 < N1 }.
Aproximemos el operador de Laplace Lu ~ Au mediante
un operador de diferencias en el mold e pentapuntual (vase
el cap. VI , 1)
Lu. - Au. = u.i,"" + u.;;,u.:,.
Sustituyamos el problema (1) - (3) por 011 problema dife-
rencial en diferencias (mtodo de las rectas):
du~t(t) = Av1 (t)+ f 1 {I) , i = (i 1, t 3), v 1 (0) = U0 (X),
o bien T
2
< ( 7if + hf
2 } -'
. (12)
~ - 0 (j/&fZ+i:1) + ( T - +) 0 (i:).
Para el error i = y - u tunemos un problema
donde
llzlf~ = lfall~. + /f.1111~ A111 = -11; 1,.1 , AsY "'" -11;1.,,,
~ N;I N~I
llzll~""' 2j ,lj h, (1;1 (i1t is))l + 2j X
i1-\ l.-t i1-i
>< 2 h,
1,...1
(z;, (ih la)).
De aqu proviene la estabilidad incondicional de la conve1-
goncia del esquema (5) en 11 A con la velocidad O (t + 1h I')
para a =/= 112, a ;;;;. t /2, y O (~ +
1h I') para o = 112.
La investigacin realizada ms arriba debe ser comple-
mentada con las condiciones de estabilidad asinttica.
Por cuanto las condiciones citadas -e ~ -c 0 se han obtenido
para un esquema operacional en diferencias con pesos y opo-
rndor ar bil.rario
A = A> 0,
pueden ser aplicadas, pue., tambin para nuestro esquema
(5). Haciendo uso de las expresiones (10) p11ra 6 0 y 6 0 , obte-
nemos las condiciones de estabilidad asinttica -e ~ -e~",
'"'~" = 2 ( rr4 + ;4 ) .1
para un esquema explicito (a = O)
-e~ -c~1 ', T~' = ./!-.-
,. oeAe
,6 0, 6 0 de (10), para un esquema
simtrico (a = 112).
En particular, para h1 = h1 = h, l1 == l 1 "" l tenemos
8 1 nh
u,= V een 21
'" nll ) - /al
o = TJa (
sen
1
,::::; 2n" .
2. Problemo1 m ullldlmendonolu de conducllblltdod 283
" ,v =;-
' [a, ce,, + 1)'' J.
'z'"Z; t)
".+1. ,,"- ..,,,,
11
-a (lh
' P
ih
1 2
t) ''-hs111 1- . J,
284 Cap. VII. Ecuac16r. de conducllbllldad trmica
Ay= - Ay=- (a yi 1 1
). 1 -(azYi,)x,.
Tenien do eo cuenta que para el caso unidimensional del
operador A: Ay =
-(ayi)"
e, (y, y)~ (Ay, y)~ e, (y, y), Ay = -u;x.
O < c 1 ~a~ c 1 ,
no es difcil convencerse de que las desigualdnde!! semejan~es
se verHican tambin parn el operador bidimensional (15):
o o o
c1 A~A~ CzA, Ay= - Y; 1x 1 - Y;;2 1,
De aqu se ve que 6E ~A ~ t::.E. 6 c,6 0 , 6 =
c,t::. 0 , =
donde 60 y 6 0 se determinan segn las frmulas (10). Para
hallar y
= yi+ :;obr" la capa nueva obtenemos el proble-
ma (6), en el que A se determine de (15). En el caso de un
esq uema explicito y
se determina en cada nodo x E w 1, por
la frmula
y = y + (1 - o) -cAy + i:tp .
Para los esquemas impllcitos (o =:/=O) se debe resolver una
ecuacin en diferencias pentapuntual con coefi c ientes varia -
.bles. Aqu se emplean los mtodos iterativos, de los cuales
el mtodo 11lternado t riongu lar resul la ser ms ecoum ico
(vase el cap. V, 5); el nmero de iteraciones pa ra dicho
mtodo se expresa por la magnitud O (V~ In~) , siemp1-e(1uc
O (h). La clescripcin del mtodo alte rna tl o lriangulnr
(i: :::::
para las ecuaciones en diferencias co11 coofii:io11tes variRble!I
se hu dado en el cap. V I; si se nplirn a In ecu a c in (6) ron
el operador A del lipo (15) , debo sur un lanlu mudificodo.
S. E1quoma econ6mlco1
3. Esquemas econmicos
t. Mtodo de direcciones variables. Al comparar los es-
quemas explicitos e implicitos (5), detendrmonos en dos
caracterstir.as: el volumen de Jos clculo!! para hallar yl+ 1
y las restricciones que se imponen sobre el paso "t.
ESQUBau sxPLfc1TO: para qeterminar yi+ 1 sobre la red roh
hay que realizar un nmero de operacioneii proporcional
al nmero de nodos, es decir, el nmero de operaciones co-
rrespondientes a un nodo no depende de la red oo,.. Sin em-
bargo, el paso 't est rigidamente limitedo superiormente por
la condicin 't ~ 'to (h): 't ~ h'/4 con h 1 = h, = h para
el esquema (i3).
BSQUEKA 1t.1PLtc1TO (a ;;;;. 1/2); para determinar yi+ se
rlebe resolver un sistema de (N1 - 1) (N 2 - 'I ) ecuaciones
en diferencias pentapuntuales; con este objeto, por lo menos
en el caso de coeficientes variables, se requiere un nmero
de operaciones (correspondientes a un nodo de Ja red c.>h)
que crece cuando 1h -O.
Surge un problema: construir los esquemas en que se com-
binen las mejores cualidades de esquemas explicitos e im-
plcitos, es decir, los esquemas a construir deben ser incondi-
cionalmente estables con un nmero de operaciones en cada
ca pa proporciona l al nmero de nodos de la red 11> 11 Los es-
quemas de este gnero suelen llamarse econmicos. Por su-
puesto, hemos de hacer una especificacin de que los esque-
mas incondiciona lmente estables en el sentido habitual han
e.le ser asintt icamente estables, lo que conduce a una restric-
cin del paso mucho ms dbil (por ejemplo, 't .;; lh/(2n)
para a = 1/2, h 1 = h 2 = h, l 1 = l, = l) que la condicin
de estabilidad ('t ~ h 2 /4) para el esquema explicito. Ade-
ms, la condicin 't = O (h) es natural para el esquema
+
o ('t2 1h 12).
Los primeros esquemas econmicos han aparecido en
:os aos 1955-'1956 y se han denominado mtodos de direc-
ciones variables. La idea algortmica principal (en la que
se radica In economia) co nsiste en que para pasar de uofl
capa t 1 a otra ca pfl l 1 + 1 se deben resolver, por el mtodo de
factor izacin, las ec uaciones en diferencias tri puntuales,
primero a lo largo de las filas y, a continuacin , n lo largo
de las columnas de la red 11>h .
286 Ca.p. VTI. Ecaacl6n ele conclucltblllclacl tlrmlca
= eJ+r+i - ~ A, (i - I).
Para ballar y1+1/I e yl+ 1 tenemos problemas de contorno en
diferencias
11-A 1yl sti-yisti= -Fi,
(8)
tiende a cero, cuando tien de e cero el paso T segn la varia-
ble t.
La idea del mtodo de aproximacin sumaria puede ex-
plicarse con un ejemplo del problema de Cauchy para la
ecuacin d iferencial ordinaria
du
y,+au=f (t), t>O. u (O)= u 0 , (9)
slH-siU/I J I
't +42Z +I _ _ ,.,., J - 0, f , .,
~=0, 1'1{- uJi;-uJ + a,ul112-f{ ,
,., =
11Jl- u J+l/t
't + zUJ - f l
Introduciendo aqul
'
obtenemos
iP. = (~2 + 41U - fi)i+l/S + 0 (T), (t5)
i'{ .. (,;2 + a,u - / 1)1+1/S +O (T).
De aqu[ se ve que i'{ = O (t), i'{ = O (t), sin embargo
+
""1l "1t = O (T) -+ O cuando T -+ O. (t6)
Todos los razonamientos aducidos ms arriba, a partir
de (tO), (U), (t4), quedan en vigor si "1 y a, representan las
matrices o los operadores, y u. /, 11 son los vectores.
De este modo, el esquema (t t ), (12) aproxima el problema
(9) en el sentido sumario (16) (tales esquemas se denominan
adittv1) .
Para demostrar le convergencia del esquema (t t ), (t2)
es menester obtener la estimacin pare el error zl~ =
= yl+ - ul+ en que se toma en consideracin la propie-
dad (16) de la aproximacin sumaria.
Pongamos
11>.. -= ~.. +11>;.
~.. =(~12+a,.u-j,.)' 1 12, ,,,;""O(T), c:.i=t, 2,
i11a = '11+1t' + ''tti zlt = '11+1 + 'lh
donde TlJ1> i 1., son las soluciones de los problemas
o
'11+112""' TIJ + Til>u '11+1""" '11+112 + "f111i,
J=o. t, ... TJo=O, (17)
(f + a,T) t1+112-i1 + tj;h (t + 42't) ')I == i1+1J1+ T~,
/=0, t, . . . (18)
~=0,
De (16) obtenemoa
1~J+l/2f<1~J1+-r1~1.
f eJ+ 1 f< 1~+ 1/2 f+ T f ~~ f < f ~J f + 'f ( 1~{ 1+ f ~ 1),
de modo que resulta lcita la estimacin
lsJ+I f<
4-1
T(f~f+I~), (21)
Lu=b.u - L 1u + L 1u,
ou
1.. u ... ---, a = 1, 2,
......
293
(26)
""<Jl L 1 v<11 + f 2,
--=- t,~ t~ti+I v{21 = v{t, 1 vJ+l..,. v(t/
De resultas obtenemos un esquema aditivo, un esquema
unidimensional local o bien un esquema de fisin
11l+llZ-11i
"'
11 J+t- 111+11a
"' = J\ 2 ( u~yJ+ 1 + (1. _ u2) yl+ 112) + q~ ,
zEwh> j=O, 1, ... , (27)
y0=u0 (z), zEw,.,
yJ+ s12 1v,.=J+112, yi+I lv,.=i+l.
las ecuaciones
J+l/2 .
T - ' = A t zJ+1/ 2 + 'I>{
,J+ 1 ,J+ 1/2
- "f =A2zJ+ I + 11>~.
donde u es la solucin del problema de partida (25), 'l>i y
'!>t son los re11iduos
1 u+;;. 1 ;. _ ,, 1 ;._u
"=A, -2- - 2 -T-+<i>1 ~ = A2u -2--r-+Cl'I,
l== u;t, u = ul.
De aqui se ve que 11>1 = O (1). 'i>s = O (t ), es decir, cada
una de las ecuaciones (27) no a proxima, tomada separada-
mente, la ecuacin (25). Tomemos la suma de residuos
u+ ;;. ;._u
'1>=11>1 +'1>z~ A 1- 2 -+ A2 u - --r- +<i>s+<i>z=
- a;
= (L1 + Lz)u-8t+cp1 +<i>i+ O ('f+ I h.I),
donde = ul+ 1ft . TomaDdo en consideracin la ecuacin
(25) para t = 11 11, , obtendremos
1 2
11> - <i>1+Cl'l- / 1+ '+ 0('f + 1hl 2) = 0('f+1 h1 1 )
1h JZ=h: +h:,
siempre que
4Jl1 + Cj)2= J+l/2 +0 ('fZ).
Esto puede ser conseguido suponiendo, por ejemplo,
cp1 - 0, Cl'l -f1+ 112 o bien qi1='Pz=/;/2.
Se puede mostrar que el esquema (27) converge unlforme-
mente con la velocidad
O ("'C + 1h I'), es deci r , ll 11J1 - ui1 lle = O (i: + 1hj').
De los ejemplos aducidos se ve que el mtodo de aproxi-
macin sumaria permite realiza r la particin de los proble-
mas complejos en unn sucesin de problemas ms sencillos
y aimpllficar considerablemente la resolucin de los pro-
blemas multidimensionalcs de lo fsica ma temtica.
Anexo
clones rKurrentea
ll1+1 = Cv1 - 111 - 1 - F., 1;;;.. t, lle - O, (2)
no aeri dlllcU notar que todas Ju incgnilu 111 pueden eer encontradas
suc:esivamente por la frmula (2), si calculamos v1 de tal o cual modo.
En este caao, cualqui er 111 ae oxpresari linealmente eo trminos de llo
e lli Todo lo dlc:ho nos permite eacribir para cualquier 1 ;;.. t una co-
rreJaci6o
111+1 - 1111 - P1 -1llt - PI (3)
con los eoellcleotee e&, p1, p 1 que por abora quedan Indeterminados.
SI pooemoa
e&0 - t, P-1 =- O, Po - O, (4)
eotoneea (3) ae varifiea tamblb iiara 1 - O. As puee, la aoluci6o del
problema ae buaeari en la forma (3) para cualquier 1 ;;. O.
Anotando (1) en forma de Ju relaelon.. rewn-enl419
111-1 - C111 -111+1 - P1, 1 ~ N - 1, llN - O (S)
y ruooaodo ao'1ogamente1 lleramoa a que la aolucln dal proble-
ma (1) para cualquer 1 ~/Y se puede buacar en la forma
1/1-1 - ~-1llN-1 - 111t-1 - 1l1N - 911 - 1 (6)
11 pooemo1
to .. t. 11- 1 - q, - o. (7)
Observemot que al VH-t queda determinado, todos l os v1 ae pueden
calcular aucealvamente aern la frmula (5).
Hallemot lit e 111'-t C01I este fin determinemos loa coeficientes,.
P1. ' ' ' 111 PI 91 Al comparar (2) Y (3) para 1 - 1, y (5) y (8), para
1 - N - t , obtendremoe
Hallemos ahora llL e l/N-i Con eate llD poqamoa en (13) 1 - Ji,
y 8D (14) 1 - k + 2. Tlllliendo prM8nte que N - 2" + t, teodremoa
Rutalldo la primere igualdad de la eegunda, obteDdremoe una ecuacl6D
rupecto de 111 llN-1:
1'-ll/1'-1 - "lit+ 11-lllt - 1'-tllN - 91'-1 - P (18)
Obtenramos una eouaol6n m para 111 e VN-a. aupoll.lendo 1 - k - l
en (13) y 1 - le + t en (14), y aualrayendo la ..,unda eouaol6n de la
primera
-czlll'l'-1 + CZ.,-11/1 - -U't + -1111' - P1'-t - '11' (19)
Teniendo prueole que ro .. IN - O, 8WD91DOI y aualtayamoe (18)
y (19). Obteodremot un 111q11ema equlvaleole
+ ru -
(a.,_, - ) 111'-&
(a.,_,+ a,.) N-t - 11.> - l- -- P -+ P-1 + '";
'f1'1 P1' Pla-& - '11' (20)
1-1
(22)
-IT
a11 +ca11_1 (c-2coe;~ 1 s).
1-1
Recurrlndo a (22) y (2S), collAruyamoa el atruleut.e alorl\mo para
ballar lla lltt-a:
" - P11 - 911-1 - P11-1 + q,., "' "" 9-a - P11 - P11-1 + 911
(c-2coa ci-;~J'I )
111- 111-1.
(c-icoa 211~t
(%3)
s) ...1-..,1-1o 1-t, 2, ... , lt,
111 - 0,S(v,.-w,.), llH-1 - 0,S (1111+"'>
Por cue.nlo cada uno de los sist.emu (23) cu.anta con una matrls tri-
diagonal (ol omero de lales alst.emu es 'lit) y puede ser resuelto por el
m6todo de ractoriuci6n roaliza.odo O IM) operaciones., entoncee para
ooconlrar 111 e11H- ae oecesltar6o O (llM) operaclooes arltm6tlcu.
A.al puea, para reaolver el sistema (t ) con l a matrls trlagooal ae
Ita construido uo m6todo en el que el nmero de operaciones arltm6-
llcu es proporcional al nmaro de ine6goitas.
Fij6moooa en que el algoritmo de marcha construido puede 98r
oumrleament.e inestable. En efecto, .al el omero C satisface la coo-
dicl6n 1 C 1> 2, entonces para ol algoritmo resul ta earact.eriateo
el crecimiento del error, expooeoclal aogo N, puesto que entre lu
ralees de la ecuaol6n ouacteriatloe q* - Cq +
t - O se tiene una que
en m6dulo es su~rior a la unidad. La ineet.ablUdad del mlano tipo
tiene lugar laml116n cuando la malrls C tiene valoree propios que aon
superion11 en mdulo a 2. Para los problemu da eata [ndole esti cons-
truida actualmente una variante del algoritmo de mucha estable en
aquel seD\ido que el error crece, al oroeer N, segn la ley potencial.
2. M6tedo de red uccl6n. En algunos casos, al resolver loa alst.emas
de ecuaciones algebraicas llnealas con una matrl1 tridiagooal, as de
mucha Importancia la exactitud de la aolucl6n obtenida. 8 1 aoilis
de las r6rmulu del m6todo de ractori11cl6n que ae emplea para re-
aolver los elstemas citados muestra que la fuente del error puede ra
dicar en las 6rmulu pua calcular los eoeflclentes da factoriaaci6n.
Estas frmuJu contienen la operacii6n de dlvlsl6n por una diferencia de
las mago ilude.a que son prximas en au valor. Mb abajo se dar' a cono-
cer al m6todo de redueca6n para resolver dichos sistemas, privado de la
deficiencia meociooad1.
A.al puea, supongamos que ae requiere ballar la solucin de un
problema de diferenelaa trlpuotual
-01v1- + c1111 - b 1111+a - J,. t < 1< N - t, 24)
' - VN - O.
Cr. 11) -
-,_,
~
N
BB >I 514&
H3AATBJILCTBO M.HP
129820, Moon . U- 11 0, rcn, t-1 PK>t<cKU nep., 2.
0pAea TpyAo-oro Kpacoro 3aU1ewt
M OC-KOacK&.A T8norp " .>t 1 B ctcpa peeonlOtlM
COt01nomrrpanpoM& l'ocyaapcnenooro koMM'l'O'H CCCP
DO tlt.'1IJll IAITUllCT'I. non:arpat101 ktiJUtto f'Oproan"
103001, Moc...a , TpUJIPJAtud nep., .
1702070000-024
e OOIJ(Ol) -86 146-86, ... 2
GolovJn L.
ALGEBRA LINEAL
Y SUS APLICACIONES
(3 edicin)
No obstante su pequelio volumen, el lib10 ()()atiene Jos pro.bJemu
fuodainentales del eurao de lgebra lineal, aal como sua dlaUntaa
aplicaciones, incluyendo la investigacin de las ourvu y llUpe!'-
flcies de segundo orden, la nocin sobre tenaores y ooa pro-
blemu.
En el libro se a:ponen l" ccnceptos primordllllee refezentee a lOll mpa-
cios lineales y euclidianos, y a l as transformaciones llnealee; adem
se estudian problemas sobre vectores y ae obtleno la forma cannica
de la.a matrices de las transformaciones autoconjugada y otrogooal en
el espacio euclidano, dndose ejemplos hlcoa de la teoria de las
form113 cuadrticas.
Como suplemento a las formu cuadrticaa, ae examina la teoria
general de las Hneu y auperficies de segundo orden. Dos capltulos
estn consagrados a las trailaformaciooes de Lorentz y nociones funda-
mentales de la teorla especial de la relatividad . En el capit ulo sobre
grupoa, adems de las de(iniciooes principales, se incluye una seleccin
de ejem>los.
Un m6r1to evidente del libro es la acertada elec<>l6n del mal.erial,
en el cual ae han examinado problemas que no entran en el programa
de los estudiantes de especialidades no matemiticas, pero que son de
cierto inters para stos. Con ello, las nociones iodispeoaables previas
y el nivel de la exposicin son tales que, al leer el texto, los estudian-
tes no encuentran ninguoa dificul tad.
La obra est6 destioada a estudiantes y profesoras de centros de enseau-
za superior. T~m bin les es de gran utilidad a los ingenieros que
deseen conocer las nociones lundamenl.ales del lgebra lineal, emple-
ando una (ueo~c que no exige inonnacio previa de matemticas
superiores.
Kag6n V.
LOBACHEVSKJ
Eo este libro ae oarra de la vida y de lu acUvidades eoclales y clentlfi-
cu del emlneo\e matemitlco ruao N. l. LobacbevslO.
Se relata en detalle sobre la nlJlez y adol"5(:encia del K' sabio, sobre
su famllia, profesores del liceo que influyeron en la formaci6n del
car,cler y U. concepc6n del mundo del cieolifico. El autor pres-
t.a especial atenci6o a loa alloa estudiantiles de Lobachevaki en
la Uo.ivereidad de Kuio, y en parUcular, a sus estudios de las male-
m.hicu.
Gran parte de la obra es" dedicada a la actividad clenUfica de Loba-
cbevekl . El \ec\or conoceri el manual de estudio \ilulado "Geometrla",
donde por pri111era vei la geometra abeolut.a fue separada de la eucli-
diana. S4l trata en delalle la creacin de la geometrla no euclidiana,
ae a oallu el trabajo "Geometrla, investigacin de la uoria de laa
llneu paralelu". lfn apartado importante del libro ae destina a loe
aos ms fructlferoa de la obra de Lobachcvski (t827-t846, eo. ~ue
fuera reot.or de la Uolveraldad de Kaztn . .En este periodo public
sua trabajos ms no la.bles: "Sobre los principios de la geometra",
"Nuevos rinclpioa de la P'!llletrla con la ~ria completa de las
paralew , e~ .. que son aoali1adoe. El siguiente apartado ha sido
dedicado a la relacin de eus contemporneos hacia lea ideaa de
Lobacbeveki, al desarrollo de estas ideaa en los afios ulteriores, a la
aplicaci6n de la geometrla de Lobachevald a olraa partes d.e lu matem-
tlcu, uf como a la meoica, tsica y cosmologa. _
Eele libro reeultari ein duda de Inters a 4!jltudiantes, maestros, prc>-
fe90rea y a toda persona aficionada a las roalem&ticu.