Practica 10 - Método de Gauss

Instituto Politcnico Nacional
Escuela Superior de Ingeniera Mecnica y Elctrica Unidad Zacatenco
Curso
Analisis Numrico
Practica #4: Solucin de sistemas de ecuaciones lineales simultaneas por el Mtodo de Gauss Jordan
Profesor
Miguel Jimnez Guzmn
Mxico D.F., Agosto de 2013
Practica #1:
Solucin de sistemas de ecuaciones lineales simultaneas por el Mtodo de Gauss-Jordan

1. Objetivo: Que el alumno comprenda y aplique el mtodo de Gauss-Jordan para la solucin de ecuaciones lineales simultneas, usando programacin en MatLab y/o C. 2. Introduccin: Un sistema de ecuaciones lineales de la forma:
donde , son los coeficientes constantes, son los trmins independientes constantes y son las incgnitas. Se dice que es un sistema que tiene n ecuaciones con n incgnitas o simplemente que es de orden nxn. Mtodo de Eliminacin Gaussiana Este mtodo se aplica para resolver sistemas lineales de la forma:
El mtodo de eliminacin Gaussiana (simple), consiste en escalonar la matriz aumentada del sistema: [ ]
para obtener un sistema equivalente de la forma:
donde la notacin se usa simplemente para denotar que el elemento cambi. Se despejan las incgnitas comenzando con la ltima ecuacin y hacia arriba. Por esta razn, muchas veces se dice que el mtodo de eliminacin Gaussiana consiste en la eliminacin hacia adelante y la sustitucin hacia atrs. Eliminacin Hacia Adelante: En esta fase reduce el conjunto de ecuaciones a un sistema triangular Superior: 1) Consiste en dividir la primera ecuacin por el coeficiente de la primera incgnita (coeficiente pivote). A este procedimiento se le conoce como normalizacin. 2) Despus se multiplica la ecuacin normalizada por el primer coeficiente de la segunda ecuacin. 3) Ntese que el primer trmino de la primera ecuacin es idntico al primer trmino de la segunda. Por lo tanto, se puede eliminar, la primera incgnita de la segunda ecuacin restando la primera a la segunda. 4) Se repite el paso 2 y 3 hasta eliminar la primera incgnita de todas las ecuaciones restantes. 5) Estos 4 pasos se repiten tomando como pivotes las ecuaciones restantes hasta convertir el sistema en la matriz triangular superior siguiente.
donde , coeficientes modificados i = reglon, j= columna, k = pivote Sustitucin Hacia Atras:
Una ves que se obtiene el sistema equivalente que es un sistema triangular superior, este se resuerlve despejando primero la variable Xn, este valor se utiliza para obtener despejando la incgnita Xn-1 hasta obtener el resultado completo del sistema. ( )
Ejemplo de Aplicacin 1: Solucin de Circuitos Elctricos de Corriente Directa mediantea Analisis Nodal. La ley de corrientes de Kirchhoff establece que para cualquier circuito elctrico, la suma algebraica de todas las corrientes en cualquier nodo en el circuito es igual a cero. En el anlisis nodal, si hay n nodos en un circuito, y se selecciona un nodo de referencia, los otros nodos pueden ser numerados desde V1 a Vn-1. Con un nodo seleccionado como el nodo de referencia, habr n-1 ecuaciones independientes. Si se supone que la admitancia entre los nodos i y j esta da como Yij, se pueden escribir las ecuaciones nodales:
(4.1) donde m=n-1 V1, V2 y Vm son los voltajes de los nodos 1, 2 y asi sucesivamente hasta..., n con respecto al nodo de referencia. Ix es la suma algebraica de las fuentes de corriente en el nodo x. La ecuacin (4.1) puede ser expresada en forma matricial como [Y][V] = [I] La solucin de la ecuacin anterior es [V] = [Y] -1 [I] donde [Y] -1 es la inversa de [Y]. (4.3) (4.2)
Ejemplo 1 For the circuit shown below, find the nodal voltages V V 1 2 , and V3 .
Figure 4.1 Circuit with Nodal Voltages
Solution Using KCL and assuming that the currents leaving a node are positive, we have Para el nodo 1,
(4.7) At node 2,
(4.8) At node 3,
(4.9) In matrix form, we have [ ] [ ] (4.10)
The MATLAB program for solving the nodal voltages is

5
MATLAB Script diary ex4_1.dat % program computes the nodal voltages given the admittance matrix Y and current vector I % Y is the admittance matrix and I is the current vector % initialize matrix y and vector I using YV=I form Y = [ 0.15 -0.1 -0.05; -0.1 0.145 -0.025; -0.05 -0.025 0.075]; I = [5; 0; 2]; % solve for the voltage fprintf('Nodal voltages V1, V2 and V3 are \n') v = inv(Y)*I diary The results obtained from MATLAB are Nodal voltages V1, V2 and V3, v= 404.2857 350.0000 412.8571 1. Desarrollo de la Prctica: a) Realizar el Diagrama de Flujo para el Mtodo de Gauss a partir del siguiente pseudocodigo siguiente:
// Eliminacin Hacia Adelante DO FOR k = 1, n 1 DO FOR i = k + 1, n factor = a(i,k) / a(k,k) DO FOR j = k + 1 to n a(i,j) = a(i,j) factor a(k,j) END DO b(i) = b(i) factor b(k) END DO END DO // Eliminacin Hacia Atras X(n) = b(n) / a(n,n) DO FOR i = n 1, 1, 1
sum = b(i) DO FOR j = i + 1, n sum = sum a(i,j) x(j) END DO X(i) = sum / a(i,i) END DO
b) A partir del diagrama de flujo escribir el cdigo del programa en MatLab y/o C. c) Ejecutar el programa y resolver el circuito elctrico del ejemplo 1, encontrando las corrientes en cada elemento resistivo: d) Escribir los resultados del proceso de eliminacin hacia adelante y del proceso de sustitucin hacia atrs para cada uno de los siguientes ejemplos:
e) El mtodo de eliminacin Gaussiana (simple) puede presentar un problema cuando uno de los elementos que se usan para hacer ceros, es cero o muy cercano acero. Este problema se puede resolver fcilmente intercambiando los renglones. A este procedimiento se le llama pivoteo parcial y aplicado a la eliminacin Gaussiana, nos d el llamado mtodo de eliminacin Gaussiana con pivoteo (parcial). Para elegir el elemento pivote en la columna se escoge el elemento con valor absoluto.
Usar eliminacin Gaussiana con pivoteo para resolver el siguiente sistema:
f) Escriba sus Conclusiones g) Guarde sus programas en una memoria USB h) Cierre el MatLab y/o el compilador de C utilizado y apague la computadora siguiendo el procedimiento correcto.

Practica 10 - Método de Gauss

Cargado por

Copyright:

Formatos disponibles

Practica 10 - Método de Gauss

Cargado por

Información del documento

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Practica 10 - Método de Gauss

Cargado por

Copyright:

Formatos disponibles

Instituto Politcnico Nacional

Escuela Superior de Ingeniera Mecnica y Elctrica Unidad Zacatenco

Miguel Jimnez Guzmn

Mxico D.F., Agosto de 2013

Solucin de sistemas de ecuaciones lineales simultaneas por el Mtodo de Gauss-Jordan

para obtener un sistema equivalente de la forma:

donde , coeficientes modificados i = reglon, j= columna, k = pivote Sustitucin Hacia Atras:

Figure 4.1 Circuit with Nodal Voltages

(4.9) In matrix form, we have [ ] [ ] (4.10)

The MATLAB program for solving the nodal voltages is

Usar eliminacin Gaussiana con pivoteo para resolver el siguiente sistema:

También podría gustarte