|
|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
2 g/ k* h/ Z0 _5 S: i请问,lsqcurvefit拟合后如何显示相关系数R^2,感谢!7 l, y8 N3 c- G& }' U" K9 p0 V$ k
clear;clc;: _8 s* U: | {! b4 S& G# }
T=[293.15 298.15 303.15 308.15 313.15];% x# F/ g8 v2 a: k* N. C
c=[0.000378156 0.000425433 0.000562474 0.000684178 0.000981096];
6 {; t6 H; `3 t0 C3 b/ _TT=T';6 V3 X2 m8 B* o# t
cc=c';
! Z. r6 j- g8 z' d- a& ux=[TT,cc];, \4 F% u- S4 }: C& {6 q% w
Hm=33509;
3 o% J1 k0 o% M6 p! MTm=408.65;
|" S4 ]6 [7 J) s* |' y+ W. z# `a=0.3;
9 r9 y5 w, D! yfun=@(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));
3 U) S9 R) A$ O& j. ]q0=[0,0];# ^' m& R$ ~ o% \* Y' g
q=lsqcurvefit(fun,q0,x,zeros(size(x,1),1))+ U" T; N' n# p/ R* q
plot(x(:,1),x(:,2),'r^');hold on;7 }' I2 z. I- d( `
ezplot(@(x,c)fun(q,[x,c]),[290,314,0.0003,0.001]);$ N/ d; Z5 H, p1 F0 j
h=ezplot(@(x,c)fun(q,[x,c]),[290,314,0.0003,0.001]);# @! D3 y6 I$ q3 Q/ n, |
|
|