《matlab建模方法》由会员分享,可在线阅读,更多相关《matlab建模方法(7页珍藏版)》请在金锄头文库上搜索。
1、matlabmatlab建模方法建模方法nfunction y=ans_f_monod(x)n%row1:t,row2:x,row3:sndata=x;nLb=0.05 0.05 0.05 0;nUb=0.5 10 0.4 1;noptions=optimset(largescale,off,display,final,tolx,1e-16,tolfun,1e-16,MaxFunEvals,6000,MaxIter,1000);np,fval,exitflag,output=fmincon(x) obj_f_monod(x,data),Lb,Lb,Ub,options);ndisp(um=);
2、ndisp(p(1);ndisp(Ks=);ndisp(p(2);ndisp(Yx/s=);ndisp(p(3);ndisp(m=);ndisp(p(4);ndisp(fval=);ndisp(fval);ndisp(exitflag=);ndisp(exitflag);ndisp(output=);ndisp(output);ndisp(revelent error=);nerror=result_f_monod(p,data);ndisp(error);function y=obj_f_monod(x,data)tt=data(:,1);xx=data(:,2:3);x0=data(1,2
3、: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);nfunction y=f_monod(t,x,a)n%a1:um;a2:Ks;a3:Yx/s;a4:m;x1:cx;x2:csny=a(1)*x(2)/(a(2)+x(2)*x(1);n -a(1)*x(2)*x(1)/(a(2)+x(2)/a(3)-a(4)*x(1);nfunction y=result_f_monod(p,data)ntt=data(:,1);nxx=data(:,2:3);nx0=data(1,2:3);nt,xj=ode45(t,x) f_monod(t,x,p),tt,x0);ndx=xj-xx;ndx1=dx./xx;nn=length(data(:,1);ny=sum(sum(abs(dx1)/n/2;nplot(tt,xx,tt,xj,*);ntitle(细胞生长及底物消耗动力学曲线拟合);nxlabel(t(h);nylabel(Cx(g/L) Cs(g/L);结束!结束!