miércoles, 14 de agosto de 2013

orderseek--Encontrar el orden de una reacción Química en Matlab

Este pequeño programa hecho en Matlab nos sirve para hallar el orden de una reacción química, cabe señalar que este programa no es creación mía, pero lo estoy usando como punto de partida  para crear un propio programa ya que en el caso de reacciones químicas tenemos muchos casos como son, volumen constante, volumen Variable, Reversibles, etc. 

jueves, 8 de agosto de 2013

Python para Bioinformatica — Introducción al Comando Seq (parte 2)

Python para BioInformatica introducción al comando Seq (parte 1)

- Cortando una Secuencia.
>>> from Bio.Seq import Seq
>>> from Bio.Seq import IUPAC
>>> secuencia = Seq ("ACGCTAGCTAGCTCGCAACGCTCGAACGCTCG")
>>> len(secuencia)
32
>>> secuenciacorta = secuencia[5:10]
>>> secuenciacorta
Seq('AGCTA', Alphabet())
-Convertir de String a Secuencia
 este tipo de conversión es muy útil cuando tenemos un archivo de texto plano o introducimos    mediante un programa nuestra secuencia.
>>> secuencia = input('introduce tu secuencia: ')
introduce tu secuencia: AGCTGCGTG
>>> secuencia
'AGCTGCGTG'
>>> type(secuencia)
<class 'str'>
>>> Secuencia2 = Seq(secuencia)
>>> type(Secuencia2)
<class 'Bio.Seq.Seq'>
>>> ## se puede hacer en la misma linea
>>> SECUENCIA = Seq(input('introduce tu secuencia: '))
introduce tu secuencia: AGCTCGACGCT
>>> type(SECUENCIA)
<class 'Bio.Seq.Seq'>

-Transcripción y Translación
>>> secuencia = Seq ("ACGCTAGCTAGCTCGCAACGCTCGAACGCTCG")
>>> transcripcion =secuencia.transcribe()
>>> transcripcion
Seq('ACGCUAGCUAGCUCGCAACGCUCGAACGCUCG', RNAAlphabet())

>>> traslacion = secuencia.translate()
>>> traslacion
Seq('TLASSQRSNA', ExtendedIUPACProtein())

lunes, 5 de agosto de 2013

Diagrama Liquido- Vapor en un proceso Isobarico-- usando Matlab

El modelo de Wilson permite que el programa para tener en cuenta la desviación del comportamiento ideal en la fase líquida. El usuario debe introducir bajo presiones moderadas ya que el término de corrección a alta presión no está incluido.

martes, 30 de julio de 2013

"Matlab for Genius", Busqueda de la solucion de una Ecuacion--- Método Grafico





Hace unos días atrás me enviaron una consulta de como encontrar la solución de  una ecuación en matlab.
supongamos que tenemos la siguiente ecuación.
2log (x) = 2 - (x)^(1/2)
en la que necesitaros encontrar el punto en el que se encuentra nuestra solución, para ello usaremos  en esta entrega el método gráfico.




>> %declaramos un vector x donde acotaremos nuestra solucion
>> x = 1:0.01:10;
>> % formaremos dos variables en las que  estaran nuestras ecuaciones
>> f1 = 2*log(x);
>> f2 = 2 -x.^(1/2);
>> plot(x,f1)
>> hold on
>> plot(x,f2)
>> hold off
el gráfico nos queda de la siguiente manera, donde el punto de intersección  es la solución.


Multiples estados estacionarios, estudio de un Catalizador

El problema catalizador adiabática se presenta en el campo de la ingeniería química. Si cinética de las reacciones químicas dentro de un cuerpo de catalizador sólido se ven influidas por la transferencia de masa múltiples estados estacionarios son posibles para una cierta gama de los parámetros que rigen: número de Arrhenius, número Prater y módulo de Thiele.