%programa;
clc;clf;
ks=0.05;
kd=0.1;
ka=0.6;
Os=input('ingrese valor de oxigeno disuelto');
while Os<0;
Os=input('ingrese valor de oxigeno disuelto');
end
Lo=input('ingrese valor de concentración BOD en el punto mezclado');
while Lo<0;
Lo=input('ingrese valor de concentración BOD en el punto mezclado');
end
Sb=input('ingrese valor de demanda de oxigeno sedimentado');
while Sb<0;
Sb=input('ingrese valor de demanda de oxigeno sedimentado');
end
tf=input('ingrese el tiempo de travesia final');
for t=1:tf
Ox=Os-kd*Lo*(exp(-ka*t)-exp(-(kd+ks)*t))/(kd+ks-ka)-Sb*(1-exp(-kd*t))/ka;
O=double(subs(Ox,t));
hold on
plot(t,Ox,'-og')
disp('********************************')
disp('tiempo concentración de oxígeno disuelto')
disp([t' O']);
end
%programa;
clc;clf;
ks=0.05;
kd=0.1;
ka=0.6;
Os=input('ingrese valor de oxigeno disuelto');
while Os<0;
Os=input('ingrese valor de oxigeno disuelto');
end
Lo=input('ingrese valor de concentración BOD en el punto mezclado');
while Lo<0;
Lo=input('ingrese valor de concentración BOD en el punto mezclado');
end
Sb=input('ingrese valor de demanda de oxigeno sedimentado');
while Sb<0;
Sb=input('ingrese valor de demanda de oxigeno sedimentado');
end
tf=input('ingrese el tiempo de travesia final');
i=1;
Ocr=10000;
tcr=0;
for t=1:tf
OX=Os-kd*Lo*(exp(-ka*t)-exp(-(kd+ks)*t))/(kd+ks-ka)-Sb*(1-exp(-kd*t))/ka;
if Ocr>OX
Ocr=OX;
tcr=t;
end
i=i+1;
O=double(subs(OX,t))
hold on
plot(t,O,'-og')
disp('********************************')
disp('tiempo concentración de oxígeno disuelto')
disp([t' O']);
end
clear all;
clc;
% Datos
Os=10;
ks=0.05;
kd=0.1;
Lo=50;
ka=0.6;
Sb=1;
% Proceso para Graficar
i=1;
Ocr=10000;
tcr=0;
% Método Streeter-Phelps
for t=0:0.2:20
O(i)=Os-kd*Lo*(exp(-ka*t)-exp(-(kd+ks)*t))/(kd+ks-ka)-Sb*(1-exp(-kd*t))/ka;
if Ocr>O(i)
Ocr=O(i);
tcr=t;
end
i=i+1;
end
t=0:0.2:20;
% Graficando la función
plot(t,O,'-b')
hold on
plot(t,O,'ro')
grid on
xlabel('t')
ylabel('O(t)')
title('Funcion O(t)')
disp('Valores Críticos:')
disp('Tiempo Crítico:')
disp(tcr)
disp('Concentración Crítica:')
disp(Ocr)
