Sin 3 RG 3 Tica
Sin 3 RG 3 Tica
Sin 3 RG 3 Tica
G. G. Malinietski
Caos,
uras y sirnu
por ordenador
URSS
M a f l u ~ e y ~ ureopzu
fi
re~~adbe6uv
M~T~M~TAW
OCHOBbI
C K ICHHepreTWKH:
~~
XaOC, CTpyKTypbI, B ~ I . ~ A c J I A T ~ J I ~ H ~ I ~
JKCII~~HM~HT
Cepea
.
C u ~ e p r e ~ iOT
~ ~npomnoro
a:
K Gyj~y~.qe~y
Malinietski Gueorgui Guendievich
Fundamentos matemticos de la sinergtica: caos, estructuras y simulacin por ordenador. Serie asinergtica: del pasado al futuro
Editor
Director financiero
Director de sistemas
Director de produccin
Revisin cientfico-tcnica
Traduccin
L ~ ~ ~ ~ ~ T ~( (JEI, ~
q ~c~To B
p~
YPCC.,
Qa n
117312, r. M O C K Bnp-T
~ , ~ Q - ~ ~ T OKTR~PR,
MH
9.
J l ~ q e ~ s i / l ~ No 05175 OT 25.06.2001 r. n o a n ~ c a ~K oneqarM 11.01.2005 r.
Q > o p ~ 60
a ~x 90/16. ner. n, 20.
O ~ n e r a ~ aB ~000
o
3noxau. 107370, r. Moc~sa,6-p Mapwana P o ~ o c c o ~ c ~ 5,
o rc"~~p~1.
.~.
Ma
Editorial URS S
Libros de ciencia
TelJfax: 7 (095) 135-42-16
Tel./fax: 7 (095) 135-42-46
E-rnail: URSS@URSS.m
Catlogo en Internet: http:l/URSS.ni
Reservados todos los derechos en lengua espaola en todos los pases del mundo. Quedan
rigurosamente prohibidas, sin la autorizacin escrita de los titulares del "Copyright", bajo
las sanciones establecidas en las leyes, la reproduccin total o parcial de esta obra por
cualquier medio o procedimiento (comprendidos la reprografa, el tratamiento inforrntico
y la exposicin en Internet) y la distribucin de ejemplares de ella mediante alquiler o
prstamo pblico.
ndice
Introduccin
.......................................
Captulo 1
La simulacin matemtica actualy los fenmenos no-lineales
.......
Experimento computacional. La terna modelo-algoritmo-programa . . . . .
Jerarqua de los modelos simplificados . . . . . . . . . . . . . . . . . . . . . .
Simulacin de La dinmica de poblaciones . . . . . . . . . . . . . . . . . . . .
Bibliografa . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Captulo 2
Modelos matemticos lineales
............................
Preguntas y ejercicios . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Bibliugrafia . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Captulo 3
Modelos matemticos no-lineales elementales
..................
E l pndulo matemtico . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
E l programa de investigacin de Henri Poincar . . . . . . . . . . . . . . . . .
Preguntas y ejercicios . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Bibliografa . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
.
Captulo 4
Atractores de la ecuacin 5 = v(x) . . . . . . . . . . . . . . . . . . . . . . . .
Teora cualitativa de \a ecuacin k = v(x) . . . . . . . . . . . . . . . . . . . .
Teoremas de comparacin . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Problemas bien planteados y modelos de fenmenos no-lineales . . . . . .
Preguntas y ejercicios . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
B i b l i o g r a ~.~. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Capitula 5
Elementos de la t e o h de las bifurcaciones
....................
Desarrollo de la teora de las bifurcaciones . . . . . . . . . . . . . . . . . . . .
Preguntas y ejercicios . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Bibliografa . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
fndice
Captulo 6
Modelos matemticos de la teora de las catstrofes
. . . . . . . . . . . . . . 128
Estabilidad estructural e ideas de la teora de las catstrofes . . . . . . . . 144
Preguntas y ejercicios . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 172
Bibliografa . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 174
Captulo 7
Sistemas elementales con tiempo discreto . . . . . . . . . . . . . . . . . . . . . 176
Captulo 8
A utooscilaciones y ciclos lmites
...........................
Perturbaciones singulares. Patos. Sistemas rgidos . . . . . . . . . . . . . .
Preguntas y ejercicios . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Bibliografa . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
180
188
190
193
200
202
204
231
242
243
Captulo 9
Mtodos topolgicos en la investigacin de los sistemas no-lineales
. . . 245
Experimento computacional diseo molecular y mtodos topolgicos . . . 265
Preguntasy ejercicios . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 269
Bibliografa . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 271
Captulo 10
Redes neuronales
.....................................
.......................................
Neurociencia
Nociones elementales sobre e l funcionamiento del cerebro . . . . . . . . . .
Modelo de Hopfield . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
E l sentido del caos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Redes neuronales multicapas . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Algoritmo de retropropagacin . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Preguntas y ejercicios . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Bibliografia . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
273
273
276
279
288
295
305
309
310
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 311
ndice de materias . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 313
ndice de autores
Introduccin
El obstculo principal para la amplia utilizacin de la simulacin matemtica
y del experimento computacional en la ciencia, la tcnica y el control, es
el dficit de especialistas calificados. La resolucin de toda una serie de
problemas importantes sewe frenada no por la falta de computadoras, sino
por la insuficiencia de colectivos con un nivel de preparacin adecuado.
Las exigencias impuestas a los especialistas en simulacin matemtica
son bastante elevadas y, al mismo tiempo, contradictorias. Por una parte,
se requiere que sean profesionales que comprendan profundamente un
rea de investigacin concreta y suficientemente estrecha. Por otra, estos
especialistas generalmente tienen que intervenir no como ejecutores, sino
como cientficos que ven el problema como un todo y estn capacitados
para precisar y, a veces, cambiar radicalmente el planteamiento del problema
propuesto por fsicos, qumicos o bilogos. El trabajo en el campo de la
simulacin matemtica presupone un estilo de pensamiento peculiar, en
el cual la profundidad y la concrecin se combinan con la amplitud y la
comprensin de las ideas generales.
El xito en la simulacin matemtica de la mayora de los problemas
serios se apoya en la terna modelo-algoritmo-programa, como en las tres
ballenas mitolgicas de los griegos.
Por eso, para trazar correctamente la estrategia de investigacin es
necesario comprender bien el instrumental disponible y los logros fundamentales en cada uno de los campos relacionados con el proyecto de
investigacin. Esto es muy importante, pues, por ejemplo, hoy la matemtica
computacional interviene cada vez con mayor frecuencia no slo como un
instrumento, sino tambin como una fuente de nuevas ideas en la simulacin
y las ciencias naturales. En este libro se destaca constantemente este aspecto.
El anlisis de fenmenos no-lineales constituye uno de los momento
claves en la resolucin de muchos problemas importantes, desde el perfeccionamiento de las tecnologas qumicas y los problemas de la fusin
termonuclear controlada, hasta los problemas actuales de la teora cuntica
de campos y la creacin de nuevas generaciones de computadoras. En el
vasto campo llamado actualmente ciencia no-lineal o dinmica no-lineal han
surgido sus propios enfoques originales, ideas y mtodos nuevos. El objetivo
de este libro es presentar esta disciplina a los estudiantes de los cursos
superiores, de postgrado y a los especialistas interesados.
Introduccin
Introduccin
El libro contiene muchos ejemplos e ilustraciones. Donde ha sido posible, se ha prestado especial atencin a las ideas claves y modelos concretos,
para cuyo estudio se crearon unos u otros instrumentos. En muchos casos
hubo que sacrificar generalizacin, rigor y detalles tcnicos en favor de la
claridad.
El libro est destinado a los lectores interesados y activos. Por consiguiente, en cada captulo se proponen varias fuentes bibliogrficas, que
van desde los textos de estudio estndares hasta monografias, resmenes
y artculos originales, los cuales ayudarn a conocer ms profundamente
el problema analizado. La bibliografa no pretende ser completa, se dan
las fuentes ms simples y accesibles. Para comodidad de los lectores la
bibliografa va acompaada de comentarios que muestran qu se puede
hallar y dnde.
Es muy importante que el lector no slo sepa mucho, sino que tambin
sea capaz de resolver problemas elementales relacionados con el anlisis
de los fenmenos no-lineales. Por eso, en el libro se propone cerca de una
centena de problemas utilizados en las clases.
El volumen del libro corresponde aproximadamente a un curso introductorio semestral acompaado de seminarios, dirigido a estudiantes
de cursos superiores de matemtica pura, matemtica aplicada y fsica.
Estos estudiantes ya han pasado los cursos de anlisis matemtico, lgebra
lineal, fsica matemtica, mtodos numricos y fsica terica. Sin embargo,
por el curso de simulacin de fenmenos no-lineales tambin se interesan
los qumicos, bilogos e ingenieros, quienes no poseen una preparacin
matemtica tan slida. Sus intereses, en cierta medida, tambin han sido
considerados. En el libro se enfatizan ciertos resultados fundamentales del
curso de matemtica, en los cuales se apoya la exposicin.
Solamente resta decir que el material propuesto refleja una pequea
parte de la ciencia no-lineal, relacionada estrechamente con el experimento
computacional, y representa apenas un fragmento de un campo vasto de
investigaciones. La eleccin del material est determinada, en parte, por las
pasiones subjetivas del autor.
El desarrollo de la ciencia del siglo XX mostr la necesidad de la
construccin de muchos modelos distintos para la descripcin de un mismo
fenmeno u objeto, la necesidad de elaborar cuadros alternativos de la
realidad. Nosotros estamos obligados a vivir no en un mundo de leyes
absolutas, verdades en ltima instancia y concepciones universales. Por eso,
los investigadores deben tratar con modelos y enfoques que fijan una cosa
e ignoran muchas otras. Nosotros debemos jugar a crear mundos en los
cuales se toma muy poco de nuestro presentecomplejo y confuso.
La avalancha de informacin provocada por las computadoras nos
ha obligado a reflexionar sobre cmo entenderla, aceptarla, ordenarla
Introduccin
CAPITULO 1
nos proporcionar el mtodo para la comprensin del contenido cualitativo de las ecuaciones. Hoy todava no somos
capaces de ello.
R. Feynman
Sin temor a exagerar, podemos afirmar que actualmente la simulacin matemtica (CM) influye en toda la actividad humana: desde la poltica estatal,
la cual se apoya en los modelos matemticos de la economa y la estabilidad
estratgica, hasta los nuevos medicamentos proyectados mediante computadoras; desde las nuevas mercancas que utilizan toda una serie de tecnologas
originales, hasta los gigantescos sistemas computarizados que venden billetes areos o aseguran las comunicaciones o la edicin de peridicos.
La construccin y el anlisis de los modelos matemticos de los
fenmenos no-lineales desempean un papel cada vez ms importante en
la esfera de la actividad relacionada con la CM. Detengmonos en las causas
de este hecho.
En los aos 50 del siglo XX tuvo inicio una revolucin cientfica, la cual
cambi radicalmente el papel desempeado en la ciencia y la tecnologa
por la matemtica en general, y la matemtica computacional en particular.
Este fenmeno estuvo relacionado con la realizacin de dos proyectos
colosales: la creacin de armas nucleares y la conquista del cosmos. Ambos
proyectos exigan gastos enormes de recursos y la eleccin escrupulosa de
las mejores variantes. La complejidad de los problemas que surgan los haca
inaccesibles a los procedimientos estndares de la fsica terica. La necesidad
de resolver estos problemas estimulaba la creacin de computadoras y de
nuevos campos de la matemtica. En ambos casos se formaron colectivos
cientficos, cuya experiencia permiti obtener resultados notables en una
serie de campos de la ciencia. En Rusia puede servir de ejemplo el Instituto
de Matemtica Aplicada M. V. Kldysh, y en los EE.UU., los centros de
investigacin de Los lamos y Livermore.
10
Captulo 1
E l problema energtico y la fusin termonuclear controlada. La prosperidad de la humanidad est estrechamente ligada con su potencial energtico.
Por ejemplo, en los pases desarrollados los gastos de energa para la obtencin de una unidad convencional de produccin agrcola han aumentado ms
de 100 veces en los ltimos 50 aos, mientras que la productividad de las cosechas de cereales durante ese mismo perodo ha aumentado slo tres veces.
Actualmente en el mundo viven cerca de seis mil millones de personas
y el consumo medio de potencia per cpita constituye un poco ms de 2 kW.
El pronstico demogrfico de la ONU predice un rpido crecimiento de la
poblacin del globo terrqueo en los prximos 60-70 aos, y la estabilizacii~
hacia el 2100 en el nivel de los 12 mil millones de personas. La utilizacin
de los desechos, el enriquecimiento de los minerales, la desalinizacin del
agua y las escalas crecientes de la proteccin del medio ambiente exigen el
aumento de la potencia de consumo per cpita hasta unos 10-20 kW.
El agotamiento de las reservas de portadores de energa y los grandes
problemas relacionados con la energa atmica exigen la bsqueda de
fuentes de energa alternativas.
Una de las posibilidades est relacionada con la utilizacin de la energa
de las reacciones termonucleares. Existen tres istopos de hidrgeno, los
cuales se diferencian por el nmero de neutrones: el hidrgeno habitual H:,
el deuterio D: y el tritio T;. Es posible la reaccin
neutrn
,
n y se liberan
durante la cual se produce un tomo de helio ~ eun ~
17,6 MeV de energa. Recordemos que un megaelectronvoltio es igual a lo6
electronvoltios, y un electronvoltio es la energa que adquiere un electrn
al desplazarse en un campo elctrico en el vaco entre dos puntos cuya
diferencia de potencial es 1V. Un electronvoltio es igual a 1,60207
J.
A modo de comparacin, podemos decir que un gramo de una mezcla
deuterio-tritio puede producir una energa superior a la de una tonelada
de combustible convencional.
11
Sin embargo, para que tal reaccin tenga lugar hay que acercar los
ncleos hasta una distancia de 10-13 cm, donde comienzan a actuar las
fuerzas nucleares.
Para esto se necesita superar las fuerzas de repulsin electrosttica, lo
cual exige temperaturas gigantescas, superiores a los 100 millones de grados.
El fsico estadounidense J. Lawson obtuvo un criterio que indica la
regin de los parmetros en la que los gastos energticos para la creacin de
las condiciones para una reaccin termonuclear son menores que la energa
generada por la reaccin misma:
Captulo 1
a la explosin de la bomba atmica se logra obtener valores de la concentracin muy grandes antes de que los productos de la explosin se dispersen.
El papel de carga atmica, la cual comprime y calienta el combustible
termonuclear, es desempeado en este caso por lseres de gran potencia.
El anlisis y la realizacin de este proyecto ya ha exigido durante muchos
aos el esfuerzo colectivo de los fsicos y los especialistas en simulacin
matemtica.
El experimento computacional desempea un papel clave en Ia resolucin no slo de muchos problemas aplicados, sino tambin de problemas
fundamentales. La astrofsica moderna contiene toda una serie de problemas
de esta ndole.
13
Experimento computacional.
La terna modelo-algoritmo-programa
La estructura de la nueva tecnologa de investigacin cientfica, el experimento computacional, se puede representar esquemticamente. A partir del
planteamiento del problema y de las condiciones impuestas a la precisin
de las predicciones, se construye el modelo matemtico. Partiendo de las
particularidades ms importantes del problema analizado, se elaboran los
algoritmos en la clase concreta dada de sistemas.
Despus de estimar cun amplio ser el crculo de investigadores que
estudiarn el problema y qu computadoras estn disponibles, se elabora
el software. El software puede consistir en programas nicos o programas
simples para computadoras personales, paquetes de programas aplicados,
si se trata de clculos masivos, o bien sistemas que incluyen la creacin
de lenguajes algortmicos especializados. Despus de esto se realizan los
clculos. Los resultados se comparan con los datos de las observaciones o
investigaciones naturales. Luego, habitualmente, se regresa al primer paso,
a la construccin de los modelos. Dependiendo de los resultados de la
comparacin de las predicciones de la teora con la realidad, el modelo se
modifica, se precisa o se construye nuevamente.
En la actualidad, un eslabn clave que determina la efectividad de
la utilizacin de las tcnicas cornputacionales en la ciencia y la tecnologa
lo constituyen los modelos matemticos. Esto se ha hecho especialmente
evidente despus de la aparicin de una cantidad enorme de computadoras personales, y durante la creacin de proyectos de cupercomputadoras.
En el primer caso, la cantidad de computadoras personales que se utilizan
para resolver problemas cientficos, en comparacin con las utilizadas para
14
.-
Captulo 1
15
Etapa moderna. En los problemas tecnolgicos surge la necesidad del clculo de una cantidad enorme de construcciones tridimensionales complejas.
Aparicin de computadoras con velocidades del orden de l o l l - 1 0 ~ ~operaciones por segundo accesibles a cientficos e ingenieros. Esquemas en
diferencias adaptados a la solucin. Renacimiento de los mtodos directos
de resolucin de ecuaciones lineales con matrices ralas. Aparicin de una
generacin nueva de algoritmos destinados a la realizacin de clculos
paralelos mediante sistemas multiprocesador y transputadoras. Creacin de
paquetes de programas aplicados, accesibles mediante redes computacionales a un crculo amplio de ingenieros. Posibilidad de realizar clculos de
construcciones nicas teniendo en cuenta un conjunto de procesos fsicos
(entrada del Shuttlea la atmsfera, clculo calorljFico de los motores de los
cohetes, etctera). Grfica computacional en colores; creacin de pelculas
en computadoras, las cuales permiten seguir la dinmica de un proceso.
Futuro cercano. Posibilidad de construir sistemas computacionales ~intelectuales., de realizar clculos calorficos que, a partir del problema planteado
y de los clculos anteriores, elijan el modelo matemtico y el algoritmo
numrico. Utilizacin de sistemas computacionales con alto grado de paralelismo como, por ejemplo, autmatas celulares que simulan los procesos
calorficos, as como sistemas neuronales capaces de sacar conclusiones a partir de la experiencia almacenada. Posibilidad de encargar a la computadora
el problema de la representacin adecuada de los resultados de los clculos.
Realizacin, en una serie de casos, de los clculos calorficos en tiempo real,
con el fin de optimizar los procesos tecnolgicos. Grandes posibilidades de
resolucin de problemas inversos, que permiten sintetizar campos calorrjTcos
con la configuracin deseada.
16
Captulo 1
Por eso, se debe proceder de otro modo. Primero se eligen los procesos
claves que desempean el papel principal en el fenmeno estudiado y en
las escalas espacio-temporales dadas. Despus se construye un modelo ms
simple, el cual tiene una regin menor de aplicabilidad y considera un
nmero menor de factores. Este procedimiento contina hasta que no se
obtenga el modelo elemental, cuyo comportamiento ya es comprensible.
Slo despus de estudiar y comprender el modelo del nivel inferior, es que
se puede pasar al nivel inmediatamente superior.
Se puede decir que el logro principal y el objetivo fundamental de
las investigaciones durante la resolucin de los problemas complejos es la
construccin de la jerarqua de modelos simplificados. En particular, se debe
establecer cul nivel del modelo es se debe utilizar en uno un otro caso. Por
ahora, casi todas las jerarquas construidas estn relacionadas con sistemas fsicos. Se estn construyendo jerarquas en una serie de campos de la qumica
y la economa matemtica. Este problema se plantea tambin en la biologa.
Un rasgo notable de la jerarqua de modelos simplificados es la
existencia de modelos matemticos bsicos, es decir, de objetos matemticos
cuya investigacin permite construir de manera efectiva y estudiar grandes
clases de modelos diferentes. Se puede decir que, a menudo, los modelos
bsicos se parecen a los cubos de juguete: con su ayuda se construye
la descripcin de los fenmenos concretos. Es importante subrayar dos
hechos fundamentales aclarados en los ltimos veinte aos. El primero
es que existen pocos modelos matemticos bsicos. Se pueden construir
modelos matemticos no-lineales extremadamente simples, pero profundos
y con un gran contenido. El segundo hecho es que con la ayuda de stos
modelos fue posible predecir fenmenos de la naturaleza decconocidos, sin
pasar por todos los escalones de la jerarqua relacionados con el
tener
detallamiento y la complejizacin de la descripcin matemtica.
La profundidad y universalidad de los modelos matemticos no-lineales
ms simples se han convertido en la base para la creacin de enfoques
interdisciplinarios.
17
.-
O2
O2
+
-, describen la movilidad o capacidad de migracin de
az2 a g 2
Laplace -
Captulo 1
18
...
..
P
P
19
--
Tanto este modelo como muchos otros objetos de la misma clase permiten
un estudio analtico detallado. Ellos sern tratados en los captulos 4, 5 y 6. .
La ecuacin (1.5) no puede describir, por ejemplo, las oscilaciones
peridicas del tamao poblacional. Si tales oscilaciones tienen lugar y son
fundamentales en el anlisis, entonces es necesario volverse atrs y precisar
el modelo. La ecuacin (1.5) puede describir cmo con el transcurso del
tiempo el tamao poblacional alcanza un valor constante, determinado por
los recursos de la regin. Sin embargo, puede ser que nos interese cmo
crece la poblacin de una especie que ocupa un nuevo nicho ecolgico,
Captulo 1
20
-
= au.
(1.6)
21
.las aplicaciones unidimensionales -analizadas en este libro- al describir
diversos fenmenos en la economa, la medicina, y al analizar las catstrofes
naturales y tecnognicas.
Adems, los modelos no-lineales conocidos surgidos en un campo,
a veces se pueden utilizar como una especie de bloques (((cubosde juguete.)
en otras disciplinas.
Veamos un ejemplo tpico. Uno de los logros experimentales ms
importantes en la ciencia del siglo XX fue el descubrimiento de B. P. Belosov
de las reacciones qumicas oscilantes. Esto trajo consigo la construccin de
los modelos matemticos correspondientes, cuyo anlisis condujo a la
hiptesis del caos qumico, es decir, de las reacciones qumicas oscilantes
con comportamiento catico. El. profesor Otto Rossler propuso modelos
hipotticos de reacciones en los que es posible el caos.
Ms tarde se demostr que, con una pequea modificacin, estos
modelos permiten describir las epidemias de una serie de enfermedades.
Y hace poco tiempo A. Yu. Andryev y M. 1. Levandovski propusieron utilizar
un sistema similar para describir los movimientos huelgusticos. El modelo
propuesto tiene la forma
La simulacin matemtica actual y los fenmenos no-lineales
x = m ( ~ x- ) - ~ X Z ,
22
Captulo 1
23
'
25
r En efecto,
veamos con los ojos de un matemtico la ley de Ohm, la cual
(-i,n,!+$~
relaciona la fension '"V aplicada a un conductor, la corriente I que circula
por ste y su resistencia R. Es natural suponer que exista cierta dependencia
I = f (V), f (O) = O. La ltima igualdad es evidente: si no se aplica tensin,
:B
Cuando AV
dI(0)
-+ O se puede dejar el primer termino AV, y llamar al
dV
1
dw)
nmero - conductividad - y al coeficiente R, resistencia. En qub
dV
R'
radica, entonces, la esencia de la ley de Ohm, si todo lo qiie hemos hecho
es vlido para cualquier funcin diferenciable f (V)? La esencia de esta ley
radica en que para un intervalo muy grande AV y para muchos conductores, el
dI(0)
trmino -AV es mucho mayor que el resto de los trminos de la serie.
dV
o Otro ejemplo es la ecuacin de Cchrodinger, la cual determina la
variacin de la densidad de probabilidad de hallar una partcula en el
tienipo y el espacio. En el mundo de los procesos atmicos, en algunos
casos es cmodo considerar la luz como un flujo de partculas (cuantos) y
en otros casos, como una onda electromagntica. Como dicen los fsicos,
tiene lugar la dualidad onda-partcula)). En 1924, Louis De Broglie enunci
u i ~ ahiptesis que permite considerar en el marco de un enfoque nico las
propiedades ondulatorias y corpusculares de los objetos del micromundo.
A finales del siglo XIX, Hamilton prest atencin a la analoga existente
entre la ptica geomtrica y la mecnica clsica. Resulta que las leyes
fundamentales de estas teoras se pueden representar matemticamente
de formas idnticas, con la particularidad de que al movimiento de una
partcula material en un campo V(x, y, t.) le corresponde la propagacin
de un rayo de luz en un medio no-homogneo con cierto ndice de
refraccin p(x, y, t.). Sin embargo, cuando la longitud de la onda d e luz no
26
Captulo 2
d
Aqu el valor propio del operador lineal ih- es la energa E, y el del
dt
li d
operador 7 - es la cantidad de movimiento p. Como funcin propia
z dx
interviene en ambos casos la funcin de onda @(x,t). Los creadores de
la mecnica cuntica supusieron que en el caso general se puede hacer lo
mismo: asignarle a las magnitudes fsicas observables los valores propios
de operadores lineales.
Consideremos una partcula material de masa m, energa E y cantidad
de movimiento p que se mueve en un campo con potencial V(X, y, Z) que
depende de las coordenadas espaciales X, y, z . La energa potencial de la
partcula es igual a V(x, y, z ) y su energa cintica es igual a p2/2m.
De acuerdo con la ley de conservacin de la energa
+ +p i
27
a@= E$.
ih-
dt
La solucin de esta ecuacin lineal es la funcin exponencial
Captulo 2
28
Los valores propios En determinan cules son los valores de la energa que
puede tomar una partcula que se mueve en este potencial. Las funciones
propias p, correspondientes indican la probabilidad de descubrir una
partcula en diferentes puntos del espacio. Por cuanto
2. Para una serie de problemas en los que no se requiere una precisin muy
elevada, o bien Las acciones sobre el sistema estudiado son pequeas,
29
Esta ecuacin simula no slo ondas que se propagan por una cuerda, sino
tambin la propagacin del sonido en un medio (aproximacin acstica).
En particular, de esta ecuacin se deduce que, independientemente de cun
grande sea una perturbacin, sta se propaga con velocidad constante c. Evidentemente, esto es as slo si las perturbaciones son pequeas. Para analizar una
perturbacin de amplitud suficientemente grande, en el primer caso se debe considerar el estiramiento y la variacin de la tensin de la cuerda. Esto conduce a la
aparicin de trminos del tipo (u,)*.Para el caso de las ondas de sonido de amplitud
grande, es conocido el efecto fsico de surgimiento de otidas de clzoque. La velocidad
de la ondas de choque depende ya de la intensidad de la perturbacin.
Si la accin sobre el medio es grande, entonces hay que renunciar al modelo
lineal (a las ecuaciones de la acstica), y considerar un objeto ms complejo: las
ecuaciones de la dinmica de gases.
donde A es una matriz y 5,b son vectores. Elijamos una nueva bace
Z = Sf, b = SE En esta bace,
ASij = SC
tal que
es decir, del mismo tipo que la ecuacin az = b, la cual sabemos resolver. As pues,
un problema complejo se desintegr6 en una serie de problemas simples.
Ejemplo 2. Ecuaciones diferenciales ordinanbs. Consideremos dos modelos simples,
pero muy tiles, relacionados con las ecuaciones diferenciales ordinarias.
Captulo 2
Modelo matemtico de desintegracin radiactiva. Un conjunto de experimentos relacionados con el estudio de la desintegracin radiactiva permiti establecer
lo siguiente. El cociente del nmero de tomos que se desintegran por unidad de
tiempo y el nmero total de tomos, es una magnitud constante que depende slo
del tipo de tomo. Este cociente 7 se denomina probabilidad de desintegracin. Aqu
se supone que el nmero total de tomos es muy grande.
Denotemos mediante N (t) el nmero de tomos que an no se han desintegrado
en el instante t . En el instante t dt este nmero es igual a N(t dt). As, durante
el tiempo dt se desintegran N() - N(t dt) tomos y por tanto,
31
coinciden con el modelo de crecimiento de una poblacin. Para Xi < O las ecuaciones
obtenidas son anlogas a la ecuacin (2.6) que describe la desintegracin radiactiva.
Hemos analizado el caso cuando los valores propios Ar son reales y diferentes.
Recordemos el mtodo tradicional de. resolucin de un sistema de ecuaciones
diferenciales para este caso. Inicialmente se resuelve la ecuacin caracterstica
det 1A - XE 1 = O y se hallan sus races Al, .,A,. Despus, a partir del, sistema de
ecuaciones algebraicas lineales se hallan los vectores propios
correspondientes a
10s valores propios Xk :
A& =
A continuacin, se desarrolla la condicin inicial Z(0) respecto a los vectores
~ropios lo cual conduce nuevamente a la resolucin de un sistema de ecuaciones
algebraicas lineales para Ck:
..
&
&,
k=l
Ejemplo 1. Estabilidad de los puntos singulares de los sistemas dinmicos. Supongamos que nos interesa la estabilidad de la solucin estacionaria 2' de un sistema
de ecuaciones diferenciales no-lineales
$3 = f@).
(2.8) .
El vector 5' se denomina punto de equilibrio o punto singular si f ( Z * ) =: O. Analicemos
una desviacin pequea del punto singular Z(t) = S*+AZ(t).Despus de este cambio
de variables y de desarrollar en serie de Taylor, obtenemos
Al igual que en el ejemplo anterior, tomemos S de tal manera que la matriz D sea
diagonal. Esto conduce a las ecuaciones
32
Captulo 2
para las componentes del vector y' (yl, ...,y,). Para hacer ms ilustrativo el
mtodo, nos hemos limitado al caso cuando los valores propios Ak son reales y
diferentes. Sin embargo, en el caso de valores propios complejos (si la matriz no
tiene valores propios mltiples o nulos) tambin llegamos al objeto del ejemplo
anterior: u n sistema de ecuaciones diferenciales lineales A$ = AAZ. Es decir, aqu
tambin el problema se redujo a uno ms simple. Est claro que si la parte real de
al menos un valor propio Ai es mayor que cero, entonces y#) -t co si ~ ~ (#0 O,)
y la solucin Z* es inestable. Si Re Xk < O para k = 1,. . . ,n, entonces y k -t O para
k = 1, .. . ,n, y el estado de equilibrio es estable.
+ uu, = vu,,,
u(x, O) = uo(x),
-00
< x < m.
cP
De Vries
+ uuz + u,,,
33
T = e a x eBt ,
lo cual proporciona la misma respuesta,
T,(x,t) = exp
{ -n2m2t
z, ) (i),
nmx
cos
m = O. 1.2,...,
T(X,t) =
CmTm(z, t);
00
00
CmTm(x,O) =
To(s)=
Cmcos
(y).
(2.12)
m=O
m=O
En la ltima igualdad se ha tenido en cuenta que para t = O la funcin T,(x, 0) coincide con cos
(y)
.
cos
(y)
).
~ q ulos
i coeficientes C, figuran
34
por cos
Captulo 2
xnx
(T)
para m # n;
dx = 1;
1 2 ( ~ 7 ~ jJ+cos
)
COS
dx =
d~ = -I
2'
entonces
1
Co = 7 /Ta(x)dx,
C,(x)=
f /To(z)cos
nnx
(--j-)
dx,
n = 1,2,... . (2.13)
T,,, ,) =
,{ F) (y)
ex,
m=o
donde
c,
O;
L a fuaci- G(xi t ) 6 denomriadaAqu-& (O) = 111y -a (m)-=211 para
frecuentemente fitncin de Green o funcin de influencia. Las soluciones con estructura
anloga a (2.15) son caractersticas de la mayora de los problemas lineales de la
fsica matemtica. Su anlisis est relacionado con la bsqueda de la funcin de
Green y la integracin posterior.
-
Ejemplo 1. Funciones generalizadas. Para ilustrar esta hecho, veamos una de las
construcciones ms elegantes de la fsica matemtica, relacionada con la construccin
de las denominadas funciones generalizadas.
Consideremos la solucin, analizada anteriormente, de la ecuacin de conduccin de calor con condiciones iniciales To(x)del tipo
T.(%)=
[f
para
(O
$0
< x < xo + e,
35
En efecto, por cuanto G(x, f , t ) es una funcin continua de sus argumentos, entonces
para todo A tan pequeo como se quiera existe un E tal que si Izo- f 1 < E , entonces
1/
& - G(X,XO, t ) ) E
( ~ ( 3 ,t)
xo
Zo
El smbolo (6, rp) denota el nmero rp(0) (el valor del funcional 6 en la funcin y ) .
Para poder definir la operacin de derivacin de las funciones generalizadas,
las funciones p(x) deben ser suaves y estar definidas en una regin acotada del
Captulo 2
36
espacio. En efecto, si g(x) es una funcin suave, igual a cero fuera de una regin
limitada G con frontera T, entonces
,1
para todos los puntos donde las funciones son diferentes de cero, tienden al lmite
igual a la derivada de la funcin cp(x).
Las bases de las teora de funciones generalizadas fueron sentadas por S. L. Sboliev y L. Schwartz. Este enfoque permiti generalizar el concepto mismo de
solucin. Sea
Lu = f(x).
Se denomina solucin generalizada de esta ecuacin en la regin G a la funcin
generalizada u que satisface la igualdad
para cualquier funcin de prueba cp de D diferente de cero en la regin G .
Aqu, el esquema general sigue siendo exactamente el mismo que en el ejemplo
analizado anteriormente de la ecuacin de conduccion de calor.
Inicialmente se busca G(x), es decir, la solucin de la ecuacin
LG = d(x).
37
(-2).
@(t)
G(x, t ) = 2 a G exP
donde O(t) = O para t < O, y O(t) = 1 para t 2 O. Aqu d(x, t) = d(x) 4 t ) .
T(x, 0) = 0,
Gt = a * ~ , , 6(x, t),
Aqu,
2
lzl)
G = ~ ( o-t 2a
< oo,
Captulo 2
-
Por cuanto la funcin p y sus derivadas son iguales a cero fuera de cierto
intervalo (a, b), integrando por partes se puede hallar la transformada de
Fourier de la derivada
aav(x),
axa
')Los que se conocen por primera vez con el curso pueden omitir este material. A los
estudiantes de fsica y matemtica que estn comenzando a estudiar la ciencia no-lineal, este
material les ayudar a recordar los enfoques tradicionales de la fsica matemtica .lineal.
39
J F ( f (C)).(O
dt =
/ [/ f
De este modo,
( F ( f)/ 9)= (f,F(cp))*
Precisamente consideraremos la ltima igualdad como la definicin de
transformada de Fourier F ( f ) de la funcin generalizada f . A partir de la
igualdad (2.21)podemos introducir la transformada inversa de Fourier
Consecuentemente,
F(1) = 2n6(x).
Integremos la funcin 6 ( t ) respecto al tiempo:
Captulo 2
40
6?(~,t ) de
la funcin G,
G(e,t ):
G(t,t) de la fun-
41
dz = J.Demos-
-00
trmosla:
I~ =
eeZ2dz
J e-y2dy = J J e-" -Y
2
dz dy.
"" J . "
2 -
+ y2 = r2,
II
de la integral de I:
Demostremos que X
P
'
= 1:
42
Captulo 2
As pues, resulta que este enfoque tan notable est orientado fundamentalmente a los problemas lineales. Y esta misma situacin caracteriza a muchos
otros enfoques analticos.
Preguntas y ejercicios
1. Construir un modelo matemtico que permita predecir la intensidad de la
radiacin emitida por las sustancias radioactivas contenidas en cierto depsito.
Suponer que se conocen las cantidades iniciales de todas las sustancias, sus
perodos de semidesintegracin y las leyes de su ingreso al depsito con el
transcurso del tiempo.
2. Un pndu.10 realiza oscilaciones pequeas. Hallar la ley de movimiento del
pndulo suponiendo que sobre l acta una fuerza exterior A sen wt . Cmo
vara este resultado si se considera un rozamiento pequeo proporcional a la
velocidad del pndulo?
3. Cul es la ley de movimiento de un pndulo sobre el que acta una fuerza
perturbadora peridica?
4. Cul es la ley de movimiento de un pndulo sobre el que acta una fuerza
que vara con el tiempo segn una ley dada?
5. Cmo vara respecto al tiempo la temperatura de una barra inicialmente
fra cuyos extremos estn aislados trmicamente? Suponer que se conoce la
ley f (x, t) de calentamiento de la barra.
6. En los extremos de una barra inicialmente fra la temperatura vara segn las
leyes L(t) y R(t). Cmo cambia la distribucin de la temperatura de la barra
con el tiempo?
-+
7. E1 momento angular L E (Lx,L,, L,) y la velocidad angular ti7 S (w,, w,, w,)
de un cuerpo slido arbitrario estn relacionado mediante la igualdad
i = Fd,
Ixz Ixv
donde
I es el tensor de inercia
I,, Iyy
I,, . La matriz i es simtrica:
h
IZX
Izy = I y x r
IXZ
Izy
12%
Ixr =Izxr
1 . y =Iy,z*
+ + + + +
43
'x
x'
y'
y'
a2
b2
a2 b2
y
Cx2 (parbola). Dada la ecuacin de segundo orden con coeficientes
buzy myy duz + euv f = 0, es posible reducirla
constantes au,,
mediante un cambio de variables de las variables independientes a la forma
cannica u,, +uYY= O (ecuacin de Laplace), u, = u, (ecuacin de conduccin
de calor), u,, - uyy= O (ecuacin de oscilaciones de una cuerda)?
10. Ciertos procesos en determinado medio se describen mediante la ecuacin
Bibliografa
Los mtodos de investigacin de los modelos matemticos lineales se analizan detalladamente en los cursos fundamentales de lgebra lineal, ecuaciones diferenciales
ordinarias y de la fsica matemtica. Por ejemplo, en los libros siguientes:
Gantmacher E R. The theory of matrices. Chelsea, 1960;
Ikhonov A. N., Vasileva A. B., Sveshnikov A. G. Differential equations. SpringerVerlag, 1985;
Arnold V. 1. Ecuaciones diferenciales ordinarias. ~ubios,1994;
Tfjonov A. N,, Samarsky A. A. Ecuaciones de la fsica matemtica.
. , .
MIR 1983;
Svishnikov A. G., Bogolibov A. N,,Kravtsov l? V. Conferencias de fsica matemtica (en ruso);
Courant R. Partial differential equations. Wiley 1962.
Las funciones generalizadas se analizan detalladamente en los libros:
0 Vladimirov I/: S. Equations of mathematical physics. Editorial URSS, 1996;
Vladimirov V: S. Generalized Functions in Mathematical Physics. Editorial
URSS, 1994.
Sobre el papel fundamental de los modelos matemticos lineales se puede juzgar
leyendo los tomos del curso de fsica terica de
0 Landu L. D., Lqsshits Ye. M , Teora del campo, Mecnica cuntica. Teora
no-relativista y <<Teora
de la elasticidad.
Los modelos matemticos lineales ayudan en muchos casos a crear nociones intuitivas
valiosas acerca de las ondas, las oscilaciones y otros procesos fsicos, nociones que
luego resultan ser muy itiles tambien en el mundo no-lineal. Sobre esto se puede
leer en el excelente libro
Trubetskov D. 1. Introduccin a la sinergtica. Oscilaciones y ondas. Editorial
URSS, 2003 (de inminente edicin en espaol en nuestra editorial).
CAPITULO 3
Modelos matemticos
no-lineales elementales
Todo conocimiento humano comienza con
la contemplacin, pasa de sta a los conceptos
y termina con las ideas.
1. Kant
Durante ms de trescientos aos los cientficos se han dedicado a la simulacin matemtica de los fenmenos no-lineales utilizando las ecuaciones diferenciales. Esto ha conducido al surgimiento de varios enfoques alternativos,
los cuales ilustraremos en el ejemplo de los modelos elementales. Aqu estudiaremos algunos esquemas de anlisis simplificados y estilizados, que nos
permitirn formarnos una idea clara de las posibilidades de cada enfoque.
E l pndulo matemtico
Analicemos el problema del pndulo matemtico. Consideremos el sistema
mostrado en la figura 3.1. Durante la simulacin hay que recurrir inevitablemente a algunas simplificaciones.Supongamos que las dimensiones del peso
son despreciables y ste se puede considerar como un punto material de
masa m. Adems, asumiremos que el hilo
es imponderable e inextensible, y que en
su movimiento el pndulo no experimenta
resistencia por parte del aire.
m ? = F.
Proyectemos los vectores de ambos miembros de la igualdad en la direccin perpendicular al hilo (fig. 3.1). Obtenemos
Fig. 3.1
mlYj
+ mgsen (o = O.
(3.1)
45
m,
1=-27,
a=O,
@+?=O,
46
Captulo 3
fi.
[ E ] = M L ~ T - ~ ,[ p ] = M L - ~ , [ r ] = L,
[ t ] = T.
.=C
(;)
115
t2~5
C2 ( E ) 115 t-315.
5 P
Este problema fue analizado detalladamente por el acadmico L. 1. Sedov
en 1945. La dependencia obtenida fue confirmada brillantemente ese mismo ao por
los datos de los experimentos nucleares realizados en Nuevo Mxico. De acuerdo
con las memorias de los testigos presenciales, durante el primer experimento nuclear
Enrico Fermi lanzaba papelillos y observaba a qu distancia eran arrastrados, para
as determinar la fuerza de la explosin. Por lo visto, l tambin utiliz las ideas de
la semejanza y las dimensiones. La estimacin de E obtenida de este modo result
ser muy aproximada a la magnitud real de la energa de la explosin.
47
1
I
60"
120"
180" cp,
Aqu hemos considerado la siguiente consecuencia de (3.4). Para @ = O tenemos p = (00; por tanto, cos po= -C/2w;. La integral (3.5) es precisamente
la solucin. La ecuacin diferencial se integra en cuadraturas. Sin embargo, no
es posible calcularla mediante funciones elementales. Formalmente, hemos
obtenido la respuesta, pero extraer de ella alguna informacin concreta y
en particular, la dependencia mostrada en la figura 3.2, no es fcil.
48
Captulo 3
cp = cp, sen wt
(3.7)
2
w,cp = wocpo
sen
wt
2-
--
--
...
--
49
la frmula (3.7) nos limitamos slo a dos trminos. Por eso recurriremos al
procedimiento estndar que se utiliza habitualmente en la construccin de
soluciones asintticas. Concretamente, consideraremos slo los trminos ms
importantes, en los que el parmetro pequeo E figura con el exponente mnimo. Asimismo, despreciaremos los armnicos superiores. Conforme a esta
ltima condicin, eliminamos todos los trminos que contienen a sen 5wt y
sen 7wt. En los parntesis que multiplican a sen wt y sen 3wt tendremos en
cuenta solamente los trminos que no contienen a 6. De este modo, quedan
slo dos trminos, los cuales estn subrayados en la frmula (3.8).
Agrupando los trminos que contienen a sen wt e igualando a cero,
obtenemos
-w 2 $wo2 - ~3W o2 ' P2 o-- o.
Por consiguiente,
De la frmula (3.9) queda claro por qu las oscilaciones, incluso para ngulos
suficientemente grandes, son aproximadamente armnicas. Esto se debe a
que el parmetro (pi/192 es pequefio. La frmula (3.9) da una idea de una
parte de la curva T = T ( E ) ,pero no da la posibilidad de analizar lo que
sucede para valores de la energa cercanos a Bo.
En la figura 3.2 se muestra la dependencia (3.9). La recta horizontal
inferior corresponde al problema lineal y la curva 2 corresponde a la
dependencia (3.9).
Comentario. A finales del siglo XVIII, despus de los trabajos de P. S. Laplace,
los mtodos asintticos se convirtieron en el instrumento principal de anlisis de
muchos modelos matemticos y, en particular, de los modelos de la mecnica celeste.
Por cuanto ya se conoca la solucin del problema de los dos cuerpos, y la masa de los
50 .
Captulo 3
cuerpos del sistema solar son incomparablemente menores que la masa del Sol, la
influencia de stos en el movimiento de un planeta dado se poda considerar como
una perturbacin pequea y utilizar los mtodos asintticos. Era muy grande la
seguridad de que se haban comprendido las leyes fundamentales de la naturaleza
y construido un aparato matemtico efectivo. En particular, Laplace respondi a
Napolen, quien se interesaba por la matemtica y haba ledo su libro sobre
mecnica celeste, que l no necesitaba de la hiptesis de un creador todopoderoso.
En efecto, por este camino se haban logrado xitos notables en la simulacin de
fenmenos no-lineales, haba sido predicha la existencia de Urano y Neptuno. En el
siglo XX, las posibilidades de utilizacin de los mtodos asintticos para el anlisis
de los modelos matemticos crecieron inconmensurablemente.Sin embargo, tambin
se comprendi cun- limitado es el campo de aplicacin de estos mtodos entre
todos los modelos no-lineales que era necesario investigar.
primera especie
4 2
I(K) =
J (1o
di,
sen2+)'12
'
sen (542)
sen (9,/2)
51
Esta frmula permite predecir un segmento un poco mayor de la dependencia T = T(p,) (fig.3.2, curva 3). Adems, existen tablas de las funciones
y utilizndolas se puede obtener toda la dependencia T = T ( E ) .
comentario. En el siglo XIX y primer tercio del siglo XX, la teora de las funciones
se consideraba una de las reas ms importantes de la matemtica
aplicada. Lamentablemente, muchas esperanzas cifradas en este aparato matemtico
resultaron en alto grado injustificadas. La cuestin es que en muchos casos la
construccin de las tablas de funciones especiales exige un gran volumen de clculos
y la utilizacin de mtodos numricos, y operar con stos no es algo fcil o evidente.
Adems, a diferencia de las funciones exponencial y trigonomtricas, las funciones
especiales no slo son ms complejas sino que tienen un campo de aplicacin
significativamente menor.
-.
52
Captulo 3
.-
para
x(t) = 5o= const.
Prestemos atencin a la profundidad y el contenido de los conceptos
de masa, fuerza y aceleracin. Para acoger la ley de Newton como una
afirmacin llena de contenido (y no como una definicin de fuerza),
es necesario tener la posibilidad de medir de manera independiente la
aceleracin, la fuerza y la masa. Y los fsicos saben cmo hacerlo.
Segn la hiptesis de Newton, la forma de la rbita est determinada
por la atraccin que el Sol ejerce sobre el planeta. De acuerdo con la ley de
la gravitacin universal, dicha fuerza es
#(O) = O,
Proyectando los vectores sobre los ejes del sistema de coordenadas cartesianas, = (x, y, z), 8 = (v,, vy, v ~ ) la
, ecuacin (3.14) se puede escribir en
forma escalar:
53
.-
+ 1 ) ~ )En. efecto,
El mtodo que permite pasar de la ecuacin diferencial (3.17) a la aplicacin (3.18) recibi el nombre de mtodo de Euler.
La figura 3.4 ilustra el mtodo de Euler. La curva continua x(t) es
sustituida por una serie de quebradas. El error de la aproximacin se debe
a que sustituimos la curva O A por la recta OB y en el paso siguiente
, en el
calculamos el segundo miembro no en el punto A ( @ ( X ( T ) ) )sino
punto B (#($(O))
TF(X(O)))
La diferencia en el primer miembro de la expresin (3.16) tiende a la
dx
d2x
derivada - cuando T + O, y la derivada - se mantiene finita.
dt
dt2
54
Captulo 3
Fig. 3.4. Solucin numrica del problema del movimiento de u n planeta alrededor de una
estrella masiva: a ) e l error cometido a l utilizar e l mtodo de Euler est relacionado con e l
hecho de que la curva se sustituye por una serie de quebradas y la derivada no se calcula
en e l punto correcto; 6 ) ejemplo de solucin numrica del problema (3.15) mediante e l
mtodo de Euler. Se observa que esta dependencia se parece a u n arco de elipse. La elipse
corresponde a la solucin exacta de ta ecuacin diferencial
La estrategia posterior es clara. El segmento [O, TI se divide en N intervalos T = T I N . Despus se calculan las iteraciones de la aplicacin (3.18).
Se puede esperar que si el paso temporal es suficientemente pequeo,
entonces las iteraciones de la aplicacin y los valores de la funcin X(nr)
para la ecuacin inicial sern prximos. El resultado tpico, correspondiente
al clculo para el problema analizado, se muestra en la figura 3.4. Como
se observa, la sucesin de quebradas aproxima bastante bien n la solucin
exacta (una elipse) de la ecuacin diferencial.
La aplicacin-(3.18), as-como la ecuacin-inicial, describe un sistema
determinstico (el comportamiento futuro del sistema est unvocamente
determinado por el estado actual) con espacio de fase de dimensin finita
(el estado del objeto est determinado unvocamente por un conjunto finito
de nmeros). La dimensin de este espacio coincide con el nmero de
componei~tesdel vector 2. Sin embargo, aqu la coordenada temporal ya
es una magnitud discreta que puede tomar slo los valores (0, T , 27, . ..).
A esta clase interesante e importante de objetos nosotros regresaremos ms
de una vez. Los sistemas de ecuaciones diferenciales ordinarias (3.17), en
los que el tiempo es continuo, y las aplicaciones tipo (3.18), en las que el
tiempo es discreto, se denominan actualmente sictearas dinbniicc~s.
En el mtodo de Euler, la diferencia entre la derivada eii diferencias
dx(nr)
x((n 1)r)- x(nr)
es de orden O@).Tiiles mtodos
y la derivada
7
dt
se denominan mtodos de primer orden. El enfoque tradicional para mejorar los
-
55
.-
nr
f (x(n + 1)r)en la
2
7
frmula (3.20) debe aparecer el trmino - f (x(nr) + T f ( ~ ( n r ) ).)As pues,
2
en lugar de la frmula (3.20) obtenemos la frmula aproximada
ms simple y menos preciso. Entonces, en lugar de
56
Captulo 3
-t
X ( t ) = $(O)
-t
+ J P(X(t))dt,
Las condiciones de convergencia de esta sucesin se analizarn ms adelante. Para simplificar hemos supuestg que contarnos con un mtodo
suficientemente preciso para el clculo de la integral.
Comentario. La ltima frmula proporciona no slo un mtodo de construccin
iterativa de una solucin numrica, sino tambin la llave para la demostracin
del resultado fundamental de la teora de ecuaciones diferenciales ordinarias: los
teoremas de existencia y unicidad.
En efecto, esta frmula puede ser escrita. en la forma simblica
57
<
~(A~x,
< X~ (A"-'
X, A ~ X )
00
La serie
se puede indicar un
n=O
El punto
2 es un punto fijo:
~2 = A a+oo
lim An2 = lim A'+'x
n+oo
-+
=X,
IIx
"Sfl
<
(T) - ZS(~)IIJ L
11%'
Captulo 3
58
Sean X , F E R ~ Supongamos
.
que la funcin F (t, 2) es continua en el paraleleppedo R: O t a, 12- fOl b, y cumple la condicin de Lipschitz
respecto a 2. Supongamos adems que M es la cota superior de I@(t, 2) en R
y que T = min (a, blM). Entonces, el problema de Cauchy
<
< <
Pasemos ahora al mtodo de resolucin analtica del sistema de ecuaciones (3.15). Estas ecuaciones son vlidas en cualquier sistema de coordenadas
cartesianas. Para simplificar el anlisis posterior, elijamos este sistema de tal
modo que z(0) = O, v,(O) = O. Pero, entonces, a partir de la tercera ecuacin
dv
se deduce que tambin 2 = 0, por lo que el movimiento tendr lugar en
dt
el plano (x,y). Este movimiento es descrito por cuatro ecuaciones de primer
orden, las cuales se pueden reducir a la forma
59
.-
();= -)
mr2
y sustituyendo
dU
( F = - %) se denomina potencial.
Captulo 3
60
Consecuentemente,
a
r
obtenemos
dz
- arccos z
fi
+ const,
donde
M
r
z=---
m2a
, c = 2mE+
rn2a2
M2
-.
En conclusin,
cp = arccos
+const.
P
M*
-=l+ecosp,
donde p = r
rna
De la geometra analtica se sabe que para e < 1 (E < O) sta es la ecuacin
m
m
'l
'l
. La dependencia de
de una elipse con semiejes a = 7
y b=
1-e
Jc-2
la coordenadas respecto al tiempo est dada por la integral (3.24), bien
por una expresin ms elegante; pero implcita.
Comparemos los dos mtodos analizados teniendo en cuenta los
ejemplos citados.
La ventaja de la estrategia relacionada con la utilizacin de mtodos
numricos es la simplicidad. Un algoritmo no muy complicado y fcil de
programar proporciona, en principio, una solucin del problema planteado.
Su otra ventaja es la universalidad. Compliquemos un poco el problema
suponiendo que tenemos no dos, sino tres cuerpos. Por ejemplo, una
61
.-
ck
62
Captulo 3
aT
63
64
Captulo 3
b)
a)
Fig. 3.6. Ejemplo de sistema disipativo (oscilador armnico con rozamiento ,viscoso). En este caso Las
trayectorias de fase son espirales
65
.de las soluciones cuando t + m. Dado que la imagen matemtica de los
estacionarios es un conjunto de atraccin en el espacio de fase
(un atractor) entonces lo ms importante, segn Poincar, son las particularidades cualitativas de los atractores. En particular, es importante saber si el
atractor describe regmenes estacionarios variables, regmenes peridicos u
otros ms complejos. Son esenciales el nmero y tipo de atractores, as como
el conjunto de datos iniciales desde los cuales se llega al atractor dado. Estos
conjuntos se denominan regiones de atraccin del atractor. Precisamente la
teora cualitativa de las ecuaciones diferenciales analiza estos problemas. Los resultados de esta teora son muy generales. Ellos muestran que una cantidad
enorme de sistemas disipativos no-lineales se comportan de igual manera.
Supongamos que la ecuacin
~ o d e l o smatemticos no-lineales elementales
t+oo
X(t+T)=X(t),
-
- -
- - - - - . .
t+oo
Y(t+T)=Y(t).
-
(3.30)
66
Captulo 3
Fig. 3.7. Tipos de atractores en -un sistema autnomo de dos ecuaciones diferenciales
ordinarias. A La izquierda se muestran las curvas integrales y a la derecha, la forma
caracterstica de las trayectorias de fase: o) punto singular estable; b) ciclo lmite
Anlisis local de los sistemas dinmicos. Durante la resolucin de ecuaciones algebraicas es un hecho frecuente que un cambio de variables acertado
desempee un papel fundamental. Por eso es natural buscar tales cambios
de variables en las ecuaciones diferenciales, que permitan reducirlas a formas ms simples. Es de esperar que sea imposible hallar semejante cambio
de variables para todo el espacio de fase (globalmente). Sin embargo, si se
considera el problema localmente, en un entorno pequeo del espacio de
fase, entonces podemos tener esperanzas de xito.
En efecto, consideremos las trayectorias de fase que figuran en este
libro. Su forma tpica lejos de los puntos singulares es aproximadamente la
mostrada en la figura 3.8a. Se puede suponer la existencia de un cambio
de variables (x, y) -+ (x', y') que <<enderece
las trayectorias de fase y les
transmita la forma mostrada en la figura 3.8 b.
El ltimo campo vectorial, es decir, el conjunto de vectores P del espacio
de fase, corresponde a la ecuacin diferencial elemental
67
Fig. 3.8. Lejos de los puntos singulares e l campo vectorial puede ser enderezado Localmente mediante un cambio de variables: a) campo vectorial inicial;
6 ) campo vectorial enderezado
--
68
Captulo 3
-
69
.-
descubriremos que la magnitud de la carga total Q se conserva. Esto impone limitaciones serias a los procesos posibles en la naturaleza. Si en un
sistema nace un barin con carga +1, entonces debe nacer un antibarin
con carga -1, un electrn puede nacer slo en pareja con un positrn.
~ o s t r e m o sque para la ecuacin x = v ( x ) tiene lugar una situacin
Supongamos que O < x 2n, v(0) = v(27r). En este caso los segundos miembros son peridicos y el espacio de fase es una circunferencia.
Mostremos que aqu tambin se puede introducir un anlogo de la carga.
Diremos que un estado de equilibrio xi ( f ( x i )= O ) tiene carga topolgica qi
dv(xi)
dv(~j)
~ V ( .)
X
igual a -1 si -< 0, igual a +1 si -> O, e igual a O si A - O.
.
dx
dz
dz
La suma respecto a todos los estados de equilibrio de tales sistemas se f (")
denomina carga topolgica del sistema y
A
se denota mediante el smbolo Q =
qi.
i
Demostremos que para todos los sistemas
dinmicos con -funcin v (x) continua, la f 4
a)
carga topolgica es igual a cero.
En el caso elemental, cuando v ( x ) es
una magnitud constante, tenemos que
Q = O (fig. 3.10 a). Demostremos que este
valor no cambia al deformar continuamente la funcin v(x). Se puede razo- f (x)
nar del siguiente modo. Supongamos que
A
v ( x ) es un hilo sujeto a los puntos A y B
(fig.3.10 b). Vamos a tensar paulatinamente este hilo. Mientras no cambie el nmero
de estados de equilibrio, la magnitud Q
no vara: Para
longitud del <<hilo Fig* 3.10. AL {(tensar e l hilo>) la carga
topolgica del sistema no cambia
surge la situacin mostrada en la figura 3.10 c (curva 2), desaparecen simultneamente dos estados de equilibrio. En el instante de desaparicin, la carga
topolgica del estado de equilibrio z* es igual a cero. Despus de desaparecer los dos estados, uno con carga +1 y el otro con carga -1, la magnitud Q
sigue siendo la misma que antes de la deformacin. Actuando de la misma
manera para el caso general, cuando tenemos un ninero finito arbitrario
de estados de equilibrio, llegamos a la conclusin de que la magnitud Q no
cambia durante una deformacin continua. Por cuanto Q = O para el caso
mostrado en la figura 3.10n, esta magnitud ser la misma tambin en los
dems casos. As pues, la carga topolgica para tales sistemas se conserva.
Razonamientos ms profuridos y de mayor contenido, los cuales se
apoyan en el concepto de continuidad y en la posibilidad de analizar las
<
Capitulo 3
70
Preguntas y ejercicios
--
+S + x
71
72
Captulo 3
Atractores de La ecuacin x = v ( x )
Intente ejecutar piezas fciles correctamente. Esto es mejor que interpretar piezas
difciles mediocremente.
= =(E)
t2
74
Captulo 4
En efecto, para que la reaccin tenga lugar, las molculas de la sustancia X deben
chocar con las molculas de la sustancia Y. Es claro que la probabilidad de que esto
75
ocurra es proporcional
- al nmero de molculas de X por unidad de volumen (es decir,
a la concentracin
- X). De igual manera, esta probabilidad debe ser proporcional a la
Y . El coeficiente de proporcionalidad k depende de las dimensiones
de las molculas, de sus velocidades, etctera. Esto es lo que refleja la frmula (4.1).
Sin en la reaccin interaccionan n molculas de X con una molcula de Y, entonces
la variacin de la concentracin de la sustancia Z es proporcional a XnY.
Consideremos una reaccin ms compleja:
C1
z,.
76
Captulo 4
x ( t ) -+ x* cuando t
+ oo,
v ( x * )= O
2)); pueden
x ( t ) -t oo cuando t -t oo
(4-5)
crecer
(4-6)
x ( t ) -+ oo cuando t + t f
(como en el modelo elemental de una reaccin en cadena x = xp, /3
(4.7)
> 1).
Teoremas de comparacin
A finales del siglo XIX los enfoques tradicionales de la matemtica clsica,
relacionados con la obtencin de soluciones en cuadraturas, con el desarrollo
asinttico y con la introduccin de funciones especiales, comenzaron a
tropezar con dificultades crecientes.
Esto-exiga nuevas ideas. Uno de los mtodos de anlisis cuyo uso
no est relacionado ni con la linealidad, ni con una u otra simetra de
la ecuacin, son los teoremas de comparacin. Los ms simple de estos
teoremas se pueden ilustrar mediante dos modelos matemticos elementales.
Movimiento de un automvil durante un rally. Supongamos que al conductor de un automvil se le entrega un grfico de movimiento donde se
indica la-velocidad v(t) que debe tener el vehculo en cada instante. Entonces
su camino x est dado por la ecuacin
77
"t
e
xl ( t )
t*
Fig. 4.1. Ejemplos de curvas integrales de Las ecuaciones diferenciales que describen e l movimiento de los automviles
Esquiador de montaa. Supongamos que por una ruta corren dos esquiadores. Sus velocidades dependen del sector de la ruta en el que se encuentran:
v = v(x). A uno de los esquiadores se le proporciona el grfico de movimiento del otro esquiador, y l trata de vencer cada sector de la ruta ms
rpido, vl(x) > v,(x). Puede afirmarse que l llegar primero a la meta?
Teorema 2. Sean
>,
78
Captulo 4
As,
x1(t*- At) > x2(t*- At),
~ , ( t *=) x2(t*).
Por consiguiente,
<
At
At
Utilizando el teorema de Lagrange llegamos a la desigualdad
V, (t* - OIAt) < v2(t*- 02At).
-t
= v(x),
dt
$(O) = x,.
Supongamos, adems, que tenemos dos ecuaciones cuyas soluciones se
conocen: .
Adems, v,(y) < v(y) < v2(y) para O < y < m. En este caso, nosotros podemos encerrar la solucin en un corredor entre las funciones z,(t) y z2(t),
es decir, z1(t) x ( t ) z2(t), lo que puede proporcionar una informacin
importante sobre el comportamiento del objeto analizado.
<
<
79
(4.9)
Captulo 4
80
.-
Fig. 4.3. Soluciones de la ecuacin = a y - py2: a) vista de La parte derecha, correspondiente a La ecuacin logstica; 6 ) algunas curvas integrales correspondientes a esta
ecuacin; c) vista de la parte derecha, correspondiente a l modelo de una reaccin en cadena;
d ) algunas curvas integrales para este caso
81
obtenemos
+
d(x*+ Ax)
1 d2f (x*+ @A%)
df (x*)
= f (x*+ Ax) = f (x*)+ -n x + (ax2),
dt
dm:
O<O<l.
dx2
82
Captulo 4
donde ,O y son constantes (fig. 4.4 a). Entonces se puede utilizar el teorema
de comparacin y la solucin de las ecuaciones (4.2):
df (x*)
donde 6 es un nmero suficientemente pequeo y a = -. La forma
dx
Elijamos un nmero E. Si 6 = E
<
83
.-
<
<
<
lim V(jJn)= V ( g ) .
n+oo
positiva L(Z), continua junto con sus derivadas parciales de primer orden,
Captulo 4
t > 0,
W(2,t)<O,
Z E G,
<
lim 5(t)= O.
t+oo
--
--
--
--
--
--
<
<O
(O
85
t+oo
<
Captulo 4
86
donde
au
=
f ( x ) = -f * ( x )< 0.
dx
axL
87
f ( x ) = x @ ,,8
> 1y U =
Para que la solucin del problema (4.8) exista durante u n tiempo finito t f es
necesario y suficiente que la integral
Captulo 4
88
En el caso extremo de una viscosidad infinitamente pequea, esta expresin se transforma en la segunda ley de Newton, la cual describe el movimiento de una partcula
en un campo de potencial U($).En el otro caso extremo, esta expresin refleja
la nocin de Aristteles acerca de que para que un cuerpo se mueva es necesario
aplicarle una fuerza. Realmente, en la literatura dedicada a la simulacin matemtica
de fenmenos no-lineales, el trmino potencialse utiliza en los dos sentidos. En lo
sucesivo, si no es evidente, nosotros aclararemos de cul de ellos se trata.
existe,
es nica,
es estable.
La estabilidad puede entenderse respecto a las perturbaciones pequeas
de las condiciones iniciales, de las condiciones de contorno, de los parmetros
de la ecuacin.
Los teoremas de existencia y unicidad para un sistema de ecuaciones
diferenciales ordinarias muestran que un problema con condiciones iniciales
(problema de Cauchy) est bien planteado en cierto intervalo de tiempo
O t < T. Este intervalo puede ser finito, depender de los segundos
miembros de las ecuaciones y de las condiciones iniciales.
A primera vista puede parecer que slo los modelos matemticos bien
planteados tienen derecho a describir algn fenmeno de la naturaleza. Sin
embargo, el anlisis de una serie de modelos matemticos no-lineales oblig
a ver de un modo nuevo el concepto de problema bien planteado. Esto
<
89
Captulo 4
90
Adems, las ecuaciones tipo (4.14) surgen en la fsica del plasma como
un modelo matemtico de la denominada inestabilidad explosiva. En la fsica
de la atmsfera a veces se analizan corrientes de gas con viscosidad negativa,
cuya descripcin exige la resolucin de problemas similares.
Procederemos del mismo modo que en el caso de una ecuacin
lineal habitual de conduccin de calor, y buscaremos la solucin como un
) (kx)):
desarrollo respecto a las funciones propias ( ~ ~ (=2cos
U(X,
t) =
Ck(t)cos (kx).
<
<
91
De este modo, la solucin analizada, para la cual la serie respecto a las funciones propias converge, existe para O t < t f , pero no puede ser representada
mediante una serie convergente en todo el semiintervalo infinito O t < OO.
Adems, si hay datos iniciales uo(z) que conducen a una solucin
acotada para O < t < m, entonces nuevamente, debido al principio de
E W ~ ( X )existir
superposicin, la solucin con datos iniciales Zo(x)= U,($)
slo hasta el instante t sin importar cun pequeo sea el valor E elegido,
siempre que la funcin~,(x)posea esta propiedad. En este caso el tiempo t
tambin se puede elegir tan pequeo como se quiera. De este modo, inclusive
si exigimos que los datos iniciales en el modelo pertenezcan a la clase de funciones infinitamente diferenciables y nos limitamos a un intervalo finito de
tiempo T (O t T), el problema (4.14) no ser un problema bien planteado.
El sentido fsico de estos razonamientos es muy simple. La ecuacin usual de conduccin de calor produce, en un tiempo t, una
funcin muy suave a partir de la funcin-6 de Dirac. u
La ecuacin de conduccin
de calor con coeficiente negativo rueda la pelcula,
en la que se filmesta solucin, en sentido contrario:
a partir de un perfil suave,
al cabo de un tiempo tf
se obtiene la funcin delta.
o
7r
x
Este perfil suave se puede
tomar en calidad de fun- Fig. 4.7. Forma caracterstica de la solucin de ia ecuacin vo(x).El cuadro tpico cin de conduccin de calor con coeficiente negativo
de tal proceso se muestra
en la figura 4.7.
Hace veinte aos an se pensaba que las ecuaciones cuyas soluciones
no existen totalmente, no tienen sentido fsico. Sin embargo,'desde entonces
han aparecido argumentos de peso a favor del anlisis de los problemas mal
planteados. Estudiemos uno de los modelos matemticos fundamentales,
surgido durante la investigacin de muchos fenmenos no-lineales.
<
<
,,
<<
92
Captulo 4
Entonces,
donde
93
= L,
(4.18)
94
Captulo 4
<<
312 312
=m,
dx .=o
y no es posible elegir la constante de
Lipschitz L.
x = (213) t
donde 6(t) es una funcin aleatoria. En este caso habr que renunciar a la
descripcin dinmica determinstica y hablar de las corrientes turbulentas
slo en el lenguaje probabilstico, estadstico.
En el concepto de problema bien planteado figuran no slo las condiciones de existencia y unicidad, sino tambin de estabilidad. La investigacin
de los fenmenos no-lineales ha hecho cambiar considerablemente la idea
de estabilidad de los sistemas no-lineales. Recordemos la definicin clsica.
Definicin. La solucin tj = V(t, yo) del problema
para todo t
> 0.
96
Captulo 4
Xt
e IIAyoll,
(4.21)
< <
<<
97
Fig. 4.11. Forma de un atractor extrao que describe e l caos dinmico en un modelo de
cintica qumica: a) proyeccin del atractor sobre e l plano (t,v); b) proyeccin panormica
del atractor. Este atractor surge en e l sistema dinmico
c, = -6, 1 = r
98
Captulo 4
cae aleatoriamente ora hacia la franja derecha, ora hacia la izquierda. Esto
sugiere la idea de una descripcin simblica de semejante trayectoria como
una sucesin infinita de ceros y unos. Si en 'la k-sima vuelta el punto
se mueve por la franja izquierda, en la k-sima posicin de la sucesin
escribimos un O; si se mueve por la franja derecha, escribimos un 1. Result
que a veces esta descripcin es muy til. Al hablar de muchos atractores
extraos de los sistemas continuos es mucho ms fcil expresarse en este
lenguaje discreto y extremadamente simple. Este enfoque recibi el nombre
de dinmica simblica y se est desarrollando activamente en la actualidad.
Supongamos que lanzamos un nmero infinito de veces una moneda. Si en
el k-simo lanzamiento cay cara,escribimos un 1 en la k-sima posicin
de la sucesin que describe este experimento; en caso contrario escribimos
un O. Un hecho notable, demostrado para una serie de atractores extraos,
consiste en que a partir de esta sucesin no es posible diferenciar un sistema
dinmico de un objeto clsico de la teora de probabilidades, como lo es
una serie de lanzamientos de una moneda.
Lo extrao>>
de los atractores extraos est relacionado con su peculiar
propiedad de sensibilidad respecto a las condiciones iniciales. Elijamos dos
puntos prximos Z1(0) y Z1'(0) pertenecientes a un atractor, y veamos cmo
vara la distancia d(t) = IZ1(t) - x"(t)l respecto al tiempo. ?'(t) y ?''(t)
son las soluciones de las ecuaciones analizadas, con condiciones iniciales
Z1(0)y Z1'(0), respectivamente (fig. 4.10). Geomtricamente, d(t) corresponde
a la longitud del segmento con extremos 2(t)y Z1'(t). Si el atractor es un
punto singular, entonces d(t) = O. Si el atractor es un ciclo lmite, entonces
d(t) es una funcin peridica del tiempo. En un atrctor extrao, d(t) ext,
X > O. Para que la magnitud X caracterice al atractor es natural analizar
trayectorias infinitamente prximas y la velocidad media de dispersin en
un intervalo grande de tiempo. Formalmenti esta magnitud se introduce
del modo siguiente:
N
X ( & ( ~ ) , d ) = l i mlim
t+oo d(O)+O
99
+ + +
a = ~ ( Z Z 3 = f($)+ + 3 = f(3~ + 5)
f(2) +
dx
a?(~)
100
Captulo 4
G,
Preguntas y ejercicios
1. En un reactor tiene lugar una reaccin en cadena durante la cual la velocidad
de variacin de la concentracin n de cierta sustancia vara segn la ley pn7
(p > O, 7 > 1). En el instante inicial la cantidad de esta sustancia en el reactor
es igual a cero. Luego la sustancia se introduce en el reactor segn la ley at2.
Estimar el tiempo al cabo del cual se acaba la reaccin, o el modelo elemental
que describe la variacin de la concentracin de slo una sustancia se hace
inadmisible.
2. Los eclogos construyeron el siguiente modelo para describir la variacin del
tamao de cierta poblacin:
k = F(z).
~F(x*)
101
x = ax,
(a)
se debe reemplazar por otro modelo que concuerda mejor con la curva de
crecimiento de la poblacin en los ltimos 100 mil aos:
x = axl+&,
> 0.
(b)
102
Captulo 4
-
CAPITULO 5
Elementos de la teora
de las bifurcaciones
Qu se puede esperar de ecuaciones nzttcho ms
104
Captulo 5
105
La constante Z pueden tomar varios valores: x,, x,, etctera. Est claro
que
F(Z, A) = O.
(5.2)
F(x,
A,)
+As, A, +AA) =F(x,, A,) + aF(xO,
6' A
~F(x
X )
O'
dx
O AX+
6F(~,f
A,)
es diferente de cero,
dx
podemos determinar de manera aproximada el nuevo estado de equilibrio
(fig. 5.3). Tal estado, a juzgar por la frmula (5.4), es nico, lo que significa que en el
punto (A,, s o ) no hay bifurcacin.
En la igualdad (5.4) nosotros consideramos slo dos trminos no-nulos de la serie
de Taylor (5.3). Surge la pregunta de cun
legtimo es esto y qu relacin tiene con el
problema no-lineal (5.2). La respuesta a es4
X
ta pregunta la proporciona el teorenla de la
funcin implcita, conocido del curso de anlisis matemtico. Este teorema desempea un Fig. 5.3. Una situacin en l a que
se puede utilizar e l teorema de
papel clave en la teora de las bifurcaciones.
la funcin implcita
Enuncimoslo:
Captulo 5
106
1. La ecuacin (5.2) tiene una nica solucin x = x(E), cuando A, - a < X <
A, a, tal que x, - p < x < 2, + p .
2. La funcin x = x(A) es diferenciable con continuidad para A, - a < X i:
A,
0.
~F(x,,
ax
# o*
dx
- -
d2F(x,A,)
0z2
# O, en lugar de (5.4)
Aqu el cuadro es otro (fig. 5.4): para A. > A, aparecieron dos soluciones,
y para X < Xo no hay soluciones (si se considera que cl es positivo).
dx
'Los puntos (xo,A,) en los cuales la derivada - cambia de signo y
dX
~F(x,,
A,)
# O se denominan puntos regulares extremales.
dA
Un ejemplo de tal comportamiento nos lo da el mismo problema de
la carga de una barra. Asumamos que en el estado inicial la barra no es
idealmente recta, sino que est un poco pandeada. Entonces la dependencia
.de la flexin mxima respecto a la carga (fig. 5.2) cambia y toma la forma
indicada en la figura 5.5.
Si la barra se encuentra en un estado estable correspondiente a un
punto de la rama inferior del diagrama de bifurcaciones (fig. 5.5), y la
carga disminuye lentamente, entonces para cierto valor de P ocurre un
salto brusco hacia el otro estado de equilibrio en la rama superior. Este
fenmeno, denominado disparo, se utiliza en la tcnica. El cambio de tipo
del diagrama de bifurcaciones, por ejemplo, durante una violacin pequea
Fig. 5.4. Ejemplo de punto regular extremal. En este punto surgen de ninguna
parte dos estados de equilibrio
107
v0 # O,
AX
x - xO, FA(x0,+AO) # O.
(5.6)
Puntos singulares F, = O , F, = O
Puntos dobles
Punto singular
extrema1
Punto conjugado
Fzx ;e 0
2
2 = X(A - x2)
D<O
Notacin:
1) X < A,;
2) X = Ao;
3) X
> A,;
109
ax
o bien
Por las condiciones del problema, en el punto regular extremal tenemos que
@'(xO, A,)
dA
# O, y existe una curva X = X(x) cuya derivada cambia de
OX
ds
signo en el punto x,. Por consiguiente, en una rama, que parte del punto
(x,,A,), los estados de equilibrio sern estables (a< O), mientras que en la
otra (a > O), no lo sern.
Prestemos atencin al cambio del potencial en el punto regular extremal.
Para X < A,, la funcin V(x) representa una curva montona creciente; para
X = A, en ella aparece un punto de inflexin. Para X > Xo surgen un
mximo y un mnimo.
Llamaremos punto singular de la curva F(x, A) = O a todo punto (so,XO)
en el que
U n punto singular doble de la curva P(x, A) = O es un punto singular por
el que pasan dos, y slo dos, ramas de F(x, A) = O con tangentes diferentes.
En tal punto no todas las derivadas segundas se anulan simultneamente.
Un punto singular extremal de F(x, A) = O es un punto singular doble
en el que x, cambia de signo en una rama (fig.5.6).
Un punto conjugado es una solucin puntual singular aislada de
F(x, A) = O.
As pues, estamos analizado una situacin en la que todas la derivadas
primeras son iguales a cero y al menos una derivada segunda es diferente de
cero. Aclaremos para qu condiciones un punto es singular doble, singular
extremal o conjugado.
En un punto singular doble, la funcin F(s,A) se puede representar
como
a
Captulo 5
110
Fig. 5.6. Bifurcaciones surgidas en los casos cuando sus tipos estn determinados por los
trminos cuadrticos: o) punto singular extremal tenedor; b) bifurcacin transcrtica
((intercambio de estabilidad)); c) cambiando La direccin del parmetro X y el. signo del
segundo miembro, e l cuadro mostrado en la figura a) se convierte en la bifurcacin subcrtica
ilustrada en la figura c)
donde
Si D < O, entonces no existen tangentes que pasen por (xo,Ao) y este punto
es conjugado.
Si D > O, entonces del punto (x,, A,), como se aprecia en la frmula (5.10), parten dos ramas con tangentes diferentes. Por tanto, ste es un
punto doble (fig.5.6 b).
Sean D > O y F,, = O. Entonces FZA# O y la ecuacin inicial se puede
escribir en la forma
AA [ 2 k ~ , ~F ~=~0. ]
yk=-
X
A"
O) ,respectivamente.
2F,,(x,, A,)
A,)
= O
111
Sin embargo, G(vo,x,) = O para cada uno de los dos valores de v, que se
obtienen de la ecuacin (5.10). Adems,
Captulo 5
112
La primera direccin est relacionada con las generalizaciones a sistemas ms complejos, en particular, a sistemas espacialmente distribuidos.
El descubrimiento de las bifurcaciones elementales analizadas anteriormente, y de otras que veremos en el captulo 7, result ser muy importante en los
modelos matemticos que describen los procesos en los medios no-lineales.
de los
Esto permiti ver con nuevos ojos el mtodo de <<comprensin
procesos que tienen lugar en los sistemas complejos. Result que, frecuentemente, aqu es cmodo utilizar un anlogo del enfoque histrico, es
decir, averiguar qu sucede para los valores anteriores y posteriores de los
en el
parmetros, observar las particularidades cualitativas de los vecinos>>
espacio paramtrico. Comprenderel surgimiento de uno un otro rgimen,
hoy significa hallar la secuencia de bifurcaciones que conduce a su aparicin.
Se descubri la existencia de varios escenarios universales de transicin de
los regmenes ordenados elementales (a los cuales pueden corresponder
los puntos singulares o los ciclos lmites) al comportamiento turbulento
(catico) en los sistemas disipativos (a los cuales pueden corresponder los
atractores extraos).
Los numerosos experimentos computacionales que se han realizado
en las ltimas dcadas han mostrado que estamos ante un nuevo nivel de
unidad. A principios del siglo XX la unidad de la naturaleza, la cual se
manifiesta en la universalidad de los modelos matemticos, era vista por
los investigadores en el hecho de que un conjunto de los ms diversos
fenmenos se podan describir mediante ecuaciones lineales en derivadas
parciales de segundo orden. Hoy en muchos casos se escriben ecuaci~nes
diferentes, pero el mtodo de complejizacin del atractor, la secuencia de
bifurcaciones o, como se suele decir, el escenario de transicin al caos (o ruta al
caos), es el mismo. La unidad resulta estar relacionada no con lamaparicin
de ecuaciones anlogas, sino con un comportamiento cualitativo universal.
Uno de los superproblemas en esta direccin est relacionado con el
anlisis del surgimiento de la turbulencia hidrodinmica y en particular,
con la investigacin de la ecuacin de Navier-Ctokes. Esta ltima describe
el movimiento de un lquido viscoso incompresible:
. div = 0,
a,
v' = ij(x),
x E an.
113
En la ecuacin (ms precisamente, en el sistema de tres ecuaciones parablicas) figuran cuatro funciones incgnitas vl(x, t), v2(x,t), v3(x, t), p(x, t).
La condicin div v' = O permite eliminar una de ellas. El nmero de Reynolds desempea el papel de parmetro de bifurcacin. El trmino (l/R)A
describe la viscosidad. Frecuentemente, los valores fsicamente interesantes
del parmetro R se encuentran en el intervalo lo3-lo6. Por consiguiente,
estamos ante un ejemplo tpico de problema con parmetro pequefio en
la derivada de orden superior. Este problema es difcil de, analizar tanto
numrica como analticamente.- Las dificultades fundamentales estn relacionadas con el trmino grad . El lector puede analizar la ecuacin
vt vv, = O y verificar que su solucin puede no existir totalmente; esto
significa que el trmino viscoso de estabilizacin (1IR)A no se puede
botar simplemente.
En la investigacin de este modelo matemtico se han concentrado
grandes esfuerzos. Durante el estudio de algunas corrientes concretas
se lograron xitos impresionantes. La teora de las bifurcaciones ayud
a analizar detalladamente la prdida de la estabilidad en las corrientes
elementales. En los ltimos aos, la utilizacin de supercomputadoras ha
permitido obtener una gran cantidad de informacin sobre las propiedades
estadsticas de las corrientes turbulentas y sobre el comportamiento de los
lquidos cerca de las paredes de los tubos o canales.
No obstante, quedan sin aclarar muchos detalles fundamentales. Todava no podemos responder convincentemente a la pregunta de si juega
Dios a los dados en estos fenmenos o no. Para simular el comportamiento catico observado es suficiente el modelo propuesto u otras ecuaciones
hidrodinmicas consistentes en ecuaciones determinsticas en derivadas
parciales, cuyas soluciones son nicas y existen para O < t < m? O es
que acaso necesitamos de fuerzas aleatorias pequeas y debemos rechazar
terminantemente la descripcin determinstica de las corrientes turbulentas?
Otro campo clsico de aplicaciones de la teora de las bifurcaciones
es la teora de cscaras elsticas. Los problemas no-lineales de la teora de la
elasticidad, los cuales generalizan el problema ya visto de flexin de una
barra, resultaron ser muy importantes en muchos problemas de mecnica
de la construccin.
Un objeto mucho ms simple para el anlisis de bifurcaciones lo
constituyen los modelos tipo ((reaccin-difusin
114
Captulo 5
Fig. 5.7. Estructuras disipativas estacionarias: a) ejemplo tpico de distribucin espacial estacionaria de las sustancias en e l modelo bruselador: 6 ) ejemplo de diagrama de bifurcaciones.
La variable M puede ser u n mximo, una norma u otra magnitud que caracterice a La estructura surgida. Las lneas verticales no estn relacionadas con e l diagrama de bifurcacin. Su papel
consiste en indicar qu valores del parmetro le corresponden a Los puntos de bifurcacin
115
116
Capitulo 5
de A. J. Toynbee .ha sido lanzado un reto histrico, al cual se puede responder de diversas maneras. Por ejemplo, disminuir las necesidades, trasladar
los problemas internos al exterior y
A
comenzar a desarrollar una poltica
dura con los vecinos (Camino de Licurgo). A esta respuesta le corresponde la rama inferior en la figura 5.8 a
(punto Al). Una segunda respuesta
t
puede ser, por ejemplo, la colonizacin de territorios de ultramar, en
los que los Estados se encuentran en
A
una etapa inferior de desarrollo y no
pueden oponer una resistencia seria.
La siguiente eleccin (punto A,) puede consistir, por ejemplo, en dirigir todos los esfuerzos a convertirse en una
potencia comercial, o bien a asentarse
en ultramar en serio y por mucho
tiempo. Los trabajos de A. J. Toynbee
y L. N. Cumiliov contienen muchos
episodios donde el desarrollo march
en correspondencia con el escenario
mostrado en el diagrama 5.8 a. El diagrama de la figura 5.8 b puede corresA* A,
A3 X4 t ponder a una crisis de una sociedad
Fig. 5.8. Ejem p[os caractefisticos de b i- de COnSUmo" cuyos estndares de vidasonmuyelevados*
furcaciones que pueden surgir en problePero, tal vez, sera mucho ms
mas relacionados con e l anlisis histrico.
El punto de bifurcacin corresponde a una interesante e importante analizar y
situacin en la que en e l sistema aparecen predecir las situaciones mostradas en
vanas vas de desarrollo, O bien suceden 1, figura 5.8 c. ptecuadro corresponcambios catastrficos
de, por ejemplo, a la destruccin del
medio ambiente como consecuencia
de la utilizacin de tecnologas tradicionales de explotacin de los recursos
naturales, a la disminucin brusca de los estndares de vida y a la llegada
con el tiempo al nivel de los recursos renovables. Las dos curvas superiores aisladas (estable e inestable) corresponden a una nueva tecnologa de
explotacin de los recursos. Y aqu precisamente queda clara la utilidad
de los diagramas similares a los mostrados en la figura. Supongamos que
nosotros no nos imaginamos de ningn modo la curva de nuestro desarrollo
histrico. Entonces nos esperan catstrofes, calamidades y serias dificultades
en los puntos X3 y X4 (fig.5.8 c).
117
OF(x*,A*)
Y*
dx
Escribamos la solucin de esta ecuacin lineal en la forma y = e p t p.
Entonces, la estabilidad del estado de equilibrio est determinada por el
problema de valores propios para cierto operador lineal A:
y=
p p = Ap,
A r
OF(z*,A)
dx
Aqu p es un valor propio, cp es cierta funcin propia, y el operador se
reduce a la multiplicacin por un nmero. Hemos visto que las bifurcaciones
OF(x*,A*)
- = O (en caso contrario, se poda
tenan lugar en los puntos donde
dx
aplicar el teorema de la funcin implcita). Pero precisamehte stos son los
puntos en los que el estado de equilibrio pierde la estabilidad.
118
Captulo 5
119
.-
~ o d oesto condujo al desarrollo de una gran clase de mtodos computacionales de la teora de las bifurcaciones. En la literatura estos mtodos
frecuentemente se denominan mtodos de prolongacin respecto al parbmetro.
La idea de estos mtodos se puede ilustrar mediante 10s ejemplos siguientes.
supongamos que nos interesa la construccin numrica de un diagrama
de bifurcaciones. Asumamos, para simplificar, que este diagrama determina
la dependencia de los puntos singulares 5 del sistema dinmico = F(X, A)
respecto al parmetro A. Estos puntos estn determinados por la races de
la ecuacin algebraica no-lineal F(x, A) = O.
La estrategia tpica de accin es la siguiente. Utilizando
analticos o numricos, primero se halla algn punto (x*, A*) perteneciente
al diagrama (fig. 5.9).
r'
A*
A,
'
(A*
+ AA) = ''(A*) +
dx*(X*)
dA AA.
+ AA) = O
Captulo 5
120
.-
l a. I
1 1,
121
.-
En el problema
dado esta ecuacin tiene solucin 5 para cualesquiera
valores de a y b. Sin embargo, el nmero de soluciones puede cambiar si
se imvone la condicin
122 .
Capitulo 5
Fig. 5.10. Esquema del dinamo de dos discos y var-iacibn del atractor en el modelo
generalizado de Rikitake al variar e l parmetro p : v, = 0,004, v, = 0,002
123
124
Captulo 5
.-
125
burg-Landu
Preguntas y ejercicios
1. Construir el diagrama de bifurcaciones para el sistema dinmico
126
Captulo 5
4. La ecuacin diferencial
- y),
Bibliografa
Entre los textos de teora de las bifurcaciones ms accesibles y simples se pueden
mencionar:
looss G., Joseph D. Elementary stability and bifurcation theory Springer-Verlag,
1980;
Bautin N. N.,Leontvich E. A. Mtodos y procedimientos de investigacin
cualitativa de sistemas dinmicos en el plano (en ruso).
A los lectores ms preparados se recomiendan los libros:
Guckenheimer J, Holmes D. Nonlinear Oscillations, Dynamical Systems, and
Bifurcations of Vector Fields. Springer-Verlag, 1990;
Arnold V 1. Chapitres supplmentaires de la thorie des quations diffrentielles
ordinaires. MIR, 1980;
'Z861 '%epa~-~aSupds
*s~o~seyqv
aSuers pue 'soeq3 'suoyes~npg:suoyenba zuarol a u 3 mouvds Q
:orqg la ua ezpue as '(zuarol ap euraqs~)sepeuIpro
sale!suamaf!p sauopensa ap oqarsuo~euraJs!s un ap oldura[a la ua :os~~?urnu
s~sgeue
la u03 ugpeu~qurosua SauopesJnjIq se1 ap Toa4 el ap sopoq?~sol ap ugpesqd~e7
'LL61:
'alty
~ o [ *sura)sAsmnpq!l!nbauo~ u! uo~$?zye%y3-~la~
*laq808pd '-3
s!lo~i~
! ~ 8 6 1' 8 e l s ~ - ~ a % u p*ssga%~au&
ds
uayw Q
:sorqy sol ua uaqureq X sorqg sour;~lpsop sol ua ezgeue as 'sa(e~sred
SepeApap ua sauopen3a aumpaur uaqpssap as sapns sol 'sopp~qys~p
seuIa$s;s
sol ap s!syyue le sauopesmnj~qse1 ap ej~oaqq ap s o p o ~ ?sol
~ ap ugpeqlde e7
'6961 ' u I U I ~ ! U ~*SUIalq
~
-0rj anpua%!z reaug-uo~pue Lroau uogesrnjtg (spa) *S uazuluv nallg
!9~61% ~ I J ~ A
-~a%upds
*suogesgddvsq pue uogesmpg 3 d o a~ u *tyuayjar2 3~ '-1uapssvty Q
LZG
sauopo~~n.Jqsol ap y o a l ol ap soquawaly
129
.-
dinmico de la forma
W x ,A)
Entonces, como vimos, las soluciones acotadas tienden a los atractores x*,
los cuales estn determinados por la igualdad
130
Captulo 6
Fig. 6.1. Cuando e l centro de gravedad corta la lnea en forma de pico, La mecedora
puede voltearse de un costado a l otro
- 4
.-.-
[--a,
a].
+ t2)2,resulta
132
--
- .- - . -
- ---
Captulo b
-- -----
Derivaida e~ti~.-expwsi6n
e igualando la derivada a cero, obtenemos h
condicln dd &twmo del potencial:
Consiguieoitemmttc,los estados de equilibrio de la mecedora estn determinados por las races del polinomio
Fig. 6.2
133
Fig. 6.3. Potencial V ( t )en diferentes puntos del espacio de control (a, b): a ) nmero de
soluciones de La ecuacin U ( t )= O; b) forma caractedstica del potencial
Excluyendo la variable
t obtenemos
27a2 = 2(2b - I ) ~ .
(6-4)
Captulo 6
134
.-
135
"d
- Lanormalalaparbolaenelpunto (t,t2)
est dada por la ecuacin
2t3
+ t(1-
2 ~-) = 0.
136
Captulo 6
tA
V(a, b, t) = V(a, b, t
+ At).
137
plano cuando la luz de una lmpara atraviesa un vaso que contiene t y la azcar est
mal mezclada (para que las propiedades pticas se diferencien considerablemente en la
diferentes partes del vaso).
Las lneas custicas surgen tambin cuando los rayos de luz atraviesan un vidrio grueso con superficie compleja. Esto se utiliza
frecuentemente en el diseo moderno.
Ntese que la evoluta de una curva suave, la parbola, tiene-una singularidad en el
eje de simetra. Esto es algo tpico en los iige
El principio de Huyproblemas de ptica. El surgimiento de sin- gens permite
hallar la posicin
gularidades se puede observar claramente en del frente
de onda en instantes
el ejemplo de la propagacin de un frente de
consecutivos
onda. Recordemos el principio de Huygens: para hallar la posicin S(t At) de un frente en el instante t At, es necesario
tomar su posicin en el instante t, trazar la familia de normales y desplazar
el punto del frente a una distancia cAt a lo largo de la normal, donde c es la
velocidad de propagacin de las perturbaciones (fig. 6.8). Supongamos que
el frente de onda es una elipse y los rayos de luz se propagan al interior.
El cuadro caracterstico que nosotros veremos se muestra en la figura 6.9.
Aqu surgen singularidades notables.
La clasificacin de las lneas custicas tpicas y de las singularidades
de los frentes de onda constituye una de las aplicaciones ms importantes
e interesantes de la teora de las catstrofes. Es curioso. que el aparato
matemtico creado para el anlisis de fenmenos no-lineales tpicos, las
bifurcaciones y las-catstrofes, es muy til en las teoras lineales clsicas,
por ejemplo, en la ptica ondulatoria y
geomtrica.
En los aos 70 del siglo XX el problema de la distribucii~a gran escala de la
materia en el Universo despert la atencin de los astrofsicos. Result que esta
distribucin tiene singularidades geomtricas caractersticas. Frecuentemente las
Fig. 6.9. Posicin del frente de onda en diferentes instantes. Con el aglomeraciones de galaxias tienen forma
transcurso de tiempo, en ei frente
de hojuelas o platillos.
aparecen singularidades
Un grupo de investigadores del Instituto de Matemtica Aplicada M. V. Kldysh. de la Academia de Ciencias de la URCS, bajo la direccin de
Y. B. Zeldvich, realiz una investigacin terica de los procesos que pue-
138
Captulo 6
139
1
a/2 (lo cual se deduce del
Utilizando el hecho de que JG
desarrollo en serie de Taylor del radical), resulta
140
Captulo 6
141
.-
Sin detenernos en el clculo de las constantes a,, a,, a,, a,, a,, sealemos
que en el punto P se tiene a = ,6 = O, y
+ - permite
4a4
eliminar el trmino cbico. (Aqu los trminos a3Pe3 a,@* pasan a formar
parte del trmino x4.) Cambiando la escala de los parmetros de control,
ct -+ Ca, p + DP, la expresin para el potencial se puede reducir a la forma
Por cuanto en este problema nos interesan slo los puntos crticos, podemos
desplazar el origen del sistema de referencia de la funcin V y suponer que
c = O. Esto nos proporciona
. l 4
2
V,,(x)= -x
+ -ax
+ bx.
4
2
W(t) = 4
(1- 2b)t2
2
+ at.
Captulo
142
a
7
,
~ ta' =, (1 - 2 b ) a , bt =
x4 a'x2
reduce este potencial a la forma Wab(x)= - + -+ b'x, que coincide c o ~
Sin embargo, el cambio de variables x =
4
2
la frmula (6.7).
El surgimiento de unas mismas catstrofes en problemas tan diferente4
es asombroso. En calidad de ejemplo se puede citar un experimento de 12
psicologa de la percepcin.
En la base de este experimento yacen las siguientes ideas. En la teora
de control, para el seguimiento de una seal se utilizan ciertos sistemas
dinmicos, los cuales pueden describirse, por ejemplo, mediante ecuaciones
del tipo (6.1) o relaciones ms complejas, donde los parmetros Xi descri'
ben ciertas caractersticas del objeto percibido, mientras que la variable z
(o un conjunto de variables) describe las propiedades de la imagen surgida,
Se puede suponer que el cerebro acta de un modo anlogo, formando
cierto sistema dinmico que permite seguir al objeto. Cuanto ms simple sea
el problema a resolver tanto ms simple debe ser el sistema que suja. Aqui
tambin estamos ante un anlogo de la autoorganizacin: de un nmero
grande de grados de libertad que posee el cerebro se elige slo un nmero
pequeo. Pero, entonces, nosotros tendremos que observar efectos cualitativos de la percepcin que pueden ser representados por la teora de los
sistemas dinmicos y, en particular, por la teora de las catstrofes. En favor
de esta tesis se puede presentar la ilusin ptica mostrada en la figura 6.11.
.
Entre las figuras mostradas, la cuarta de la fila superior se percibe con
igual probabilidad como la cara de un hombre y como la figura de una mujer.
Si esta figura est incluida en una sucesin (fila superior en esta figura), entonces la percepcin de las figuras centrales se desplaza en dependencia del
orden en que se mire esta fila, de izquierda a derecha o viceversa. Se observa
Hombre
143
una biestabilidad, es decir, dos posibles respuestas para unos mismos valores
de las variables de control, y una histresis. Ambos fenmenos ya han sido
analizados al considerar el problema de la mecedora. Introduzcamos un
parmetro ms, que determine el grado de detallamiento de la imagen. ste
cambia a lo largo del eje de las ordenadas. Experimentando con esta figura, el
lector puede comprobar que la regin de biestabilidad, donde la figura se interpreta como la cara de un hombre o como la figura de una mujer en dependencia del orden en que se mire, se encuentra en el interior del pico.Esto
permite suponer que aqu tambin tenemos un caso de cspide, igual a la del
caso de la mecedora y de la mquina de Zeeman. Es natural considerar que el
cambio de percepcin, el cual sucede generalmente en la frontera del p'ICO,
esta relacionado con saltos de catastrficos en cierto sistema dinmico.
Sin embargo, la respuesta a la pregunta sobre cmo sucede el proceso
de autoorganizacin y qu mecanismos aseguran el surgimiento y la seleccin de los parmetros de orden, exige esfuerzos conjuntos de psiclogos y
matemticos.
La teora de las catstrofes permite comprender en cierta medida por
qu en los ms diversos problemas surgen unas mismas catstrofes. Resulta
que las catstrofes tpicas estn determinadas por el nmero de parmetros
de control en el sistema analizado. Cuando hay un solo parmetro de
control, a, la nica catstrofe posible es el pliegue, la cual est determinada
por el potencial
((
m3
U,(%) = 3
+ ax;
cuando los parmetros de control son dos, son posibles un pliegue (6.8) y
una cspide (6.7).A medida que aumenta el nmero de parmetros, es decir,
al aumentar el nivel de complejidad del sistema el nmero de posibles catstrofes crece y la superficie de respuesta se vuelve ms caprichosa. Uno de los
creadores de la teora de las catstrofes, R. Thom, mostr que si el nmero de
variables de control no es mayor que cuatro, entonces en las situaciones tpicas son posibles slo siete catstrofes elementales. Este resultado se puede
considerar como una realizacin del programa de investigacin de Poincar
para los sistemas gradientales (6.1). Limitando el tipo de modelos matemticos
a una clase pequea de sistemas dinmicos tipo (6.1),as como su anlisis
al estudio de los entornos de los puntos singulares, los investigadores
obtuvieron la posibilidad de elaborar una clasificacin de las catstrofes.
Al analizar la mquina de Zeeman nosotros tropezamos varias veces
con la necesidad de desechar los trminos superiores de los desarrollos en
series. Intuitivamente es claro que precisamente eso era lo que se deba
hacer, pero la teora de las catstrofes permite fundamentar rigurosamente
dicho procedimiento. Esto es especialmente importante si el problema es
multidimensional y muchas consideraciones evidentesresultan falsas.
'
+
+ d(u + p) + e(u + q ) + f = au2 + buv + cv2 + u(2ap + bq + d ) +
+ ~ ( ~ V + U ~ + V ~ + ~ ~ ) + C ( V ~ + ~ ~ V + ~ ~ )
+v(2cq+bep+c)+ ( f +ap2+bpq+cq2+dp+eq).
(6.1
145
(6.12)
El coeficiente g es igual a la expresin entre llaves en (6.10). Geomtricamente, la transformacin realizada significa que hemos trasladado el origen de
coordenadas al centro de la figura analizada (fig. 6.12). La expresin (6.12)
es invariante respecto a la transformacin (u + -u, v -+ -u). Reduzcamos
laforma cuadrtica au2 buv cv2 a una forma an ms simtrica. Elijamos
las nuevas coordenadas de modo que (6.12) sea simtrica respecto a los
nuevos ejes de coordenadas.
Esto se puede lograr por dos mtodos. Uno de ellos se utiliza generalmente en la teora de las catstrofes, mientras que el otro, en la geometra
analtica. Veamos ambos enfoques.
Por cuanto el determinante es diferente de cero (condicin (6.11)), entonces al menos uno de los i~merosa, b, c es diferente de cero. Supongamos
que este nmero es, por ejemplo, uno de los coeficientes de los trminos
cuadrticos, digamos, a. Completemos un cuadrado perfecto:
Captulo
146
4
-4
donde r = (&x+Jq),
eda
A=d- -
JE'
e
eL
donde s = r + - g = f 4 c .
2Je'
Puesto que # = 4ac, entonces a y c deben ser del mismo signo,
Anteriormente supusimos que a > O, c '> O. En caso contrario se pueden
dividir todos los trminos por a (si a # O). Esto reduce el problema al caso
anterior.
Si a = O, entonces tambikn b = O, y la transformacin se simplifica:
y2+dx+ey+ f = O
As pues, los cambios de variables que acabamos de describir, relacionados con la completacin de un cuadrado perfecto, permiten reducir la igualdad (6.9) a la forma cannica Ar2 Bs2 = G . Si A, B y C tienen el mismo
signo, entonces sta es la ecuacin de una elipse. Si el signo de G no coincide
con los signos de A y B, entonces. esta ecuacin no determina ninguna
curva en el plano (r, S ) . Si A y B tienen signos diferentes, entonces tenemos
la ecuacin de una hiprbola. Para b2 = 4ac la ecuacin (6.9) corresponde a
una parbola. Se puede cambiar la escala y reducir la frmula (6.9) a la forma
p2 + q2= 1,
p2 - q2= 1,
p - q 2 =o,
con lo cual se pierde informacin. Por ejemplo, todas las elipses, independientemente de la magnitud de sus ejes y de la posicin que ocupen en el
plano, se transforman en una circunferencia de radio unidad y centro en el
origen de coordenadas (fig. 6.12).
La transformaciones expuestas cambiaron la escala a lo largo de los ejes
de coordenadas y desplazaron el origen de coordenadas. Dichas transformaciones cambiaron los ngulos, convirtiendo el sistema de coordenadas rectaiigulares en un sistema de coordenadas oblicuas. El lector puede comprobar
147
+ u' sen a.
= a(ul cos a - v' sen al2 b(ul cos a - v' sen a)(vf cos a
+ u' sen a) +
Captulo B
148
q(x) =
X,xixj*
ij
Xli
..
< <
<<
= Xll
( x: + 2
n
= ~ i [i ( z l
149
'
+ ~ p i j ~ -j ~) ( p l )j 2 + no contienen a 2,
j=2
j=2
)=
Yl = X 1
kjxjf
Yi
= xi (i > l),
donde q es una forma cuadrtica de las variables y,, .. ,y,. Por eso, se
puede seguir realizando este procedimiento, aplicndolo primero a la forma
cuadrtica 4, luego a la forma de las variables y3, . .,yn, etctera. Como
resultado de esto la forma se reduce a la forma diagonal
q(y) = dly:
+ + dnyi*
* * e
Yk
Modifiquemos la escala mediante el cambio de variables zk = -
JWY
< 7-8.
150
Captulo 6
"
.-
+ v2
= 2,
= 2);
2) u2 - V 2
(r = 2,
= O);
3) - u2 - v2
(r = 2,
= -2);
4) u2
1,
= 1);
5) - u 2
(r = 1, S = -1);
6) 0
(r = O,
= O).
1) u2
(Y
151
+ py + + 6y3 = 0.
Fig. 6.14. Posibles posiciones de las rectas de races para la forma cbica
Capitulo 6
152
= uv(u - v).
Hacemos ahora el cambio de variables
C(U, v)
+
+
+
+
+
+
Entonces, t2
+ s2= u2 + v2, y
153
.-
d = uv(v - U ) ( W - (YU).
'
donde yl,y2, y3, y4 SO^ las ordenadas de los puntos de interseccin con la
recta x = l. Si una. de las rectas coincide con el eje y (y, + m), entonces
c=
- Y3
Captulo 6
154
.-
155
a%?),
156
.-
Captulo 6::
157
.-
x u n vector m-dimensional, y un vector n-dimensional y la funcin f , u n vector p-dimensional. Supongamos que para cierto punto (5,
y) E IRm x IIB" el conjunto
(donde D f l(z,F)es la aplicacin lineal (6.16) determinada por el jacobiano
tomado en el punto (Z, y)) es el grfico de la funcin y = y (x). Entonces,
tambin para la funcin no-lineal inicial la igualdad
-x 2 =
+C.
158
Capitulo 6
la igualdad (6.17) habr una funcin negativa y en el segundo miembro, una funcin
positiva. Esta contradiccin demuestra que la aplicacin y(x) no existe.
Ejemplo 2. Las funciones xk y xm no son equivalentes por la derecha para k # m.
Supongamos nuevamente que y(0) = O, y que, por ejemplo, k < m. Si estas funciones
son equivalentes, entonces se debe cumplir la igualdad
xk = (y(x))" C.
Analizando el punto x = O de igual manera que en el ejemplo anterior, llegamos a
la conclusin de que C = O. Derivemos ambos miembros de la ltima igualdad:
m-1
kxk-' = rny(x)
y (x),
m-1
rl
y (x),
B f (0)
6xi
xi dt.
O
i=l
159
obtenemos f =
k: = O
160
Captulo @*
9
k!
llk
x = 0,
x # O.
161
. . o f
ax1
03,
af-(0) - -- af (0) = o,
(3x1
axn
se denominarn puntos crticos. El valor de la funcin en un punto crtico se
llama valor crtico de la funcin,
Geomtricamente, la condicin formulada significa que en un punto crtico el plano tangente (recta tangente, para n = 1) es horizontal. En el caso unidimensional los
puntos crticos pueden ser mnimos
e . .
Af
"J
5
< O)
ox2
o puntos de injlexibn
163
canal, son aislados, es decir, en sus entomos no hay otros puntos crticos.
En el caso del canal, el punto crtico es, en realidad, toda una recta. Si el
potencial est determinado por un canal, entonces el punto del espacio de
fase que determina el estado del sistema se encuentra en equilibrio indiferente.
Son posibles tambin otros tipos ms complejos de puntos singulares.
Parece natural que en los casos tpicos que
se presentan en la simulacin de diferentes sistemas, deban surgir puntos crticos elementales
(mximos, mnimos, puntos de silla). La afirmacin denominada lema de Morse formula las condiciones suficientesde surgimiento de tales puntos crticos. De hecho, este lema generaliza las
condiciones suficientes de mximo y mnimo de
una funcin suave y = y(x) para el caso de varias variables. El teorema de Morse permite, ade- Fjg- 6-22. fk~ncincanal. Aqu
ms, reducir a la forma cannica una clase deter- 10s Puntos crticos constituyen
toda una recta
minada de funciones cerca de los puntos crticos.
Introduzcamos la definicin siguiente.
Definicin. La funcin f tiene en el punto u un punto crtico no-degenerado si D f (u) = O y la forma cuadrtica
I
o2f (u) =
2.2
a2f(u)
jaziaxjf
determinada por las derivadas segundas de la funcin f , es no-degenerada.
En otras palabras, el rango de esta forma cuadrtica es igual a n o la
matriz
i,j
* * *
< <
Captulo 6
164
Sf (0)
dx;
Por consiguiente,
f (x) =
xixjh,(x).
Procederemos de la misma manera que cuando redujimos las formas cuadrticas a formas cannicas. Cambiando la numeracin de las coordenadas
o girando el sistema de coordenadas, si es necesario, logramos
que
el
elemento diagonal superior sea diferente de cero, Hll # 0.
Tomemos la funcin g(u,, ... ,u,) = JIH~~('IL~, . ,un)I y hagamos el
cambio de variables
vi = ui para i 2 2.
165
crtico
en el origen de coordenadas y que en este punto la matriz de Hesse tiene
Capitulo 6
166
donde
r:
2
fz,
f**.f
Rn-'
2,
- ,..., x,),
+ f(z,+,
+ IR es cierta funcin
suave.
167
.-
Para
< 0, x =
Captulo 6
168
4
A
f(z,A) = const
s (Y,P) = f (~A(x),
P ( ~ )+)
= conSt
169
.-
propiedades t {picas.
Las definiciones formuladas permiten enunciar un resultado fundamental de la teora de las catstrofes, el teorema de clasificacin de Thom.
Teorema de Thom. En el caso tpico, la familia r-paramtrica de funciones
suaves Rn + Rn para todo n y todos r 5, es estructuralmente suave y
<
<<
o
o
o
o
+ Alul + ( M ) ;
A, cspide: f(u;' + A,u: + X1ul) + (M);
A, cola de milano: u: + X,u? + X2u: + Xlul + (M);
A, mariposa: f(u! + X,u; + X3u: + X2u: + X1ul) + (M);
A, tienda india (wigwam): u: + X,U; + X4ut + X3u: +
+
X,u, + ( M ) ;
A, pliegue: u:
X,U:
catstrofes umbilicales:
o
o
o
o
Ap,)
+(N);
X1ul
+(N);
+ X4u2 + X,U: +
segundo ombligo hiperblico: u:u2 + u; +
X2u2+ X,u, + ( N ) ;
E, ombligo simblico: *(u; + u: + X ~ U , U ; + X ~ U +
; X3u1u2+ X2u2+
o D:
o
X1u,)
+ (N);
Captulo 6
170
.-
< <
<
< <
< <
catstrofes
171
Fig. 6.24. En muchos problemas fsicos, los parmetros para Los cuales ocurren los saltos
catastrficos estn determinados no por el principio de retardo mximo (o), sino por e l
principio de Maxwell (6)
del conjunto de bifurcacin en la que nacen o desaparecen los mnimos locales (como en el problema de la mecedora). Se puede decir que en este caso
el sistema aminimiza localmente a la funcin potencial. El comportamiento
exacto se ilustra en la figura 6.24a. Sin embargo, si en el sistema investigado
el nivel de ruido es suficientemente alto (una descripcin detallada de las
fluctuaciones exige un modelo ms complejo que la ecuacin (6.1)), entonces
el sistema tiene la posibilidad de estimarno slo el mnimo local donde se
encuentra, sino tambin los valores del potencial cercanos donde pueden estar sus mnimos. En este caso puede ocurrir una minimizacin no-local, sino
global. El estado del sistema bajo estas condiciones estar determinado por el
mnimo global de la funcin potencial. El comportamiento tpico en este caso
se muestra en la figura 6.24 b. Esta aproximacin, utilizada frecuentemente
en la teora de las transiciones de fase, se denomina principio de Maxwell.
En la descripcin de las transiciones de fase de segundo gnero es
utilizado ampliamente el modelo Gunzburg-Landu. Este modelo est dado
por el potencial
172
Capitulo 6
el potencial
1 4
1
2
V(x,a, b) = -x
-ax
bx.
4
2
Este potencial coincide con una catstrofe cspide, y los efectos no-lineales en tal sistema sern los mismo que en el problema de la mecedora
o de la mquina de Zeeman. En otras palabras, el problema se redujo a lo
estudiado anteriormente. Algunas diferencias pueden surgir si, en virtud de
las particularidades fsicas del problema, se debe partir no del principio de
retardo mximo, como sucede habitualmente, sino del principio de Maxwell.
Al describir las transiciones de fase de primer orden a veces es necesario
considerar un potencial an ms general
1
1
1
V(x,a, b, c,d) = -a: 6 -ax4
-bx3
-a2
dx.
6
4
3
2
Mas este potencial determina simplemente una catstrofe mariposa As, cuya
descripcin detallada se expone en los textos de estudio de teora de las
catstrofes.
De esta manera, en el primer nivel se logra hacer corresponder los
modelos matemticos de diversos sistemas no-lineales con los modelos ya
estudiados de la teora de las catstrofes. En una serie de casos esto permite
descubrir efectos nuevos en muchos campos.
Se debe anotar que en muchas disciplinas, en particular en la termodinmica, la teora de la elasticidad, la fsica de los lseres y la teora de
estabilidad de las naves, los problemas fundamentales fueron resueltos sin
utilizar la teora de las catstrofes.
Aun as, la teora general permite formular las regularidades descubiertas en una forma ms simple y clara. Adems, surge la posibilidad de
destacar los rasgos cualitativos generales de diferentes sistemas no-lineales.
Este tipo de trabajos se puede atribuir al segundo nivel.
El tercer nivel, el ms profundo e interesante, est relacionado con el
anlisis de objetos ms complejos, los cuales obligan a salirse del marco
de la teora elemental de las catstrofes, exigiendo frecuentemente tanto la
creacin de nuevos modelos como el desarrollo del aparato matemtico.
Preguntas y ejercicios
1. Analizar el conjunto de bifurcaciones para la catstrofe cola de milano
173
2. Para describir las transiciones de fase Van der Waals propuso la ecuacin de
estado fenomenolgica
rpidamente oscilantes
u(x, r) =
eiTl(xn)a(x,
a, r) da,
174
Captulo 6
7. En ciertos casos las abejas salvajes no forman colonias. En otros casos forman
colonias en las que viven decenas de miles de ejemplares. Prcticamente no
hay variantes intermedias. Explicar este fenmeno suponiendo que la densidad
de abejas se debe mantener en un nivel de B abejas/krn2. Considerar que
la colonia vive en una regin circular de radio R. El nmero de abejas es
N z ~ R ~ElBtiempo
.
que gasta una abeja en un viaje de ida y vuelta desde el
A
R. Asumir que la productividad de una abeja es P(R)= c - BR,
dS
y que la ganancia S ( R ) debida a la cooperacin es una funcin tal que - > O
dR
dS
para R = R,, y lim - = O cuando R + OO. Analizar la funcin de utilidad
dR
F ( R ) = P(R) S(R) y determinar el tipo de catstrofe.
centro es
u = o (x2
Bibliografa
Una introduccin elemental a la teora de las catstrofes y sus aplicaciones, con
muchos ejemplos e ilustraciones, se presenta en el libro:
Posto~zT., Stewart 1. Catastrophe theory and its applications. Pitman, 1978.
La exposicin en este libro comienza con las nociones fundamentales de lgebra
lineal, geometra analtica y anlisis matemtico. Se consideran las aplicaciones a
la teora de estabilidad de los barcos, ias corrientes de lquidos bidimensionales, la
ptica geomtrica y-ondulatoria, la teora de la elasticidad, la fsica de los lseres y
un crculo grande de problemas relacionados con la simulacin blanda.
Un libro de texto ms complejo, el cual revela las particularidades matemticas de
la teora y aplicaciones importantes a problemas de la aerodinmica, la climatologa
y la mecnica cuntica, es
Gilrnore R. Catastrophe Theory for Scientists and Engineers. Wiley 1981.
La particularidad de este libro consiste en que en l se analiza la transicin de la
teora de las catstrofes a la teora de las bifurcaciones, la cual permite investigar
modelos matemticos ms complejos, que no se pueden describir mediante una
funcin potencial. En el libro
Arnold V; 1. Teora de catstrofes. Alianza Universidad, 1989.
Se consideran las aplicaciones de la teora de las catstrofes al problema del anlisis
de la distribucin a gran escala de la materia en el Universo, los problemas de
optimizacin y el problema de evasin de las dificultades. En el libro tambin se analizan
las relaciones profundas de la teora con otras reas de la matemtica moderna.
Enfoques matemticos ms rigurosos se proponen en los libros:
175
.-
Sistemas elementales
con tiempo discreto
cada poca hay que darle lo que se merece! Los tiempos nuevos tambin han alcanzado
logros brillantes.
jA
R. Shumann
177
La investigacin de las aplicaciones unidimensionales permiti introducir nuevos conceptos, aplicables a una clase grande de sistemas disipativos,
descubrir una serie de fenmenos nuevos y responder a varias preguntas
fundamentales. Cmo ocurre la transicin de los regmenes ordenados
ms simples a los regmenes caticos al variar los parmetros? Cmo se
suceden en el espacio de parmetros las regiones en las que se observan el
orden y el caos? Cmo ocurre la complejizacin del orden cuando varan
los parmetros? Cmo son los tipos elementales de caos en los sistemas
no-lineales abiertos y cules son los mtodos para describirlos?
Veamos algunas situaciones tpicas donde las aplicaciones unidimensionales participan como modelos matemticos, o bien surgen en el proceso
de resolucin numrica de las ecuaciones;
Aplicaciones como modelos de procesos con tiempo discreto. Anteriormente analizamos las propiedades de la ecuacin diferencial no-lineal
denominada ecuacin logstica. En este modelo el tiempo t se considera
continuo y puede tomar cualquier valor real en el intervalo O t < OO.
El tamao poblacional x(t) se considera una funcin diferenciable y continua
de este argumento. En realidad, nosotros asumimos implcitamente que
tenemos la posibilidad de medir la magnitud x(t) en cualquier instante y
compararla con las predicciones del modelo.
Pero es posible una situacin completamente diferente. Supongamos
que existe la posibilidad de estimar el tamao poblacional slo una vez al
ao. Entonces, es natural considerar a la variable temporal como una variable
discreta, no-continua. Por ejemplo, ella puede tomar slo valores enteros
n = 1,2,. . . . Entonces, el tamao poblacional se expresar mediante una
funcin de argumento discreto x(n), o utilizando la notacin habitual, x,.
Para simplificar denotaremos la sucesin x,, x2, ... mediante {x,).
Se puede suponer que x, es el tamao poblacional en el ao nmero n.
Por lo visto, entre los nmeros x, hay 'cierta regularidad. Es natural esperar
que el tamao poblacional en el ao x,+, dependa de cuntos individuos
haba un ao atrs, es decir, de la magnitud x,. De este modo, en el caso ms
simple (cuando x,+~depende slo del tamao poblacional del ao anterior y
no de x,-,, xnm2,
etctera.), obtenemos un modelo matemtico del tipo (7.1).
En este modelo f es una funcin unvoca continua de sus argumentos, X es
un parmetro que depende del problema concreto analizado, y x' (valor
inicial) es el primer trmino de la sucesin {x,). Frecuentemente se utiliza
una funcin f de la forma Xx(N - z):
<
178
-
Captulo 7
.-
2,
o < X;
< l.
(7.3)
De esta manera, hemos 0btenid.o nuevamente la aplicacin unidimensional (7.3). Como veremos ms adelante, sus propiedades son similares a
las propiedades de la ecuacin diferencial slo para valores pequeos de 7.
Cuando el paso temporal T supera cierto valor crtico, el comportamiento
de estos dos objetos comienza a diferenciarse cualitativamente.
179
.-
Aplicaciones elementales que surgen al resolver numricamente ecuaciones algebraicas no-lineales. Un mtodo estndar de resolucin numrica
de las ecuaciones algebraicas no-lineales de la forma
es la construccin de una sucesin (2,) convergente a su raz z*. Esto se
puede hacer utilizando el mtodo de iteracin simple (fig.7.1 a)
o bien el mtodo de Newton (fig. 7.1 b)
180
Captulo 7
Fig. 7.2. Una aplicacin unidimensional puede surgir a l analizar datos experimentales: o) determinacin de los mximos locales en la dependencia z(t);b) forma tpica d e la aplicacin
surgida durante e l anlisis de las reacciones qumicas oscilantes con comportamiento catico
181
Construyamos el grfico de
la recta y = x y de la curva
y = f (x) para el valor elegido
de X (fig. 7.3). Sealemos x, en el
eje de las abscisas y tracemos una
vertical desde este punto hasta
cortar la curva y = f(x) (punto A); luego tracemos una horizontaldesdeelpunto Ahastala Fig.7.3. Representacin grfica de La sucesin
interseccin con la recta y = x {x,) para La aplicacin Logstica en e l caso
x, -+ 0 cuando n -+ m
(punto B). Ahora tracemos una
vertical desde el punto B hasta el
eje x y denotemos el punto de interseccin mediante x2. NO es difcil verificar que x2 = f (x,) Tomando el punto x2 como punto inicial y repitiendo el
procedimiento anterior obtendremos x,, despus obtendremos x,, etctera.
Este procedimiento, denominado construccin de la escalera de Lamerey, permite hallar grficamente los trminos de la sucesin {x,). En la figura 7.3 se
observa que x, -+ O cuando n -t m. A partir de la frmula (7.3) se deduce
que la funcin f (x) transforma el segmento [O, 11 en el segmento [O, X/4].
Si X 4, entonces todos los valores x, pertenecen al segmento [O, 11 siempre
que O x, 1. Precisamente por esto se dice que la frmula (7.3) define
una aplicacin de un segmento en s mismo.
Supongamos ahora que X es un poco mayor que la unidad. En este
caso, la sucesin {x,) se comporta de otro modo (fig. 7.4): la sucesin {x,)
tiende a un valor constante x* > O. Aplicado al problema biolgico inicial, esto significa que el-tamao
poblacional de la especie se estabiliza al cabo de varios aos, y con
el tiempo deja de cambiar.
'El valor de x* se puede hallar
a partir de la ecuacin
<
< <
x* = f (x*,A).
(7.4)
Todos los puntos que satisfaFig. 7.4. Funcin f (3) y primeros trminos de
cen esta ecuacin se denominan
La sucesin {xnlfX n -f Z* cuando n
m
puntos fijos de la aplicacin, pues
x, = x* para todo n si 2, = x*.
de las iteraciones de la
En la figura 7.4 se muestra un cuadro tpico
aplicacin unidimensional en un entorno del punto fijo x*. Para n grandes,
en el entorno de este punto las magnitudes x, - x* tienden a una progresin
geomtrica:
x,
0,2
0,4
0,6
0,8
1 x
-f
182
Captulo 2
.-
Sealemos que todos los puntos fijos yacen en la interseccin del grfico de
la funcin x,+~= f (x,, A) con la bisectriz del ngulo del primer cuadrante
x,+, = x,. En efecto, el punto de interseccin de estas lneas satisface la
definicin de punto fijo, x,+, = f (x,, A) = x,.
~x~erimentando'con
la escalera de Lamerey o bien analizando la aplicaO f (x*)
cin elemental x,+~= TX,, en la que el cero es el punto fijo x* y
dx = 7,
se puede comprobar que si fx(x*,A) > O, entonces los elementos de la
sucesin {x, ) convergen montonamente a x* para 1 > f x (x*,A) > 0, divergen para f,(x*, A) > 1; y si f,(x*, A) < O, entonces los elementos de la
sucesin (2,) oscilan alrededor del punto fijo (fig. 7.1 a), convergiendo si
O > f (x*,A) > -1 y divergiendo si fx(x*,A) < -1.
Para A < 1 la ecuacin cuadrtica x* = Ax*(l - x*) tiene una raz
no-negativa x* = O. Para X > 1 tiene dos races no-negativas: X* = O y
A-1
x* = -. Para A = 1 el punto fijo x* = O pierde la estabilidad, y el
X
punto que aparece se convierte en el punto estable.
No es difcil determinar la estabilidad del punto fijo x* de la aplicacin
f (x, A). Sea x, = x* Ax,, donde Ax, es un nmero pequeo. Si el punto es
estable, entonces al crecer n la magnitud IAx,l debe disminuir. Escribamos
la frmula (7.1) en la forma
x*
A)
Ax,.
+ af (x*,
dx
Pero esto no es otra cosa que una progresin geomtrica habitual. Por
consiguiente, lAx,l -+ 0 cuando se cumple la desigualdad
183
superes fables
Continuemos aumentando el parmetro A. El comportamiento del sistema cambia nuevamente: a partir de valores de n suficientemente grandes,
en la sucesin {x,) se alternarn dos nmeros: a, y a,. (Hablando ms
exactamente, la sucesin {x,) es tal manera que x,,+, + a,, x,, -+ a,
cuando n -+ co.) Estos nmeros estn relacionados mediante las expresiones
a, = f (a,), a, = f (a,). En este caso se dice que la aplicacin (7.3) tiene un
ciclo estable de perodo dos y se denota con S2,y los nmeros al y a, se denominan elementos del ciclo. En la figura 7.5 se muestra la forma de un ciclo S'.
La transicin del punto fijo (se puede considerar como un ciclo S')
al ciclo S2 tuvo lugar como resultado de una bifurcacin denominada
bifurcacin de duplicacin del perodo. En este caso, el punto x* no desapareci,
- Of (x*,A)
se hizo menor que -1.
pero la magnitud
ax
Con el aumento ulterior de X la sucesin {x,) cambia nuevamente.
Surge un ciclo S4 (fig. 7.6):
x~~ -+ a,,
x4,+,
-+
a,,
x4,+, -+ a,,
x4,+,
-+
a,
cuando m -+
s2'+'
184
Captulo 7
.-
Fig. 7.8. Forma de la funcin f2(x) para la aplicacin logstica cuando f (z)determina
ciclos superestables: a) ciclo S2;b) ciclo S4
185
+
p veces
1 1
Captulo 7
Fig. 7.9. Complejizacin de Los ciclos estables de la aplicacin x,+, = Ax,(lcomo resultado de Las bifurcaciones de duplicacin del perodo
3,)
con A,, A,, A,, . . . los valores del parmetro X para los cuales ocurren las
duplicaciones, y con Al, A,, A,, ... los valores del parmetro para los cuales
x = 112 es uno de los elemento del ciclo S2,S4,SS, etctera. Tales ciclos
se denominan superestables (el nombre queda claro de la desigualdad (7.8)).
Introduzcamos tambin las magnitudes d,, d,, ... ,d,, .. ., iguales a la distancia entre x = 112 y el elemento del ciclo s2'ms cercano a l para
A = A2'. Todas estas notaciones se aclaran en la figura 7.9. De esta figura se
ve que dn = f:"31/2) - 112.
Los clculos realizados en computadoras muestran que para valores
grandes de n los nmeros A, y A, se comportan como progresiones geomtricas. Su razn est dada por la constante 6 = 4,6692016.. . . En otras
palabras,
187
Existe una analoga palpable entre los puntos fijos de las aplicacio
unidimensionales y los objetos de la teora de la universalidad.
En efecto,
(-
1
= niim
+o(-a)(-c~)~jf:+~
(-a)m
Para demostrar que existe una constante universal 6 y que los valores A,
forman una progresin geomtrica, se puede razonar del modo siguiente.
Supongamos que despus de desplazar la coordenada x en 112,
la coordenada del vrtice coincide con x = O. Consecuentemente, los
valores A,, para los que hay ciclos superestables, satisfacen la ecuacin
')
189
f (x,
d f (G A,)
donde 6 j ( x ) z
dX
igualdad proporciona
+ ( A - X,)6f
= f (x,A),
(x),
+ ( A - X,)L(
T f (x,A) = T f (x,A,)
(2,A,))6
+ O ((6f
)2).
+ ( A - A,)LT~-'J(X,A,).
.. +
+ Lf (2,A,)6f + 0 ( ( 6 f )2)
Sin embargo, por cuanto
i+oo
n 2"
(-a)f
A,)
(-,x
An+i) z
= T" f ( x ,A,)
para n B 1,
para
T n f(2,A,)
g(x)+ - A,)p:(O,(x)c,.
Recordemos que g(0) = 1, y que para los valores A, correspondientes al
ciclo superestable
s2",tenemos que
Tnf(O,An)= ( - ( y ) n f
2n
(O,An) = O.
190
Captulo 7
-
Consiguientemente,
De este modo, los nmeros A, forman una progresin geomtrica con una
misma razn p-', independientemente de la funcin f . De aqu queda claro
que el valor propio pl coincide con la constante universal de Feigenbaum 6.
191
192
Captulo 7
2 = XY - bt,
193
194
Captulo 7
anterior del parmetro), entonces puede resultar que en el sistema hay varios
atractores y un cambio del parmetro X traslada el punto xl de la regin
de atraccin de un atractor a la regin de atraccin & otro. Para aceptar o
rechazar esta hiptesis es necesario saber cuntos atractores puede tener el
sistema investigado y cules son sus regiones de atraccin. Incluso para ias
aplicaciones bidimensionales es muy difcil dar respuesta a estas preguntas.
Afortunadamente, para la aplicacin logstica la situacin es suficientemente simple. La pregunta sobre el nmero de atractores fue respondida
por D. Singer para la denominadas aplicaciones S-unirnodales. El modelo (7.3)
tambin puede ser reducido a estas aplicaciones mediante un cambio
de variables.
Una funcin se denomina S-unirnodal si:
1) es continua;
2) f(0) = 1;
3) f es estrictamente decreciente en [O, 11 y estrictamente creciente
en [-1, O];
4 ) fl(0) = O, f" < 0;
5) f aplica el segmento [f (l),11 sobre s mismo;
6) la magnitud
195
.-
< <
N+oo
N+oo
&
P(X)
= - ):6(x - x,),
,=1
donde x,, 1 1
Captulo 2
196
f(x)= 1-22
f(z)=l-244
f ( x )= 1 - M / ( l + ( x - a ) * )
a = 0,52
a)
b)
c)
Fig. 7.12. Ejemplos de medidas invariantes de las aplicaciones unidimensionales: a ) La medida tiene singularidades en s = -1 y x = 1; b) La medida es una funcin acotada, continua
y diferente de cero en e l interior de todo el intervalo; c) la medida esta concentrada en dos
islas, lo que corresponde a un ciclo ruidoso
197
dx
X,
198
Captulo 7
En la figura 7.13 b estos puntos se denotan mediante las letras M,, M,, M3
y N,, N,, N,. Todos ellos son puntos
. los punfijos de la aplicacin f 3 ( ~ )En
tos M,, M2,M3 la derivada de f 3 ( ~es
)
la misma y no supera en mdulo a la
unidad (vase la condicin de estabilidad de un ciclo). Precisamente, estos
puntos determinan el ciclo estable S3.
La pendiente del grfico en los puntos
Y
N,, N,, N, tambin es la misma, pero
0,8
b)
all se cumple la desigualdad inversa.
En la aplicacin f a ellos les corres0,6
ponde el ciclo inestable S3que apareci
O,4
= 3,83
simultneamente con el estable.
El surgimiento simultneo de pun02
tos singulares estables e inestables recio
0,2 0,4 0,6 0,8 1 x
bi el nombre de bifurcacin tangencial.
Este nombre se debe a que en el punto
Fig. 7.13. Forma de La funcin f3 ande bifurcacin la curva f (x) es tangente
tes y despus del surgimiento del ciclo
a la diagonal y = x. Correspondienteestable S3: a) X < Ac; 6) A > A,
mente, la derivada 0f /0x en este punto
es igual a la unidad. Subrayemos la analoga de esta bifurcacin con el
punto regular extremal, analizado en el captulo 5, donde los estados de
equilibrio estable e inestable aparecieron simultneamente de la nada.
Tomemos un valor inicial cualquiera y veamos cmo acta la aplicacin f 3(x) para X > X y para X < X. En otras palabras, vamos a hacer
un seguimiento de cada tercer elemento de la sucesin {x,). Los clculos
muestran que, en el primer caso, despus de largo un proceso de transicin
los puntos se acercan al ciclo S3 (fig. 7.14 a). En el segundo caso tiene
lugar un movimiento lento hacia el punto M (fig. 7.14 b), pero despus
los elementos de la sucesin se alejan de este punto. Posteriormente, ellos
comienzan nuevamente a aproximarse a l. En la figura 7.15a se ilustra tal
comportamiento para una de las aplicaciones con vrtice agudo. Como se
observa, los intervalos de movimiento hacia el punto M, cuando la solucin
se parece a una solucin regular, se alternan en la sucesin con saltos caticos rpidos. En la dependencia de z, respecto a n los destellos turbulentos
cortos se alternan con intervalos largos de fase laminar (fig. 7.15 b). En este
modelo elemental hay intermitencia.
El anlisis de los sistemas reales con intermitencia cerca de un punto
de bifurcacin tangencial, encierra dificultades considerables, puesto que
durante un largo tiempo se puede observar solamente la fase laminar.
199
.-
s3.
Fig. 7.15. Cuadro tpico de intermitencia: o) algunas de Las primeras iteraciones de la aplica.
cion xn+, = 1 - Ixn - ~ 1 " ~ [ 1(+x n- X ) 2 ] - 1 f X = 0,769; b) dependencia de La sucesin respecto al nmero de La iteracin. Para una mejor ilustracin, Los puntos x , y xn+, estn unidos
200
Captulo 7
Preguntas y ejercicios
1. Supongamos que se conocen la aplicacin unidimensional xl+, = f (x,) y el
k-simo trmino de la sucesin { x , ) . Es posible hallar a partir de estos datos
los trminos 2,-,, zk-,, etctera? Esto siempre es posible? Se requiere para
esto alguna informacin adicional?
2. Qu modelos matemticos de tiempo discreto pueden ser propuestos para los
sistemas tipo depredador-presa, parsito-hospedadom, para la descripcin
de la dinmica de dos especies que compiten por recursos comunes?
+ +
COS x1
= x2,
COS x2
= x3/
201
.7. Un plano inclinado infinito forma un ngulo a con el horizonte. Una bola
absolutamente elstica cae sobre el plano con un ngulo de cada 7r/2 - S y
una velocidad de cada v . Rebotando en el plano y cayendo nuevamente, la
por el plano. Hallar la coordenada del n-simo rebote
bola comienza a <<subir>>
y estimar la distancia mxima que avanza la bola a lo largo del plano desde el
punto donde cay por primera vez.
8. Consideremos un tringulo ABC. Tracemos sus alturas. Denotemos sus bases
(es decir, los puntos donde ellas cortan los lados correspondientes) mediante
A,, B,, C,.Ahora tomemos el tringulo AIB,C, y denotemos las bases de
las alturas mediante A,, B,, C,,etctera. Denotemos los ngulos del tringulo
A,, B,, Cl mediante a,, B, ,y, ; los ngulos del tringulo A, BnC, mediante
u,, B,, y,. Construir la aplicacin
< <
202
13.
14.
15.
16.
17.
Captulo 7
= f (S,, Y,),
= 9(xn,Y,)?
18. Supongamos que x, = 1, x, = a > O, y
Yn+l
- 22.;
xn+l -
cJ2;
Bibliografa
Una visin general de la teora de las aplicaciones unidimensionales y un anlisis
de sus resultados fundamentales se exponen en los libros:
Collet P., Eckmann J. P. Iterated Maps on the Interval as Dynamical Systems.
Birkhauser, 1980;
Guckenlzeimer J., H~lmesP. Nonlinear Oscillations, Dynamical Systems and
Bifurcations of Vector Fieldes. Springer-Verlag, 1990;
Schuster H. G. Deterministic Chaos. An Introduction. Physik Verlag, 1984.
La teora de la universalidad se discute tanto en estos libros como en los artculos:
Feigenbaum M. J. Universal Behaviour in Nonlinear.Systems. Los Alamos Science 1, 4-27, 1980;
Kuznetsov A. P., Kuznetsov S. P. Dinmica crtica de las aplicaciones unidimensionales. Parte 1. Escenarios de Feigenbaum (en ruso).
203
La preparacin matemtica debe ser 'ms seria para leer los artculos originales
* Feigenbaum M. J. Quantitative Universality for a Class of Nonlinear Transformations // J. Statist. Phys. 1978. Vol. 19, No 1. P. 25-52;
Feigenbaum M. J. The Universal Metric Properties of Nolinear Transformations //
Ibid. 1979. Vol. 21, No 6. P. 669-706;.
Vul Ye. B., Sinay Y. G., Janin K.M. Universalidad de Feigenbaum y formalismo
termodinmico (en ruso).
Un resumen de una serie de resultados claves en este campo, orientado a los
matemticos, se presenta en el libro
Sharkovski A. N.,Maistrenko Yu. L., Romanenko E. Yu. Difference Equations and
their Applications. Kluwer Academic Publishers, 1993.
CAPITULO 8
-. ---
--
--
205
206
Captulo 8
.-
+ +
+ +
+ +
c (x,+ a,, y, + b,) -+ (xl + a, + At f (x, + a,, y, + b,),
y, + b, + Atg(x, + a,, y, + 4)) c'.
.--
207
+ a,, b1 +y,)
=f
+ fxal + jybl +
+ At(a,f, + b,fy),
b; = b1 + At(alg, + blgo),
a; 2 a, + At(a,f, + 4f,X
a;
a,
4 4 + At(a2gx+ 49,)*
La ltimas relaciones permiten hallar el rea del tringulo A'B'C' (se pueden
despreciar las alteraciones pequeas de las aristas):
SABc = V(tl),
SAlBjc = V(tl
+ At)-
donde
a = f,+g,.
(8-4)
a f i ( x , , ,xN)
.
axi
i=l
La frmula (8.5) permite precisar el concepto de sistema disipativo, Denominaremos sistemas disipGivos a los sistemas dinmicos en los que ' 0 < O
v = OV,
donde S2 =
208
Captulo 8
4
209
.-
de Stokes
Supongamos que A = (-9, f ); dF = (dx, dy) = (f dt, g dt). Precisamente en la ltima igualdad se considera el hecho de que estamos tratando con
el campo vectorial correspondiente al sistema dinmico (8.1). Est claro que
la integral del vector A por el ciclo lmite es igual a cero. A partir de la
frmula del rotacional de un campo vectorial,
210
Captulo 8
?
. Y 3
Esto contradice las condiciones del teorema, puesto que la integral por
la regin G debe ser o bien positiva, o bien negativa. Esto demuestra la
afirmacin. Veamos en calidad de ejemplo uno de los modelos tipo (8.1)
ms conocidos.
Ejemplo. Modelo de Lotka-Volterra.
211
Por tanto,
r2Nl -
log N,
er(2N~ ~ - ' 2
= ce-71N2 ~
2 '
En el plano de fase (N,, N,) a cada valor de C le corresponde una curva cerrada
(fig. 8.3).En el interior de estas curvas est
\
Si se considera que en el instante inicial el
O ) prximo
)
al punto
punto ( N , ( O ) , N ~ ~est
singular, entonces se puede linealizar el
sistema dinmico (8.8) en un entorno de
este estado de equilibrio:
N2 = - + A N 2 .
71
Despus de linealizar el sistema (8.8) obtenemos
Nl = 2 + A N l ;
'72
Ni
Fig. 8.3. Trayectorias de fase tpicas en
e l modelo de Lotka-Volterra
-- --
--
--
'72W
AN, = -AN2(0)
sen (wt).
&,'y1
212
Captulo 8
a las armnicas. Lejos de este punto las oscilaciones tienen una forma ms compleja.
Estas oscilaciones se pueden interpretar del modo siguiente. Mientras el tamao
poblacional de las presas sea suficientemente grande, el tamaiio poblacional de los
depredadores crecer (sector AB en la figura 8.4). Sin embargo, tiene lugar una
retroalimentacin negativa: el aumento del nmero de depredadores conduce a la
disminucin de la poblacin de presas. Esto hace que comience a disminuir el
nmero de depredadores (sector BC en la figura 8.4). Pero cuando N,(t) se hace
suficientemente pequeo, la poblacin de presas comienza a crecer nuevamente,
Esto provoca nuevamente el crecimiento de la poblacin de depredadores, etctera.
Come~itandolas consecuencias
biolgicas del modelo matemtico
construido, Volterra seala: A aquellos que se ocupen de la comprobacin experimental de los resultados
obtenidos, es decir, de la aplicacin de
la teora, les corresponder realizar el
anlisis detallado de las suposiciones
1
iniciales y de la legitimidad biolgica
A
B
C
de las conclusiones basadas en los
Fig. 8.4. Forma caracterstica de La curva
experimentos, las observaciones y su
integral en e l modelo (8.7)
anlisis estadstico.Se debe subrayar
que en muchos problemas relacionados con la simulacin matemtica, esta parte del trabajo es la ms complicada y
laboriosa. Precisamente por esto es que en el estadio de la investigacin matemtica
del modelo se debe tratar de establecer sus posibles defectos.
En el caso del sistema dinmico (8.7),como uno de tales defectos, por lo visto,
se debe considerar la existencia de la integral (8.8). En efecto, es natural suponer
que si en cierta rea del hbitat de los animales se observan oscilaciones del tamao
poblacional, entonces su amplitud y perodo deben depender de las propiedades del
rea y de las especies que en ella habitan. Al mismo tiempo, se puede suponer que la
amplitud y el perodo no dependen de cambios pequeos de las condiciones iniciales
N,(O) y- N2(0).En otras palabras, en el sistema deben ocurrir autooscilaciones.
Sin embargo, el modelo (8.7) no posee esta propiedad. En efecto, las magnitudes
N,(O) y N2(0) determinan unvocamente la constante C de acuerdo con la frmula
(8.8). A su vez, a diferentes valores de esta constante les corresponden diferentes
curvas cerradas en el plano de fase y, por consiguiente, oscilacion~scon amplitudes
diferentes, es decir, en tal sistema las condiciones iniciales nunca son olvidadas.
Ntese que esto sucede a pesar de que en este sistema i2 # 0.
Regresemos a las preguntas tradicionales que se intentan aclarar al investigar
los sistemas de dos ecuaciones diferenciales ordinarias.
213
.-
Ay = g(x* Ax, y*
es decir,
P* - (a1, + ~
_-_
__
_-
__
- a,,a,,) = 0,
p2 a l p a. = o,
2 2 ) (alla,
~
-_
-__
donde
al = - tr A,
a0 = det A,
214
Captulo 8
215
.fase del sistema es una elipse. Para una desviacin pequea respecto al
estado de equilibrio, el sistema comienza a describir una elipse alrededor del punto (x*,y*). La posicin de equilibrio es estable, pero no
asintticamente estable.
Como se deduce de la ecuacin (8.12), el'tipo
del punto singular est determinado por la traza y
el determinante de la matriz del sistema linealizado
(por los nmeros a, y a,). Si se tiene en cuenta que
el sistema (8.1) depende del parrnetro A, y consiguientemente, a, = cu,(A), a, = a2(A),entonces
para comprender cmo puede cambiar el tipo del
punto singular al variar el parmetro, resulta muy
til la figura 8.6.
Esta figura ha sido construida partiendo del Fig. 8.6. Divisin del esteorema de Viite y de las propiedades elementales de pacio de parmetros sela ecuacin cuadrtica (8.12). En particular, en esta gn el tipo de los estafigura se observa que para una variacin pequea dos de equilibrio: 1) nodo
un centro se convierte, usualmente, en un foco estable; 2) nodo inestable;
estable o inestable. En esta figura tambin se ve que 3) foco estable; 4) foco inen la lnea a, = 4 1 4 un foco se convierte en nodo. estable; 5) silla; 6) centro
Existen ciclos lmites en e l sistema (8.1)? Cmo son? La existencia de
ciclos lmites es la principal diferencia entre un sistema de dos ecuaciones
diferenciales y una sola ecuacin. Analicemos inicialmente el modelo ms
simple, en el que el nico atractor es un ciclo lmite, y en el que se puede
hallar analticamente todo el conjunto de curvas integrales. Este sistema dinmico se presenta en la cintica qumica, en la hidrodinmica y en la teora de
ondas, durante la descripcin asinttica de un conjunto de modelos tipo (8.1)
cuyos parmetros se encuentran cerca de la lnea al = O, a, > O (fig. 8.6).
Es cmodo escribir este sistema de dos ecuaciones respecto a x(t) e y(t)
como una ecuacin respecto a una funcin compleja W ( t ) :
W(t) = x(t) iy(t),
216
Captulo 8
217
Fig. 8.7. Bifurcacin de nacimiento de un ciclo lmite: o ) dependencia del atractor del
sistema (8.13) respecto al parmetro A; b) nacimiento de un ciclo estable a l aumentar e l
parmetro A; c) ((aniquilacin)) de un foco estable y de un ciclo Lmite inestable
Captulo 8,
218
-
lCI
lC
El sistema dinmico (8.19) se deymina bruselador. Sus puntos singulares son las
soluciones de las ecuaciones algebraicas
Los valores X = O y X = B/Y satisfacen la segunda ecuacin. Sin embargo, X = O
no verifica la primera ecuacin. De la segunda ecuacin se obtiene que Y # O, pues
de otro modo X = O. Por consiguiente, tomando en consideracin que XY = B,
hallamos las coordenadas del punto singular:
X'
eI
'
(B+~YY
4A
de donde a, = A2, al = A2 1 - B .
El parmetro que se puede cambiar en el
sistema es la concentracin B. La figura 8.6
nuestra qu sucede con el punto singular
de este sistema al aumentar el valor del
parmetro B. Por cuanto a, = A2 > 0,
el punto que caracteriza al sistema en
el espacio de parmetros pertenece al semiplano superior. Cuando B es pequeo
( B < A2-2A+l, lnea a, = a:/4),el punto
singular es un nodo. Despus se convierte en un foco estable. Finalmente, para
B = A2 1 el foco estable se vuelve inestable. Aqu ocurre una bifurcacin de Hopf.
El ciclo lmite se convierte en un atractor
del sistema. En la figura 8.8 se muestra un
ejemplo de ciclo de este tipo para valores
suficientemente grandes del parmetro B.
Hasta comienzos de los aos 70 del
siglo pasado, la mayora de los qumicos
219
.-
220
(+,)e
Capitulo 8
.el cual pertenece al interior del ciclo lmite. Introduzcamos una coordenada r a
lo largo de este rayo. Consideremos una trayectoria que parte de un punto A
perteneciente al rayo (fig. 8.9 a). Supongamos que esta trayectoria corta por primera
vez al rayo en el punto B . Introduzcamos la funcin
la cual asigna a cada punto de coordenada rA la coordenada r, del punto B .
Supongamos tambin que el tiempo que tarda una trayectoria en llegar desde el
punto A hasta el punto B es finito tanto para el ciclo como para todas las trayectorias
cercanas a l.
La funcin f es continua y unvoca. Supongamos lo contrario. Si la funcin f
es discontinua en el punto C,
lim f (re
&+O
&-+O
r* = f (r*).
En efecto, la igualdad (8.22) indica que partiendo del punto de coordenada r* en el
eje OA, al cabo de un tiempo finito la trayectoria regresa a este mismo punto.
Evidentemente, si r, + r* para todo r, perteneciente a un &-entornodel punto r*, entonces elgciclo lmite, al cual corresponde el valor r*, ser asintticamente
estable.
La condicin suficiente de estabilidad coincide con la condicin de convergencia
del mtodo de iteracin simple (8.21) de resolucin de la ecuacin algebraica r = f (r):
221
('lq). (T)
si
S,
(F)3
(F)
<
= S, ql = f
S,
S,
(y).
q1 =
e,
< N,
as que qn + r*,
2"
N
N
dh(u)
dU
cosa que el jacobiano calculado a lo largo de la solucin ( t ) . Si escribimos el
sistema (8.25) en forma coordenada
Esta ecuacin es lineal respecto a AG. La derivada
entonces
En calidad de A(0) habra que tomar 63.Sin embargo, por cuanto la ecuacin (8.26)
es lineal, es ms cmodo tomar el vector unitario W en la direccin que nos interesa.
En nuestro caso de dos variables, el \rector W est dirigido a lo largo del rayo OA,
223
. ..
Oscilador de Van der Pol. La ecuacin de Van der Po1 es uno de los modelos
matemticos bsicos. Este sistema dinmico es uno de los sistemas de
ecuaciones diferenciales ms investigados.
Dicho modelo surgi en la electrnica durante la descripcin de
la generacin de oscilaciones en esquemas electrnicos. La presencia de
&"
7
i = a,v,
a, = const
< 0.
-_
-_
--
--
225
Fig. 8.10. Modelo de Van der Pol: a) esquema del oscilador; 6) caracterstica
tensin-corriente del elemento activo
fi
Derivando esta igualdad respecto al tiempo y teniendo presente la dependencia i = i(v), obtenemos la ecuacin de Van der Po1
G/C es
Para E = O la ecuacin de Van der describe oscilaciones cercanas a [as siPo1 se convierte en la ecuacin de las
nusoidales
oscilaciones de un pndulo matemtico
con frecuencia wo = 1.
Con ayuda de computadoras no es difcil verificar que para valores
pequeos de E el sistema posee un nico ciclo lmite y que las oscilaciones
estacionarias son cercanas, por su forma, a las sinusoidales (fig. 8.11). Por
Captulo 8
226
+ u0 = O,
1
2
E : l + u1 = u0 - ~ ~
en ambos miembros.
,
(8.33a)
: o
EO
i i ~ ,
(8.33b)
.,
'
S
ii
uO= A sen (t
+ po),
po = const.
A = const,
+
= A cos(t +
= A cos (t
cos (t
>
227
.-
1 - COS 2a
3
U#)
+ q(t) + r(t),
+ p = 0,
(8.35a)
(8.35b)
(8.35~)
q + q = Ccos(t+po),
F r = D cos (3t 3p0),
donde
+ p,),
p, = const,
Q=
Ct
tos (t + loo),
3
D
r = - - cos (3t + 39,).
8
Es especialmente importante la dependencia q(t), en la cual figura
explcitamente la coordenada temporal. (La ecuacin (8.35b) se puede
resolver utilizando el mtodo de variacin de la constante.)
De este modo, el desarrollo directo proporciona la solucin asinttica
+
A = A, +&Al+ E ~ A ~...,
$0,
+ + e2cp2+
= qo +cpl
228
Captulo 8
tos (37
~ E W sen
~ T (37
+ 3p0) + ... .
u = A sen (T
+ po)+ E
229
Por cuanto no nos interesa la solucin trivial, sino el ciclo lmite, tomamos
A = 2.
Consiguientemente, el ciclo lmite en la ecuacin de Van der Po1 est
determinado por la igualdad r = t
u = 2 sen (t
(8.42)
En lugar de una funcin incgnita u(t), la cual satisface la ecuacin (8.31), hemos obtenido tres funciones u(t), ~ ( t )P(t),
,
relacionadas por
la igualdad (8.42). Esta igualdad proporciona cierta libertad, que podemos
utilizar para hallar la solucin en una forma ms simple. Esto se puede hacer
por analoga con el caso lineal. De (8.41) se deduce que u = -U sen (t P).
Asumamos que esto es vlido tambin para la ecuacin no-lineal:
+ P(t)).
(8.43)
230
Captulo 8
Y
(8.45)
El sistema de dos ecuaciones diferenciales ordinarias de primer orden (8.44), (8.45) para las funciones a y p, se puede resolver respecto a
a y B. Multiplic.ando la igualdad (8.44) por cos (t P), la igualdad - (8.45)
por sen (t ,O) y sumando los resultados, obtenemos
+ p)).
cos (t
+ /3)
y de
(8.46b)
A ,
Por cuanto el seno y el coseno de cualquier argumento real pertenecen al segmento [-1,1], para toda funcin a acotada obtenemos = O(&)
y 1= O(&).
Por consiguiente, para un valor suficientemente pequeo del parmetro E
las funciones a(t) y P(t) son funciones de variacin lenta dependientes del
tiempo. Por eso, en el segmento [O, 2ir] ellas cambian poco, y en una primera
aproximacin se las puede considerar constantes en este intervalo.
Esto da la posibilidad de promediar en el segmento [O, 2?r] los miembros
primeros y segundos de las ecuaciones (8.47):
231
Aqu hemos utilizado el hecho de que (o x t (de lo contrario habramos tenido que elegir de otro modo los lmites de integracin y promediar respecto al
perodo) y de que los valores medios de las funciones sen (np) y cos (mp)
(n, m = 1,2,3, ...) respecto al perodo, son iguales a cero. Promediando de
la misma manera la ecuacin (8.47), para la funcin P(t) obtenemos
As pues, las ecuaciones (8.48) proporcionan la ley del variacin de las variables a(t) y P(t). Junto con la igualdad (8.43) ellas determinan la solucin
aproximada de la ecuacin de Van der Pol. En particular, en el marco del
mtodo de promedio se logra describir no slo los parmetros del ciclo
lmite, sino tambin el rgimen de transicin.
232
Captula~
.-
5%
Para
Ntese que los espacios de fase de los sistemas (8.49) y (8.50) son diferentes.
En efecto, sea z un vector de m componentes y sea y iin vector de n
componentes. Para determinar unvocamente el estado del sistema (8.49)'
es necesario definir m n nmeros. La trayectoria est determinada por2
m n condiciones iniciales. En el caso del sistema (8.50) el espacio de fase
es n-dimensional y slo se necesitan n condiciones iniciales.
Respecto a ( 8 5 0 ) el sistema (8.49) se denomina frecuentemente sistema
perturbado, y respecto al sistema (8.49) el sistema (8.50) se denomina
sistema degenerado. Cuando en presencia de una perturbacin el tipo de
un sistema cambia respecto al sistema no-perturbado de tal manera que'
para determinar las soluciones del sistema perturbado se requieren ms
condiciones complementarias que para el no-perturbado, la perturbacin se
denomina singular.
Se puede mostrar adnde pueden conducir las perturbaciones singulares en el ejemplo del modelo lineal elemental
Comparemos las frmulas (8.51) y (8.52). Ntese que si t > & / a ,entonces el
primer trmino en la frmula (8.52) se vuelve pequeo y x(t) -+ -b/a cuando
t + m, es decir, la solucin del sistema perturbado tiende a la solucin
del sistema no-perturbado. Sin embargo, en la regin O < t < / a estas
soluciones se diferencian sustancialmente, sin importar cun pequeo sea
el parmetro E. Este fenmeno, el cual consiste en que en presencia de
una perturbacin singular la solucin de la ecuacin perturbada (8.49) se
233
234
Captulo 8'
.El ejemplo lineal considerado muestra que no siempre se puede pasar
del sistema (8.49) al sistema (8.50). En efecto, si a < O en la frmula (8.52),
entonces no existe ningn valor de t para el cual la funcin z(t) tiende
a la solucin del sistema degenerado (8.51). Un resultado clsico que da
las condiciones suficientes para que la solucin del sistema degenerado
est cerca de la solucin del sistema inicial cuando E -+ O, es el teorema de
A. N. Tjonov. Analicemos este teorema y su demostracin.
Regresemos al sistema degenerado (8.50). Supongamos que se logr
hallar la solucin 2 = p(y, t) de la ecuacin F(z, y, t) = 0, que Z es una
funcin continua de los argumentos y y t, y que k es una raz aislada de la
ecuacin F ( z , y, t) = O (es decir, incrementando z en un vector pequeo u,
obtenemos F(z u, y, t) # O. Despus de sustituir la solucin hallada en la
ecuacin para y, resulta
--
para
para
&+O
&+O
< t < T,
O < t < T,
O
235
..
inicial. Sin embargo, en una serie de casos es importante tener una idea
de cmo cambian las variables en la capa lmite. En este caso tambin se
utilizan mtodos asintticos. Las series asintticas correspondientes se diferencian sustancialmente de las que aparecen en los problemas regularmente
perturbados. En los ltimos la solucin aproximada x,(t, E ) a la solucin
exacta ~ ( tE ,) , se busca frecuentemente en forma de una serie de potencias:
Generalmente esta serie tiene una precisin Ilx(t, E ) - z,(t, E)II de oruniformemente en el segmento O t T . En los problemas
den o(~"+'),
singularmente perturbados que hemos analizado, la serie contiene'tambin
trminos denominadosfunciones de capa lmite nk(7):
< <
236
Captulo pi
I
(x2- l ) x + x = O.
(8.56)
+ x = O se puede integrar:
r,
<
Fig. 8.13. Ecuacin de Van der Po1 con e l parmetro pequefio en la derivada superior:
a) forma caracterstica del ciclo lmite en e l plano de Linard; 6 ) curvas integrales de La
ecuacin para E = O; C) forma tpica de las oscilaciones de relajacin
237
.-
,
n
m3
238
Captulo 4
239
/-j
-****e.
i
a =0399*7404512
Fig. 8.14. Ecuacin de Van der Po1 con trmino independiente singularmente perturbada. Crecimiento
superrpido de La amplitud del ciclo lmite despus
de La bifurcacin de Hopf. Aqu e = 0,01
Cuello
Fig. 8.15. Forma caracterstica de las resonancias denominadas patos. Las trayectorias
patos siguen durante algn tiempo un tramo repelente de La curva Lenta
En el anlisis-d e las-ecuaciones diferenciales no-lineales singularmente perturbadas se utiliza ampliamente la simulacin computacional. Sin
embargo, aqu es necesario prestar una especial atencin a los algoritmos computacionales que se aplican, y utilizar frecuentemente mtodos
numricos especiales.
La esencia de las dificultades surgidas se puede ilustrar mediante
el ejemplo de un problema lineal de capa lmite. Consideremos un peso
pequeo de masa E unido a un resorte dentro de un medio con rozamiento
viscoso. Supongamos que su dinmica es descrita por la ecuacin diferencial
x ( t ) = C, exlt + c2eA2',
donde A,, X2 son en este caso los valores propios de la matriz de Jacobi A,
y C,,C2 son constantes determinadas por las condiciones iniciales. Aqu
(o
241
< <
Zn+1
242
Captulo 8
.-
Preguntas y ejercicios
<
x = Xx , y = by?
4. Para la descripcin de la reaccin qumica oscilante de Belosov-Zhabotinski,
R. Field y R. Noyes propusieron el esquema de reaccin siguiente:
lo-',
lo2; w es el
243
S + g(t)x
+ w2(t)x= o,
donde w(t) es una funcin montona positiva. Obtener las condiciones suficientes de acotacin de las oscilaciones en este sistema. (Indicacin. Realizar
t
el cambio de variables
u($)
ds y analizar la variacin de energa del
i
sistema.)
10. Investigue el comportamiento cualitativo de los sistemas dinmicos
Bibliografa
Las nociones bsicas de la teora clsica de oscilaciones, en la cual los ciclos lmites
desempean un papel clave, se presenta en los libros:
Andronov A. A., Vitt A. A,, Khaikin S. E. Theory of Oscillators. Dover, 1989;
Rabinvich M. 1,Trubetskov D. 1. Introduccin a la teora de oscilaciones y ondas
(en ruso);
244
Captulo 8
.-
Mtodos topoLgicos
en la investigacin
de Los sistemas no-lineales
En cierto sentido la topologa es la ciencia que
estudia la continuidad: partiendo de la continuidad
del espacio o de las formas, ella pasa a las generalizaciones, las cuales despus conducen a un nuevo
concepto de continuidad, y el espacio habitual,, tal
como nos lo imaginamos, se queda muy atrs, lejos.
S. Barr
246
Captulo 9
x' = ax + by + p ,
y' = cx
+ dy + q,
ad - bc # O.
247
.-
248
Captulo 9
< <
1
'
249
250
Captulo 9
Fig. 9.5
) - hl (1 (a +90') ) = -g (a + 90),
251
conjunto de todos los puntos de X cuyas distancias hasta x son menores que r, se .
denomina entorno de radio 7 del punto x en el conjunto X y se denota mediante
el smbolo N ( r , x, X).
Definicin. Un subconjunto U de un conjunto X (U C X) se denomina
conjunto abierto en X si para todo punto x E U existe un nmero r > O
tal que N($,
T , U)C U.
El lector puede verificar que todos los entornos son conjuntos abiertos,
252
Captulo 9
253
----
i,j
254
Captulo 9
-------
255
.-
y {DiK; DiK; DiK;), y que los vrtices DIKl D2K2D3K3son los mismos para ambas configuraciones, pero las quebradas pueden ser diferentes (fig.9.8).
Supongamos que en el primer caso D I K , es la quebrada a, y en el segundo,
Di K; es la quebrada
(fig.
i 9.8). Las quebradas ri forman un ciclo. Denotmoslo
mediante a&.Pero D2K,D3 K3D2 tambin
forma un ciclo. Sin embargo, el grado de
interseccin de estos ciclos es igual a ce- a l a
ro, es decir, sustituyendo el sendero DIKl
por la quebrada DiK;, no cambiamos la
paridad de la interseccin. Razonando del
Fig. 9.8
mismo modo para los dems senderos
concluimos que los senderos restantes se
pueden deformar, convirtiendo {DIKl D2K2D3K 3 ) en {DiK; D;K;Di K j } .
Supongamos que existe una configuracin de los senderos en la cual no
hay intersecciones y, por tanto, el grado de interseccin In de los senderos
no-adyacentes es igual a cero.
Sea esta configuracin {DiKiDiKiDiKi). Pero el grado de interseccin In de esta configuracin es exactamente el mismo que el de la
configuracin mostrada en la figura 9.6, donde In = 1. Esto significa que
es imposible trazar senderos desde las casas hasta los pozos sin que dichos
senderos se corten.
Problemas similares, aunque ms complejos, relacionados con la bsqueda de invariantes topolgicos, surgen en la teora de nudos. En esta
teora se determina si es o no posible transformar un nudo de una cuerda
cerrada en otro nudo sin cortar la cuerda. En los ltimos tiempos la teora
'de nudos ha comenzado a aplicarse activamente en la fsica estadstica.
Probablemente, aparecern otros campos de aplicacin, sobre todo en la
qumica, donde se trata no de deshacer nudos gordianos, sino de analizar
sus configuraciones.
(9.1)
Captulo 9
256
.-
Para todo poliedro cuya superficie es horneomorfa a una esfera y cada cara es
horneomorfa a u n crculo, se verijica la igualdad (9.1).
Fig. 9.9
257
Fig. 9.10
V-A+C
.-.
~PP
>O
(p - 2)(q - 2) < 4.
El producto de dos nmeros positivos enteros puede ser menor que cuatro
si los factores son iguales, respectivamente, a
Esto conduce a los siguientes pares {p, q): {3,3) - tetraedro, {4,3) cubo, {3,4) - octaedro, {5,3) - dodecaedro, {3,5) - icosaedro.
258
-
Captulo 9
.-
A D
D'A'
Fig. 9.J 1
259
.-
260
Captulo 9
.-
261
Ca ptu16 !#
262
263
.-
Tomemos una arista del polgono, por ejemplo, r,, la cual une los
vrtices A y B. Supongamos que al movernos de A hacia B seguimos cierto
vector W que coincide con El(x, y), y que al movernos de B hacia A, este
vector coincide con v2(x,y).
Al moverse desde A hasta B y, despus, desde B hasta A, el vector W
realiza un nmero entero de vueltas. Esto se deduce de la construccin
del campo vectorial Zl(x, y): 5,(A) = v2(A), %(B) = v2(B). Denotemos
con d(rl) el nmero de vueltas que da el vector G al moverse a lo largo de
la arista rl. En la figura 9.15, por ejemplo, d(r,) = 1, d(r2)= O, d(r,) = -1.
Sealemos que si recorremos la arista r1 en sentido contrario, desde B
hasta A, y regresamos a B, entonces la magnitud d(rl) cambia de signo.
Recorramos uno de estos polgonos M en sentido positivo. Supongamos que, en este caso, el vector G1(x, y) realiza z,(M) vueltas. Recorramos
ahora el contorno en sentido contrario y supongamos que el vector v2(x,y)
realiza z2(M)vueltas. Naturalmente, z, y z2 son nmeros enteros, pues los
vectores considerados, despus de recorrer el contorno cerrado, regresan a
la posicin inicial.
Como resultado, el vector w' realiza zl(M) - z2(M) vueltas. Al mismo
tiempo, observando el movimiento a lo largo de cada una de las aristas, se
puede concluir que
Ezi(Mk) E
z2(Mk)= O*
(904)
Sin embargo,
k
264
Captulo 9:
.-
Supongamos que en el interior del poligono M hay un punto singular xo (fig. 9.17),
Construyamos un sistema de lneas cerradascercanas la una de la otra, como se muestra
n o suna de ellas
en la figura. ~ u ~ o n ~ a r que
coincide con una circunferencia pequea coit
centro en el punto x0, y que otra coincide
con el contorno del polgono M. Al pasar
de una lnea a otra cercana a ella, como se
deduce de las consideraciones de continufdad, el nmero de vueltas de! vector G,(X,9)
Fig. 9.17; Construyendo un siste- debe cambiar POCO. Sin embargo, el nmero
por tanto, no puede
ma de lneas cerradas se puede de vueltas es entero
comprobar que el ngulo de giro cambiar en una magnitud pequea, es dedel vector 3 al recorrer e l contar- cir, se mantiene constante. Pero al recorrer la
est determinado por jndice circunferencia pequea con centro en el pundel punto singular en el interior
to x, obtendremos el ndice de este punto,
del polgono
y al recorrer el contorno del polgono obtenemos la magnitud zl(M).Consiguientemente,
estas dos magnitudes son iguales. Entonces, de la frmula (9.4) se deduce
que la suma de los ndices de los campos vectoriales vl(x, y) y v,(x, y)
(y de otros campos en la superficie G) es la misma.
Fig. 9.18. puntos singulares de un campo
vectorial cuyo ndice no es difcil calcular:
a ) nodo situado en un vrtice; b) punto de
silla situado en el punto medio de una arista;
c) nodo situado en el centro de una cara
b)
a)
c)
265
De este modo, existe una relacin notable entre las propiedades topolgicas del espacio de fase G y los campos vectoriales que se pueden
construir en l.
266
Capitulo Q
a
.-
267
268
Captulo 9
.relacionadas con las diferentes tecnologas que utilizan las ciclodextrinas. Las
investigaciones de Usnov mostraron que la participacin de la ciclodextr~.
en las reacciones habituales del libro de texto escolar puede ejercer una:
influencia cardinal en la marcha de dichas reacciones. En los ltimos ao$
han aparecido razones de peso para considerar que la ciclodextrina y sus
anlogos no son tan exticos como parece, y que la naturaleza los utiliza en
muchas de sus estructuras.
La sntesis orgnica directa de las ciclodextrinas es compleja y costosa:
Por eso, actualmente se trabaja en la bsqueda de microorganismos que puedan convertir el almidn habitual en ciclodextrina. Para qu necesitan las
bacterias cambiar la topologa del almidn produciendo ciclodextrina? Existen dos hiptesis. Conforme a .la primera, ciertas cepas de bacterias realizan
este proceso para almacenar alimento, creando una especie de conservas
que no representan ningn inters para los competidores. De acuerdo con la
segunda, las bacterias utilizan las ciclodextrinas para cambiar la solubilidad
de una serie de sustancias que participan en su ciclo vital.
Sin embargo, la bsqueda de microorganismos productores de anlogos
de las ciclodextrinas y de objetos afines, exige grandes esfuerzos. Por eso
se tiene que acudir nuevamente a la simulacin matemtica. Estimar con
ayuda de computadoras la estabilidad de las molculas cuyos productores
se busca, y modelar las reacciones qumicas en las cuales los compuestos
creados deben participar en el marco de una u otra tecnologa.
Hay otro campo de las investigaciones aplicadas donde los mtodos
topolgicos conducen a resultados importantes. Hasta el da de hoy se
han descubierto y descrito ms de siete millones de compuestos qumicos;
Cmo predecir a partir de los compuestos conocidos las propiedades de
las sustancias an no creadas? Cmo pronosticar las caractersticas de los
objetos ya conocidos sin tener que recurrir a experimentos costosos?
Se puede proponer el siguiente enfoque. Suponer que la molcula es
cierto grafo, ignorar su forma tridimensional, la magnitud de los ngulos
entre los enlaces qumicos y otras caractersticas. Introducir ciertos ndices
topolgicos que reflejen, por ejemplo, con cuntos tomos est enlazado
cada tomo de una molcula dada, y si stos forman cadenas lineales o
estructuras ramificadas.
Hace ya bastante tiempo se conoce un ndice elemental, relacionado
con un invariante topolgico de un <<grafoqumico. Si se analizan los
hidrocarburos, es decir, los compuestos de carbono e hidrgeno, entonces
se puede considerar que el grafo est compuesto slo por los tomos de
carbono (vrtices) y los enlaces entre ellos (aristas). El invariante elemental
es el nmero de vrtices. Este invariante, conocido en la qumica como
nmero de carbono, est relacionado con toda una serie de caractersticas
fsicas de los hidrocarburos.
269
Sin embargo, el anlisis de las molculas ramificadas requiere parmetros topolgicos ms complejos. Uno de ellos es el ndice de Wiener. A veces
en la teora de grafos se utiliza el concepto de distancia topolgica entre dos
vrtices, que es el nmero de aristas que intervienen en el camino ms corto
entre estos dos vrtices. Para el caso de una molcula, el ndice de Wiener es
igual a la suma de las distancias topolgicas entre todos los pares de tomos
de su esqueleto de carbono. Resulta que conociendo el ndice de Wiener se
puede determinar con mucha precisin la temperatura de ebullicin de una
amplia clase de compuestos.
Existe toda una serie de ndices topolgicos que permiten pronosticar tales caractersticas como la temperatura de fusin, el octanaje y el rendimiento
de holln en la combustin incompleta de los hidrocarburos. Se ha logrado
relacionar algunos ndices con la manifestacin de actividad fisiolgica en las
molculas: con su aroma, su toxicidad, su capacidad de detener el crecimiento de una serie de bacterias, y su carcinogenicidad. La topologa molecular
ha resultado ser un camino efectivo de simplificacin del anlisis de los
objetos estudiados, lo cual es muy importante en la simulacin matemtica.
Preguntas y ejercicios
1. Demostrar que toda aplicacin continua f de una circunferencia en una recta,
transforma cierto par de puntos diametralmente opuestos en un punto.
2. Demostrar que la ecuacin algebraica no-lineal de grado impar
donde
f
tiene al menos una raz real.
3. Se pide construir un sistema dinmico con las propiedades siguientes:
\
3 = f ($1 y),
Y = 9(x,y),
f (x, Y) = f (3 a, y), g(x, y) = 9 ( +~ a, y),
f(x,Y) = f(S, y b), g(x, y) = g(x, y + b).
.Es posible conectar diez ciudades entre s mediante caminos que no se corten,
y de modo que de cada ciudad salgan cinco caminos que lleguen a otras cinco,
ciudades?
Demostrar que en todo poliedro convexo existen dos caras con el micm0
nmero de aristas.
Demostrar que no existe un poliedro, en cada vrtice y en cada cara del cu4
inciden no menos de cuatro aristas.
En un plato yacen dos hojuelas de forma irregular. Es posible cortarlas con
solo golpe de un cuchillo en dos partes equivalentes (de reas iguales)?
Se puede cortar en cuatro partes iguales una hojuela en un plano mediante
dos rectas perpendiculares?
Uno de los primeros resultados de la geometra proyectiva, obtenido en el
siglo XVII, 'es la afirmacin siguiente. Supongamos que dos tringulos ABC
y A'B'C estn dispuestos en un plano de tal manera que las rectas AA', BB'
y CC' se cortan en un punto. Denotemos mediante P, Q y R los puntos de
interseccin de las rectas AB y A'B', AC y A'C', BC y B'C', respectivamente;
Demostrar que los puntos P, Q y R pertenecen a una misma recta.
Cada uno de los puntos de una circunferencia es negro o blanco. Demostrarque existe un tringulo issceles inscrito en esta circunferencia cuyos vrtices
son de un solo color.
12. Tres torres de un castillo estn ubicadas en los vrtices de un tringulo ABC.
La cuarta torre est en un punto O del interior del tringulo. Las murallas del
castillo se extienden a lo largo de los segmentos AB, BC, AC, AO, BO, CO.
Es posible visitar todo el castillo, es decir, recorrer todas las murallas y visitar
todas las torres, sin pasar dos veces por una misma muralla?
13. Se puede cubrir todo el plano con un nmero finito de interiores de parbolas?
14. En un plano se han trazado n lneas rectas. Demostrar que las regiones, en las
que estas rectas dividen el plano, se pueden pintar con dos colores, de modo
que cualesquiera dos regiones vecinas (que tienen en comn un segmento de
recta) poseen diferente color.
-- .Imtryamos
uEZSnadeMo5iusy despus cortmosla a lo largo de la lnea
central. Es no-orientable la superficie obtenida?
16. Demostrar que toda aplicacin continua de un crculo en sf mismo, tiene al
menos un punto fijo. (Indicacin. Considerar la analoga de este problema con
el problema del peinadode una esfera o un toro.)
17. Dado un tringulo acutngulo ABC, hallar en su interior un punto O tal
que la suma de las longitudes de los segmentos OA, OB, OC sea mnima.
(Indicacin. Prstese atencin a los problemas similares que surgen en relacin
con las pelculas de jabbn, las superficies mnimas, etctera.)
18. Los poblados A, B, C,D se encuentran en los vrtices de un cuadrado. Cul
es el sistema de caminos de longitud mnima que une los poblados A, B, C, D
de tal manera que desde cada uno se puede llegar a cada uno de los restantes?
Se puede generalizar este resultado a un sistema de n poblados A,, . ,A,?
Qu algoritmos de clculo propondra usted en este caso? (Indicacin. Utilice
el resultado del problema anterior.)
--
- -
..
.-
271
Bibliografa
Las nociones fundamentales de la topologa a nivel elemental se pueden ver en los
libros:
272
Captulo g
Redes neuronales
Hoy no podemos decir con certeza si la ecuacin de Schrodinger contiene ranas, compositores
o alguna moral, o si ella no puede contener nada
semejante. Nosotros no podemos decir si se necesita
o no algo superior a la ecuacin, alguna especie
de dioses. Por eso cada cual puede tener su propia
opinin al respecto.
R. Feynman
Neurociencia
Las computadoras han ayudado mucho a la dinmica no-lineal. Ellas han
permitido investigar modelos matemticos que surgen en diferentes reas
y descubrir todo un conjunto de efectos no-lineales interesantes. Hoy la
dinmica no-lineal est pagando las deudas a las ciencias y tcnicas
computacionales.
La simulacin computacional desempea un papel importante en
el enfoque interdisciplinario denominado neurociencia. Este enfoque est
llamado a responder cmo trabaja el cerebro, cules sucesos psquicos
estn determinados por patrones de impulsos nerviosos en el cerebro.
La neurociencia surgi a caballo entre varias disciplinas cientficas, entre las
que se cuentan la neurofisiologa, la bioqumica, la psicologa cognitiva, la
matemtica y la tcnica computacional.
El desarrollo de la neurociencia mostr que la construccin e investigacin de modelos matemticos son el mtodo principal de comprensin de
los procesos que transcurren en el cerebro, de interpretacin de los datos
. experimentales disponibles y de planteamiento de problemas nuevos. Naturalmente, hay que reconocer que el camino desde el estado actual de este
campo hasta la comprensin profunda de los principios de funcionamiento
del cerebro, por lo visto, es muy largo. Los modelos ms bien responden al
cmo podran funcionar unos u otros sistemas que concuerdan en algunos
aspectos con los datos sobre la arquitectura, las funciones y las particularidades del cerebro. No obstante, las investigaciones en la neurociencia ya han
abierto el camino para la creacin de nuevas arquitecturas de computadoras
274
Captulo 10
Redes neuronales
275
.-
Si <<olvidamos
los valores de x en algunos puntos, perdemos suficientemente rpido
la idea de toda la figura (fig. 10.1 b). Si olvidamos los coeficientes superiores de Fourier,
lo nico que ocurrir es que el contorno se
276
Captulo 10
hace ms difuso (fig. 10.1c). Por tanto, al trabajar con elementos poco fiables
es preferible utilizar los coeficientes de Fourier en lugar del contorno inicial:
es mejor guardar todos los huevos en una misma cesta.
En este captulo estudiaremos primero las nociones elementales sobre
el funcionamiento del cerebro, que de una u otra forma se reflejan en los
diferentes modelos, y tambin analizaremos algunos modelos matemticos
fundamentales de este campo.
;i
*c
277
Redes neuronales
.-
278
Captulo 10
279
Redes neuronales
Modelo de Hopfield
Por lo visto, la red neurona1 de Hopfield es hoy el modelo matemtico ms
popular en la neurociencia. Esto est condicionado por su simplicidad y claridad. Adems, este modelo result cercano por su descripcin matemtica a
un objeto investigado detalladamente en la fsica estadstica: el vidrio de espn.
La red de Hopfield muestra de qu modo, en principio, se puede
organizar una memoria en una red de elementos que no son muy fiables.
Prestemos atencin a algunas particularidades importantes de la memoria de los animales y el hombre, las cuales sera interesante expresar
mediante un modelo.
{(F;
{(F}
{(rt),
<rt
280
Captulo 10
~5
Ejemplo.
6: - imagen memorizada por el dispositivo de memoria,
- algn fragmento de esta imagen o la imagen con un ruido superpuesto,
<yt- imagen no distorsionada, almacenada en la memoria;
rr
' :e
fr
' :e
llrr
c:Il
Carcter distribuido del almacenamiento de la informacin. Alta fiabilidad. Partiremos de las ideas ya expuestas sobre el funcionamiento del
cerebro. Por cuanto en el cerebro no se logra indicar un elemento concreto
responsable del almacenamiento de un patrn dado, consideraremos que la
informacin de un objeto es almacenada no por una neurona elemental
de la caja negra, sino por todo un grupo. Adems, sera conveniente
que el dispositivo creado, como en el
caso del cerebro, fuera capaz de funcionar
tambin en el caso de prdida o deterioro
de una parte de sus elementos. Tambin
sera bueno que para la red neurona1
\
creada, la probabilidad de reconocimiento
\
,
\
correcto se comportara, en dependencia
O
50
100 d, %
del porcentaje de elementos deteriorados,
como se muestra en la figura 10.4.
m A
Redes neuronales
281
IlCb ["l.
--
282
Captulo 10
(F,
cuando t -+
v((;) = 0.
En otras palabras, si tanto los patrones de entrada como los patrones
claves fueran nmeros reales, entonces el sistema computacional que resuelve la ecuacin (10.1) podra jugar el papel de caja negra con memoria. Al entrenamiento de semejante
red correspondera la construccin de
la funcin v(x) con puntos singulares estables en los puntos $' ( p = 1,...,M)
del eje numrico; a estos corresponden
los mnimos de la funcin potencial.
La funcin v(x) se puede construir de
tal manera que los puntos singulares
6
52
6
inestables xk (v(xk)= O) tomen valores
Fig. 10.5. La ecuacin (10.1) permite dados (fig.10.5). Ellos definen las regiones de atraccin de los patrones claves
realizar una variante elemental de memona asociativa
y reflejan cules nmeros se parecen
entre s y al patrn clave. Destaquemos un detalle importante: aqu podemos definir por nuestra propia cuenta
no slo los patrones claves, sino tambin sus regiones de atraccin.
Las situaciones que se presentan en la realidad y en las que se aplican
las redes neuronales son, evidentemente, mucho ms complejas. Ante todo,
porque los mostrados y los patrones memorizados t; son vectores de
grandes dimensiones. Por ejemplo, stos pueden ser imgenes, y en 5 estar
grabada la informacin del color de cada pxel de la pantalla.
Sin embargo, la idea general y el planteamiento del problema siguen
siendo los mismos. Si se considera que la investigacin de 10s atractores
del sistema dinmico dado es un problema directo de la teora cualitativa,
entonces aqu estamos tratando con un problema inverso. Para su resolucin
se debe construir un sistema dinmico que tenga una coleccin dada de
atractores (y en el caso ideal, que tambin tenga las regiones de atraccin de
estos atractores). Este problema es un problema mal planteado. Son muchos
y muy diversos los sistemas dinmicos que pueden tener los mismos
atractores (de hecho, son muchas y muy diversas las funciones U ( x ) que
pueden tener extremos en los mismos puntos). Por eso se puede simplificar
el problema considerando, por ejemplo, que estamos tratando con sistemas
dinmicos de un tipo especial, para los cuales existe adems la funcin de
Liapunov. Este enfoque se utiliza en el modelo de Hopfield y en un conjunto
de otras redes neuronales. En este campo de investigaciones, las redes
neuronales que poseen funcin de Liapunov o sus anlogos se denominan
frecuentemente redes convergentes.
Para describir una red neurona1 es necesario:
N
283
Redes neuronales
.-
h(t)> O,
h(t)< O
= -1,
(el caso h(t)= O no se considera). Evidentemente, sign (h) es simplemente
el signo de la funcin h. Se considera que el tiempo es discreto: t = 1,2, .;-:y
Tendremos en cuenta que el potencial de membrana de la neurona i en
el instante t depende del estado de otras neuronas. Para describir esta
circunstancia en el modelo, consideraremos que la accin sobre una neurona
dada por parte de otras neuronas est determinada por la suma ponderada
S(t
1,
284
Captulo 10
= AS,
(-
285
Redes neuronales
, :<
ciI,
< <
aprendizaje de Hebb:
J
t j P para
=-
i # j;
Jii = 0.
p=l
Captulo 10
286
frmula (10.7) para Oi = 0:
donde m p =
&
entonces m, = 1, m, = O, m, = O,
-
N- M
a~
amP
a2
@mpam,)
El valor de la energa
en este caso es aproximadamente igual a cero
-
(- 5 $) ,
-
As pues, las expresiones (10.4), (10.5) y (10.9) determinan completamente la dinmica de la red de Hopfield, en la cual, bajo ciertas condiciones,
a los patrones claves les corresponden atractores.
El test preferido para el modelo de Hopfield es el reconocimiento de
imgenes que representan distintas letras del alfabeto. En la figura 10.6 se
muestran tres patrones de entrada: las imgenes de las letras A, Y y K, con
287
Redes neuronales
.-
step :
step
8I
8l
81
811
U # i @ II
I@
W8
II
lil
a:
II
Ir
1
1
1
am
1:'
BI
UI
91
88
111111
11188
e:
m8
88
R#1188188111818118
P8
II
88
II
1
1
1
aaeu
.1211E+06
UDlb@1061181W110
1111
1116
#8
88
II
8
al
II
Ia
U8
II
88
10
111
e:
.1740E+06
111 r
@M
I D
II 88
1I
II
I1
88
m
r
1)1
48
II
88li
m1
11881818
1118
llllllllli
P I
11111 81
a1
M I
umn
',a
step:
8B
ri
DI
8li
88
11i
811
8U
a#
88
88
I
II
88
18~B0181018188DDIB
i88l#8~l98#0D88I88
il
111
RI
ID
U@
DI
QI
UI
06
II
e:
.1861E+06
au
8I
Fig. 10.6. Ejemplo de reconocimiento de un patr6n por la red de Hopfield. En la fila superior
se muestran Los patrones claves. En la inferior tenemos e l reconocimiento de una letra A
a la cual se le suQerpuso un ruido. E l parmetro step corresponde a l nmero de paso T ;
e son las energas de la configuracin mostrada S ( t ) . Al cabo de vanos pasos la red reconoce
.correctamente el patrn
Captulo 10
288
.-
289
Redes neuronales
.-
step :
M
m+
.
:"".:
.
.
.
:
:
..!
.
m :
1
llLN
Da 8
m .
11 111 11
1 111 111
mmm .B
8 111
.m
m
.
mm m m m
.
1
11
1 Ia :
1
a
BID B
18 D 1
I
18 111 8
m. m 8 1
m
.m
umm
111
111
D
mi1 118 1 1C 11
DIIDIIB
II
11 111
m1 amm
1 1 1
I1
1I
11
m1
m
D
1 11.
88,
1 E
1 1
1 1
m
1 1
e:
m m
'
1
11 m1
1 8 1
e : -.3692E+04
step :
step:
'D
e 1
:
LI
I
D H
- .7544E+05
M
\
e:
- .1530E*06
Fig. 10.7. Efecto de memoria falsa. En La fila superior se muestran los patrones claves conocidos por la red. La fila inferior ilustra la Llegada de la red a una configuracin estacionaria,
independiente del tiempo. Esta configuracin no corresponde a ninguna imagen clave
290
Captulo 10
Fig. 10.8. Cada neurona se modela mediante una aplicacin S(t 1) = F(S(t),V(t))
unidimensional dependiente de la accin exterior V(t) como de u n parmetro: a) para
valores grandes de V la neurona analizada es equivalente a una neurona de McCulloch-Pitts;
b) para valores no muy grandes V > 1, e l estado de excitacin +1 todava es un punto
fijo estable; c) para V < 1 la aplicacin se convierte en una aplicacin de dilatacin y la
dinmica de la neurona se torna catica
Supongamos ahora que el extremo derecho del segmento est fijo, y que el
extremo izquierdo puede girar alrededor de este punto como alrededor de un
eje. Es evidente que si el ngulo entre el segmento y la horizontal no es igual a
cero, la transicin al estado de excitacin ya no ser instantnea (fig. 10.8 b).
En este caso la magnitud S(t) se aproximar exponencialmente a +1,
y cuanto mayor sea e ngulo tanto ms lenta ser la aproximacin.
Denotando mediante K(V) la dependencia del ngulo respecto a la
accin exterior V, la aplicacin F(S,V) en la frmula (10.11) para V > O se
escribe como
F(S,V) = 1 (S - l)K(V).
Respectivamente, para
W < 0,
F(S, V) = sign V
+ (S - sign V ) K ( V ) ,
292
Captulo 10
Fig. 10.9. Forma de la funcin K(V) que define la pendiente de la aplicacin unidimensional
continua a trozos y lineal a trozos que determina la dinmica de una neurona solitaria.
En los cuadrados de arriba se muestra la funcin F(S, V) para los valores de V indicados
en el eje de abscisas. Para IVI > 1 la aplicacin unidimensional es de contraccin, y para
IVI < 1 es de dilatacin
F(S, V) = (sign V
(10.12)
2
l + lvl Y la
aplicacin F(S,V) para valores diferentes de V.
As, la aplicacin unidimensional construida, la cual describe la dinmica de la neurona, se ha elegido de tal manera que para una accin fuerte positiva (V > +1)6 negativa (V < -1) la neurona pasa correspondientemente
al estado de excitacin o de inhibicin,
V=i,l
mientras que para una accin exterior
dbil la actividad de la neurona es ca-\
tica. En la figura 10.10 se muestran diferentes tipos de dinmica de una neurona solitaria para una entrada fija V.
V=-0,7
El entrenamiento, como en el moi9i
\-- delo clsico de Hopfield, se realiza seFig. 10.10. Actividad de una neurona en gn la regla de aprendizaje de Hebb:
u/"--
a-s......
- IZJ
v=- 123
para
i # j;
Jii= 0.
293
Redes neuronales
<<
/K
'
+ +
294
Captulo 10
cp,
295
Redes neuronales
sectores peridicos grandes en los electroencefalogramas son, habitualmente, evidencias de ataques de epilepsia. En muchos casos, el ordenamiento
excesivo es testimonio de enfermedades serias. En la literatura existe ms de
una decena de hiptesis sobre los mecanismos del cerebro que utilizanel
caos dinmico. La generalizacin del modelo de Hopfield permite promover
otra hiptesis ms. El caos puede ser uno de los medios para la lucha contra
la memoria falsa. l puede corresponder a los estados no s, en los
cuales la informacin necesaria para tomar una decisin, o bien el tiempo,
es insuficiente.
296
Captulo 10
CI
< <
a,+, = f (am)*
(10.13)
En este caso, la estrategia de resolucin del problema es la siguiente.
Para simplificar, supongamos que todos los elementos de la sucesin {a,)
pertenecen al intervalo j , es decir, la funcin f aplica este intervalo en s
Redes neuronales
297
a,
AI
a,+,=
f(a,),
a1a2
m=1,
..., N - 1 .
(10.14)
f,hallamos el
< <
298
Captulo 10
< <
7.
299
Redes neuronales
.satisfactorias. Sobre su base se logran crear notables sistemas de reconocimiento, diagnstico y prediccin. Sin embargo, la teora de estos sistemas
que permita estimar a priori qu se puede esperar en cada problema, se
encuentra en las primeras etapas de su desarrollo.
Fig. 10.15. Esquema de un perceptrn elemental que debe distinguir las consonantes
de las vocales
(4,
300
Captulo 10
Dividamos en cuadrados elementales el dispositivo sensible (por ejemplo, un escner) al cual se le muestra la imagen (fig. 10.15). Supongamos
que si el m-simo cuadrado es oscuro (en l hay ms pxeles negros que
blancos), entonces la magnitud x, es igual a la unidad; en caso contrario, es
+
igual a cero. Es decir, t1= (x:, ... ,xh) ,donde el 1en los ndices del vector
corresponden a la primera letra mostrada, t2= (x:, .. . ,
etctera.
La idea del entrenamiento (es decir, el ajuste de los pesos w,) se
puede explicar mediante la parbola siguiente. Supongamos que en cierto
reino, el rey tiene que tomar decisiones responsables con mucha frecuencia,
para lo cual el rey tiene que aconsejarse con sus ministros o expertos. l
no tiene posibilidades de contar con otros expertos, y debe apoyarse en
los que tiene. Cmo aumentar la efectividad de las decisiones tomadas?
Es necesario ver qu consejos han dado los expertos y si fueron correctas la
decisiones tomadas. Sobre la base de esta experiencia el rey puede aclarar
cmo debe valorar a los expertos. Naturalmente, resolver los problemas por
votacin es irracional. Unos expertos deben ser odos y se debe proceder
como ellos aconsejan, mientras que otros deben ser escuchados, pero se
debe proceder al contrario. Las magnitudes {w,) definen precisamente las
diferentes posiciones del rey para con los distintos expertos.
El algoritmo de aprendizaje se puede describir del modo siguiente: si
en el (i 1)-simo paso del entrenamiento el perceptrn dio una respuesta
correcta, entonces los pesos sinpticos no cambian:
wi+l - a
- w,.
4
$N),
<
301
Redes neuronales
.-
regla delta.
-+
(61,
302
Captulo 10
las redes interpolantes, lo cual estuvo relacionado con el paso a las redes
neuronales multicapas y a los nuevos algoritmos de aprendizaje.
El ejemplo clsico de Minsky est relacionado con la funcin XOR
(disyuncin exclusiva). sta es una funcin elemental de dos argumentos,
cada uno de los cuales puede ser cero o uno. Ella es igual a la unidad
nicamente si un solo argumento es igual a la unidad. Denotemos una
entrada mediante x y la otra mediante y. Correspondientemente, los
vectores de entrada ( tendrn dos componentes (&, ty).La muestra de
aprendizaje para esta funcin se presenta en la tabla siguiente:
A estos patrones les corresponden los puntos A*, B,, A,, B1del plano
(ex,tg) (fig*10*16)*
Consideremos ahora un perceptrn elemental. La salida del elemento
sumador es la funcin S = wl& w2tg. La conmutacin del elemento
no-lineal tendr lugar si se cumple la condicin
"
1 Bi
A,*
303
Redes neuronales
(ex,
304
Captulo 10
shpticos entre los elementos de la primera capa iguales a la unidad, y estableciendo el umbral del elemento de la segunda capa un poco inferior a tres.
As pues, el paso de una red neuronal monocapa a una de dos capas
permite en este caso construir el sistema de diagnstico. En la parbola del
rey esto:significa que cada ministro (neurona de la primera capa) necesita
un personal de expertos. stos estn definidos por los pesos {u;) de las
neuronas de la primera capa. Sin embargo, el sistema construido tiene un
t3)
defecto muy serio: las regiones que se pueden realizar en el plano (tZ,
deben yacer a un mismo lado de las rectas que la forman, es decir, las regiones
deben ser convexas. Las regiones mostradas en la figura 10.12, las cuales
realmente pueden surgir en los problemas de diagnstico computacional,
no son accesibles a las redes de dos capas (vase tambin la figura 10.19).
Fig. 10.20. Una red neuronal de tres capas puede resolver e l problema que se ilustra en la
figura 10.19
305
Redes neuronales
Algoritmo de retropropagacin
La idea de este algoritmo es simple. El rey los ministros y los ayudantes
deben, a su nivel, estimar quin realiz la mayor contribucin en la toma
de una decisin errnea. La actitud hacia el culpable se debe cambiar en
primer lugar. dem ms, este algoritmo se puede considerar como una especie
de generalizacin de la regla delta.
dy
= ~ (-1y).
dx
Captulo 10
306
Paso inverso.
Paso 1. Despus de determinar la actividad de los elementos de salida yj,
como resultado del paso directo se calcula el error
--
donde yj y iRj son las componentes j-simas de los vectores correspondientes (omitiremos el argumento para simplificar). Si E = O, los pesos no
cambian. En caso contrario, debemos realizar el paso siguiente.-
Redes neuronales
307
.Paso 4. Estimamos cun rpido cambia el error al cambiar los pesos sinpticos wii :
Paso 6. Ajustamos los pesos wij mediante una variante de la regla delta:
donde n es el nmero del paso (se supone que todas las magnitudes en
las frmulas (10.23)-(10.27) pertenecan al paso n), y 71 es la velocidad de
aprendizaje, un coeficiente que yace en el intervalo 1-0,Ol.
dz
308
Captulo 10
309
Redes neuronales
Preguntas y ejercicios
<<
< <
Capitulo 10
310
Bibliografa
1. Izhikivich E. M., Malinietski G. G. Sobre el posible papel del caos en los sistemas
neuronales (en ruso).
2. Kostyliov 1.A., Malinietski G. G., Potpov A. B. Parmetros de orden en la red
neuronal de Hopfield (en ruso).
ndice de autores
Adams J. C. 62
Andryev A. Yu. 21
Andrnov A. A. 67
Augusti G. 174
Barr S. 245
Bsov N. G. 11
Belosov B. P. 21, 219
Bessel E W. 50
Burgers J. M. 32
Butzov V. E 235
Darwin ch. R. 18
De Broglie L. 25, 26
De Vries C. 32
Dorodnitsyn A. A. 238
Duffing G. 70
Ebbinghaus H. 309
Eigen M. 75
Euler L. 104
Feigenbaum M. J. 121,187,197
Fermi E. 46
Feynman R. Ph. 9
Fibonacci L. 201
Field R. J. 242
Freeman W. J. 290
Fuller R. B. 267
Gardner C. S. 32
Golgi C. 276
Greene J. M. 32
Gunzburg V. L. 171
Gumiliov L. N. 116
Hadamard J. S. 88
Haken H. 233
Hamilton W. R. 25
Heaviside 0. 40
Hebb D. 0.285
Hnon M. 191
Hermite Ch. 50
Hopfield J. J. 279
Hopf E. 32
Hube1 D. H. 278
Hutchinson G. E. 18
Jacobi C. 50,104
Jordan M. 260
Klein F. 245
Korteweg D. J. 32
Krojin O. N. 11
Kruskal M. D. 32
Kurdimov S. P. 8
Lagrange J. L. 104
Laguerre E. N. 50
Landu L. D. 171
aplace P. S. 49
Lawson J. D. 11
Le Verrier U. 62
Lefever R. 217
Leray J. 93
Levandovski M. 1. 21
Liapunov A. M. 117
Lipschitz R. 57
Lorenz E. N. 97, 126, 192
Malthus T. R. 20, 30
Manneville P. 197
Mathieu E. L. 50
May R. 20
McCulloch W. S. 283
ndice de autores
312
.Minkowski H. 245
Minsky M. 301
Miura R. M. 32
Mobius A. E 258, 260
Newton I. 51,204
Noyes R. M. 242
Nuckolls J. 11
Ohm G.S. 25
Oseliedets V. 1. 98
Osgood W. E 87
Pitts W. H. 283
Poincar J. H. 63,104
Pomeau Y. 197
Pontriaguin L. S. 67
Poston T. 128
Prandtl L. 126, 192
Prigogine 1. R. 114, 217
Purkinje J. E. 276
Ramn y Caja1 S. 276
Rayleigh J. 126
Reynolds 0 . 1 1 3
RikitakeT.121
Rosenblatt F. 301
Rossler, O. E. 21
Sjarov A. D. 11
Samarski A. A. 8
Schmidt E. 117
Schrodinger E. 25, 27
Schwartz L. 36,41
Schwarz H. A. 194
Sedov L. 1.46
Shandarin S. E 138
Shumann R. 176
Singer D. 194
Smale S. 154
Sboliev S. L. 36
Stewart 1. 128
Tamm I.E. 11
Teller E. 11
Thom R. 128,143
Tjonov A. N. 89,234
Toynbee A. J. 115,116
Turing A. M. 114
ndice de materias
Aislamiento 120
algoritmo de aprendizaje 283
-- retropropagacin 305
algoritmos de interpolacin 297
Andrnov-Hopf, bifurcacin 216
aplicacin S-unimodal 194
- de contraccin 57
-- Hnon 191
-- un segmento en s mismo 181
- homeomorfa 246
- logstica 20, 178
aplicaciones unidimensionales 176
arista de un grafo 252
armnico de Fourier 192
atractor 65
- extrao 96
autmata celular 15
autooscilaciones 204
axn 276
Bendixson, criterio 209
bifurcacin 104
- de Andrnov-Hopf 216
-- duplicacin del perodo 183
-- nacimiento de un ciclo lmite 216
- tangencia1 198 ,
- transcrtica 110
biocenosis 15
biunivocidad 246
botella de Klein 260
brucelador 114,217,218
buckminsterfullereno 267
bulbo olfatorio 289
Cadena trfica 18
clculo calorfico 15
campo calorfico 15
- vectorial 74
-- rpido-lento 237
campos locales 284
-- interna 236
capacidad de aprendizaje 301
carcter umbral de un proceso 81
caracterstica de Euler 258
- tensin-corriente 224
carcinogenicidad 269
carga topolgica 69
-- de un sistema 69
catstrofe 128
- cuspidal 169
- elemental 170
clula de Purkinje 276
- piramidal276
centro 214
ciclo 253
- estable de una aplicacin 183
- lmite estable 65
- ruidoso 196
- superestable 183, 186
ciclodextrina 267
cinta de Mobius 259
coeficiente de Malthus 30
cola de milano 169
complejidad organizada 309
cornpletitud asinttica 124
componente cuasiestacionaria 233
componentes conexas 252
comportamiento adaptativo 310
coneccionismo 283
confinamiento magntico 11
conjunto abierto 251
- de bifurcacin 134
---de la funcin V 133
--entrenamiento 299
-- informacin 265
- discriminante de la funcin U 133
- vaco 251
constante de Lipcchitz 57
-- una reaccin 74
ndice de materias
- reducida de Planck 26
- de umbral 299
-- condticcin de calor 28
-- dispersin 118
-- Duffing 70
-- Hutchinson 18
-- Korteweg-De Vries 32
-- reposo 277
- refractaria 277
--
-- Malthus 20
-- onda 28
-- ramificacin 118
-- Schrodinger 27
-- Van der Po1 223,225
- en variaciones 222
- generalizada de
Gunzburg-Landu
125
- logstica 177
- modlica 123
- rgida 240
- singularmente perturbada 233
ecuaciones de Navier-Stokes
efecto de la capa T 11
93
ndice de materias
315
.-
-- prueba 36
-- seguimiento 219
- definida positiva 83
- delta de Dirac 35
- generalizada 36
Jacobiano 31
- de uni aplicacin lineal 155
jerarqua de modelos simplificados 15
-, ecuacin 20
mariposa 169
matrices ralas 15
matriz de conexiones 284
--Hesse 163
-- Jacobi 31
-- una forma cuadrtica 148
- mal condicionada 240
May, aplicacin logstica 20
mecedora 131
medida invariante 196
memoria falsa 288
mtodo de aproximaciones sucesivas
56
-- complejizacin 112
-- iteracin simple 179
-- la capa lmite 235
-- Lagrange 298
-- Newton 179
--orden N 55
-- primer orden 54
'
316
fndice de materias
.-
-- promedio 229
- simblico 169
-- retropropagacin 307
operador de Laplace 27
oregonador 242
orientacin del espn 284
oscilaciones glicolticas 243
oscilador anarmnico 48
- de bloqueo 277
papirotazos catastrficos 174
parlisis de la red 308
parmetros de orden 124,288
patrn 273
- de actividad 289
-- salida 279
- espacial de amplitud 290
patrones claves 279
perceptrn 299
perodo de semidesintegracin 30
-- una solucin 65
perturbacin 170
- regular 231
- singular 232
pesos sinpticos 283
plano de Linard 237
-- Poincar 223
pliegue 136, 143,169
poliedro convexo 257
potencial 59, 86, 129
- de Coulomb 28
-- membrana 283
- energtico 10
potenciales de acci6n 276
principio de Huygens 137
-- Maxwell 171
-- retardo mximo 134
- -- subordinacin de los modos de
corta vida a los de larga vida 233
-- superposicin 24
probabilidad de desintegracin 30
problema bien planteado 88
- de Cauchy 88
-- correspondencia 61
-- evasin de las dificultades 174
-- la explosin fuerte 46
-- las casas y los pozos 253
-- los cuatro colores 271
ndice de materias
ndice de materias
.sucesin de Cauchy 57
suma ponderada 283
sumidero lineal 80
superficie bilateral 259
- cerrada no-orientable 260
-- orientable 260
- de respuesta 135
- no-orientable 259
- orientable 259
- unilateral 259
Suzuki-Kuramoto, ecuacin 125
Tamao poblacional 17
tecnicas computacionales 13
tempo-mundos 23
tenedor 110
teorema de clasificacin de Thom 169
-- enderezamiento de un campo
vectorial 67
-- la funcin implcita 105, 156
---- inversa 156
-- Picard-Lindelof 57
-- Stokes 209
-- Tjonov 234
--Vikte 215
teora de cscaras elsticas 113
-- estabilidad hidrodinmica 120
-- imperfecciones 107
-- interpolacin 298
-- la universalidad 187
-- las bifurcaciones 103
--- formas normales 123
-- oscilaciones 66
-- semejanza y dimensiones 45
-- variedades centrales 124
--- inerciales 124
trmino de correccin de orden 227
- secular 227
Thom, teorema de clasificacin 169
tiempo de vida medio 30
tienda india 169
tokamak 11
transformacin de duplicacin 187
-- Fourier de una funcin de
prueba 38
- inversa 38
transformada de Fourier de una
funcin generalizada 39
transparencia de la red 308
transputadora 15
trayectoria heteroclnica 243
- semiperidica 196
Umbral local 283
unilateral 259
Valor crtico de una funcin 162
- principal de la integral de l / x 41
Van der Pol, ecuacin 225
variable de control 131
-- estado 131
variables de fase 73
variedad central 124
- estable 214
- inestable 214
velocidad 74
- de aprendizaje 301
-- fase 74
--reloj 278
ventanas de periodicidad 193
vrtices de un grafo 252
vidrio de espn 279
volumen de fase 206
volvox 269