videotutoriales de matlab, octave, numpy, python, POA, algoritmia, simulacion de procesos en ingenieria quimica, y algo de ingenieria quimica
miércoles, 5 de agosto de 2015
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);
martes, 30 de junio de 2015
martes, 23 de junio de 2015
martes, 17 de marzo de 2015
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;
Suscribirse a:
Entradas (Atom)