EDA365电子论坛网
标题:
请问，lsqcurvefit拟合后如何显示相关系数R^2？
[打印本页]
作者:
limopcb
时间:
2020-4-28 14:15
标题:
请问，lsqcurvefit拟合后如何显示相关系数R^2？
) q0 y: K3 E( L" z9 X b
请问，lsqcurvefit拟合后如何显示相关系数R^2，感谢！
6 y4 j& w5 o& V7 R* V7 H
clear;clc;
1 F. L' b l) b: q/ B$ C
T=[293.15 298.15 303.15 308.15 313.15];
% L" h1 O: `( }; h
c=[0.000378156 0.000425433 0.000562474 0.000684178 0.000981096];
1 ?6 A0 ~% u+ b
TT=T';
1 }+ ]7 O. a, }
cc=c';
2 V0 x0 M7 i9 O: V
x=[TT,cc];
1 v# `" V0 Z) {" o
Hm=33509;
/ b' e# [5 B- L, w C
Tm=408.65;
( ^8 H9 W! Y$ ?& l( w# |( `
a=0.3;
. x! A* J# h* I9 S3 @7 w2 I5 T$ v5 ?
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" N% z5 o& G; Z
q0=[0,0];
* O+ _, T% y2 C/ z
q=lsqcurvefit(fun,q0,x,zeros(size(x,1),1))
) D+ a. w* H4 Y+ S4 o- d
plot(x(:,1),x(:,2),'r^');hold on;
7 ^! o# z) Y' S
ezplot(@(x,c)fun(q,[x,c]),[290,314,0.0003,0.001]);
( e9 X6 C; g9 z k9 o$ y: ]! f
h=ezplot(@(x,c)fun(q,[x,c]),[290,314,0.0003,0.001]);
+ X3 c% X+ v' k% T( y) {7 G7 l+ H
作者:
Allevi
时间:
2020-4-28 15:08
非线性拟合，R2没有意义的
! d( c6 n- J. ]+ G. u1 t
参考
8 u2 e2 i/ f! u2 {: a
statisticsbyjim.com/regression/r-squared-invalid-nonlinear-regression
4 i, S( H; L3 q3 O7 j& d1 w5 f
blog.minitab.com/blog/adventures-in-statistics-2/why-is-there-no-r-squared-for-nonlinear-regression
欢迎光临 EDA365电子论坛网 (https://bbs.eda365.com/)
Powered by Discuz! X3.2