| 
 | 
	
    
 
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册  
 
x
 
% u1 G( q) O7 }! P7 _ 
请问,lsqcurvefit拟合后如何显示相关系数R^2,感谢! 
9 w+ T; @7 k) p9 V2 N9 g" m& Uclear;clc; 
: F) c$ l. @1 c+ V% f3 p/ `T=[293.15        298.15        303.15        308.15        313.15];# \/ B; B+ _! p+ U& { 
c=[0.000378156        0.000425433        0.000562474        0.000684178        0.000981096];8 m+ m0 |  ~: a& M! j% A% J$ F6 T 
TT=T';' ]7 Z9 j/ b# t& K8 F 
cc=c';4 u0 r$ ^& o5 ~% J% B  {) F; ? 
x=[TT,cc];1 f0 h1 N' J0 D1 k* | 
Hm=33509; 
! W" [" I" S) ~  N4 DTm=408.65; 
+ f5 t# r" p( Xa=0.3;& Z  u$ f7 z  i) a9 }+ | 
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));' e! e: r- l1 _4 p 
q0=[0,0]; 
$ |/ i+ F4 y' g' l7 _: i- tq=lsqcurvefit(fun,q0,x,zeros(size(x,1),1)) 
4 X: O+ V8 c/ M& ]5 D1 c& h$ Splot(x(:,1),x(:,2),'r^');hold on; 
4 H# Z7 q7 K  I2 Kezplot(@(x,c)fun(q,[x,c]),[290,314,0.0003,0.001]); 
, g" _0 k, S3 P  H5 l9 \h=ezplot(@(x,c)fun(q,[x,c]),[290,314,0.0003,0.001]);* x0 ?, X/ N4 N6 k! n 
 |   
 
 
 
 |