| 查看: 2081 | 回复: 0 | ||
[求助]
matlab 拟合优度的判断语句？
|
|
用matlab软件中的lsqnonlin 拟合了几条曲线，想判断拟合的好坏程度，可是lsqnonlin输出的是误差平方和及残差等，加上什么样的语句可以显示类似拟合优度R2，及平均相对误差这类的数据呢？谢谢～ 主程序： format short global Umax Ks Kp Y a Uo Ko D Ux Kx Ku Ki; Ki=0.0005;Ku=0.0001; k0=[0.2 0.5 0.1 10 0.01 0.01 0.01 0.1 0.01 0.01]; %参数的初始值 x0=[0.0514 0.250 0]; %菌体浓度、底物浓度和产物浓度的初始值 t1=[0 6 12 18 24 30 36 42 48 54 60 66]';%发酵时间的实验数据 tspan=[0 6 12 18 24 30 36 42 48 54 60 66 ]';%时间间隔 %yexp实验数据[菌体浓度x1、底物浓度x2和产物浓度x3] yexp=[[0.0514 0.0711 0.0997 0.1887 0.3478 0.6634 0.8711 1.1023 1.2336 1.4239 1.5849 1.6739];[0.25 0.2447 0.2366 0.2258 0.2196 0.203 0.1872 0.1698 0.1453 0.1124 0.0855 0.0608];[0 0.0012 0.0021 0.0039 0.0088 0.0137 0.0183 0.0268 0.0361 0.0438 0.0557 0.0643]]'; lb=[0 0 0 0 0.0001 0.0116 0.001 1 0.001 0.001];ub=[1 10 13 50 50 2 1 1.001 1 1];%参数的下、上限 [k,resnorm,residual,exitflag,output,lambda,jacobian]=lsqnonlin(@ObjFunc4LNL83,k0,lb,ub,[],x0,yexp);%非线性拟合 y1=[yexp(:,1)];y2=[yexp(:,2)];y3=[yexp(:,3)]'; [t4plot,x4plot]=ode45(@kineticseqs83,[0 100],x0,[],k); plot(t1,y1,'b*',t4plot,x4plot(:,1),'k-'),xlabel('T(h)'),ylabel('Cx/（g/L）'); figure plot(t1,y2,'g*',t4plot,x4plot(:,2),'k-'),xlabel('T(h)'),ylabel('Cs/(mmol/L)'); figure plot(t1,y3,'r*',t4plot,x4plot(:,3),'k-'),xlabel('T(h)'),ylabel('Cp/(mmol/L)'); disp(k) ======================================== function dxdt=kineticseqs83(t,x,k) %模型方程 global Umax Ks Kp Y a Uo Ko D Ux Kx Ku Ki Umax=k(1); Ks=k(2); Kp=k(3); Y=k(4); a=k(5); Uo=k(6); Ko=k(7); D=k(8); Ux=k(9); Kx=k(10); dxdt1=Umax*x(2)*x(1)*(1-Kp*x(3))/(Ks+x(2)); if dxdt1<0 dxdt1=0; end dxdt3=a*dxdt1+Uo*x(2)*x(1)/(Ko+x(2))*(Ki*x(1)/(Ki+x(3))); if dxdt3<0 dxdt3=0; end dxdt2=-1/Y*dxdt1-Ux*x(2)*x(1)/(Kx+x(2))*(Ku*x(1)/(Ku+x(3)))-D*dxdt3; dxdt=[dxdt1;dxdt2;dxdt3]; ========================================== function f=ObjFunc4LNL83(k,x0,yexp) tspan=[0 6 12 18 24 30 36 42 48 54 60 66 ]; [t1,x]=ode45(@kineticseqs83,tspan,x0,[],k); y(:,1)=x(:,1);y(:,2)=x(:,2);y(:,3)=x(:,3); f1=y(: ,1)-yexp(: ,1);f2=y(: ,2)-yexp(: ,2);f3=y(: ,3)-yexp(: ,3); f=[f1 f2*5 f3*10]; |
回复此楼» 猜你喜欢
如何从铁电相到顺电相。
已经有1人回复
-
各位大佬，求一份最新Wien2k版本作为学习用途
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有193人回复
-
南科大活性流体和软物质课题组诚招2027级博士、硕士研究生和博后
已经有42人回复
-
瑞典林雪平大学博士后招聘｜PEC Water Splitting 方向
已经有24人回复
-
PRB投稿，Acknowledgment sent to author状态十天了
已经有3人回复
-
PVA溶解
已经有2人回复
-
基金申请
已经有44人回复
-
CSC与新西兰维多利亚大学PhD奖学金项目
已经有0人回复
-
新西兰Robinson研究所 招聘CSC公派访问人员
已经有0人回复
-
帮我的英语口语老师找学生
已经有0人回复
找到一些相关的精华帖子，希望有用哦~
递推最小二乘法在什么书本中有详细的讲解和matlab仿真介绍？
已经有5人回复
- Matlab数学手册
已经有88人回复
- 关于Matlab作图的问题
已经有14人回复
- 联苯类手性的判断
已经有22人回复
- 在网络层如何判断HTTP协议的请求
已经有7人回复
- matlab 求指点 动力学方程拟合过程中导数的获取
已经有10人回复
- matlab作业哪位大神麻烦给做一下,十万火急,周五之间要交
已经有3人回复
- 热力学第二定律的疑问——自发过程的判断？
已经有7人回复
- MATLAB 数据导入
已经有13人回复
- 【求助】matlab 中，几个m文件调用
已经有14人回复
- MATLAB安装后出现一堆java的问题
已经有12人回复
- 【求助】帮我编个程啊~C语言，C++，Matlab，VB...都行~
已经有14人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
100