Código de Python - Algoritmo de eliminación de Gauss-Jordan

Imágen de perfil

Algoritmo de eliminación de Gauss-Jordangráfica de visualizaciones


Python

estrellaestrellaestrellaestrellaestrella(1)
Publicado el 29 de Octubre del 2018 por Administrador
8.726 visualizaciones desde el 29 de Octubre del 2018
Es un algoritmo del álgebra lineal que se usa para determinar las soluciones de un sistema de ecuaciones lineales, para encontrar matrices e inversas.
https://es.wikipedia.org/wiki/Eliminación_de_Gauss-Jordan

jordan-gauss

Versión 1
estrellaestrellaestrellaestrellaestrella(1)

Publicado el 29 de Octubre del 2018gráfica de visualizaciones de la versión: Versión 1
8.728 visualizaciones desde el 29 de Octubre del 2018
estrellaestrellaestrellaestrellaestrella
estrellaestrellaestrellaestrella
estrellaestrellaestrella
estrellaestrella
estrella


Forma parte de Python Fiddle
 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
#      --- Importando librerias ---      #
import fpformat
 
#      --- Se definen variables ---      #
dimx = 0
dimy = 0
mini = 0
 
#  --- Se definen funciones/metodos ---  #
 
# Revision del dato ingresado (flotante)
def error_test2(x):
    try:
        int(x)
        return False
    except:
        return True
 
# Revision del dato ingresado (entero)
def error_test3(x):
    try:
        float(x)
        return False
    except:
        return True
 
# Revision general
def error_test(y, z=1):
    if z == 1:
        while error_test2(y):
            print "\n\n\t\tERROR!!!\n\t\tSe han ingresado valores o caracteres no validos."
            y = raw_input("\n\nPor favor, reingrese el valor deseado: ")
        return int(y)
    elif z == 2:
        while error_test3(y):
            print "\n\n\t\tERROR!!!\n\t\tSe han ingresado valores o caracteres no validos."
            y = raw_input("\n\nPor favor, reingrese el valor deseado: ")
        return float(y)
 
# Generado de la matriz
def generador(y,x):
    mat = []
    for k in range(y):
        mat.append([])
    for i in range(y):
        for j in range(x):
            print
            mat[i].append(error_test(raw_input("a_" + str(i+1) + str(j+1) + " = "), 2))
    return mat
 
# Despliega matriz
def impri(m):
    print "\n"
    print "\t||" + ("\t" * (len(m[0]) + 1)) + "||"
    for i in range(len(m)):
        print "\t||\t",
        for j in range(len(m[0])):
            print fpformat.fix(m[i][j],2),"\t",
        print "||"
        print "\t||" + ("\t" * (len(m[0]) + 1)) + "||"
    print "\n"
 
# Escoger valor minimo
def peque(a,b):
    if a > b:
        return b
    elif a < b:
        return a
    else:
        return a
 
# Sumador de listas
def adic(l1, l2):
    ln = [0]*len(l1)
    for i in range(len(l1)):
        ln[i] = l1[i] + l2[i]
    return ln
 
# Multiplicador por constante
def multi(l1, k):
    ln = [0]*len(l1)
    for i in range(len(l1)):
        ln[i] = k * l1[i]
    return ln
 
# Buscador y cambiador de filas
def swap_finder(m, l, x):
    if l == (x-1):
        print "La matriz es una matriz singular."
        print "Eso implica que no se puede resolver.\n\n"
        return m, False
    else:
        for i in range(l, x):
            if m[i][l] != 0:
                m[l], m[i] = m[i], m[l]
                return m, True
        print "La matriz es una matriz singular."
        print "Eso implica que no se puede resolver.\n\n"
        return m, False
 
# Sistema de resolucion Gauss-Jordan
def solver(m,p):
    for i in range(p):
        boo = True
        if m[i][i] == 0:
            m, boo = swap_finder(m, i, p)
        if boo:
            m[i] = multi(m[i], 1/float(m[i][i]))
        else:
            return m
        for j in range(i+1,len(m)):
            m[j] = adic(m[j], multi(m[i], -1*float(m[j][i])))
    for k in range(p-1, -1, -1):
        for l in range(k-1, -1, -1):
            m[l] = adic(m[l], multi(m[k], -1*float(m[l][k])))
    return m
 
#       --- Programa Principal ---       #
 
print "\n"*5
 
dimy = error_test(raw_input("Ingrese la cantidad de filas de la matriz: "), 1)
dimx = error_test(raw_input("Ingrese la cantidad de columnas de la matriz: "), 1)
 
matriz = generador(dimy, dimx)
 
impri(matriz)
 
mini = peque(dimx, dimy)
 
solucion = solver(matriz, mini)
 
print "\t\t*** Solucion ***"
 
impri(matriz)



Comentarios sobre la versión: Versión 1 (1)

Diego
13 de Noviembre del 2019
estrellaestrellaestrellaestrellaestrella
De que es la libreria que importas? No entiendo esa parte
Responder

Comentar la versión: Versión 1

Nombre
Correo (no se visualiza en la web)
Valoración
Comentarios
Es necesario revisar y aceptar las políticas de privacidad

http://lwp-l.com/s4902