matlab建模方法(1).ppt

文档介绍：模型程序
以MONOD方程为例
.
1
function y=ans_f_monod(x)
%row1:t,row2:x,row3:s
data=x;
Lb=[ 0];
Ub=[ 10 1];
options=optimset('largescale','off','display','final','tolx',1e-16,'tolfun',1e-16,'MaxFunEvals',6000,'MaxIter',1000);
[p,fval,exitflag,output]=fmincon(@(x) obj_f_monod(x,data),Lb,[],[],[],[],Lb,Ub,[],options);
disp('um=');
disp(p(1));
disp('Ks=');
disp(p(2));
.
2
disp('Yx/s=');
disp(p(3));
disp('m=');
disp(p(4));
disp('fval=');
disp(fval);
disp('exitflag=');
disp(exitflag);
disp('output=');
disp(output);
disp('revelent error=');
error=result_f_monod(p,data);
disp(error);
.
3
function y=obj_f_monod(x,data)
tt=data(:,1);
xx=data(:,2:3);
x0=data(1,2:3);
a=x;
[t,xj]=ode45(@(t,x) f_monod(t,x,a),tt,x0);
dx=xj-xx;
dx1=dx./xx;
y=sum(sum(dx1.^2));
.
4
function y=f_monod(t,x,a)
%a1:um;a2:Ks;a3:Yx/s;a4:m;x1:cx;x2:cs
y=[a(1)*x