Informe de Curso
Informe de Curso
Informe de Curso
FACULTAD DE INGENIERIA
ESCUELA PROFESIONAL DE
INGENIERIA AMBIENTAL
TEMA: ANALISIS NUMERICO
DOCENTE: ING.HAYDEE SISA
YATACO
CURSO: METODOS NUMERICOS
ESTUDIANTE:
GIAN CARLOS CAIRO GONZALEZ 2018062427
CLAUDIA ADELINA CCAPERA MOROCHARA 2020066190
DARIO ULISES SILVA CRUZ 2016056359
ARIAN HUAMAN COAQUIRA 2018060835
TACNA-PERÚ
2020
INDICE
INTRODUCION
ANALISIS NUMERICO
Aplicación
El análisis numérico lo podemos aplicar en distintas áreas y procedimientos
matemáticos, tales como:
- Cálculo de derivadas
- Integrales
- Ecuaciones diferenciales
- Operaciones con matrices
- Interpolaciones
- Ajuste de curvas
- Polinomios
- Los métodos numéricos se aplican en áreas como: Ingeniería Industrial,
Ingeniería Química, Ingeniería Civil, Ingeniería Mecánica, Ingeniería eléctrica,
etc...
Aproximación
Una aproximación es una representación inexacta que, es suficientemente fiel como
para ser útil. Esta aproximación nunca es utilizada en ciencias exactas a grado
profesional debido a la pérdida de información. Aunque en matemáticas la
aproximación típicamente se aplica a números, también puede aplicarse a objetos tales
como las funciones matemáticas, figuras geométricas o leyes físicas.
Ejemplo
3.1416 es una aproximación numérica de π
2.7183 es una aproximación numérica de e
1.4142 es una aproximación numérica de 2
0.333333 es una aproximación numérica de 1/3
Error
Cada vez que realizamos una medición tenemos una gran probabilidad de cometer
algún tipo de error, el cual nos arrojara un resultado más o menos alejado del que
realmente deberíamos obtener. Y es que medir, es más bien un proceso aproximado que
exacto. De entre los errores más comunes podemos distinguir dos grandes grupos:
Errores sistemáticos. Son errores relacionados con la forma, en la que su utiliza
el instrumento de medida.
a. Error de calibrado. Se trata de uno de los errores más frecuentes y está
ligado directamente al instrumento. Muchos de ellos deben ser calibrados
antes de ser usados.
b. Error de paralaje. Es propio de instrumentos de medida analógicos como
por ejemplo aquellos que poseen agujas para marcar los valores. Dos
observadores situados en posiciones oblicuas a la aguja pueden leer valores
diferentes.
Errores aleatorios o accidentales. Se tratan de errores que se producen debido
a causas que no se pueden controlar. Para intentar reducir el efecto de este tipo
de errores se suele medir varias veces en las mismas condiciones y se considera
como valor final más probable la media aritmética de los datos obtenidos.
Dado que todas las medidas están afectadas por un error experimental, en el mundo
científico es común hacer constar cada resultado obtenido en una medición junto con
la incertidumbre sobre esa medida.
Incertidumbre
La incertidumbre es un valor numérico que se obtiene por medio de dos
nuevos conceptos denominados error absoluto y error relativo.
Error absoluto
El error absoluto de una medida es el valor absoluto de la diferencia entre el valor
verdadero y el valor experimental de la magnitud que ha sido objeto de la medición.
Error relativo
El error relativo de una medida es el cociente entre el error absoluto y el valor
verdadero de la magnitud que ha sido objeto de medición.
Error absoluto
Error relativo=
Valor verdadero
Fuentes de error
Error de redondeo
Surgen debido a la limitación del espacio, ya que los números reales presentan una
infinidad de dígitos, pero en la maquina solo se pueden representar una cantidad finita.
Error de truncamiento
Surge debido a la limitación en tiempo, ya que se deben realizar un numero finito de
accione
Ejemplo
- Evaluar funciones con la serie de Taylor
- Proceso iterativo convergente
- Evaluar por intervalos
EJEMPLOS
1. 3,5 m como longitud de un terreno que mide realmente 3,59 m, 60 m como la
distancia entre dos postes que están situados a 59,91 m.
1. 2: En la medida de 1 m se ha cometido un error de 1 mm, y en 300 Km, 300
m. ¿Qué error relativo es mayor?
1. Ejercicio propuesto
Al medir la distancia entre dos pueblos, sabemos que el valor real es de 5,478 Km.
a) Aproxima la distancia hasta las décimas tanto por redondeo y por truncamiento.
b) Calcula el error absoluto cometido en ambos casos.
c) Calcula el error relativo cometido en ambos casos.
d) Justifica, de acuerdo con los resultados anteriores, qué aproximación es mejor
Estructura selectiva
Se caracteriza por la presencia de condiciones, que se evalúan mediante la sentencia
if, else, elseif, en sentencia else puede ser omitido sino que puede presentarse en caso,
en observaciones la condiciones la condición if una condición matricial del tipo A==B
son del mismo tamaño, para que cumpla la condición necesario de dos elementos de los
matrices A Y B.
Ejemplo 1
Una función cuyo nombre sea f, definida por
Ejemplo 2
Estructura repetiva o iterativa
Se caracteriza por la presencia de iteraciones o repetición de un conjunto de
instrucciones. Básicamente se implementan mediante dos bucles: for, while.
Bucle de for
Este bucle usa una variable llamada contador, la primera vez la variable contador
Toma un valor igual a inicio y ejecuta bloque A incrementa en un valor igual a paso (si
se omite, se le asigna automáticamente el valor 1). Cuando contador llega al valor final,
el ciclo repetitivo se acaba y el programa continúa con las sentencias que siguen a end.
Ejemplo
Una función cuyo nombre sea fac, que reciba un entero no negativo n y calcule su
factorial n!. Probar con n = 5 y n = 0.
Bucle while
En ese caso es apropiado el empleo del bucle while. El bucle while, se inicia y
continúa mientras la condición se cumple (sea verdadera) y finaliza cuando la condición
no se cumple (sea falsa). Se acostumbra colocar antes del ciclo la inialización de la
variable de la condición, para que esta se cumpla al menos la primera vez.
Ejemplo
Una función numsumandos (tope), que retorne el máximo valor n de manera que la
suma 1+2+3+4+: : :+n no supere el valor tope.
METODOS ITERATIVOS
METODOS DIRECTOS
Son procedimientos para obtener resultados realizando una secuencia finita de
operaciones aritméticas. La cantidad de cálculos aritméticos depende del tamaño del
problema. El resultado obtenido será exacto siempre que se puedan conservar en forma
exacta los valores calculados en las operaciones aritméticas, caso contrario se
introducirán los errores de redondeo.
Modelo matemático
a1,1x1 = b1
a2,1x1 + a2,2x2 = b2
a3,1x1 + a3,2x2 + a3,3x3 = b3
....
an,1x1 + an,2x2 + ax,3x3 + . . . + an,nxn = bn
METODO DE LA BISECCION
Recordemos que, si f es continua sobre [a;b] y f (a) f (b) < 0, entonces hay al menos
una raíz real de f entre a y b. El método de bisección, es un tipo de búsqueda en el que
se inicia con un intervalo que contiene a la raíz, luego este intervalo se divide por la
mitad y se escoge el subintervalo que contiene a la raíz, con el subintervalo escogido se
efectúa el mismo proceso, así sucesivamente. El proceso se repite hasta obtener una
mejor aproximación, considerando como aproximación a la raíz el punto medio de cada
intervalo escogido.
METODO DE NEWTON
Si f es una función definida en un intervalo [a;b] ; tal que f 2 C2 [a;b] (es decir, que
la derivada de f es continua en [a;b]) y si x0 2 [a;b] es el punto que inicia el proceso de
convergencia; entonces, la fórmula del método de Newton .
METODO DE MATRICIAL
Resolver una ecuación matricial de la forma:
AX +B =C;
Asumiendo compatibilidad en las dimensiones de las matrices. En este caso la
Solución es:
X = A=1 (C=B):
Siempre y cuando A sea una matriz inversible. Lo importante aquí es que la inversa
multiplica a por la izquierda a la matriz C=B, de lo contrario, si fuera por la derecha, el
resultado sería completamente diferente, debido a la no conmutatividad de la
multiplicación de matrices.
ARREGLOS Y MATRICES
Arreglos
Corresponde a un conjunto de números, que están ordenados en filas y columnas, los
elementos de las filas y columnas, los elementos que se encuentran en las filas, están
separados por comas y los que se encuentran dentro de la columna, están separados por
un punto u coma.
Ejemplo
Arreglos unidimensionales
Conjunto de elementos que solo están ordenados en una columna o una fila
Length
Este código calcula el numero de elementos encontrados en un arreglo
Ejemplo
Arreglos bidimensionales
Como su nombre lo indica, bidimensional se refiere a un conjunto de números, los
cuales se encuentran ordenados en mas de una fila y columna
Ejemplo
Operadores aritméticos
Adición
Es la suma de los elementos de un arreglo
Ejemplo
Sustracción
Es la resta de los elementos de un arreglo
Ejemplo
Multiplicación
Es el producto de los elementos de un arreglo
Ejemplo
División
Es el cociente de cada elemento de los arreglos
Ejemplo
Matriz
Corresponde a un arreglo rectandular de números reales o complejos ordenados en
filas o columnas
Ejemplo
Operadores aritméticos
Adición
Es la suma de los elementos de las matrices según su posición
Ejemplo
Multiplicación
Se realiza el producto de cada elemento de la matriz por el numero dado
Ejemplo
Matriz simétrica
Una matriz es simétrica si la traspuesta de la matriz es la misma matriz
Ejemplo
Matriz identidad
Es una matriz cuadrada donde sus elementos de la diagonal son uno y los demás son
0
Ejemplo
Matriz Inversa
Una matriz A de orden n (n filas y n columnas) tiene inversa cuando su rango es n, es
decir, cuando el rango de dicha matriz coincide con su orden, o también, cuando su
determinante sea distinto de cero.
Ejemplo
Determinante de una matriz
El determinante de una matriz cuadrada A es un número real
Ejemplo
eye
Genera una matriz identidad de orden mxm
Ejemplo
Zeros
Genera una matriz nula de orden mxn
Ejemplo
Ones
Genera una matriz de orden mxn, donde cada elemento es 1
Magic
Genera una matriz cuya suma de elementos de sus filas, columnas y diagonales da
como resultado un mismo valor
Ejemplo
Funciones en matrices
Diag
Devuelve los elementos de la diagonal de una matriz cuadrada
Ejemplo
Sum
Devuelve la suma obtenida en cada una de las columnas de la matriz
Ejemplo
Rand
Genera una matriz de orden nxn en forma aleatoria, donde sus elementos son
uniformemente distribuidos en el intervalo 0,1
Ejemplo
Rref
Devuelve de forma canonica la matriz
Ejemplo
Find
Devuelve los índices de los elementos de la matriz que no son cero, enmerandolos de
forma vertical
Ejemplo
Norm
Calcula la norma de un vector o matriz
Ejemplo
Ejemplo 2
Cross
Calcula el producto vectorial entre 2 vectores
Ejemplo
Eig
Calcula los valores y vectores propios de la matriz
Ejemplo
Length
Calcula la dimensión del vector
Ejemplo
Size
Calcula la dimensión de la matriz
Ejemplo
Vertcat
Concatena matrices verticalmente
Ejemplo
Horzcat
Concatena matrices horizontalmente
Ejemplo
Repmat
Crea una matriz de copias de una matriz existente
Ejemplo
Blkdiag
Crea una matriz diagonal de bloques, en el orden que se ubican las matrices, los
demnas elementos de la matriz recién creada son ceros
Ejemplo
METODOLOGIA
4. Un lápiz mide 8,0cm. El profesor dejó como tarea a Juan que lo midiera y le
reportara los errores de medición. Al medir el lápiz, Juan obtuvo que la longitud era
8,1cm.
Lapiz Juan
cm 8 8.1
5. Determinar el volumen del agua que tomará. Midió 3 veces y obtuvo 10,0mL;
10,2mL y 10,1mL. ¿Cuál fue el error absoluto y relativo?
datos (ml) Error absoluto (ml) Error relativo
10.0 0.1
10.2 0.1
10.1 0.0
promedio 10.1 0.1 1%
6. Determinar la masa de harina que usará para la torta. Midió 3 veces y obtuvo 1,2kg;
1,5kg y 1,3kg. ¿Cuál fue el error absoluto y relativo de la medida 1,5kg?
Datos (kg) Error abs (kg) Error relativo
1.2
1.5 0.2 13 %
1.3
Promedio 1.3
Ejercicios de estructuras
1. Leer dos números y calcular su suma, resta y multiplicación.
2. Ingresar un número e indicar si e número es par y múltiplo de 5
>> A=[1,4;5,2];
>> A=[1,4; 5,2];
>> B=[2,1; 3,4];
>> C=A+B
C=
3 5
8 6
>> A=[1,4;5,2];
>> B=A+2
B=
3 6
7 4
5. Sustracción de arreglos
>> A=[1,4;5,2];
>> B=[2,1;3,4];
>> D=A-B
D=
-1 3
2 -2
6. Multiplicación de arreglos
>> A=[1,4;5,2;3,7];
>> B=[2,1;3,4;2,6];
>> E=A.*B
E=
2 4
15 8
6 42
>> D=[3-27;654;549];
>> E=-4.*D
E=
96
-2616
-2196
>> E=-4*D
E=
96
-2616
-2196
8. División de arreglos
>> A=[1,4;5,2];
>> B=[2,1;3,4];
>> F=A/B
F=
-1.6000 1.4000
2.8000 -0.2000
>> A=[7 14 8 3; 6 19 7 2; 3 13 4 1]
A=
7 14 8 3
6 19 7 2
3 13 4 1
>> B=[1 3 4 7; 3 4 5 6; 2 1 13 5]
B=
1 3 4 7
3 4 5 6
2 1 13 5
>> C=A+B
C=
8 17 12 10
9 23 12 8
5 14 17 6
>> C=A-B
C=
6 11 4 -4
3 15 2 -4
1 12 -9 -4
>> C=A.*B
C=
7 42 32 21
18 76 35 12
6 13 52 5
>> C=A./B
C=
7.0000 4.6667 2.0000 0.4286
2.0000 4.7500 1.4000 0.3333
1.5000 13.0000 0.3077 0.2000
>> A=[7 14 8 3; 6 19 7 2; 3 13 4 1]
A=
7 14 8 3
6 19 7 2
3 13 4 1
>> A'
ans =
7 6 3
14 19 13
8 7 4
3 2 1
4. HALLAR LA INVERSA DE
LA SIGUIENTE MATRIZ.
>> B=[1 3 4 7; 3 4 5 6; 2 1 13
5; 34 6 8 10]
B=
1 3 4 7
3 4 5 6
2 1 13 5
34 6 8 10
>> y=inv(B)
y=
-0.0117 -0.0411 -0.0012 0.0335
-0.4507 0.6725 -0.0951 -0.0405
-0.1197 0.0810 0.0880 -0.0088
0.4061 -0 .3286 -0.0094 0.0176
>> det(A)
ans =
0
6. GENERE UNA MATRIZ NULA DE (3,6).
Código:
Did you mean:
>> B=zeros(3,6)
B=
0 0 0 0 0 0
0 0 0 0 0 0
0 0 0 0 0 0
7. GENERE UNA MATRIZ DE ORDEN (6*6) EN FORMA ALEATORIA.
Código:
>> X=rand(6)
X=
0.8147 0.2785 0.9572 0.7922 0.6787 0.7060
0.9058 0.5469 0.4854 0.9595 0.7577 0.0318
0.1270 0.9575 0.8003 0.6557 0.7431 0.2769
0.9134 0.9649 0.1419 0.0357 0.3922 0.0462
0.6324 0.1576 0.4218 0.8491 0.6555 0.0971
0.0975 0.9706 0.9157 0.9340 0.1712 0.8235
8. CALCULAR LA NORMA DE LA SIGUIENTE MATRIZ.
Código:
>> A=[50 5 27 400 7; 0 67 90 6 97; 30 14 23 251 490]
A= 50 5 27 40 7
50 5 27 400 7 0
0 67 90 6 97
0 67 90 6 97 30 14 23 25 490
30 14 23 251 490 1
>> norm(A,1)
ans =
657
>> norm(A,2)
ans =
606.8897
>> norm(A,inf)
ans =
808
Código
>> u=[0 1 2 3 5 15]
u=
0 1 2 3 5 15
>> length(u)
ans =
Código
>> C=[1 3;5 4;7 9]
C=
1 3
5 4
7 9
>> D=[3 5 7;4 1 16;9 61 13]
D=
3 5 7
4 1 16
9 61 13
>> E=[C,D]
E=
1 3 3 5 7
5 4 4 1 16
7 9 9 61 13
12. CREE UNA MATRIZ DIAGONAL DE BLOQUES CONSIDERANDO EL
PRIMER BLOQUE (2*2)
15.
x=[3,4,90,- 4;5,5,34,0;6,-
53,7,9;4,5,6,5]
Sum(x(:,2))
17. GENERE UNA MATRIZ CUYA SUMA DE ELEMENTOS DE SUS FILAS,
COLUMNAS Y DIAGONALES DA COMO RESULTALDO UN MISMO
VALOR, DICHA MATRIZ ES DE ORDEN(6)
X=magic(6)
c=a+b
19. GENERE UNA MATRIZ SOLO CON EL VALOR 1 DE (8.9). GENERE UNA
MATRIZ NULA DE (9*8) Y OBTENGA SU MULTIPLICACIÓN.
A= ones(8,9)
B=zeros(9,8)
C=A*B (No aplico. *, debido que son matrices diferentes)
X=[3,4,90,-4,-2;5,5,34,0,4;6,-53,7,9,7;4,5,6,5,11]
Surf(X)
EJERCICIOS PARA APLICACION
s=
p=
2.9580
Observar que si ejecutamos
>> calculo(vector_a,vector_
b) ans =
Ejercicios EA y ER
b) Ea=300 m= 0,3 Km, por lo que el error relativo en el segundo caso es:
Er=0,3/300=0,001, es decir, el 0,1%.
Vemos que, en ambos casos, el error relativo es el mismo por lo que las dos
mediciones, aunque no son iguales, tienen comparativamente la misma precisión.
Ejercicio 3. Como medida de un radio de 7 dm hemos obtenido 70.7 cm. Calcula
el error absoluto y el relativo.
Solución:
Antes de operar, tenemos que expresar el valor del radio en cm, para
poder realizar las operaciones: 7dm=70cm. De aquí tenemos:
a) Ea = |70 – 70,7| = 0,7 cm.
b) Er= 0,7/70 = 0,01
Ejercicio 4. Al medir la distancia entre dos pueblos, sabemos que el valor real es
de 5,478 Km.
a) Aproxima la distancia hasta las décimas tanto por redondeo y por
truncamiento.
b) Calcula el error absoluto cometido en ambos casos.
c) Calcula el error relativo cometido en ambos casos.
d) Justifica, de acuerdo con los resultados anteriores, qué aproximación es
mejor.
Solución:
a) La distancia, aproximada por redondeo hasta las décimas, es de
5,5Km y, aproximada por truncamiento, es de 5,4Km.
b) El error absoluto en el primer caso viene dado por: Ea = |5,478 – 5,5| =
0,022
Km y, en el segundo caso: Ea = |5,478-5,4| = 0,078 cm.
c) El error relativo en el primer caso es de Er=0,022/5,478= 0,00402 y,
en el segundo caso de Er= 0,014.
d) La mejor aproximación es la que hacemos por redondeo, ya que el
error absoluto es más pequeño. De ahí también que el error relativo
también sea más pequeño.
B
.
'
a
n
s
=
1 + 2i 2 -i
3 + 4i 4 -i
(
1
:
4
)'
a
n
s
=
1
2
3
4
v1 = 1:4;
v2 = [v1, 5; 0.1:0.1:0.5]
v3 = [v2', [11,12,13,14,15]']
0 1
12345 2 0:2 12
1 0:1 11
v1 = v2 = 0:1 0:2 0:3 0:4 0:5 v3 = B 4 0:4 14 C
B C
1 2 3 4 B 3 0:3 13 C
B
4 0:5 15
@ A
Estructuras repetitivas while y for
E j e m p l o 1 ( S u m a Im p a r e s ) ( U s o d e w h ile )
E s c r ib ir u n a M - f u n c ió n q u e c a lc u le y d e v u e lv a e l m a y o r n ú m e r o im p a r n p a r a e l c u a l la s u m a
d e t o d o s lo s n ú m e r o s im p a r e s e n t r e 1 y n e s m e n o r o ig u a l q u e 5 0 0 0 .
fu n c tio n [n ] = S u m a Im p a r e s
%
% n = SumaImpares es e l mayor numero impar t a l que l a suma de todos
% los impares desde 1 hasta n es < = 5000
%
n = 1;
suma = 1; inicio
suma = suma + n;
false
end suma <= 5000
n = n - 2;
true
n = n + 2
end
suma = suma + n
n = n - 2
Salida: n
fin
E je m p lo 2 (U s o d e fo r)
E s c r ib e u n a M - f u n c ió n
fu n c tio n [s u m a ] = S u m a v ( v )
fu n c tio n [s u m a ] = S u m a v ( v )
%
% Sumav(v) es l a suma de las componentes del vector v
%
suma = 0;
fo r k = 1 : le n g th ( v ) Entrada: v
suma = suma + v ( k ) ;
end suma = 0
k = 1
false
k <= longitud(v)
+ v(k)
k = k + 1
Salida: v
O b s e r v a c i ó n : L o q u e e s i m p o r t a n t e o b s e r v a r y c o m p r e n d e r a q u í e s q u e e n e s t e ti p o d e b u c l e , n o
h a y q u e o c u p a r s e e x p líc it a m e n t e d e la in ic ia liz a c ió n d e la v a r ia b le - ín d ic e k n i d e s u in c r e m e n t o e n
c a d a it e r a c ió n . L a in s t r u c c ió n f o r s e o c u p a d e t o d o c u a n d o s e p r o p o r c io n a la li s t a d e v a lo r e s q u e
d e b e to m a r k .
Programación de Matlab
1. ) Escribir un script que calcule el coste de enviar un paquete de correos, en base a la siguiente
tabla de precios:
T ip o d e e n v ío P e s o (0 -1 k g ) P e s o (1 -5 k g ) P e s o (5 -2 5 k g )
T ie r r a 1.50 e 1 .5 0 e + 1 e a d i- 5 .5 0 e + 0 .6 0 e
c io n a l p o r c a d a k g o adicionales por cada kg
fra c c ió n d e k g , a p a r - o fracción de kg, a
ti r d e 1 k g d e p e s o . partir de 5 kg de peso.
A ir e 3 e 3 e + 1 e a d ic io - n a l 1 0 .2 0 e + 1 .2 0 e
por c a d a k g o fra c - a d ic io n a le s p o r c a d a
c ió n d e k g , a p a r ti r d e k g o fr a c c ió n d e k g , a
1 kg de peso. p a r ti r d e 5 k g d e p e s o .
Urgente 18 e 1 8 e + 6 e a d i- No se realizan envíos
c io n a le s p o r c a d a k g o urgentes de paquetes
f r a c c i ó n d e k g , a p a r ti r con un peso superior a
d e 1 k g d e p e so . 5 kg.
E l p r o g r a m a p e d ir á p o r p a n t a lla lo s d a t o s p e s o y t i p o d e e n v í o . E n e l c a s o d e u n e n v ío p o r
ti e r r a o a i r e d e m á s d e 2 5 k g , e l p r o g r a m a e s c r i b i r á e n p a n t a l l a e l m e n s a j e “ N o s e r e a l i z a n
re p a r to s d e m á s d e 2 5 k g ” . E n e l c a s o d e e n v ío s u r g e n t e s d e m á s d e 5 k g , e l p r o g ra m a
e s c r ib ir á e n p a n t a lla e l m e n s a je “ N o s e r e a liz a n r e p a r t o s u r g e n t e s d e m á s d e 5 k g ”
Arreglo y matrices
Crea un arreglo o array multidimensional que contenga 3 columnas y las filas que tu quieras,
las dos primeras columnas tendrán números y en la 3 columna sera el resultado de sumar el
número de la primera y segunda columna. Muestra el array o arreglo de la siguiente forma:
3+5=8
4+4=7
Construya un algoritmo que al recibir los montos de ventas mensuales de cinco departamentos
de una fábrica, proporciona la siguiente información:
Realizar un subproceso para obtener el monto del departamento que tuvo la mayor
venta en el mes de julio.
BIBLIOGRAFIA
https://www.fisicalab.com/apartado/errores-absoluto-relativos
https://aula.upt.edu.pe/pluginfile.php/171822/mod_resource/content/0/Matrice
s-y-Arreglos.pdf
https://es.slideshare.net/tysonLlanos/programacion-matlab-37386312