Algebra II: Tema 10.
1
TEMA 10.- NORMAS DE MATRICES Y
N uMERO DE CONDICIoN
Indice
1. Introduccion 1
2. Norma vectorial y norma matricial. 2
2.1. Norma matricial inducida por normas vectoriales. . . . . . . . . 4
2.2. Algunos ejemplos de normas matriciales inducidas. . . . . . . . 7
3. N umero de condicion de una matriz 8
3.1. Caso de matrices normales. . . . . . . . . . . . . . . . . . . . . . 10
4. Inversas de matrices perturbadas. 11
1. Introduccion
Frecuentemente, el estudio de un sistema fsico pasa por la resolucion de un
sistema de ecuaciones lineales Ax = b (que se supone compatible determinado
en todo este tema). Sin embargo, incluso asumiendo que este modelo lineal
representa perfectamente la realidad, las matrices
A y
b de las que se dispone
no son identicas a las A y b reales, principalmente debido a errores numericos
de redondeo o errores en la medicion de parametros fsicos. As, en lugar de
obtener la solucion exacta x
0
del sistema Ax = b, en realidad se obtiene la
solucion x
0
del sistema
Ax =
b. Naturalmente, interesa tener una aproximacion
de la distancia (en un sentido a un por precisar) entre x
0
y x
0
. El objetivo de
este tema es precisamente profundizar en esta idea. En las proximas secciones
se vera que dicha distancia depende esencialmente de una caracterstica de la
matriz A a la que se denomina condicionamiento o n umero de condicion.
Ejemplo 1. El sistema de ecuaciones
_
10
4
1
1 1
_ _
x
1
x
2
_
=
_
1
2
_
tiene como solucion (con 7 dgitos signicativos)
_
x
1
x
2
_
=
_
1,00010001
0,99989999
_
.
Sin embargo, una resolucion mediante el metodo de Gauss (con solo 3 dgitos
signicativos) proporcionara la solucion
Algebra II: Tema 10. 2
_
x
1
x
2
_
=
_
0
1
_
,
que no parece una aproximacion suciente a la solucion anterior. Este error se
puede reducir mediante pivoteo parcial o total.
Ejemplo 2. Un ejemplo clasico de matriz mal condicionada es la matriz
de Hilbert:
H =
_
_
1
1
1
2
1
3
1
n
1
2
1
3
1
4
1
n+1
1
3
1
4
1
5
1
n+2
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
1
n
1
n+1
1
n+2
1
2n1
_
_
,
que es muy sensible a errores numericos.
2. Norma vectorial y norma matricial.
Denicion 1 Sea (E, K, +, ) un espacio vectorial, con K = R o K = C. Una
norma en E es cualquier aplicacion
| | : E R
que verique K y z, w E las siguientes propiedades:
N1) |z| 0, y |z| = 0 z = 0.
N2) |z| = [[|z|
N3) |z +w| |z| +|w| (Desigualdad triangular).
Ejemplo 3. Tres ejemplos clasicos de norma en E = C
n
son:
norma 1:
|z|
1
= [z
1
[ + +[z
n
[
norma 2 o eucldea:
|z|
2
=
_
[z
1
[
2
+ +[z
n
[
2
norma o norma del supremo:
|z|
= max [z
1
[, , [z
n
[ .
Algebra II: Tema 10. 3
Por supuesto, se debe vericar que cada una de las expresiones anteriores
satisfacen las tres condiciones de la denicion de norma. Detallamos la com-
probacion para la norma 1 a continuacion y dejamos la norma al lector. La
norma 2 se discute un poco mas abajo en el contexto de normas inducidas por
un producto escalar.
N1)
|z|
1
= [z
1
[ + +[z
n
[ 0 y ademas
0 = |z|
1
z
1
= = z
n
= 0 z = 0.
N2)
|z|
1
= [z
1
[ + +[z
n
[ = [[ ([z
1
[ + +[z
n
[) = [[|z|
1
.
N3)
|z +w|
1
=
n
i=1
[z
i
+w
i
[
i=1
[z
i
[ +
n
i=1
[w
i
[
= |z|
1
+|w|
1
.
Ademas, las normas 1 y 2 son casos particulares de la norma p, denida
p N 0 como
|z|
p
= ([z
1
[
p
+ +[z
n
[
p
)
1
p
.
Observacion.- Toda norma || en E induce una distancia d(z, w) := |zw|
en E. En particular, la norma eucldea en R
n
induce la distancia habitual.
Denicion 2 Se dice que dos normas | |
, | |
en E son equivalentes si
existen , > 0 tales que
u E, |u|
|u|
|u|
.
Proposicion.- En dimension nita, todas las normas son equivalentes.
Algebra II: Tema 10. 4
Observacion.- Todo producto escalar , ) en E induce una norma en E,
denida como
|z| := +
_
z, z).
Es trivial comprobar que esta denicion satisface las condiciones N1 y N2.
N3 se demuestra empleando la desigualdad de Schwarz.
En particular, la norma 2 anteriormente mencionada coincide con la norma
inducida por el producto escalar estandar en C
n
, es decir:
|z|
2
=
z
h
z.
Ejemplo 4. Es conocido que la forma sesquilineal (es decir, lineal en la
primera componente y antilineal en la segunda)
, ) : C
mn
C
mn
C
(A, B) tr(B
h
A)
constituye un producto escalar en C
mn
. La norma inducida por este producto
escalar es
|A|
F
=
_
tr(A
h
A) =
_
m
i=1
n
j=1
[a
ij
[
2
,
denominada norma de Frobenius.
2.1. Norma matricial inducida por normas vectoriales.
Denicion 3 Sean | |
, | |
dos normas en C
m
, C
n
respectivamente. Se
denomina norma matricial | | en C
mn
inducida por dichas normas
vectoriales a
| | : C
mn
R
A |A| := sup
x=0
|Ax|
|x|
.
La denicion anterior verica las propiedades de norma. Efectivamente:
N1)
_
x C
n
0,
|Ax|
|x|
0
_
|A| 0,
y
|A| = 0 x C
n
0, |Ax|
= 0 x C
n
0, Ax = 0 A = [0].
N2)
|A| = sup
x=0
|Ax|
|x|
= [[ sup
x=0
|Ax|
|x|
= [[|A|.
Algebra II: Tema 10. 5
N3)
|A +B| = sup
x=0
|Ax +Bx|
|x|
sup
x=0
_
|Ax|
|x|
+
|Bx|
|x|
_
sup
x=0
|Ax|
|x|
+ sup
x=0
|Bx|
|x|
= |A| +|B|.
Nota: Hasta el momento se ha insistido en la existencia de dos normas
vectoriales | |
, | |
(diferentes, en general) y una norma matricial | |.
En adelante, frecuentemente se omitiran los subndices y para evitar una
notacion recargada. El lector debera deducir del contexto que norma es la
aplicada en cada caso (y, por supuesto, si es vectorial o matricial).
Ademas, sera tambien frecuente que ambas normas | |
y | |
sean
del mismo tipo (es decir, por ejemplo, ambas la norma p) pero en C
m
y C
n
respectivamente.
Observacion.- Como consecuencia trivial de la denicion,
x C
n
, |Ax| |A||x|.
Observacion.- Se verica
sup
x=0
|Ax|
|x|
= sup
x=0
_
_
_
_
1
|x|
Ax
_
_
_
_
= sup
x=1
|Ax| .
Observacion.- El cociente
Ax
x
alcanza su supremo en S = x C
n
: |x| = 1,
es decir,
x
1
C
n
: |x
1
| = 1, |Ax
1
| = |A|.
Demostracion:
Por un lado, un subconjunto de C
n
es compacto si y solo si es cerrado y
acotado. El conjunto S es cerrado (por ser complementario de un abierto) y
acotado, luego compacto.
Por otro lado, las aplicaciones
Algebra II: Tema 10. 6
C
n
C
m
x Ax
y
C
m
R
y |y|
son ambas continuas. Puesto que la composicion de dos aplicaciones continuas
es continua, la aplicacion
C
n
R
x |Ax|
tambien lo es. Finalmente, toda funcion continua denida sobre un compac-
to alcanza sus extremos en puntos del compacto, quedando as probada la
observacion.
Corolario.- El cociente
Ax
x
tambien alcanza su supremo en C
n
0.
La observacion y el corolario anterior permiten sustituir el termino supremo
por el termino maximo y escribir
|A| = max
x=0
|Ax|
|x|
= max
x=1
|Ax|,
expresion que se empleara en el resto del tema.
Proposicion.- Si || es una norma en C
nn
inducida por normas vectoriales,
entonces
A, B C
nn
, |AB| |A||B|
Demostracion:
x C
n
, |ABx| |A||Bx| |A||B||x|
|AB| = max
x=1
|ABx| max
x=1
|A||B||x| = |A||B|.
A cualquier norma que verique la anterior propiedad se le denomina sub-
multiplicativa.
Observacion.- Si m = n y las normas | |
y | |
coinciden, la norma | |
en C
nn
inducida por ella(s) verica
|I| = max
x=0
|Ix|
|x|
= 1.
Algebra II: Tema 10. 7
De esto se deduce que la norma de Frobenius no es inducida por ninguna
norma vectorial, pues |I|
F
=
n. (Nota: se puede demostrar que tampoco es
inducida por dos normas vectoriales, aunque estas sean distintas).
2.2. Algunos ejemplos de normas matriciales inducidas.
Ejemplos relevantes de normas matriciales inducidas por normas vectoriales
son los siguientes:
Norma 1: Si las normas | |
y | |
son la norma 1 en C
m
y C
n
respectivamente,
|A|
1
= max
x
1
=1
|Ax|
1
= max
1kn
_
m
j=1
[a
jk
[
_
.
Norma 2: Si las normas | |
y | |
son la norma 2 en C
m
y C
n
respectivamente,
|A|
2
= max
x
2
=1
|Ax|
2
,
tambien denominada norma espectral, y tratada en la proxima proposi-
cion.
Norma : Si las normas | |
y | |
son la norma en C
m
y C
n
respectivamente,
|A|
= max
x
=1
|Ax|
= max
1jm
_
n
k=1
[a
jk
[
_
.
Proposicion.- Sea A C
mn
. Entonces, |A|
2
=
_
(A
h
A), donde repre-
senta el radio espectral.
Demostracion : Recordamos que A
h
A es hermtica, y por tanto diagonali-
zable unitariamente con autovalores reales
1
. . .
n
. Ademas, dichos
autovalores son no negativos: si es autovalor de A
h
A asociado al autovector
u,
|Au|
2
2
= (Au)
h
Au = u
h
A
h
Au = u
h
u = |u|
2
2
=
|Au|
2
2
|u|
2
2
0.
Tambien sabemos que v ,= 0 se verica:
Algebra II: Tema 10. 8
1
= mn
u=0
u
h
A
h
Au
u
h
u
v
h
A
h
Av
v
h
v
max
u=0
u
h
A
h
Au
u
h
u
=
n
.
De este modo,
|A|
2
= max
u=0
|Au|
2
|u|
2
=
max
u=0
u
h
A
h
Au
u
h
u
=
_
n
.
Proposicion.- Cualquier norma || en C
nn
inducida por normas vectoriales
verica
|A| (A).
Demostracion: Sean
1
, . . . ,
n
los valores propios de A, con [
1
[ [
n
[.
Sea u un vector propio de A asociado a
n
. Entonces:
|A|
|Au|
|u|
=
|
n
u|
|u|
= [
n
[ = (A).
3. N umero de condicion de una matriz
Sea A C
mn
, con rg(A) = n m. Se dene el n umero de condicion de A
asociado a una norma | | como
c(A) =
max
x=1
|Ax|
mn
x=1
|Ax|
.
Comentario: El cociente anterior esta bien denido, porque
mn
x=1
|Ax| = 0 x ,= 0, Ax = 0 rg(A) < n.
En el caso rg(A) < n, se dene c(A) = .
El n umero de condicion proporciona una cota superior para el error en la
resolucion de un sistema de ecuaciones. Veamoslo en dos casos.
Perturbacion del termino independiente b.
Se considera el sistema Ax = b, con A C
mn
y rg(A) = n m. Su-
pongamos que tiene solucion unica x
0
,= 0. Se desea estudiar la variacion
del vector solucion x
0
ante variaciones del vector b. De este modo, se
considera el sistema Ax = b +b, y se asume que tambien tiene solucion
unica x
0
+ x
0
. Se pretende estimar |x
0
|. Observamos
Algebra II: Tema 10. 9
Ax
0
= b
A(x
0
+ x
0
) = b + b.
Restando ambas ecuaciones,
A(x
0
) = b.
Denotando M = max
x=1
|Ax| y m = mn
x=1
|Ax|, se tiene:
m = mn
x=0
|Ax|
|x|
|Ax
0
|
|x
0
|
=
|b|
|x
0
|
|x
0
|
|b|
m
M = max
x=0
|Ax|
|x|
|Ax
0
|
|x
0
|
=
|b|
|x
0
|
1
|x
0
|
M
|b|
.
Multiplicando ambas inecuaciones,
|x
0
|
|x
0
|
M
m
|b|
|b|
= c(A)
|b|
|b|
.
En conclusion, el n umero de condicion proporciona una cota superior
para la amplicacion del error relativo. Ademas, se puede demostrar que
esta cota es optima, es decir, en la anterior inecuacion se alcanza la
igualdad para alg un x
0
.
Perturbacion de la matriz A.
De nuevo se considera el sistema Ax = b, que (se asume) tiene solucion
unica x
0
,= 0. Ahora se desea estudiar la variacion del vector solucion x
0
ante variaciones de la matriz A. De este modo, se considera el sistema
(A + A)x = b. Supongamos que este nuevo sistema sigue teniendo
solucion unica x
0
+ x
0
, con x
0
,= 0. En este caso se verica
|x
0
|
|x
0
+ x
0
|
c(A)
|A|
|A|
.
La demostracion se omite. De nuevo, el n umero de condicion proporciona
una cota superior para la amplicacion del error relativo.
Algebra II: Tema 10. 10
Observacion.- Si A C
nn
es regular y || es una norma matricial inducida
por una norma vectorial, entonces
c(A) = |A||A
1
|.
Demostracion: Si A es regular,
|A
1
| = max
x=0
|A
1
x|
|x|
= max
y=0
|y|
|Ay|
= max
y=0
_
|Ay|
|y|
_
1
=
_
mn
y=0
|Ay|
|y|
_
1
=
_
mn
y=1
|Ay|
_
1
,
de donde se obtiene inmediatamente la expresion anterior.
3.1. Caso de matrices normales.
Si A C
nn
es normal e invertible, U C
nn
unitaria tal que
U
h
AU = D = diag(
1
, . . . ,
n
),
con 0 < [
1
[ . . . [
n
[. Ademas
U
h
A
h
AU = (U
h
A
h
U)U
h
AU = D
h
D = diag([
1
[
2
, . . . , [
n
[
2
)
[
1
[
2
= mn
u=0
u
h
A
h
Au
u
h
u
v
h
A
h
Av
v
h
v
max
u=0
u
h
A
h
Au
u
h
u
= [
n
[
2
,
para todo v no nulo. De aqu se deduce
|A|
2
= [
n
[ y
c
2
(A) =
[
n
[
[
1
[
,
donde c
2
(A) representa el n umero de condicion asociado a la norma 2, que
tambien se denomina n umero de condicion espectral.
Algebra II: Tema 10. 11
4. Inversas de matrices perturbadas.
Sea A C
nn
regular, A una perturbacion de A tal que A+A tambien
es regular. Entonces,
|(A + A)
1
A
1
|
|(A + A)
1
|
c(A)
|A|
|A|
,
donde | | es una norma matricial inducida por una norma en C
n
.
Demostracion:
A
1
= IA
1
= (A + A)
1
(A + A)A
1
= (A + A)
1
_
I + AA
1
_
(A + A)
1
A
1
= (A + A)
1
(A + A)
1
_
I + AA
1
_
= (A + A)
1
_
I I AA
1
_
= (A + A)
1
AA
1
.
Tomando normas,
|(A + A)
1
A
1
| |(A + A)
1
||A||A
1
|
|(A + A)
1
A
1
|
|(A + A)
1
|
|A||A
1
| = c(A)
|A|
|A|
.
