Crear cuenta / - Bitacoras.comAgregador → Enlace permanente ¬

0puntos votar

Resolver Sistema Tridiagonal Disperso Con Scipy, Python

Hola , como paso previo para poder resolver una E.D.P del calor se necesita saber resolver sistemas tridiagonales lineales dispersos. Con Python es relativamente sencillo si utilizamos las librerías Scipy, con las cuales podemos implementar matrices dispersas con funciones específicas como dia_matrix. Una vez implementada la matriz dispersa de nuestro sistema tan sólo hay que elegir el método para resolverlo: descomposición LU, Choleski, etc. En mi caso he utilizado la descomposición LU (también conocida como Algoritmo de Crout) y que viene implementada con Scipy. Éste método es eficiente para resolver la ecuación del calor, además; nuestra matriz de coeficientes dispersa de la ecuación del calor es tridiagonal y cumpliendo las condiciones necesarias para que funcione el algoritmo de Crout: semidefinida positiva, etc. He realizado un pequeño algoritmo en Python 3 y Scipy en el que vemos cómo resolver primero un sistema lineal cualquiera 5 x 5 con Crout y luego cómo aplicarlo con un sistema 7 x 7 en el que nue...
tags
Continuar leyendo

Recomienda esta anotación por e-mail

Ayúdanos a hacer de Bitacoras.com un servicio mejor para todos. Lee nuestros consejos.

Ningún usuario registrado ha votado aún.