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

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

[复制链接]

该用户从未签到

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

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
  • TA的每日心情
    开心
    2019-11-19 15:19
  • 签到天数: 1 天

    [LV.1]初来乍到

    2#
    发表于 2020-4-28 15:08 | 只看该作者
    非线性拟合,R2没有意义的
    * g, k& |, Q" \参考
    9 z- t1 u- \) y7 J9 ?! Rstatisticsbyjim.com/regression/r-squared-invalid-nonlinear-regression
    5 I3 n9 m* r/ mblog.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-8-18 16:27 , Processed in 0.109375 second(s), 23 queries , Gzip On.

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

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

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