|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
0 o ^# U, P% [" {请问,lsqcurvefit拟合后如何显示相关系数R^2,感谢!% P) W0 V7 a# d$ d/ @0 D
clear;clc;, S& Y, [0 Q3 O' [# j- w% f* Q
T=[293.15 298.15 303.15 308.15 313.15];
$ {& K2 `( V& J* T6 Yc=[0.000378156 0.000425433 0.000562474 0.000684178 0.000981096];
1 ?5 G; Y) z2 d8 }TT=T';
% _5 N8 `" M/ acc=c';
1 ^9 t( O$ z6 B. D4 Dx=[TT,cc];
% g% y5 @0 h, f6 b7 SHm=33509;
7 }; t( }7 i0 E( P+ H! a; ^Tm=408.65; `0 }$ v, i: y$ S# s
a=0.3;% |0 p, V# B2 J! C3 ?; S
fun=@(q,x)(-Hm.*(x(:,1)./Tm-1)./(8.314.*x(:,1))+log(x(:,2))+(1-x(:,2)).^2.*((q(1)./(8.314.*x(:,1))).*(exp(-a.*q(1)./(8.314.*x(:,1)))./(x(:,2)+(1-x(:,2)).*exp(-a.*q(1)./(8.314.*x(:,1))))).^2+(q(2)./(8.314.*x(:,1))).*exp(-a.*q(2)./(8.314.*x(:,1)))./(1-x(:,2)+x(:,2).*exp(-a.*q(2)./(8.314.*x(:,1)))).^2));1 q7 ^' k( X4 ]2 x8 F% f
q0=[0,0];+ U8 Y! x+ t/ A
q=lsqcurvefit(fun,q0,x,zeros(size(x,1),1)) q @# t3 {' |' }, J9 B
plot(x(:,1),x(:,2),'r^');hold on;1 D" @6 B2 I1 R% l$ R8 A
ezplot(@(x,c)fun(q,[x,c]),[290,314,0.0003,0.001]);
% M0 U7 N) _/ I; j2 h, |h=ezplot(@(x,c)fun(q,[x,c]),[290,314,0.0003,0.001]);- U+ k. f( R$ g( V8 \9 Y; K
|
|