|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
2 ^+ ]% p* v7 s1 f; t请问,lsqcurvefit拟合后如何显示相关系数R^2,感谢!( n; s9 d9 B1 R" e2 X
clear;clc;
% @$ O% o' h3 \9 I, u* {T=[293.15 298.15 303.15 308.15 313.15];9 I! T! {& U. ~) ?) { m+ r
c=[0.000378156 0.000425433 0.000562474 0.000684178 0.000981096];9 {, x8 u1 v) v9 d }/ i& S- g
TT=T';
7 ]$ K+ ?2 Q3 D3 L. C. |2 Ucc=c';
9 R }0 E8 {9 L7 k. [x=[TT,cc];. r1 a, [2 o! m2 p; B' {$ y
Hm=33509;
3 G% c. M& y4 u" M6 B; V$ r$ Z4 I0 _* uTm=408.65;" x1 Q+ Y5 u7 `5 i/ u$ `' l
a=0.3;
! P" s) M( ?; }- I9 a( }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));
# s8 w m( s) t$ C h! c7 ? ?q0=[0,0];3 Q3 @7 b% G' I. H$ d9 I# {+ M3 {
q=lsqcurvefit(fun,q0,x,zeros(size(x,1),1))
( Y' |0 H$ X$ z/ P; Jplot(x(:,1),x(:,2),'r^');hold on;
4 @( m0 Q0 Y) w9 wezplot(@(x,c)fun(q,[x,c]),[290,314,0.0003,0.001]);/ _$ i6 q5 [" t
h=ezplot(@(x,c)fun(q,[x,c]),[290,314,0.0003,0.001]);4 N. B: g0 h6 q% u
|
|