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

Publicado el 29 de Octubre del 2018 por Administrador (703 códigos)
17.267 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(2)

Publicado el 29 de Octubre del 2018gráfica de visualizaciones de la versión: Versión 1
17.268 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 (2)

Diego
13 de Noviembre del 2019
estrellaestrellaestrellaestrellaestrella
De que es la libreria que importas? No entiendo esa parte
Responder
luisa
26 de Febrero del 2020
estrellaestrellaestrellaestrellaestrella
muy bien hecho
Responder

Comentar la versión: Versión 1

Nombre
Correo (no se visualiza en la web)
Valoración
Comentarios...
CerrarCerrar
CerrarCerrar
Cerrar

Tienes que ser un usuario registrado para poder insertar imágenes, archivos y/o videos.

Puedes registrarte o validarte desde aquí.

Codigo
Negrita
Subrayado
Tachado
Cursiva
Insertar enlace
Imagen externa
Emoticon
Tabular
Centrar
Titulo
Linea
Disminuir
Aumentar
Vista preliminar
sonreir
dientes
lengua
guiño
enfadado
confundido
llorar
avergonzado
sorprendido
triste
sol
estrella
jarra
camara
taza de cafe
email
beso
bombilla
amor
mal
bien
Es necesario revisar y aceptar las políticas de privacidad

http://lwp-l.com/s4902