[已解决] 跪求1stopt代跑,代数方程组参数拟合问题,跪求,感谢

[复制链接]
化工坑中人 发表于 2021-1-14 16:38:24
用matlab求出的置信区间特别大,只能麻烦论坛里的大神用 1stopt代跑一下,感激不尽
代数方程组求参数k1,k2,k3

x1=1.0e-05 *[0.0045    0.0600    0.1525    0.3396    0.5740];
x2=[-0.0046   -0.0046   -0.0058   -0.0049   -0.0051];
x3=1.0e-05 *[0.0023    0.0317    0.0818    0.1944    0.3517];

y1=[0.0020    0.0963    1.0882    2.1391    3.0805  ];
y2=[0.0071    0.5481    1.2768    1.9309    2.4807  ];

lb=[0 0 0]; ub=[+inf +inf +inf]; k0=[100000 10 100000];

y1=k(1)*x1+k(3)*x3;
y2=k(1)*x1-k(2)*x2;

最佳答案


shihe 发表于 2021-1-14 23:04:11
参考下:
均方差(RMSE): 0.217145961509
残差平方和(SSR): 0.471523685996682
相关系数(R): 0.97608884075042
相关系数之平方(R^2): 0.952749425037498
修正R平方(Adj. R^2): 1.18867892327847
确定系数(DC): 0.952572882212911
F统计(F-Statistic): 33.0783234436198

参数                  最佳估算
--------------------        -------------
k1        408415.139351425
k2        67.4393102054562
k3        259151.422070523
回复此楼

4 条回复


shihe 发表于 2021-1-14 23:04:11
参考下:
均方差(RMSE): 0.217145961509
残差平方和(SSR): 0.471523685996682
相关系数(R): 0.97608884075042
相关系数之平方(R^2): 0.952749425037498
修正R平方(Adj. R^2): 1.18867892327847
确定系数(DC): 0.952572882212911
F统计(F-Statistic): 33.0783234436198

参数                  最佳估算
--------------------        -------------
k1        408415.139351425
k2        67.4393102054562
k3        259151.422070523
回复此楼

化工坑中人 发表于 2021-1-19 16:12:47
shihe 发表于 2021-1-14 23:04
参考下:
均方差(RMSE): 0.217145961509
残差平方和(SSR): 0.471523685996682

大神,求下代码,因为我发现我还有很多相同的问题需要处理,本人1stopt比较小白,求代码,感谢感谢

shihe 发表于 2021-1-19 19:22:08
很简单的
  1. Parameter k(3)>0;
  2. SharedModel;
  3. Variable x(3),y(2);
  4. Function
  5. y1=k1*x1+k3*x3;
  6. y2=k1*x1-k2*x2;
  7. Data;
  8. x1=1.0e-05*[0.0045    0.0600    0.1525    0.3396    0.5740];
  9. x2=[-0.0046   -0.0046   -0.0058   -0.0049   -0.0051];
  10. x3=1.0e-05*[0.0023    0.0317    0.0818    0.1944    0.3517];
  11. y1=[0.0020    0.0963    1.0882    2.1391    3.0805];
  12. y2=[0.0071    0.5481    1.2768    1.9309    2.4807];
复制代码

化工坑中人 发表于 2021-1-20 00:58:08

感谢感谢,我用你的代码又拟合了问题,不过拟合效果太差了,不知道是我数据的原因还是我的版本太低,劳烦你再帮我跑一下吧,感激不尽,谢谢
Parameter k(3)>0;
SharedModel;
Variable x(3),y(2);
Function
y1=k1*x1+k3*x3;
y2=k1*x1-k2*x2;
Data;
x1=[4.18653754388550e-16        7.67213235342095e-16        1.18230989742350e-15        1.49564708515217e-15        1.88831063901827e-15        7.21448218553861e-16        1.03152390812815e-15        1.33737191124301e-15        1.74941751407285e-15        2.18579585000466e-15        8.19436210215603e-16        1.12336594927635e-15        1.44308924838643e-15        1.85684403295872e-15        2.73979125460571e-15        3.62941403454823e-16        5.44711707676504e-16        7.07258537897752e-16        1.05504854630030e-15        2.14404880388626e-15];
x2=[-1.09763659492295e-05        -7.63540479112841e-06        -4.18228552495061e-06        -1.94884053501297e-06        -1.27421093256162e-06        -2.89183714557005e-06        1.06486102128333e-06        2.51521952710162e-06        1.80088172734856e-06        5.43983305443430e-07        1.15669773037338e-06        3.43749349619594e-06        1.37169872490826e-06        1.09004014721339e-06        2.38892533332713e-07        6.35474248951947e-06        7.56163939155184e-06        4.25349252104060e-06        2.31560576722522e-06        5.03511450434112e-07];
x3=[1.77412085829895e-16        4.06981811895747e-16        8.32577523305119e-16        1.31959756417481e-15        2.18261128836886e-15        4.58778017971032e-16        9.00702390759861e-16        1.58766887357866e-15        3.02253826847468e-15        5.33122866988718e-15        8.52970721966541e-16        1.76929468664514e-15        3.26620080671926e-15        6.43720325532972e-15        1.91457857286791e-14        5.09951388380118e-16        1.22808819051568e-15        2.24832343952675e-15        6.17976895233189e-15        4.40120833804113e-14];
y1=[ 0.0228    0.0345    0.0423    0.0477    0.0517  0.0152    0.0203    0.0273    0.0336    0.0390  0.0092    0.0114    0.0139    0.0161    0.0179  0.0040    0.0048    0.0064    0.0080    0.0093];
y2=[0.0024    0.0070    0.0149    0.0216    0.0271 0.0087    0.0167    0.0226    0.0271    0.0305   0.0115    0.0182    0.0236    0.0277    0.0310  0.0116    0.0185    0.0239    0.0281    0.0313];
您需要登录后才可以回帖 登录 | 注册

本版积分规则

相关帖子
热门教程
站长推荐
快速回复 返回顶部 返回列表