Examen Sustitutorio - Model - 20183 - Luz Marlene Alfaro Asto (2014105100) - Ayacucho
Examen Sustitutorio - Model - 20183 - Luz Marlene Alfaro Asto (2014105100) - Ayacucho
Examen Sustitutorio - Model - 20183 - Luz Marlene Alfaro Asto (2014105100) - Ayacucho
MODELAMIENTO AMBIENTAL
EXAMEN
SUSTITUTORIO En Números
En Letras
x y
1 148.00
2 1487.00
3 3807.00
4 10498.00
5 17551.00
6 34057.00
7 48905.00
8 76987.00
9 109193.00
M-ES20172
10 147413.00
clear
x=[1,2,3,4,5,6,7,8,9,10]
y=[148,1487,3807,10498,17551,34057,48905,76987,109193,147413]
hold on
p=polyfit(x,y,2);
x1=linspace(0.0,10,200);
y1=polyval(p,x1);
plot(x1,y1,'k.-')
hold on
p1=polyfit(x,y,3);
x2=linspace(0.0,10,200);
y2=polyval(p1,x2);
plot(x2,y2,'r.-')
grid minor
hold on
p2=polyfit(x,y,4);
x3=linspace(0.0,10,200);
y3=polyval(p2,x3);
plot(x3,y3,'c.-')
grid minor
hold on
p3=polyfit(x,y,5);
x4=linspace(0.0,10,200);
y4=polyval(p3,x4);
plot(x4,y4,'y.-')
grid minor
M-ES20172
xlabel('x')
ylabel('y')
legend('2°','3°','4°','5°')
hold on
PREGUNTAS 2:
En todos los procesos en este ejemplo se asumen mezcla completa y es entonces aplicable
a cuerpos pequeños de agua tales como Estanques de Control de la Contaminación del
Agua, o como una primera aproximación a cuerpos grandes como lagos.
M-ES20172
Considerando solamente la DBO y el proceso de re-aeración, el modelo matemático
será:
dDBO
K 1 DBO
dt
clear
%tiempo final
tf=20;
%variables
k1=0.1;
k2=0.6;
od=10;
%valores iniciales
X(1)=50;
Y(1)=7;
h=0.1;
n=tf/h;
t=0:h:tf;
%desarrollo por euler
for i=1:n
dX(i)=-k1-X(i);
X(i+1)=X(i)+h*dX(1);
dY(i)=k2*(od-Y(i));
Y(i+1)=Y(i)+h*dY(1);
end
%graficos
hold on
plot(t,X,'r.--')
M-ES20172
plot(t,Y,'y.--')
grid minor
xlabel('tiempo')
ylabel('DBO y reareación')
legend('DBO','Reaireación')
M-ES20172