Universidad Nacional de Ingenieria
Universidad Nacional de Ingenieria
Universidad Nacional de Ingenieria
FAClJLTAD DE CIENCIAS
1876
INFORME DE SUFICIENCIA
LICENCIADO EN MATEMÁTICA
PRESENTADO POR
LI�IA - PERlJ
2003
Dedicatoria
A mis padres que me dieron la posibilidad del conocimiento, a mis hijos que con
su sonrisa me hacen creativo� a mi esposa por ser el amor de mi vida y a mis
amigos que ayuda.ron a iniciar y te1n1inar esta. tarea..
Contenido
Introducción ··················································�·-················································· l
1.1. Preliminares.................................................................................................. 8
Anexo.................................................................................................................... 57
Bibliografía .......................................................................................................... 62
Introducción
Los problemas de valores propios son una clase especia] de problemas con
valores en la frontera que son comunes en el contexto de problemas en la
ingeniería que provocan vibraciones� elasticidad y otros sistema oscilantes.
Además se usan en una amplia variedad de contextos en ingeniería que van más
alla de problemas de valores de frontera. Antes de describir un método nútnerico
para resolver estos problemas� presentareinos alguna info11nación general como
antecedente, esta incluye el análisis de la importancia tanto matemática como
ingenieril de los valores propios.
Dado una matriz cuadrada A, en el presente inforrne se estudia una técnica para
el cá1cu1o de los vectores no nulos x y los escalares l que cumplen la siguiente
propiedad
A X-AX (1)
k(A x)=k(lx)=l(kx)
A( k x ) == l ( k x )
!·
Ai = A x ; .�' •i = l (2�)
du(t)
dt
= 2u(t) + 6v(t)
dv(t)- -- (3)
-2u(t) - 5v(t)
dt
2
(4)
Áe ux1=2e x1 + 6e x1 u ll
2x1 + 6x2 == h1
(5)
-2Tt - 5x2 == h2
2 6 ...fJ ...l'J
=Á
-?- -5 ...t2 ..t2
2 6
..4 ==
-2 -5
3
- - 1 , .x = 2 -3
y l = -2 , x -
-] 2
i,(t) 2 -3
(6)
v(t) -] 2
4
dada la condición inicial
7 -== C¡. 2 -3
-4 -1 2
luego:
2 -3 CJ 7
-1 2 -4
u(t) 2 -3
v(t) -1 2
5
depende de los signos de los \lalores propios.
Los valores propios están relacionados con las propiedades de las estructuras.
Cuai1do una fuerza o una carga es aplicada a lo largo de una viga cuyos soportes
estan en sus extremos, uno fijo y e1 otro con desplazamiento.
• • • •
Figura 01
C11ando la carga llega a un valor critico l,. A medida que la carga se incrementa,
Figura02
_ c/2x(t) = .AX(I)
O<t<I
- - (9)
dt 2
�
6
·'
d2u(t)
t 2
-
�- X;+J - 2x; +
2
X;-J
d (Af)
2 -l O ... ... O
-l 2 -1 O O O
A=
o ••
•
•
••
o o o (10)
o ••
•
••
o o
•
o o o ••
2 -·1
O O O · .. -l 2
7
Capitulo l.
Definiciones Previas
I. l Preliminares
r = a - ifl
2 + p2
!YI = Ja
-·-
Observe que r == r Y r r == lrl 1
.
X = (X 1, .. ., X D
)t
L x; Y;
n
Observe que
9
vector columna), entonces x* =(x)' es una matriz de 1 x n (o vector fila), y
n
y•x = {x,y} = Lx; .Y;
i=I
n n
2
x*x = {x,x} = llxll�= Lx; x; = ....... 1
iz: 1 i=I ;
Tv =l v (l.l)
(1.2)
10
cor1esponde al vector propio v y, viceversa. Entonces para todo ro= a v con a *
O, tiene A m == A OJ esto es todo múltiplo de un vector propio también es un
vector propio.
801
11
multiplicidad.
Para deter111inar los vectores .propios asociados a un valor propio hay que sustituir
en el sistema de ecuaciones lineales (1.3) por el valor propio concreto Á;. Puesto
que el siste1na es compatible indete1111inado lo cual esta garantizado por el hecl10
de que det(A - U) = O las ecuaciones no son line.almente independ.ientes (L.I.) y
la for1na de resolver el sistema, en la práctica, consiste en dar valores arbitrarios
para tantas componentes de los valores propios como indique el valor de la
multiplicidad.
Por último, 11na vez determinados los vectores propios estos pueden
no1111alizarse, ·por ejernplo, tomando la no11na euclídea:
N
�V�=]
L.J 1
i=]
-
0= QJ VJ + ·· · + an V,,
n
Teorema 1.4. Si { v1 , ••• ,vn } es un conjimto de n valores LJ. en R entonces
11
Demostración
12
Supongamos que A es la matriz cuyas columnas son { v 1 , ••• , vn } , entonces el
Teorema 1.5 Si A es una matriz y Á¡,...;,n son distintos valores propios de A (es
decir l ¡ * li si i=1=-j) con los vectores propios asociados { v1 , ... , Vn } entonces {
Demostración
-
Sean VJ , •.• , Vn los vectores propios co11espondientes a los valores propios
! 1 , ••• ,1 ,. :
-
/J i V¡ + .... + fJn V n =O (*)
-
A(/J1 \1J + ·· · + /Jn Vn ) =O
13
·,
-
P1A Vt + .... + /JnA Vn =O
-
/J1A1V1 +· .. ·+fJnAnVn ·=O **
( )
-
(A1-An)/l1 VJ +· .. ·+{An-1-ln)Pn-l Vn =O
pero A¡- A n * O , para todo i= 1 7 ..... ,.n-1 � /J¡ = O para todo i = 1 �---,n-1, luego de (*)
se puede afirmar que /Jo = O por consiguiente el conjunto { v1 , •••, Vn } es
linealmente independiente .
Demostración
- /J1 - /Jn
=> \1¡ =- VJ
- ··· -
Vn
/J; /J¡
14
ortogonal :=)
- - /Jn -
V,· •V·I =- /J1 -
VJ
-
•V¡ � ... - -
Vn • V¡ =O
/J¡ p¡
lo cual es imposible por ser cualquier vector del conjunto { v 1 , ••• , vn } diferente
de cero. Por consiguiente el comjunto {111 , .•• , Vn } es linealmente independiente.
nxn
Definición J .8. Se dice que dos matrices A , B E R son similares, si existe
DXD
Demostración
-
Supongamos que x * O es tal que
S-'BSi-Ax-lx
B(S x) = l(S x )
-
Puesto que x * O y S es no singular se tiene Si =1= O. Por tanto Sx es un vector
propio de B asociado al valor propio l.
15
det{A - U)= n (a¡¡
n
i=J
- l)
Teorema 1.10. Sea A una matriz de orden nxn. La matriz A es simétrica si, y
sólo si, A es diagonalizable mediante una matriz ortogonal. Es decir, si existe
una matriz ortogonal Q tal que Q'AQ es diagonal
Demostración
D = Q'AQ
A 1 = QD'Q' == QDQ' = A
16
Luego la matriz T'AT se puede escribir .por bloques como
l1Z
T'AT= (1)
O B
T'AT= (T'A 1) 1
=
lt o
(2)
Z' B'
R=
l o
o Q,
)�1 0
R'T'ATR=R' R
O B
-
1 O l1 0 1 O
o Q\ O B 0 Q1
-
-
17
Por consiguiente A es semejante a una matriz diagonal, por lo tanto es
diagonaliz.able.
Puesto que se puede construir una base de vectores propios ortono11nales de una
matriz simétric� diremos que toda matriz simetrica es diagonalizable por una
matriz ortogonal. Es decir si A es simétrica entonces existe una matriz ortogonal
Q tal que Q'AQ =Des diagonal__
Teorema 1.11. Sea A una matriz simetrica de orden nxn, entonces los vectores
Demostración
(A u , v} = (A u)' \1 = U
1
A' V
= u' (A v) = {u � v} = µ(u , v)
18
Lema 1.12. Toda matriz simétrica tiene., al menos un valor propio real.
Corolario 1.13. Si A es una molriz simetrica de orden nxn, entonces los valores
propios de A son números reales., y existen n vectores propios de A que forma,1
un conjunto ortonormal.
Demostración
D = P- AP � AP
1
= PD
x Ax>O, V x �R -{O}.
t
Teorema 1.15 Una matriz simétrica A es definida positiva si y sólo si todos Los
19
valores propios de A son positivos.
Demostración
1
o <x A X
Así que l>O. En consecuenc� todo valor propio de una matriz definida positiva
es positivo.
-
linealmente independiente. Asi para todo x =1= O existe un conjunto í1nico de
constantes fJ 1 , • • • ,fJn no todos ceros tal que
X ¿p ¡ Vi
i=I
n n n n
i i j i
LP¡ A
1 v LP¡l¡ v
x' A x - x =x• = LLP; P¡l¡{v )'v
i=) i=J j=I i=J
1
Pero los vectores v , ••• ,v" fo111,an un conjunto ortonor1nal y, por tanto
o, i * j
l, i = j
20
y como l¡>O =>
n n n
= L L /J /j¡l¡ (v
j=l i=I
j
j
)t v i
= L
i=I
l ¡p f >o
Hemos dicho que dos mat1·ices A y B son semejantes entre si� ct1ando existe una
matriz no singular P tal que B = PAP-1 • La impc,1 tancia de este concepto surge
de) teorema que establece qt1e dos matrices que representan a tina mis1na
t,·ansfo1111ación lineal, con respecto a dos bases distintas, son semejantes entre sí.
Teorema 1.16 Todas las matrices semejantes tienen los mismos valores propios.
Demostración
B = PA.P-·1
21
=det P det (A - ll) det p- 1
=det(A-ll)
Hemos incluido dos hechos básicos en esta prueba: el detenninante del producto
de dos matrices es el producto de sus dete1111inantes; el dete1111inante de la inversa
de una matriz es el reciproco de su dete11ninante.
El teorema 1.16 sugiere una estrategia para encontrar los valores propios de A:
convirtiendo A en una matriz B por medio de una transf01·111ación de semej�
B=PAP-1, y calculando los valores propios de B. Si B es más simple que A, el
cálculo de sus valores propios quizá sea más simple. Específicamente si B es
triangular los valores propios de B ( y los de A) son simplemente los elementos
diagonales de B. Esto nos lleva de un modo natural al teorema de Schur, seg1'10 el
cual la estrategia precedente siempre será posible (al menos teóricamente·).
Recuerde que una matriz U es unitaria si uu· = l. Aqui, u* denota a la
conjugada transpuesta de U: (U*) ij = ¼í. El siguiente es el teorema de Schur.
triangular
Demostracion
22
un vector propio correspondiente . No habrá pérdida de generalidad si suponemos
que Uxll, == l. Como es usual� denotamos con e == (1, O, ...,O) . Sea fJ == ,x 1, si
1 T
k
IXI 1
/J x == u* e de modo que
1
1
Esto demuestra que la primera columna de UAU* es l e • Sea A la matriz que
se obtiene al eliminar en UAu·· la primera fila y la primera columna. Por
hipótesis de inducció� hay una matriz unitaria Q de orden n-1, tal que QAQ* es
triangular. La matriz unitaria con que se reduce A a una fo11na triangular es
V=
1 O
u
o Q
pues
1 O -
1 O l. W 1 O
VAV*= UAU* -
o Q o Q O A o Q*
1 O 1 WQ* WQ*
- -
o Q O A Q* O QA Q
En esta ecuación, w es ttn vector fila de orden n-1, y los ceros representan
vectores nulos de orden n-1. Como la matriz
23
1 O
o º
es unitaria.
Demostración
Si A es htrnlltiana, entonces A= A•. Sea U una matriz unitaria tal que UAU* es
triangular superior. En tal caso (UAU*)* es triangular infe1ior. Pero
Así� la matriz UAU* es ambas cosas, triangular superior y triangular interior; por
consiguiente, se trata de una matriz diagonal.
Demostn1ción
.. - . .
Para que la matriz U = I - vv sea un1tana se requiere
• * •
I = UU = (I - VV )(I - vv )
24
Una condición necesaria y suficiente para que v es que
-·v .
v = 2 o vv = O.
Lema 1.21 Sean x y y dos vectores tales que llx 11 2 == IIY 11 2 y (x , y) es real. En
este caso hay una matriz unitaria U de la forma l - vv • tal que Ux == y.
Demostración
.fi
Hx-y R 2
- --* -
Ux - y = (I - vv* ) x - y - X -VV X - y 2
y X - -a ( X - y) (X* - y* ) X
=
= (X - y) [ 1 -a2 (X* X - y* X) ]
Esta últirn.a expresión es igual con cero� pues el factor entre paréntesis es nulo.
,. - -* -
Para ver esto, apoyémonos en las hipótesis x * x = y* y y y X-X y� pues
<x,y> es real para calcular:
--·
x+x
--·
x-y
- -·
x-y x)
25
= 1 - 1 a-"(X-* - y_$ )(X- - y-)
2
Teorema 1.22 El espec'tro de una matriz A den x n (es decir el conjunto de sus
valores propios) está contenido en la unión de los siguientes n discos, D¡, en el
plano complejo:
e
D n
D¡= zE : lz-a;; 1 :S I)ay 1 (l:Si<n)
j=I
- -
p:t
Demostración
Sea A un elemento del espectro de A. Elija un vector x tal qL1e llxll 00=l y Ax=l x.
Sea i un indice para el cual lx¡1=1.. Como (A x) ¡=l X¡ , tenemos que
26
D
l X¡ = :Eaij Xj
j-=1
En consecuencia
o
(l-a¡¡) X¡ = :Eaij Xj
- .
j=I
J*I
n n
ll-a¡¡ 1 :S :Elaij Ux¡ 1 :S :Elaij 1
j• =i• - -
j=I
J*' J�
!·
de modo que l � O¡
-1 O J_
4
A== .l.
4
l J_
4
1 l 3
se muestran en la siguiente figura, donde podemos deducir que todos los valores
propios de A satisfacen la desigualdad½ :S lll :S 5.
27
Eje imaginario 2
-
l
Eje real
-2 J 4
Discos de Gershgoñn
Demostración
28
Aplicando el teorema de Ge1shgorin a D + C� con C = p-1 BP., concluimos que el
espectro de A+B se encuentra en la 11oión de los discos de Gershgorin de D+C,
los cuales son
n o n
Este último teorema indica que c11aodo se perturba 11na matriz A, sus valores
propios se ven perturbados por 11na cantidad que no excede de Xoo(P)IIBll 00 ,
donde B es la matriz de perturbación y xoo(P)es el número de condición de P con
la DOI 1118 11 • II oo •
Para las matrices hermitianas (es decir� aquellas para las que A·=A) se puede
elegir una matriz P unitaria; se trata en realidad del corolario 1.18 del teorema de
Suchr. por lo tanto, en este caso las ftlas de P son vectores que satisfacen 11 x li 2 =
1. De lo anterior se obtiene que IIP 11 00 ::; /n. Esto tambié11 es cierto para p-1, de
modo que xoo (P) � n. En consecuencia, para toda matriz B, los valores propios de
A + B están situados en la t1oíón de discos
M C : ll-l¡I � x�(P)IIB 11 oo
29
donde l 1 ., .....ln son los valores propios de la matriz he1111ítiana A.
30
Capitulo II
Método de la Potencia
Es un método iterativo para aproximar el valor propio (dominante). Supongamos
que la matriz A de orden nxn tiene n va1ores propios A.1,·· ·,Á n con un conjunto
asociado de vectores propios linealmente independientes {v >
< 1 . n)
, •. , v } .
<
L
D
x<º> = /J¡ v
i
(/Ji * O) (2.1)
i=J
x(l} = A x<0> -
•
•
•
luego:
D
::)X (t)
= I':JJ¡Ak(v i )
i=J
= LP¡ lf v <i)
i==l
t k
V(2) + ... + Rn
An
----
-n
p A1 V
(2.2)
Sea
k
k .,,c2) Án -(n)
+ ... + fJ
n At
V (2.3)
(2.4)
k
Lim =0 (2.5)
k-+00
n
<1> ) O.
Sea 'P una funcional lineal sobre C tal que 'P ( v =t-
32
(2.7)
Sea
'P(x<A:+1) ) _ 1
11+ ['P ( v( I> ) + 'P ( t < A:+J) ) ]
- (2.8)
(
'P ) )
x < A:
11['1'( v (I)
) + 'P(.: <A:)
) J
de (2.6) y (2.8) se obtiene
de (2.4)
(2.10)
de (2.1 O) y (2.6)
1,(I)
33
Entrada: A€ 9l ffXII
xcmn
M : ní1mero de iteraciones
Para k = I, 2., .. .. .,M hacer
y:= Ax
r ·•= 'P(Y)
\J'(..T)
X :=y
Salida: r (aproximación del valor propio)
x (aproximación del vector propio)
1
Para asegurar la convergencia de x en la práctica, esto es para evitar que tienda
al vector cero o crezca sin limite, en el metodo de la potencia se introduce una
no11nalin1ción de los vectores x<.t)
Entrada: A€R
XE R
M : 01'1roero de iteraciones
Para k == 1,2, ...,M hacer
y:== Ax
r : == 'P(y)
'P(x)
X :=
IIY 1111)
Salida: r (aproximación del valor propio)
x (aproximación del vector propio)
34
•·
(l+I)
r<lc) = Xp¡
(i)
Xpi
lim r >
< 1 = l1
r-+«:;
Para estabilizar el cambio de indices PA. y para k grande y evitar overflow, en cada
paso se elige
·- •
p1c .. - m 1n V
f -€ { l , 2 , ..... , n x
} 1 (i)
j 1 -
- 11-x 11 � }
(k)
35
n
Entrada: AER
n
xER
M : número de iteraciones
to/ : tolerancia
p : = min{;E{l,2, ... ,n} : l:rj l = llxll 00 };
-
X .·= X .
X ,
p
Para k == 1,2, .. . ,M hacer
y:== A i;
r :== yp;
p : = min{íE{l,2, ... ,n} : lYj l = II.Yll 00 }
Si _Yp == O, entonces escribir(' criterio falla')
escribir(' correr nuevamente')
Parar!
..
SlDO
-
X:= Y .
YP,
Er r : = llx - ) y ILxi
P
¡-1 es un valor propio de A-1• Esto sugiere un mecanismo, para calcular el valor
propio más pequeño de A como sigue:
36
Luego los valores propios de A-1 son tales que
/\
x(k+J > =A -(k)
X
1
Evitando el cálculo de A- , en cada paso se debe dete11ninar x < k+t)
solucionando
el sistema lineal
37
2.2.1. ALGORITMO DEL METODO DE LA POTENCIA INVERSA
Entrada: A€R
n
X€ R
M : número de iteraciones
to/ : tolerancia;
p := min{;E{l,2, ... ,n} : lxi l = llxll 00 };
X .·= X .
Xp ,
X- 1-11
.X •
2 ,
x-=y;
E" := 11 x -i 2
k+1 1
1
Si Err < tol entonces escribir('vector propio aprox', x )
escribir('valor propio aproximado', r)
Parar!
fin_si
fin_para
Sea Ák el valor propio de A tal que O< IAt - µI < E, donde µE C ha sido
deter1ninado con anterioridad y
IA.j - µI > E, para las otros valores propios Áj * At de A. Como los valores
38
propios de A - µ/ son los ní1meros lj - µ � el método de la potencia inversa se
1
puede aplicar a (A - µI) para obtener 11na aproximación de /3 = (Ak - µ)- • Para
dete11ninar x < k+I)
de
Entrada: AER
X€ R
M : número de iteraciones
to/ : tolerancia.;
•·=
X X .
Xp ,
39
2.4. METODO DE LA POTENCIA CON DESPLAZAMIENTO
40
2.4.1. ALGORITMO DE LA POTENCIA CON DESPLAZAMIENTO
Entrada: AER
xER
M : número de iteraciones
tol : tolerancia
µ : para deter111inar un valor más proximo
p : = min{i€{l,2, ... ,n}: lx¡I = llxll 00 };
X.·= X .
X ,
p
Para k == 1,2, . . . ,M hacer
y:= (A - µ/) x;
,. := yp;
p : = min{iE{l,2, ... ,n} : lvj l = IIYII C(J }
Si YP = O, entonces escribir ("criterio falla')
SIDO
x:= Y .
Xp,
41
lim
:
entonces {P m} =0 convergera a P, oon orden a con 11na constante de error.
DEFINICIÓN 2.6. Supongamos que {/J 0 }: 1es una sucesión que converge a
cero y que {a n }:: 1 converge en un número a. Si 3 una constante positiva K con
1 a n- a 1 :S Kl /Jo I para un valor mas alto de n, entonces diremos que {aº}� 1
converge en a con un rápidez de convergencia 0(/Jn)
(a 0 =a+O(/J ó ªº - a con rapidez de convergencia 0(/Jn))
0)
(*)
ª
2.5.l. PROCEDIMIENTO DE Al'l K'4:N A 2
42
es suficientemente grande como para que
entonces
por tanto
=>P=------
Po+z- 2Pn+I+ P n
P=----------------
(Pn+ 1- p D) 2
P= .Pn - -------
Pn+2- 2P nt-1+ p n
43
A
El metodo ll 2 de Aitken se fundamenta en que )a sucesión P n defmida por
/\ 00
44
2.5.2. ALGORITMO ACELERADO DEL MÉTODO DE LA POTENCIA
Enbada: AE R
X€ R
M : número de iteraciones
to/ : tolerancia
k-== I·,
µ0=0;
µ1=0;
p := min{;E{l.2•...• n} :
-
lxj l = Hxll 00 };
X•·= XXp ,.
Err := 11 X - )P Y 00 ll
Si Err < tol entonces escribir(_"vector propio aprox", x )
escribir(_'valor propio aproximado', r)
Parar!
jin_si
jin_si
µo = µ1;
µ1 = µ;
jin_pa,ra
45
Nota: ver programa en el anexo
2
volver a aplicar el procedimiento 6. de Aitken.
46
Entrada: A€R
i ER
M : número de ite.raciones
to/ : tolerancia
k = 1 ,·
-
x:= 11f112;
Para k = 2, .. . .,M hacer
y: = Ax;
-T -
µ =x •y;
Si IIY 11 2 == O, entonces escribir('criterio falla')
Parar!
•
SIDO
-
X:= y .
11.Y 11 2 ,
Err := tlx - ) y 11P «>
47
TEOREMA 2.9. Si A es una matriz simétrica de nxn con los valores propios
l 1 , ••• ,ln y 11 A x - l x 11 2 < E para algún vector x con II x 11 2 = 1 y con el número
rea) A, entonces
1 n
Demostración: supongamos que v ,... ,v fo1,oan un conj11nto ortogonal de
vectores propios de A asociados, respectivamente, a los valores propios l1�-··),n.
e11tonces para algún conjunto único de las constantes P 1,...,/Jn, podemos expresar
- .
x como sigue:
x= LP.i j
v
j=I
Por tanto,
2
LI
n n n
2
/Jj (Aj - l) vj = Ll/Jj 1 llj- l 1 � mio l)�j-ll LI/Jj 1
2 2 2
pero
o
LIPj1 =
2
nx 11 � = 1
j=J
,
as1 que
48
-
Capítulo 111.
3 .1. Introducción
Describimos el siguiente sistena fisico que puede servir como contexto para
examinar este tipo de problema.
p p
y
(0..0)
•
t
y •
-------
M
X
._ t•
"'
P
(LO)
b)
X o)
La curvatura de una columna esbelta sujeta a una carga axial P se puede modelar
por
(1)
___ _
d 2y 2
�+P y = O
dx2 (2)
donde
(3)
y(O) = O (4)
y(L) = O (5)
donde A y B son constantes arbitrarias que serán evaluadas por medio de las
condiciones de frontera. De acuerdo con la primera condición (4)
50
De acuerdo con la segunda condición (5)
O = Asen(pL) + B cos(pL
. )
Pero, puesto que B = O, Asen (pL) ==O. Ya que A = O representa una solución
trivial, concluimos que sen(pL )==O • Para que esta igualdad se cumpla
Así existe un número finito de valores que cumplen con las condiciones a la
front� la ecuación (7) se _puede resolver para
D1C
p= L (8)
l - l 1 l
o) n = l b) tF=2
La figura muestra la solución para los primeros cuatro valores propios, puede
proporcionar el significado fisico de los resultados. Cada valor propio
cu11esponde a una manera en la que se pandea la columna. Combinando las
ecuaciones (3) y (8), se obtiene
51
n 2n:2EI
p = L2 paran=l,2,3, ... (9)
ésta se pueden tomar como cargas de pandeo, pues representan los niveles a los
cuales se mueve la columna en cada configuración de pandeo sucesiva. En un
sentido práctico, usualmente el primer valor es el interés, ya que en general, las
fallas ocWTen cuando primero se pandea la columna Así, una carga crítica se
puede definir como
Yi-1- 'l:y¡+Yi-1
-ly¡ =0
h2
-1
h2
o
-1 -1 =0
b2 h2
(11)
o -)
h2
52
Ax == lx (12)
ésta es la base para una técnica de solución iterativa que fmalmente da e) valor
propio más grande y su vector propio asociado.
3.556x1 -1.778x2 o Xt
53
potenci_ acelerada. nb
,',
3.556 -1.778
..... ri 7- _;.;._... =
.:..:.:it -1.778 3 .�::>o
... - ,,. -i."778
-1.778 3.SSó
Q'JC,2. '"l 6}
r
t
(Q
l • -· L. ..J .J "'- ' {l. / } U(S]=S.90843 �.00338
::o
'..\
0911 ..... ,.."i}
, _.,, --,Q:)¿_ ,
{ -
1- .
, ,.._--,.,..
-iUjjQj' {1.}} U(8)=6.03767 6.0704
•'-" -M<l'.
{-1 • Jr, ...,.�, {1.}} U(9]=6.05117 6.J7068
' U(iOJ=ó.05913 6.07058
·�- l .¼,1 J..v
' '"' -t
1,- .
:, {1.}}
, .. .. .. "" ,,.
6.07047
. . U[- :::.1-"" 0 9 9
l,..¿ O
"' ... • 4 - - ,.. ..
1 - 1 . -i.J....;:1J · , -0 . O
-�mor002c1cr. "'O
,... .
.... __ ... -
v r u..-;,o.:; ....
¡ ____ , o· n . • • •
..�J< - L � .• ) = '· { o . o
,...
- •......, •¡
, 0 r
V ? ---,1}}
.. ,-: "
-.J - '- - � ...)
Conclusiones
o El método de la potencia también puede ser aplicado para hacer el calculo del
valor propio má s peq ueñ o, par a ell o cal cul am os la ma triz inversa de A y (A- 1), y
· e -- 1 va 1o r pr o .
pio 1
.
el m et od o de la po te nc ia co nv rg era a T
© También podemos hacer el cálculo de los valores propios aproximados
inte11nedios, reemplazando la matriz original por 11na que incluya sólo los valores
propios restantes
56
Bibliografia.
...coHEN ALAN M.; cu·rrs J.F; FIELDER R; JONES O.E.; RIBBANS J.;
STUART E. Análisis Numérico. Reverté S.A., Barcelona (1977).
'Título:
Decano o su representante:
" (firma)
..
Profesor Especialista:
(firma)
Dado en la ciudad de Lima en la fecha: 0-6 l_éJ //<!J 's a horas: ¡9-: 15