comment
IRC Chat
play_arrow
Este sitio utiliza cookies propias y de terceros. Si continúa navegando consideramos que acepta el uso de cookies. OK Más Información.

[Uso de matplotlib y numpy] Potencial de Morse

  • 1 Respuestas
  • 1166 Vistas

0 Usuarios y 1 Visitante están viendo este tema.

Desconectado fortil

  • *
  • Underc0der
  • Mensajes: 175
  • Actividad:
    0%
  • Reputación 0
    • Ver Perfil
    • Email
« en: Octubre 05, 2012, 01:02:29 am »
Hola a todxs

Como yo estudio química me incursioné un poco en python para científicos (matplotlib y numpy, entre otros) y he echo una programita sobre el potencial del Morse, que sirve para describir la energía de las funciones de onda etc. xD

acá se los dejo, además puede servir de ejemplo para muchas cosas mas, como transformar nuestra cámara normal en una "infrarojo" (entre comillas porque no mas es efecto de imagen) luego le investigo mas sobre el tema y les paso como.

Código: Python
  1. ##POTENCIAL DE MORSE BY FORTIL##
  2.  
  3. import numpy as np
  4. import matplotlib.pyplot as plt
  5.  
  6. def rango(x,x1):
  7.     return np.arange(x,x1,0.1)
  8.  
  9. def f(x,v0):
  10.     return v0*(1-x)**2
  11.  
  12. def ex(x, a=0.5, x0 = 2.0, v0 = 0.05):
  13.     return f(np.exp(-a*((x)-x0)),v0)
  14.  
  15. def coseno(x,a):
  16.     return a+0.0005*np.cos(10*x+np.pi)
  17.    
  18. x = np.arange(0, 80.0, 0.1)
  19.  
  20. fig = plt.figure()
  21. ax1 = fig.add_subplot(111)
  22.  
  23. ax1.plot(x, ex(x), 'b-',x,ex(x,0.05),'r-')
  24. ax1.set_title("Potencial de Morse por Fortil", fontsize=16)
  25. ax1.set_ylabel(r'$V_f$')
  26. ax1.set_xlabel(r'$X$')
  27.  
  28. ax1.text(20,0.025,r"$\alpha_1 $ mayor",fontsize=12, color='r')
  29. ax1.text(7,0.035, r"$\alpha_2 $ menor", fontsize=12, color='b')
  30. ax1.text(60,0.051, r"$V_0 = 0.05$", fontsize=12, color='y')
  31. ax1.text(60,0.009, r"$X_0 = 2.0$", fontsize=12, color='black')
  32. ax1.text(20,0.06,r"$\alpha_2 = \frac{\alpha_1}{10} $",
  33.          fontsize=14, color='g')
  34. ax1.text(25,0.12,r"$V_f = V_0\times[1-e^{-\alpha(X-X_0)}]$",
  35.          fontsize=18, color='black')
  36.  
  37. ax1.plot(rango(ex(0.6),80.0),
  38.          coseno(rango(ex(0.6),80.0),0.05),'c-')
  39. ax1.acorr(x, usevlines=False, normed=True, maxlags=1, lw=2)
  40. ax1.set_ylim(0,0.16)
  41. ax1.set_xlim(0,70)
  42.  
  43. plt.show()
  44.  
  45.  

Conectado ANTRAX

  • *
  • Administrator
  • Mensajes: 5401
  • Actividad:
    41.67%
  • Reputación 31
  • ANTRAX
    • Ver Perfil
    • Underc0de
    • Email
  • Skype: underc0de.org
  • Twitter: @Underc0de
« Respuesta #1 en: Octubre 05, 2012, 08:54:42 am »
Buenisimo fortil!
viene de 10! muchas gracias


 

¿Te gustó el post? COMPARTILO!



Cifrado Morse

Iniciado por Sanko

Respuestas: 6
Vistas: 3274
Último mensaje Septiembre 08, 2013, 12:55:17 am
por arthusu
[Visión artificial: parte II] Introducción a OpenCV y Numpy

Iniciado por Once

Respuestas: 0
Vistas: 1475
Último mensaje Abril 11, 2017, 02:33:30 am
por Once
[Visión artificial: parte III] Introducción a NumPy

Iniciado por Once

Respuestas: 0
Vistas: 1172
Último mensaje Julio 12, 2017, 12:52:49 am
por Once