Loading [MathJax]/jax/output/HTML-CSS/jax.js
Saltar al contenido principal
Library homepage
 

Text Color

Text Size

 

Margin Size

 

Font Type

Enable Dyslexic Font
LibreTexts Español

9: Adecuación de soluciones

( \newcommand{\kernel}{\mathrm{null}\,}\)

Objetivos de aprendizaje

Después de leer este capítulo, deberías poder:

  1. conocer la diferencia entre sistemas de ecuaciones mal condicionados y bien acondicionados,
  2. definir y encontrar la norma de una matriz
  3. definir y evaluar el número de condición de una matriz cuadrada invertible
  4. relacionar el número de condición de una matriz de coeficientes con el mal o bien condicionamiento del sistema de ecuaciones lineales simultáneas, es decir, cuánto se puede confiar en la solución de las ecuaciones lineales simultáneas.

¿Qué quiere decir con sistema de ecuaciones mal acondicionado y bien acondicionado?

Se considera que un sistema de ecuaciones está bien acondicionado si un pequeño cambio en la matriz de coeficientes o un pequeño cambio en el lado derecho da como resultado un pequeño cambio en el vector de solución.

Un sistema de ecuaciones se considera mal condicionado si un pequeño cambio en la matriz de coeficientes o un pequeño cambio en el lado derecho da como resultado un cambio grande en el vector de solución.

Ejemplo 1

¿Este sistema de ecuaciones está bien condicionado?

[1223.999][xy]=[47.999] 

Solución

La solución al conjunto de ecuaciones anterior es

[xy]=[21]

Hacer un pequeño cambio en el vector del lado derecho de las ecuaciones

[1223.999][xy]=[4.0017.998]  

da

[xy]=[3.9994.000]

Hacer un pequeño cambio en la matriz de coeficientes de las ecuaciones

[1.0012.0012.0013.998][xy]=[47.999]  

da

[xy]=[3.9940.001388]

Este último sistema de ecuación “parece” mal condicionado porque un pequeño cambio en la matriz de coeficientes o el lado derecho resultó en un gran cambio en el vector de solución.

Ejemplo 2

¿Este sistema de ecuaciones está bien condicionado?

[1223][xy]=[47]  

Solución

La solución a las ecuaciones anteriores es

[xy]=[21]

Hacer un pequeño cambio en el vector del lado derecho de las ecuaciones.

[1223][xy]=[4.0017.001]  

da

[xy]=[1.9991.001]

Hacer un pequeño cambio en la matriz de coeficientes de las ecuaciones.

[1.0012.0012.0013.001][xy]=[47]  

da

[xy]=[2.0030.997]

Este sistema de ecuación “se ve” bien condicionado porque pequeños cambios en la matriz de coeficientes o el lado derecho dieron como resultado pequeños cambios en el vector de solución.

Entonces, ¿y si el sistema de ecuaciones está mal condicionado o bien condicionado?

Bueno, si un sistema de ecuaciones está mal condicionado, no podemos confiar tanto en la solución. Vuelva a examinar el problema de la velocidad, Ejemplo 5.1 en el Capítulo 5. Los valores en la matriz de coeficientes[A] son cuadrados de tiempo, etc. Por ejemplo, si en lugar dea11=25, lo usastea11=24.99, querrías que este pequeño cambio marcara una gran diferencia en el vector de solución. Si lo hiciera, ¿confiaría en la solución?

Posteriormente veremos cuánto (términos cuantificables) podemos confiar en la solución en un sistema de ecuaciones. Cada matriz cuadrada invertible tiene un número de condición y junto con la máquina epsilon, podemos cuantificar cuántos dígitos significativos se puede confiar en la solución.

Para calcular el número de condición de una matriz cuadrada invertible, necesito saber qué significa la norma de una matriz. ¿Cómo se define la norma de una matriz?

Al igual que el determinante, la norma de una matriz es un simple número escalar único. Sin embargo, la norma es siempre positiva y se define para todas las matrices, cuadradas o rectangulares, y matrices cuadradas invertibles o no invertibles.

Una de las definiciones populares de una norma es la norma de suma de filas (también llamada norma de matriz uniforme). Para unam×n matriz[A], la norma de suma de filas de[A] se define como

A=max1imnj=1|aij|

es decir, encontrar la suma del valor absoluto de los elementos de cada fila de la matriz[A]. El máximo dem dichos valores es la norma de suma de filas de la matriz[A].

Ejemplo 3

Encuentra la norma de suma de filas de la siguiente matriz [A].

A=[107032.0996515]

Solución

A=max1i33j=1|aij|=max[(|10|+|7|+|0|),(|3|+|2.099|+|6|),(|5|+|1|+|5|)]=max[(10+7+0),(3+2.099+6),(5+1+5)]=max[17, 11.099, 11]=17

¿Cómo se relaciona la norma con el condicionamiento de la matriz?

Empecemos a responder a esta pregunta usando un ejemplo. Volver al sistema mal condicionado de ecuaciones,

[1223.999][xy]=[47.999]

que da la solución como

[xy]=[21]

Denotando el conjunto anterior de ecuaciones como

[A][X]=[C]

X=2

C=7.999

Haciendo un pequeño cambio en el lado derecho,

[1223.999][xy]=[4.0017.998]

da

[xy]=[3.9994.000]

Denotando el conjunto anterior de ecuaciones por

[A][X]=[C]

vector lado derecho es encontrado por

[ΔC]=[C][C]

y el cambio en el vector de solución se encuentra por

[ΔX]=[X][X]

entonces

[ΔC]=[4.0017.998][47.999]=[0.0010.001]

y

[ΔX]=[3.9994.000][21]=[5.9993.000]

entonces

ΔC=0.001

ΔX=5.999

El cambio relativo en la norma del vector de solución es

ΔXX=5.9992

 =2.9995

El cambio relativo en la norma del vector del lado derecho es

ΔCC=0.0017.999

 =1.250×104

Ver el pequeño cambio relativo de1.250×104 en la norma del vector del lado derecho da como resultado un gran cambio relativo en la norma del vector de solución de 2.9995.

De hecho, la relación entre el cambio relativo en la norma del vector de solución y el cambio relativo en la norma del vector del lado derecho es

ΔX/XΔC/C=2.99951.250×104=23993

Volvamos ahora al sistema bien acondicionado de ecuaciones.

[1223] [xy]=[47]

da

[xy]=[21]

Denotando el sistema de ecuaciones por

[A][X]=[C]

X=2

C=7

Hacer un pequeño cambio en el vector del lado derecho

[1223][xy]=[4.0017.001]

da

[xy]=[1.9991.001]

Denotando el conjunto anterior de ecuaciones por

[A][X]=[C]

el cambio en el vector del lado derecho se encuentra entonces por

[ΔC]=[C][C]

y el cambio en el vector de solución es

[ΔX]=[X][X]

entonces

[ΔC]=[4.0017.001][47]=[0.0010.001]

y

[ΔX]=[1.9991.001][21]=[0.0010.001]

entonces

ΔC=0.001

ΔX=0.001

El cambio relativo en la norma del vector de solución es

ΔXX=0.0012=5×104

El cambio relativo en la norma del vector del lado derecho es

ΔCC=0.0017=1.429×104

Ver el pequeño cambio relativo en la norma del vector del lado derecho de1.429×104 los resultados en el pequeño cambio relativo en la norma del vector de solución de5×104.

De hecho, la relación entre el cambio relativo en la norma del vector de solución y el cambio relativo en la norma del vector del lado derecho es

ΔX/XΔC/C=5×1041.429×104=3.5

¿Cuáles son algunas de las propiedades de las normas?

  1. Para una matriz[A],A0
  2. Para una matriz[A] y un escalar k,kA=|k|A
  3. Para dos matrices[A] y[B] del mismo orden,A+BA+B
  4. Para dos matrices[A] y[B] eso se puede multiplicar como[A] [B],ABAB

¿Existe una relación general entreΔX/X yΔC/C o entreΔX/X yΔA/A? Si es así, podría ayudarnos a identificar sistemas de ecuaciones bien acondicionados y mal condicionados.

Si existe tal relación, ¿nos ayudará a cuantificar el condicionamiento de la matriz? Es decir, ¿nos dirá en cuántos dígitos significativos podríamos confiar en la solución de un sistema de ecuaciones lineales simultáneas?

Existe una relación que existe entre

ΔXXandΔCC

y entre

ΔXXandΔAA

Estas relaciones son

ΔXXA  A1ΔCC

y

ΔXX+ΔXA  A1ΔAA

Las dos desigualdades anteriores muestran que el cambio relativo en la norma del vector del lado derecho o la matriz de coeficientes puede amplificarse tanto comoAA1.

Este númeroAA1 se llama el número de condición de la matriz y junto con la máquina épsilon, podemos cuantificar la precisión de la solución de[A] [X]=[C].

Demostrar para

[A] [X]=[C]

que

ΔXX+ΔX  A  A1ΔAA

Prueba

Let

[A]  [X]=[C](1)

Entonces si[A] se cambia a[A], el[X] cambiará a[X], tal

que

[A]  [X]=[C](2)

De las ecuaciones (1) y (2),

[A]  [X]=[A][X]

Denotando cambio en[A] y[X] matrices como[ΔA] y[ΔX], respectivamente

[ΔA]=[A][A]

[ΔX]=[X][X]

entonces

[A][X]=  ([A]+[ΔA]) ([X]+[ΔX])

Ampliando la expresión anterior

[A][X]= [A][X]+[A][ΔX]+[ΔA][X]+[ΔA][ΔX]

[0]=[A][ΔX]+[ΔA]([X]+[ΔX])

[A][ΔX]=[ΔA]([X]+[ΔX])

[ΔX]=[A]1[ΔA]([X]+[ΔX])

Aplicando el teorema de las normas, que la norma de matrices multiplicadas es menor que la multiplicación de las normas individuales de las matrices,

ΔXA1ΔAX+ΔX

Multiplicando ambos lados porA

A ΔXAA1ΔAX+ΔX

ΔXX+ΔXA A1 ΔAA

¿Cómo utilizo los teoremas anteriores para encontrar cuántos dígitos significativos son correctos en mi vector de solución?

El error relativo en una norma de vector de solución es Cond (A) error× relativo en la norma vectorial derecha.

El posible error relativo en la norma del vector de solución esCond(A)×mach

De ahí que nosCond(A)×mach debieran dar el número de dígitos significativos, m que son al menos correctos en nuestra solución al averiguar el mayor valor de m para el cualCond(A)×mach es menor que0.5×10m.

Ejemplo 4

¿Cuántos dígitos significativos puedo confiar en la solución del siguiente sistema de ecuaciones?

[1223.999] [xy]=[24]

Solución

Para

[A]=[1223.999] 

se puede mostrar

[A]1=[3999200020001000] 

A=5.999

A1=5999

Cond(A)=AA1=5.999×5999.4=35990

Asumiendo una sola precisión con 23 brocas utilizadas en la mantisa para números reales, la máquina épsilon es

mach=223=0.119209×106

Cond(A)×mach=35990×0.119209×106=0.4290×102

Para qué valor positivo máximo de m,Cond(A)×mach sería menor o igual a0.5×10m

0.4290×1020.5×10m

0.8580×10210m

log(0.8580×102)log(10m)

2.067m

m2.067

m2

Entonces, dos dígitos significativos son al menos correctos en el vector de solución.

Ejemplo 5

¿Cuántos dígitos significativos puedo confiar en la solución del siguiente sistema de ecuaciones?

[1223] [xy]=[47]

Solución

Para

[A]=[1223] 

se puede mostrar

[A]1=[3221] 

Entonces

A=5,

A1=5.

Cond (A)=A  A1=5×5=25

Asumiendo una sola precisión con 23 brocas utilizadas en la mantisa para números reales, la máquina épsilon

mach=223=0.119209×106

Cond(A)×mach=25×0.119209×106=0.2980×105

Para qué valor positivo máximo de m,Cond(A)×mach sería menor o igual a0.5×10m

0.2980×1050.5×10m

m5

Entonces, cinco dígitos significativos son al menos correctos en el vector de solución.

Cuestionario de adecuación de soluciones

Quiz 1

La norma de suma de filas de la matriz

[A]=[673131921232941475161]es

(A)29

(B)61

(C)98

(D)200

Quiz 2

La adecuación de la solución de ecuaciones lineales simultáneas[A][X]=[C] depende de

(A) el número de condición de la matriz de coeficientes[A]

(B) la máquina epsilon

(C) el número de condición para la matriz[A] y el epsilon de la máquina

(D) norma de la matriz de coeficientes[A]

Quiz 3

Dado un conjunto de ecuaciones en forma de matriz[A][X]=[C],A=250,A1=40andεmach=0.119×106,, entonces el número de dígitos significativos en los que al menos puede confiar en la solución son

(A)1

(B)2

(C)3

(D)4

Quiz 4

La solución a un conjunto de ecuaciones lineales simultáneas

[a11a12a13a21a22a23a31a32a33][x1x2x3]=[4494138]se da como

[x1x2x3]=[247]La solución a otro conjunto de ecuaciones lineales simultáneas viene dada por (tenga en cuenta que la matriz de coeficientes es la misma que la anterior)

[a11a12a13a21a22a23a31a32a33][x1x2x3]=[43.9993.98138.03]

se da como

[x1x2x3]=[214.01208.0160]

Según la norma de suma de filas, el número de condición de la matriz de coeficientes es mayor que (elija el mayor valor posible)

(A)1

(B)138

(C)4500

(D)139320

Quiz 5

El número de condición de la matriz den×n identidad basado en la norma de suma de filas es

(A)0

(B)1

(C)n

(D)n2

Quiz 6

Vamos[A]=[12+δ2δ1]. Basado en la norma de suma de filas y dado queδ0,δ>0, el número de condición de la matriz es

(A)3δ3+δ

(B)9δ23δ2

(C)(3+δ)23δ2

(D)32δδ23δ2

Ejercicio de Adecuación de Soluciones

Ejercicio 1

La adecuación de la solución de ecuaciones lineales simultáneas depende de

  1. Número de condición
  2. Máquina epsilon
  3. Producto de número de condición y máquina epsilon
  4. Norma de la matriz.
Contestar

C

Ejercicio 2

Si un sistema de ecuaciones[A] [X]=[C] está mal condicionado, entonces

  1. det(A)=0
  2. Cond(A)=1
  3. Cond(A)es grande.
  4. Aes grande.
Contestar

C

Ejercicio 3

SiCond(A)=104 ymach = 0.119× 10 -6, entonces en[A] [X]=[C], al menos estos muchos dígitos significativos son correctos en su solución,

  1. 0
  2. 1
  3. 2
  4. 3
Contestar

C

Ejercicio 4

Hacer un pequeño cambio en la matriz de coeficientes de[1223.999][xy]=[47.999]
y encontrar
ΔX/XΔA/A

¿Es un número grande o pequeño? ¿Cómo se relaciona este número con el número de condición de la matriz de coeficientes?

Contestar

Cambiando[A] a

[1.0012.0012.0014.000]

Resultados en solución de

[59992999]

ΔXXΔAA=5999.720.0025.999

4=8.994×106 $

Ejercicio 5

Hacer un pequeño cambio en la matriz de coeficientes de
[1223][xy]=[4 7]
y encontrar
ΔX/XΔA/A. ¿Es un número grande o pequeño? Compara tus resultados con el problema anterior. ¿Cómo se relaciona este número con el número de condición de la matriz de coeficientes?

Contestar

Cambiando[A] a

[1.0012.0012.0013.000]

Resultados en solución de

[2.0030.9970]

ΔXXΔAA=0.0032.0000.0025

=3.75

Ejercicio 6

Demostrar

ΔXX  A A1 ΔCC

Contestar

Utilice el teorema de que si[A][B]=[C] entoncesABC

Ejercicio 7

Para

[A]=[107032.0996515]

da

[A]1=[0.10990.23330.27990.29990.33320.39990.049950.16666.664×105]

  1. ¿Cuál es el número de condición de[A]?
  2. ¿Cuántos dígitos significativos podemos al menos confiar en la solución de[A] [X]=[C] ifmach=0.1192×106?
  3. Sin calcular la inversa de la matriz[A], ¿se puede estimar el número de condición de[A] usar el teorema en el problema #6?
Contestar

A.A=17
A1=1.033
Cond(A)=17.56

B. 5

C. Pruebe diferentes valores del lado derecho deC=[±1  ±1  ±1]T con signos elegidos al azar. LuegoA1X se obtiene de resolver el conjunto de ecuaciones[A] [X]=[C] comoC=1.

Ejercicio 8

Demostrar que elCond(A)1.

Contestar

Sabemos que

A BAB

entonces si

[B]=[A]1,

A A1AA1

IAA1

1AA1

AA11

Cond(A)1


This page titled 9: Adecuación de soluciones is shared under a CC BY-NC-SA 4.0 license and was authored, remixed, and/or curated by Autar Kaw via source content that was edited to the style and standards of the LibreTexts platform.

Support Center

How can we help?