|
|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
7 C7 Q, D+ ^4 G6 A$ L请问,lsqcurvefit拟合后如何显示相关系数R^2,感谢!
) b2 \% R$ ^; z, o8 Uclear;clc;8 L, w* T$ a' [$ c% s1 s$ S) x. ]
T=[293.15 298.15 303.15 308.15 313.15];
9 B- G& J( h* o+ z% M; E7 Sc=[0.000378156 0.000425433 0.000562474 0.000684178 0.000981096];2 C% k4 D" o$ v7 r! p: J& S6 I
TT=T';
a5 u$ }' P1 Vcc=c';
+ U: X( }4 B5 ~+ S. Ox=[TT,cc];5 {7 H c/ r; D2 K2 F* D. w
Hm=33509;
! b- e2 \& a7 x6 ?Tm=408.65;0 q% r+ R1 B' F; n7 T9 j
a=0.3;. H: g5 K% q' _4 `8 l
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));
' J6 U( a+ B( q ? ?q0=[0,0];
B" g; I& ?5 l4 R4 g* Gq=lsqcurvefit(fun,q0,x,zeros(size(x,1),1))
7 k4 P4 e& m1 `2 M2 H8 Lplot(x(:,1),x(:,2),'r^');hold on; @' {& f' p( D6 M, s" V
ezplot(@(x,c)fun(q,[x,c]),[290,314,0.0003,0.001]);
. e' P, z7 R( {h=ezplot(@(x,c)fun(q,[x,c]),[290,314,0.0003,0.001]);( F7 Y0 }! I# ?' \
|
|