|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
& n+ M7 E5 e4 a H2 B8 z请问,lsqcurvefit拟合后如何显示相关系数R^2,感谢!
. j" r2 M8 ?- m, g# b, G( ?clear;clc;. n/ E8 p5 x% \$ ^7 d
T=[293.15 298.15 303.15 308.15 313.15];% C" L6 J" t4 G1 E* }
c=[0.000378156 0.000425433 0.000562474 0.000684178 0.000981096];
. f, l2 V) Y* z/ K5 j. n& z& pTT=T';
; [9 J* a$ R1 }# B7 rcc=c';
6 M! D5 Z1 [/ n4 a7 o6 wx=[TT,cc];
3 |7 }* E( x9 p6 Z6 \5 DHm=33509;" @: C! Z2 d' d
Tm=408.65;
0 e1 O" n. v- L: }, O- _a=0.3;4 |4 Y7 Y Z# P) @ G7 d
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));
8 ~) ?4 k0 F, g+ C: E9 O3 e" @+ y: zq0=[0,0];$ Y- ^ g4 n' |+ T7 P
q=lsqcurvefit(fun,q0,x,zeros(size(x,1),1))
/ S' K5 w; M% Q }plot(x(:,1),x(:,2),'r^');hold on;; k7 z: Y: n: s
ezplot(@(x,c)fun(q,[x,c]),[290,314,0.0003,0.001]);2 \, V) y- Q/ J/ }
h=ezplot(@(x,c)fun(q,[x,c]),[290,314,0.0003,0.001]);
' Z+ D% ~, I3 T ]& ^/ B |
|