找回密码
 注册
关于网站域名变更的通知
查看: 602|回复: 1
打印 上一主题 下一主题

请问,lsqcurvefit拟合后如何显示相关系数R^2?

[复制链接]

该用户从未签到

跳转到指定楼层
1#
发表于 2020-4-28 14:15 | 只看该作者 回帖奖励 |倒序浏览 |阅读模式

EDA365欢迎您登录!

您需要 登录 才可以下载或查看,没有帐号?注册

x

- c. t/ }2 d" B) N2 L  I请问,lsqcurvefit拟合后如何显示相关系数R^2,感谢!
( p* n! a0 j  H- Zclear;clc;
  ]" `# E$ B) |  Q  ]# x- CT=[293.15        298.15        303.15        308.15        313.15];
$ Q- a; ?( A; B& D/ F: J# W6 m3 n/ Oc=[0.000378156        0.000425433        0.000562474        0.000684178        0.000981096];1 [4 \. ^( N# q/ R
TT=T';
# ]* X. S& E  {( xcc=c';) e, Q) B( d0 V8 }# W" q
x=[TT,cc];
/ q* c  a& z- H* m- H1 K3 Z: p& h$ qHm=33509;
, A, L) S( m6 u& ?6 e/ [9 E/ OTm=408.65;
0 g3 l" @4 O* x" ]8 Ra=0.3;
. }! O" _4 S/ J8 R3 z5 yfun=@(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));5 O9 \, i1 @0 K% x  p
q0=[0,0];6 y' L9 _2 D/ S0 ^* D) \
q=lsqcurvefit(fun,q0,x,zeros(size(x,1),1))
1 R' Z7 }8 G" X  U" O- X$ Yplot(x(:,1),x(:,2),'r^');hold on;% M9 @; d. n+ K9 a$ M, I
ezplot(@(x,c)fun(q,[x,c]),[290,314,0.0003,0.001]);
" D  Q! x4 O0 z" zh=ezplot(@(x,c)fun(q,[x,c]),[290,314,0.0003,0.001]);
' `6 g4 s8 }$ W, y! e' O5 {
  • TA的每日心情
    开心
    2019-11-19 15:19
  • 签到天数: 1 天

    [LV.1]初来乍到

    2#
    发表于 2020-4-28 15:08 | 只看该作者
    非线性拟合,R2没有意义的% k0 @1 s% t; g# j
    参考
    : z- Z* X( G8 Y; s& s  F  Ostatisticsbyjim.com/regression/r-squared-invalid-nonlinear-regression. k" {1 }" _) Y) S. m9 K
    blog.minitab.com/blog/adventures-in-statistics-2/why-is-there-no-r-squared-for-nonlinear-regression
    您需要登录后才可以回帖 登录 | 注册

    本版积分规则

    关闭

    推荐内容上一条 /1 下一条

    EDA365公众号

    关于我们|手机版|EDA365电子论坛网 ( 粤ICP备18020198号-1 )

    GMT+8, 2025-7-22 22:44 , Processed in 0.109375 second(s), 24 queries , Gzip On.

    深圳市墨知创新科技有限公司

    地址:深圳市南山区科技生态园2栋A座805 电话:19926409050

    快速回复 返回顶部 返回列表