ClubEnsayos.com - Ensayos de Calidad, Tareas y Monografias
Buscar

Práctica 3 Método de Gauss-Seidel


Enviado por   •  21 de Mayo de 2019  •  Apuntes  •  2.185 Palabras (9 Páginas)  •  256 Visitas

Página 1 de 9

Práctica 3

Método de Gauss-Seidel

Objetivo: Programar y probar el Método de Gauss-Seidel para la Solución de Sistemas de Ecuaciones Lineales.

Metodología: 

1.- Implementar en el lenguaje de su elección el pseudocódigo del Método de Gauss-Seidel con Relajación del siguiente listado:

  1. n=No. de renglones de la matriz a
  2. x(n)=0
  3. for i=1 hasta n
  4.    temp=a(i,i)
  5.    Si temp==0 entonces
  6.       error=1
  7.       Sale de la rutina
  8.    for j=1 hasta n
  9.       a(i,j)=a(i,j)/temp
  10.    Termina for j
  11.    b(i)=b(i)/temp
  12. Termina for i
  13. for i=1 hasta n
  14.    suma=b(i)
  15.    for j=1 hasta n
  16.       Si i ≠ j entonces
  17.         suma=suma-a(i,j)*x(j)
  18.    Termina for j
  19.    x(i)=suma
  20. Termina for i
  21. for iteracion=2 hasta imax
  22.    marcador=0
  23.    for i=1 hasta n
  24.       anterior=x(i)
  25.       suma=b(i)
  26.       for j=1 hasta n
  27.          Si i ≠ j entonces
  28.             suma=suma-a(i,j)*x(j)
  29.       Termina for j
  30.       x(i)=λ*suma+(1- λ)*anterior
  31.       Si marcador=0 y x(i) ≠0 entonces
  32.          ea=|x(i)-anterior/x(i)|*100
  33.       Si ea>es entonces marcador=1  
  34.    Termina for i
  35.    Si marcador ==0 entonces
  36.       Sale de la rutina
  37. Termina for iteracion

Parámetros de Entrada:

a – Matriz de Coeficientes

b – Vector de Términos Independientes

es – Error permisible

imax – No. máximo de iteraciones

λ – Coeficiente de relajación

Parámetros de Salida:

x – Vector Solución

error – Código de error

             0 – No hay error

             1 – Hay un cero en la diagonal principal

2.- Probar su programa con el siguiente ejemplo usando como datos de entrada es=0.01, imax=100, λ=1 y mostrar la pantalla de salida con los resultados parciales y finales.

[pic 1]

[pic 2]

3.- Anotar sus observaciones

Al realizar esta práctica, no presente problemas a la hora de escribir el código principal del método. El único inconveniente fue en la compresión a primer contacto con el método, ya que era algo distinto a los ya realizados anteriormente. Adicionalmente a eso, no hubo alguna otra problemática que impidiera la ejecución correcta del mismo.

4.- Anotar sus conclusiones.

Los conocimientos adquiridos durante la realización de la práctica fueron los suficientes para lograr la comprensión total del funcionamiento del método, y además lograr en futuras practicas o proyectos la correcta aplicación y ejecución de las técnicas adquiridas. Gracias a esto creo que mis conocimientos en programación han crecido de forma bastante lineal, y espero en un futuro lograr realizar grandes proyectos

5.- Agregar el código del método.

/*######################################################

****                        Archivo: gauss_seidel.c                        ****

****         Este archivo contiene la funcion dek metodo                ****

****         de Gauss Seidel a del proyecto de la                        ****

****        Unidad 2 para la solucion de Sistemas de                 ****

****        Ecuaciones Lineales.                                        ****

****                                                                 ****

****        Elaboro: Luis Esteban Maritnez Salcido                        ****

****        Materia: Metodos Numericos                                ****

****        Fecha: 05/04/2019                                        ****

########################################################

*/

/*Incluye archivos de encabezados*/

#include "unidad3.h"

//Funcion del Metodo Gaus-Jordan

int gauss_seidel(float **a, float *b, int n, float es, int imax, float lambda)

{

        /*Declaracion de variables locales*/

        int i,j,iteracion,flag;

        float *x,suma,temp,ant,ea=100;

        //Aparta Memoria para el Vector x

...

Descargar como (para miembros actualizados) txt (6 Kb) pdf (200 Kb) docx (43 Kb)
Leer 8 páginas más »
Disponible sólo en Clubensayos.com