Correccion de Programacion
Publicado por Cesar (2 intervenciones) el 22/09/2016 21:57:41
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
clear all, clf, global t O DBOu O0 Os
DBOu=20; Os=8.8;
t=[0 .083 .16 .25 .5 .75 1 2 3 4 5 6 7 8 9 10]';
O=[8.3 8 7.9 7.8 7.7 7.6 7.4 7.2 7 6.8 8.3 6.1 5.6 5.2 4.8 4.2]'; O0=O(1);Opciones=optimset('TolX',1e-20,'TolFun',...
1e-20,'MaxFunEvals',10000)
Est=fminsearch('AutodepurMin',[0.4 0.9],Opciones);
k1=Est(1); k2=Est(2)
tc=linspace(t(1),t(length(t)),100);
Oc=Os-k1*DBOu/(k2-k1)*exp(-k1*tc)+(k1*...
(Os+DBOu-O0)+k2*(O0-Os))/...
(k2-k1)*exp(-k2*tc);
plot(tc,Oc,'k-','LineWidth',1), hold on, plot(t,O,'ko')
axis([0 10 3.5 9])
ylabel('\fontsize{13} Conc. de Oxígeno (mg/L)')
xlabel('\fontsize{13}Tiempo (días)')
text(3,8.65,'\fontsize{13}Conc. de oxígeno a saturación')
line([t(1) t(length(t))],[Os Os],'Color','k','LineWidth',1)
Text1=['\fontsize{13}k_1= ',num2str(round(1000*k1)/1000)];
Text2=['\fontsize{13}k_2= ',num2str(round(100*k2)/100)];
text(7,4.5,Text1), text(7,4,Text2)
%----------------------------------------------------------
%------ Funcion llamada por Autodepur.m -----------
function q=AutodepurMin(p);
global t O DBOu O0 Os
k1=p(1); k2=p(2);
Oc=Os-k1*DBOu/(k2-k1)*exp(-k1*t)+...
(k1*(Os+DBOu-O0)+k2*(O0-Os))/(k2-k1)*exp(-k2*t);
Error=O-Oc; q=sum((Error).^2);
%--------------------------------------------------Valora esta pregunta
-2
