• Author / Uploaded
• AntonioMonteolivaGarcía

Citation preview

clear clc % Ajuste potencial H=[1.50,2.00,2.50,3.00,3.50,4.00]-1; Qo=[1.70,5.00,9.00,14.00,20.00,26.00]; Y=log(Qo); X=log(H); p=polyfit(X,Y,1); d=p(1); D=exp(p(2)); subplot(3,2,1),plot(Qo,H,'o',D*[0:0.01:4].^d,[0:0.01:4]),title('Ajuste potencial'),xlabel('Qout, m3/s');ylabel('H, m'); % Datos Hinicial=1; % m A=20000; % m2 k=0.1/24/3600; % s-1 Cin=100*1000; % mg/m3 % Caudal de entrada, Qin t=[0,50,100,200]*60; Q=[0,10,0,0]; tin=[0:1:t(end)]; Qin=interp1(t,Q,tin); % m3/s subplot(3,2,2),plot(tin/60,Qin),title('Caudal de entrada'),xlabel('t, min'),ylabel('Qin, m3/s'); V(1)=A*Hinicial; C(1)=0; for i=2:length(tin)+1 Qout(i-1)=D*(((V(i-1)/A)-1)^d); V(i)=V(i-1)+Qin(i-1)-Qout(i-1); C(i)=C(i-1)+((Qin(i-1)*Cin-Qout(i-1)*C(i-1)-k*V(i-1)*C(i-1)-C(i1)*(V(i)-V(i-1)))/V(i-1)); end subplot(3,2,3),plot(tin/60,Qout),title('Caudal de salida'),xlabel('t,min'),ylabel('Qout, m3/s'); subplot(3,2,4),plot(tin/60,V(1:length(V)-1)),title('Volumen de la presa'),xlabel('t,min'),ylabel('V, m3'); subplot(3,2,5),plot(tin/60,C(1:length(V)-1)/1000),title('Concentración contaminante'),xlabel('t, min'),ylabel('Concentración, mg/L');