miércoles, 1 de julio de 2015

Ejemplo de un modelo de reactor lote. En Matlab

function lotecomplejo
[t,x]=ode45(@lote,[0 2],[16 0 0 0])
subplot(4,1,1); plot(t,x(:,1))
subplot(4,1,2); plot(t,x(:,2))
subplot(4,1,3); plot(t,x(:,3))
subplot(4,1,4); plot(t,x(:,4))


function dx=lote(t,x)

% definicion de parametros

k1=5;
k2=4;
k3=3;
k4=3;

dx=zeros(4,1);

% ecuaciones de diseño (balances de masa), x1 = Ca, x2 = Cb, x3 = Cc, x4 = Cd

dx(1) = - k1*x(1);
dx(2) = k1*x(1) - k2*x(2) + k3*x(3) - k4*x(2);
dx(3) = k2*x(2) - k3*x(3);
dx(4) = k4*x(2);

Método de McCabe Thiele para destilación fraccionada.

**** DESCARGAS ****

lunes, 12 de enero de 2015

Método de Regula Falsi en Matlab

Función principal 


function []=regu_fal(x1, x2, n);

fx1=f(x1);
fx2=f(x2);
c = x2 - ((f(x2)*(x2-x1))/(f(x2) - f(x1)))
fc=f(c)
if fx1*fx2<0;
     for i=1:n
     if fc<0
        x1=c;
        fx1=f(x1);
        c = x1 - (fx1*(x1-x2)/(fx1-fx2))
        fc = f(c)
       else
          x2=c;
          fx2=f(x2);
          c = x2 - (fx2*(x2-x1)/(fx2-fx1))
          fc = f(c)
      end
    end
else
    disp('Error en los Valores X1, X2');
end
end
Ecuación a resolver.

function fx = f(x);
	fx = x^3-9*x+1;  %% puedes cambiar por la funcion que gustes
	return;