• Author / Uploaded
• Francisco Castillo

Citation preview

Reporte

Nombre: Francisco Javier Castillo Montes Matrícula: 2764254 Nombre del curso: Nombre del profesor: Métodos Numéricos Módulo:

Santiago Banda Muñoz Actividad:

Módulo 2. Solución de sistemas de

Actividad 9

ecuaciones Fecha: 3 de octubre de 2016

Bibliografía TecMilenio. (3 de octubre de 2016). Blackboard. Obtenido de https://miscursos.tecmilenio.mx/bbcswebdav/institution/UTM/semestre/profesio nal/ti/ti13201/bb/tema8.htm

Objetivo: Encontrar la solución a un sistema de ecuaciones no lineales. Procedimiento: 1. Leer el tema 2. Realizar ejemplos 3. Realizar reporte Resultados: El método como se vio anteriormente, se puede hacer a mano si se tienen pocos puntos, pero si hay 100 puntos o más es prácticamente imposible solucionar esto a menos que se utilice una computadora. Si se tienen 20 puntos que describen la cantidad de voltaje generada por una celda solar cada hora:

Reporte

{ static void Main(string[] args) { double pivote, factor; double[,] jacobiana = new double[20, 3]; double[,] matriz = new double[3, 4];

-19.5771,

//Temperaturas de cada mes double[] voltaje = { -2.4684, 56.8038, 4.8393, 13.2878, -0.1931, 55.4146, 74.4829, 62.0632, 37.3616, 14.1876, 72.9088, 81.5205, 120.7811, 156.3607,

141.7361, 15, 16,

110.898, 129.7614, 192.8516, 253.4778 }; double[] tiempo = { 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 17, 18, 19, 20}; for (int i = 0; i < { jacobiana[i, 0] jacobiana[i, 1] jacobiana[i, 2] }

20; i = i + 1) = tiempo[i] * tiempo[i]; = Math.Sin(tiempo[i]); = Math.Exp(tiempo[i] / 10);

//Multiplicación de matrices //-----------------------------for (int i = 0; i < 3; i = i + 1) for (int j = 0; j < 3; j = j + 1) for (int k = 0; k < 20; k = k + 1) matriz[i, j] = matriz[i, j] + jacobiana[k, i] * jacobiana[k, j];

Reporte for (int i = 0; i < 3; i = i + 1) for (int j = 0; j < 1; j = j + 1) for (int k = 0; k < 20; k = k + 1) matriz[i, 3] = matriz[i, 3] - voltaje[k] * jacobiana[k,

i];

//-----------------------------//Eliminación Gaussiana //-----------------------------for (int reng = 0; reng < 3; reng = reng + 1) { pivote = matriz[reng, reng]; for (int colu = 0; colu < 4; colu = colu + 1) matriz[reng, colu] = matriz[reng, colu] / pivote; for (int reng_elimi = 0; reng_elimi < 3; reng_elimi = reng_elimi + 1)

if (reng_elimi != reng) { factor = matriz[reng_elimi, reng]; for (int colu_elimi = 0; colu_elimi < 4; colu_elimi = colu_elimi + 1) matriz[reng_elimi, colu_elimi] = matriz[reng_elimi, colu_elimi] - factor * matriz[reng, colu_elimi]; }

} //------------------------------

//Imprime los valores de las variables Console.WriteLine(matriz[0, 3] + " " + matriz[1, 3] + " " + matriz[2, 3] + " "); Console.ReadLine(); } }

}

Conclusión:

Reporte

En este ejercicio vimos la utilización del método para resolver mínimos cuadrados y por el cual se podrán resolver más métodos iguales o similares a este tipo. Podemos concluir que en esta serie de datos y otra más podremos resolver o saber sus valores dependiendo del método a utilizar.