Material Detail

Consideramos el sistema de ecuaciones lineales, M x = d, con M=[dd, -1, 0;-1, dd, -1;0, -1, dd]; d = [2,8,-6]'; donde los elementos de la diagonal principal de M, dd, es un dato a introducir, siendo dd>=2 para que no haya problemas de convergencia. El sistema anterior se resolverá usando los algoritmos de Gauss-Seidel y SOR. Se puede elegir el parámetro w del algoritmo SOR, siendo 1 menor o igual a w menor o igual a 2. Los resultados que se presentan son: Se representan (entre a y b) el módulo de los valores propios de inv(D+wL)*(wU+(w-1)D) en función de w, donde M=D+L+U. El w óptimo es el w para el cual el máximo del módulo de los valores propios es mínimo. Se proporciona el w óptimo y el número de iteraciones para Gauss-Seidel y SOR con el w elegido. El objetivo de este laboratorio es, fijada una matriz, observar la diferencia entre el número de iteraciones usadas para resolver el sistema con cierta tolerancia prefijada con el método Gauss-Seidel y con SOR. Cuando Gauss-Seidel conver

Disciplines:

• Mathematics and Statistics  / Mathematics  / Applied Mathematics