|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
( K2 K0 u! E+ o2 f
请问,lsqcurvefit拟合后如何显示相关系数R^2,感谢!9 k2 \7 b2 E& B6 U! X5 d, N; S
clear;clc;
; z! T+ t1 F0 c3 X7 YT=[293.15 298.15 303.15 308.15 313.15];( N( L/ ?/ J5 ^" h% ]0 B9 Q
c=[0.000378156 0.000425433 0.000562474 0.000684178 0.000981096];: B) ]! X! l2 i9 A+ `8 [1 n
TT=T';# g3 i+ o! D. a* s
cc=c';8 r0 h4 }3 o; S$ _( X; E$ ?6 N
x=[TT,cc];
, s! J/ B. b4 Q- M" l$ @ XHm=33509;. s; G, w0 H: a- e; Y
Tm=408.65;
4 c C) e- I9 g4 Z0 na=0.3;
; ^ E# C; N- q zfun=@(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));
9 b2 \' s% l% u3 J! f5 B4 wq0=[0,0];) a! b9 y" ], m8 Q5 _
q=lsqcurvefit(fun,q0,x,zeros(size(x,1),1))9 C: T4 a) Z! ]
plot(x(:,1),x(:,2),'r^');hold on;6 Q4 }. s* p2 r }$ ]3 M; x" M
ezplot(@(x,c)fun(q,[x,c]),[290,314,0.0003,0.001]);: B1 [) j) [+ x
h=ezplot(@(x,c)fun(q,[x,c]),[290,314,0.0003,0.001]);
7 r# j5 r) V* M7 ~! J7 g |
|