|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
$ C c& d/ d: J! d( i: A* ^" Z3 T请问,lsqcurvefit拟合后如何显示相关系数R^2,感谢!
* I; R; d* k* Z* a, zclear;clc;+ ?* {: k% i/ v
T=[293.15 298.15 303.15 308.15 313.15];
4 F( k4 A# y* f% E# Nc=[0.000378156 0.000425433 0.000562474 0.000684178 0.000981096];
: z" Y% h( y% J, s5 ^. z1 WTT=T';. _, R# I6 h- R5 |1 d
cc=c';$ m- V* F: o2 k' x3 }
x=[TT,cc];; Y o- S" C0 n! M3 I" R3 e) r1 S
Hm=33509;
& N5 k: _6 _; q* V, xTm=408.65;2 x2 m& S9 p9 r" }8 G/ ^ i" ?$ ?
a=0.3;$ s+ J. A4 v( \( S9 i/ @$ c: B
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));( u( m( ^: B- H H( Y8 S5 L
q0=[0,0];& [ a1 P! O t- ]- o; P
q=lsqcurvefit(fun,q0,x,zeros(size(x,1),1))/ u! T( h }, C* m0 a
plot(x(:,1),x(:,2),'r^');hold on;9 |/ b# y4 z: ?. s `+ u
ezplot(@(x,c)fun(q,[x,c]),[290,314,0.0003,0.001]);
- n/ }+ m5 N8 M7 Ph=ezplot(@(x,c)fun(q,[x,c]),[290,314,0.0003,0.001]);
6 Z7 Y: o2 Y7 X4 Q |
|