|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
+ h/ j0 A9 V4 `) i/ |( j
请问,lsqcurvefit拟合后如何显示相关系数R^2,感谢!1 Q% I# h4 i, o3 i# [/ Z" Z2 U
clear;clc;
e& ]9 u2 u2 t0 Y7 nT=[293.15 298.15 303.15 308.15 313.15];( K4 W3 A# s" `5 |, @9 `% c
c=[0.000378156 0.000425433 0.000562474 0.000684178 0.000981096];
7 Y2 d3 s1 S2 sTT=T';8 L1 T7 E5 S4 ^+ x5 u) h
cc=c';# u8 l; b/ ]2 w, k. N
x=[TT,cc];
' P4 [# d/ M/ m d4 N, y& U1 z/ T2 S2 BHm=33509;; h$ r) M1 B/ V2 y
Tm=408.65;
; y% n& N4 W# s/ ~a=0.3;
2 N. a& p* P5 T# O9 T* 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));' x9 y% u" Y" [8 C
q0=[0,0];
; a. B8 i' T, | |q=lsqcurvefit(fun,q0,x,zeros(size(x,1),1))
, ^; W0 g3 a6 Cplot(x(:,1),x(:,2),'r^');hold on;
% @( y' `! G8 B$ s* pezplot(@(x,c)fun(q,[x,c]),[290,314,0.0003,0.001]);) ]$ f* M& I/ M8 M) Y
h=ezplot(@(x,c)fun(q,[x,c]),[290,314,0.0003,0.001]);& {% X$ u7 r- |+ Y5 {( f: r- A
|
|