CÁLCULO NUMÉRICO Y ESTADÍSTICA
APLICADA
Cálculo Numérico y
Estadística Aplicada
EJERCICIOS VARIOS
1
CAPÍTULO VARIOS
Ejercicio 1
Utilizando el polinomio característico y la resolución numérica de la ecuación lineal
resultante, determinar los valores de la matriz real y simétrica A.
Nota: convergencia 8 decimales.
( )
4 0 1
𝐴= 0 1 1
1 1 2
Solución:
2
3 2
− 𝑥 +7 𝑥 +12 𝑥+ 3
𝑥𝑛 +1 =𝑥 𝑛 − 2
− 3 𝑥 +14 𝑥 +12
3
Calculando
4
3 2
𝑃 ( 𝑥 )=− 𝑥 + 7 𝑥 +12 𝑥+3
5
Ejercicio 2
Una matriz de energías correspondientes a un problema triatómico es proporcional a la
matriz
Determinar los valores de la matriz real y simétrica H, dando pruebas de los resultados
obtenidos.
Solución:
3 2
𝑃 ( 𝑥 )= − 𝑥 + 7 𝑥 − 12 𝑥 +5=0
6
7
Calculando
8
Ejercicio 3
En un problema de correlación se desea encontrar un valor de ρ que satisfaga
la ecuación:
Solución:
Determinamos ρ numéricamente:
Se trata de una función monótona -1<ρ<1
Análisis del cambio de signo
- - - - +
-0,5 -0,25 -0 +0,25 +0,5
9
El valor de ρ estará comprendido entre 0,25 y 0,5.
Utilizando el método de Newton-Raphson
Si sustituimos este valor en , tenemos:
10
11
12
Ejercicio 4
En la teoría de calores específicos de sólidos (Debye) a bajas temperaturas aparece la
integral
Utilizando las reglas del trapecio y de Simpson con 5 puntos, estimar el valor de I para
ZD=1,2 y comparar ambos resultados (redondeados por exceso a 6 decimales).
Integración numérica
13
Trapecio:
ZD=1,2 al ser 5 puntos y el primero es y0=0, el espaciado es h=1,2/4=0,3
x 0 0,3 0,6 0,9 1,2
y 0 0,077174 0,262736 0,499451 0,744790
La integral por trapecio viene dada por:
Entonces, tenemos que la aproximación por el método del trapecio sería:
14
Método de Simpson
15
Ejercicio 5
En la teoría de Debye para sólidos a bajas temperaturas aparece la
integral
Para el aluminio a temperaturas suficientemente elevadas
Evaluar la integral anterior para la temperatura T=190 K mediante
integración numérica utilizando.
a) Regla del trapecio (5 puntos).
b) Regla se Simpson.
16
Solución:
𝜃𝐷 /𝑇 2 3 3 2
𝑥 𝑥 3𝑥
𝐼= ∫ 𝑦 ( 𝑥 ) 𝑑𝑥=∫ 𝑥 𝑑𝑧=𝑙𝑖𝑚 𝑥 =𝑙𝑖𝑚 𝑥 =¿0;𝑛𝑜h𝑎𝑦𝑖𝑛𝑑𝑒𝑡𝑒𝑟𝑚𝑖𝑛𝑎𝑐𝑖ó𝑛¿
0 0 𝑒 −1 𝑥→0 𝑒 −1 𝑥→0 𝑒
17
a) Integral por trapecios
b) Integral por Simpson
18
Ejercicio 6
La energía libre de Helmholtz de exceso para un modelo de fluido atómico clásico, en
condiciones de densidad y temperaturas dadas, puede calcularse integrando la
ecuación de estado y es directamente proporcional a la integral.
Donde representa una densidad adimensional.
Calcular la integral I,
a) Analíticamente.
b) Mediante la regla del trapecio, utilizando 9 puntos para . Comparar con el
resultado exacto.
19
Solución:
a)
b) La integración de trapecios con 9 puntos, requiere una partición del intervalo (0,1/4) en
8 subintervalos. Esto implica un espaciado h=0,25/8.
20
21
Ejercicio 7
En muchas cuestiones de observaciones ópticas aparece la integral de Fresnel
Mediante el método de los coeficientes indeterminados obtener la siguiente aproximación
numérica
Estimar el valor de C(A=2) y discutir el resultado.
22
Solución:
Método de coeficientes indeterminados (páginas: 164, 1651 180 y 182)
Tenemos que crear un algoritmo exacto para , pues hay cuatro coeficientes que
debemos determinar (cuatro ecuaciones).
Primer caso:
23
Segundo caso
Tercer caso
24
Cuarto caso
25
Del sistema resultante se obtienen los cuatro coeficientes
Con ellos evaluamos la función pedida para A=2
26
Ejercicio 8
Calcular por mínimos cuadrados los parámetros a y b de la relación
27
Solución:
reescribimos la ecuación
28
29
Ejercicio 9
Calcular mediante el método del polinomio característico los valores y vectores propios
de la matriz siguiente. Discutir los resultados.
Solución:
Determinamos los autovalores
Las raíces son:
30
Obtenemos los vectores propios asociados a los autovalores anteriores.
31
32
Ejercicio 10
Diagonalizar mediante el método de Jacobi, obteniendo autovalores y autovectores la
matriz simétrica.
Para diagonalizar la matriz debemos “eliminar” los términos a13=a31=-3
La rotación viene definida por:
33
La transformación de semejanza es:
34
Los valores propios son x1=-1, x2=1 y x3=5
Los autovectores asociados son las correspondientes columnas de la matriz O
35
Ejercicio 11
Diagonalizar por el método de Jacobi la matriz 3x3 con elementos:
a11=a22=a33=a13=a31=1, y a12=a21=a23=a32=0.
Solución:
Construimos la matriz con los términos indicados, y así:
Para diagonalizar la matriz debemos “eliminar” los términos a13=a31=1
36
La rotación viene definida por:
La transformación de semejanza es:
37
Con un solo paso de rotación la matriz se diagonaliza. Los valores propios son:
Los valores propios son x1=0, x2=1 y x3=2
Los autovectores asociados son las correspondientes columnas de la matriz O
38
Ejercicio 12
Diagonalizar por el método de Jacobi, determinando valores y vectores propios la
matriz A
Solución:
Para diagonalizar la matriz debemos “eliminar” los términos a13=a31=-3
La rotación viene definida por:
39
La transformación de semejanza es:
40
Los valores propios son 1=0, 2=1 y 3=6
Los autovectores asociados son las correspondientes columnas de la matriz O
41
CÁLCULO NUMÉRICO Y ESTADÍSTICA
APLICADA
Cálculo Numérico y
Estadística Aplicada
EJERCICIOS VARIOS
42
