|
|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
9 `. s1 V. \# M0 \% I
请问,lsqcurvefit拟合后如何显示相关系数R^2,感谢!
0 K8 q, \& l) ?6 m, k, Qclear;clc;, ]) b) R- M: d* k: S$ R: B
T=[293.15 298.15 303.15 308.15 313.15];
" [, Q' N$ ]! G1 b, @$ J/ D0 Pc=[0.000378156 0.000425433 0.000562474 0.000684178 0.000981096];. X, v y, @5 Q: k4 D
TT=T';% G: F$ W9 m @$ t- G) {1 P- `5 G
cc=c';
) y7 r5 m; [- c" |4 B' s+ o9 i, zx=[TT,cc];, ]) m6 [) G/ U/ S( L4 N( m
Hm=33509;
7 P. [- `' l1 m$ e( `Tm=408.65;1 Q* e& i& r+ j M0 F# P4 ?( ^
a=0.3;# c/ ^. P/ ~& Z( f
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));% k; p, c- T7 n. ~- q$ E* I- g
q0=[0,0];: o/ l. i. ~- [
q=lsqcurvefit(fun,q0,x,zeros(size(x,1),1))) j+ @0 [+ j. r4 Q9 o
plot(x(:,1),x(:,2),'r^');hold on;
) v% C# L2 v" J/ e6 Tezplot(@(x,c)fun(q,[x,c]),[290,314,0.0003,0.001]);. ?- \4 s8 J, M0 L0 S" q
h=ezplot(@(x,c)fun(q,[x,c]),[290,314,0.0003,0.001]);
$ F+ A) ~! k1 j+ \" ~% k |
|