|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
- c. t/ }2 d" B) N2 L I请问,lsqcurvefit拟合后如何显示相关系数R^2,感谢!
( p* n! a0 j H- Zclear;clc;
]" `# E$ B) | Q ]# x- CT=[293.15 298.15 303.15 308.15 313.15];
$ Q- a; ?( A; B& D/ F: J# W6 m3 n/ Oc=[0.000378156 0.000425433 0.000562474 0.000684178 0.000981096];1 [4 \. ^( N# q/ R
TT=T';
# ]* X. S& E {( xcc=c';) e, Q) B( d0 V8 }# W" q
x=[TT,cc];
/ q* c a& z- H* m- H1 K3 Z: p& h$ qHm=33509;
, A, L) S( m6 u& ?6 e/ [9 E/ OTm=408.65;
0 g3 l" @4 O* x" ]8 Ra=0.3;
. }! O" _4 S/ J8 R3 z5 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));5 O9 \, i1 @0 K% x p
q0=[0,0];6 y' L9 _2 D/ S0 ^* D) \
q=lsqcurvefit(fun,q0,x,zeros(size(x,1),1))
1 R' Z7 }8 G" X U" O- X$ Yplot(x(:,1),x(:,2),'r^');hold on;% M9 @; d. n+ K9 a$ M, I
ezplot(@(x,c)fun(q,[x,c]),[290,314,0.0003,0.001]);
" D Q! x4 O0 z" zh=ezplot(@(x,c)fun(q,[x,c]),[290,314,0.0003,0.001]);
' `6 g4 s8 }$ W, y! e' O5 { |
|