|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
0 g) x6 Q& k' A5 q请问,lsqcurvefit拟合后如何显示相关系数R^2,感谢!
3 ~- |: }1 x+ k+ gclear;clc;
* N* b% B( o! P9 e4 ]% _$ x5 kT=[293.15 298.15 303.15 308.15 313.15];
2 S4 ]2 d4 H {- B6 p+ p/ Vc=[0.000378156 0.000425433 0.000562474 0.000684178 0.000981096];
- U" C5 l& r( E) zTT=T';
: \5 U+ c9 k) N: f( xcc=c';0 L" |/ X5 F3 n! J
x=[TT,cc];
% X) _/ p, B: X: \% U. [, {, cHm=33509;
: m2 K) v/ i/ s: u+ p, w4 gTm=408.65;9 s( G6 @+ C& c% Y$ A( F/ e! }* I) S
a=0.3;) P/ j! z- F2 [1 e
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));
2 Y1 c" m3 @; \" Rq0=[0,0];, ^- s( J1 g" |, J) n
q=lsqcurvefit(fun,q0,x,zeros(size(x,1),1))
0 K, z1 `; q" J# S$ Y( A+ _/ lplot(x(:,1),x(:,2),'r^');hold on;& z+ S4 h6 ~ {- ]! w2 S
ezplot(@(x,c)fun(q,[x,c]),[290,314,0.0003,0.001]);. Y& u1 P/ {1 G* v9 i) m5 n* N
h=ezplot(@(x,c)fun(q,[x,c]),[290,314,0.0003,0.001]);# G# V: H3 k6 L# C6 \9 F
|
|