Métodos Numéricos: Eliminación Gaussiana (Principiante)/Gaussian Elimination (Beginner)
ฝัง
- เผยแพร่เมื่อ 21 ก.ย. 2024
- En este vídeo programamos en Python el método de Eliminación Gaussiana, para resolver sistemas de ecuaciones lineales en Canopy.
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
Copia y pega el código desde:
bit.ly/1OQiPRC
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
Software para programar en Python
Canopy Express:
www.enthought....
SageMath:
cloud.sagemath...
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
Vídeo para dar formato de salida de tablas:
• Python: Formato de tab...
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
Nuestro blog:
cctmexico.blogs...
increible el video me ayudo mucho para un proyecto que tengo en unos dias
+Endry Castillo ¡Muchas gracias por tomarte el tiempo para dejarnos tu comentario tan amable ^_^! ¡Éxito en tu proyecto!
LO MEJOR DE LA VIDA, MIL GRACIAS
¡Nos da mucho gusto saber que el vídeo te fue de utilidad!
Gracias. El vídeo fue muy útil.
¡Muchas gracias a tí, por vernos!
Muy bueno, muchas gracias
¡Qué bueno que te gustó! ¡Muchas gracias a tí, por vernos!
Buenos días! Su código funciona muy bien y he hecho los cambios que dijo y funciona fenomenal. Sin embargo, estoy intentando poner las operaciones en funciones, con un main incluido, y me está costando mucho. Tiene alguna idea de como hacerlo? Si me pudiera ayudar con el código con las funciones le agradecería mucho. Muchas gracias!
+MChiribogaD ¡Qué bien que funcionó todo! Justamente para las funciones tenemos un vídeo que te puede ayudar a darle la estructura a tu programa: ➡th-cam.com/video/WzTzGWXtj0s/w-d-xo.html
Muchas gracias! Pude resolver el problema!
+MChiribogaD 😀 ¡Excelente! ¡Nos da mucho gusto saber que pudiste implementar las funciones en tu código! 👍💯💯
Me arroja un error:
#Esta es la matriz ingresada:
[[ 1. 2. 3.]
[ 1. 2. 3.]
[ 5. 4. 3.]]
#Este es el error que me arroja
C:/Users/Christopher/Desktop/Unipoli/Algebra lineal/Gauss-Jordan.py:19: RuntimeWarning: divide by zero encountered in double_scalars
factor=(matriz[r,k]/matriz[k,k])
-inf
Resultado matriz
[[ 1. 2. 3.]
[ 0. 0. 3.]
[ 0. 4. 3.]]
Resultado del vector
[ 27. 5. -inf]
Resultados:
[ 0. 0. -inf]
A que se refiere con el "-inf"?
Hola! el sistema de ecuaciones no se puede resolver, ese inf es infinito, puedes revisar tu resultado en WolframAlpha, para que no te quedes con dudas es muy sencillo, aquí te mostramos como th-cam.com/video/QVJfqhkC7Bg/w-d-xo.html
Muchisimas gracias por su respuesta tan rapida!
canal 100% recomendado! gracias!
Muchas gracias a tí, por dejarnos tu amable comentario ^_^ Saludos!!
Muy buen video , ya cheque varias veces el codigo pero solo tengo un error en el primer print , me dice missing parentheses in call to 'print'
Hola, en la nueva versión de Python se requieren los paréntesis para imprimir líneas, por ejemplo: print (Hola)
Hola, en clase nos piden usar trinket para realizar nuestros códigos de Python.
Siguiente las indicaciones del video, cuando quiero correr el programa me arroja un error que dice:
**ValueError: Tuple must contain values for all dimensions on line 15 in main.py**
no se que tenga que hacer.
Le agradecería su mas pronta respuesta.
Puede ser que sea un problema con las dimensiones de tus arreglos (o listas?) Aquí nos encontramos una pregunta similar en stackoverflow, si gustas puedes pegar tu codigo en repl.it para ver si encontramos el error, o puedes hacer tu pregunta directamente en Stackoverflow, también tienen una sección en español, stackoverflow.com/questions/28520366/numpy-valueerror-size-of-tuple-must-match-number-of-fields
Hola.... en la matriz resultados, la diagonal principal no esta en 1's y reemplazando los valores 10 para x, -3 para y; y 5 para z no cumple con las ecuaciones, ¿por favor en donde estoy haciendo mal?
Hola Migiuel, ¿Estás usando los mismos datos que nosotros? Si es así, revisa que no haya algún error cuando pasaste el código.
si quiero que la matriz sea 4x4 siempre debo fijar m y n igual con 4 al principio?
Si! Es una de las opciones más fáciles! Aunque también puedes pedirla al usuario.
Hola buenas noches. Cuando trato de correr el script en la terminal de linux de ubuntu 17.10 me arroja el siguiente error: AttributeError: 'module' object has no attribute 'open'; a que creen que se deba? La versiòn de python que uso es 2.7
No hemos tenido ese error, pero aquí dejamos dos sugerencias:
1) Probar el código en colaboratory, para descartar algún problema del código: colab.research.google.com
2) Hacer la pregunta stackoverflow.com responden muy rápido y está la opción de preguntar en Inglés o Español, seguro ayudarán a resolver ese problema.
una pregunta, que viene siendo el vector solucion y que numero debor poner ahi?
+Luis Pichardo el vector solución, son los resultados de las ecuaciones, por ejemplo para 3x+2y+5z=10 el 10 sería parte del vector solución 😊
Hola. Como puedo modificar el programa para que me de el resultado de la matriz en su forma escalonada reducida? Es decir que encuentre todo los ceros en la parte inferior y superior y que todo los pivotes sean "1"
Hola +edgar! Hay que modificar los ciclos de las operaciones en renglones y columnas, aquí te dejamos el código:
for k in range (0,m):
for r in range(0,m):
if (r != k):
factor=(matrix[r,k]/matrix[k,k])
vector[r]=vector[r]-(factor*vector[k])
for c in range(0,n):
matrix[r,c]=matrix[r,c]-(factor*matrix[k,c])
Y para imprimir los resultados:
print 'Resultados: '
for r in range(0,m):
x[r]=(vector[r]/matrix[r,r])
print x[r]
+cctmexico ¡Muchas gracias!Saludos ;)
De nada edgar ^_^ Saludos!!
@@cctmexico Hola tengo una pregunta, como se define factor en las variables? me da error por no definirlo
Hola, lo que pasa es que quiero ejecutar lo que muestran en el video pero no lo ejecuto por canopy, si no que lo ejecuto por pyscripter, pero al ingresar lo que muestran acá, me ejecuta error, quisiera que me ejecutara en pyscripter, podrían ayudarme por favor?
+Zuleima Escalante Pérez Hola Zuleima, no hemos utilizado Pyscripter, si te sirve aquí te dejamos el en Sage (listo para ejecutar), no tienes que instalar algún software, incluso puedes compartir el proyecto y para crear una cuenta, se puede hacer desde la cuenta de Google o Facebook, el proyecto está disponible en: cloud.sagemath.com/projects/984c76e2-1d8b-49aa-86ec-407825b4535a/files/GaussJordanYT.sagews
+Zuleima Escalante Pérez También te recomendamos, la comunidad de Python en Google+, son muy buenos y rápidos para contestar dudas de Python, y seguro te encuentras a alguien que ya esté utilizando Pyscripter (planeamos en un futuro hacerlo nosotros también) la comunidad la encuentras en: plus.google.com/communities/103393744324769547228
+cctmexico Listo, Muchas gracias!!
Excelente Zuleima Escalante Pérez!! ^_^
hola buenas tardes, es que tengo una duda en el procedimiento de la parte cuando sacan los ceros por debajo de la diagonal principal, esa parte ahí de las formulas para los ceros no la entiendo , si me pueden ayudar
Hola Daniel, aquí tenemos un vídeo que te puede ser de utilidad: th-cam.com/video/0qICNqRJAig/w-d-xo.html
hola en la parte de matrix[(r),(c)]=(raw_input("Elemento a["+str(r+1)+","+str(c+1)+"] "))
me dice expected an indented block
como puedo corregirlo?
Hola Edith, revisa que no tengas espacios (indentaciones) cuando comenzaste a escribir las instrucciones.
hola implemente todo el programa en python, pero no se ejecuta y el archivo incluso se queda oculto al intentar buscarlo, tuve que descargar e importar la libreria numpy ya que python de por si no cuenta con esta libreria almenos que se le agregue, pero aun asi al intentar de ejecutar el programa en python no hace nada, alguna forma de solucionar el problema?
hice unos arreglos pero sigue siendo el mismo programa, observa a ver si puedes encontrar el problema:
import numpy
a = int(raw_input('Introdusca el valor de A:'))
b = int(raw_input('Introdusca el Valor de B:'))
matriz = numpy.zeros((a,b))
vector = numpy.zeros((b))
x = numpy.zeros((a))
print 'Introduce la matriz de coeficientes y el vector solución'
for r in range(0, a): #r = renglones
for c in range(0, n): #c = columnas
matriz[(r),(c)] = (raw_input("Elemento a[" +str(r + 1)+ ","+str(c + 1)+"] "))
vector[(r)]=(raw_input('b['+str(r + 1)+"]: '))
print(matriz)
for k in range (0,m)
for r in range(0,m): #Renglon pivote
factor=(matriz[r, k]/matriz[k, k])
vector[r]=vector[r]-(factor*vector[k])
for c in range(0,n):
matriz[r, c]=matriz[r, c]-(factor*matriz[k, c])
#sustitución hacia atrás
x[a - 1]=vector[a - 1]/matriz[a - 1, a - 1]
print x[a - 1]
for r in range(a - 2, -1, -1): #para que valla retrocediendo de 1 en 1 hacia atras
suma = 0
for c in range(0,b):
suma = suma + matriz[r, c]*x[c]
x[r]=(vector[r] - suma)/matriz[r, r]
print 'Resultado matriz'
print(matriz)
print 'Resultado del vector'
print(vector)
print 'Resultados: '
print(x)
+Endry Castillo A ver que te parece esta opción ^_^ te dejamos el programa aquí en Sage (listo para ejecutar), no tienes que instalar algún software, incluso puedes compartir el proyecto y para crear una cuenta, se puede hacer desde la cuenta de Google o Facebook, el proyecto está disponible en: cloud.sagemath.com/projects/984c76e2-1d8b-49aa-86ec-407825b4535a/files/GaussJordanYT.sagews
+cctmexico Nota: Para correrlo en SageMath, ocupas una cuenta (se hace en dos o tres clics) y después seleccionar la opción Run
+cctmexico no bien, al final era mi pc que no quería ejecutar python, termine usando canopy por que para mi presentación debia mostrarlo así, en fin le hice unos arreglos para que el usuario tuviera una interfaz para elegir si quería este método de reducción o el Gauss Jordán al igual de que también puede hacerlo para que los valores fueran aleatorios y que el usuario supiera lo que quiere, al final yo y mi grupo tuvimos la calificación mazima, gracias
¡Qué curioso! pero así pasa, con estas cosas, si vas a seguir programando en Python, prueba Sage, a ver que tal te funciona. Saludos!
+Endry Castillo ¡Muchas felicidades por su calificación! ^_^
Como haria para que el metodo funcione ingresando la matriz y el vector ya armados? lo intente pero me arrojo un resultado diferente :/
¿Nos podrías compartir tu código para checarlo? lo puedes pegar en repl.it/sessions o aquí mismo en los comentarios :D
Vale muchisimas gracias, ya lo he podido solucionar el problema era que al ingresar los arreglos con esos valores enteros, los resultados iban a ser unicamente enteros osea que se perdia parte del resultado por no aceptar floats, la solucion era especificar que los arreglos eran float simplemente jaja, ya te he dejado el codigo corregido
repl.it/IYDq
¡Excelente! Qué bien que ya quedó resuelto, y muchísimas gracias por dejar el código, ya verás que será muy útil para todos los que andamos en esto de los métodos en Python ¡Saludos desde México!
¡Muchas gracias! ^_^
en la linea 11 que es esta vector[(r)]=(input("b["+str(r+1)+"]: ")) me da esta advertencia " unindent does not match any outer indentation level " y no me deja correrlo.......Me gustaria que me corrijieran
Hola, ese error es porque hay un espacio (tabulación) dónde no se requiere, puede ser en una línea anterior o posterior, aquí te dejamos un video donde revisamos los errores más comunes:
th-cam.com/users/edit?o=U&video_id=X46F2-OxxqY
OK gracias.... que bien que si respondan las dudas lo mas pronto posible....mil gracias
Hola me podrian ayudar cuando mando ha hacer una matriz 4X3 phyton me manda este error
Traceback (most recent call last):
File "prueba.py", line 21, in
vector[r]=float(input('b['+str(r+1)+']: '))
IndexError: index 3 is out of bounds for axis 0 with size 3
Hola! Puedes hacer tu pregunta (e incluir el código, si así lo deseas) en stackoverflow.com es completamente gratis, responden muy rápido y tienes la opción de preguntar en Inglés o Español, seguro te ayudarán con ese problema. También aquí te dejamos un vídeo, que te puede servir como guía, para agregar tu pregunta: th-cam.com/video/X46F2-OxxqY/w-d-xo.html
como seria la sustitucion hacia adelante? saludos
En este vídeo utilizamos la sustitución hacia adelante, esperamos que te sea de utilidad: th-cam.com/video/FpVeXhAQg9w/w-d-xo.html
prefecto gracias
a mi no me funciono con la ultima version de python 3.9, tuve que modificar la sintaxis, creo que python a cambiado un poco y deberian hacer una actualizacion al codigo
Hola Edgar, puedes revisar los últimos vídeos donde manejamos las nuevas sintaxis de Python, por si te son de utilidad aquí te los dejamos:
Método Format:
th-cam.com/video/mjw9n0K1Ud4/w-d-xo.html
Método F-string:
th-cam.com/video/Pxb5O0QydTA/w-d-xo.html
Hola, sus videos me están siendo de muchísima utilidad pero ya para ejecutar el código tengo problemas y me arroja los siguientes códigos:
NameError: name 'vector' is not defined
TypeError: 'str' object does not support item assignment
Para el vector que no esta definido debajo de la linea donde esta "x=numpy.zeros((m))" definí vector=" " y se soluciono, corrí el código y sin problemas pero ya cuando b[1]:7 me arroja que el str no soporte el item asignado y de ahí ya no eh podido avanzar.
Tengo la ultima versión de Canopy corriendo en windows 10 con Python 3 instalado.
Gracias.
¡Muchas gracias por tan amable comentario! Puedes poner tu código (si así lo deseas) en stackoverflow.com y hacer tu pregunta en inglés o español, la comunidad responde muy rápido y te podrán ayudar con ese problemita, saludos!
Hola este programa sirve para m y n diferentes? Por qué no me jala me podrían ayudar?
PD si ingreso la matriz y los vectores y no tiene solución como enviara la solución?
Hola! Te mandaría un mensaje de error :)
Puedes hacer tu pregunta (e incluir el código, si así lo deseas) en stackoverflow.com es completamente gratis, responden muy rápido y tienes la opción de preguntar en Inglés o Español, seguro te ayudarán con ese problema. También aquí te dejamos un vídeo, que te puede servir como guía, para agregar tu pregunta: th-cam.com/video/X46F2-OxxqY/w-d-xo.html
Como podria introducir partial pivoting con este codigo?
+Anel Ruiz Hola Anel! Aquí te dejamos un código que encontramos en la red y te puede ayudar para el pivoteo: anh.cs.luc.edu/331/code/gaussJordan2.py nosotros estamos por hacer una versión actualizada del método, te avisamos cuando este lista, por si te es de utilidad. Saludos!
Buenos días! Su código funciona muy bien y he hecho los cambios que dijo y funciona fenomenal. Sin embargo, estoy intentando poner las operaciones en funciones, con un main incluido, y me está costando mucho. Tiene alguna idea de como hacerlo? Si me pudiera ayudar con el código con las funciones le agradecería mucho. Muchas gracias!
Una pregunta como seria el código para Gauss Jordan?
Hola Jackeline! Hay que modificar los ciclos de las operaciones en renglones y columnas, aquí te dejamos el código:
for k in range (0,m):
for r in range(0,m):
if (r != k):
Muchas gracias, y tengo otra duda quisiera calcular el determinante en numpy para poder establecer si el sistema es singular. Pero no se en que parte debería colocarlo del código cctmexico
Otra pregunta, probé el código que me indicaste pero como hago para que hayan unos en la diagonal?
Este metodo es tambien conocido como Gauss simple, o son diferentes?
¡Hola! No habíamos escuchado ese nombre para el método ¿De dónde nos vés?
@@cctmexico Saludos! Desde México, específicamente desde Tijuana
:O Paisano! Parece que es el mismo checa: www.cimat.mx/~dnvr30/CursoMinas/elim_gauss.htm
@@cctmexico Gracias! También por todo sus vídeos, realmente tienen temas muy variados y útiles... Por cierto, el lenguaje de pyrhon me encanta jaja
De nada ^_^ Qué bueno que te guste el contenido del canal y si Python es literalmente encantador!!
Hola disculpa me aparece este error:
NameError: name 'raw_input' is not defined
en este caso ¿Qué solución puede haber?
Borra los raw_
estoy desesperada T-T sale del rango al querer cambiar el valor a las columnas
n=len(matriz)
for k in range(0, n):
for r in range(k+1,n):
factor=(matriz[r][k]/matriz[k][k])
vector[r]=vector[r]-(factor*vector[k])
for c in range(0, n):
matriz[r][c]=matriz[r][c]-(factor*matriz[k][c])
Hola Yolanda! ¿Qué tipo de error te manda? Aquí te dejamos un vídeo que hicimos para detectar y arreglar errores en Python, por si te es de utilidad:
th-cam.com/video/X46F2-OxxqY/w-d-xo.html