查看: 143|回复: 4|关注: 0

[已解决] 关于降雨后 土壤水分 响应过程 的拟合参数求值

[复制链接]

新手

12 麦片

财富积分


050


1

主题

4

帖子

0

最佳答案
本帖最后由 xidaly 于 2019-3-16 00:13 编辑

已知,降雨后土壤水分的响应变化符合一个指数方程:

突变更准确一些

突变更准确一些

y是不同时间的土壤含水量
y0是初始(降水发生前)的土壤含水量
t为时间(天)
求a和b的值
恳求了 :handshake
祝各位大神诸事皆宜

部分数据如下:
降水y0(雨前)y1(雨后1天)y2(雨后2天)y3(雨后3天)
7.5
0.156075
0.155081
0.153531
0.151986
13.6
0.143239
0.14179
0.139307
0.136419
25.9
0.098147
0.157719
0.153984
0.147144
43.6
0.110921
0.186931
0.168722
0.160657
12.5
0.032922
0.033904
0.034842
0.035059



论坛优秀回答者

中级

893 麦片

财富积分


5001500


0

主题

2279

帖子

194

最佳答案
  • 关注者: 127
发表于 7 天前 | 显示全部楼层 |此回复为最佳答案
用1stOpt求解试试,数据太少

  1. VarConstant y0=[0.156075,0.143239,0.098147,0.110921,0.032922];
  2. Variable t,y;
  3. Function y=y0+a*t*exp(-b*t);
  4. Data;
  5. 1        0.155081
  6. 2        0.153531
  7. 3        0.151986
  8. Data;
  9. 1        0.14179
  10. 2        0.139307
  11. 3        0.136419
  12. Data;
  13. 1        0.157719
  14. 2        0.153984
  15. 3        0.147144
  16. Data;
  17. 1        0.186931
  18. 2        0.168722
  19. 3        0.160657
  20. Data;
  21. 1        0.033904
  22. 2        0.034842
  23. 3        0.035059
复制代码


均方差(RMSE): 0.0296618176571017
残差平方和(SSR): 0.0131973514008472
相关系数(R): 0.594439988079594
相关系数之平方(R^2): 0.353358899428068
修正R平方(Adj. R^2): -100
确定系数(DC): -73.4028059003556
F统计(F-Statistic): 0.992133966384286

参数                  最佳估算
--------------------        -------------
a        0.0561173518324212
b        0.774038044657029

新手

12 麦片

财富积分


050


1

主题

4

帖子

0

最佳答案
 楼主| 发表于 7 天前 | 显示全部楼层
shihe 发表于 2019-3-16 11:25
用1stOpt求解试试,数据太少

谢谢您,第一次接触这个软件,感觉很神奇!
刚试用了下。
发现结果跟您运算有差异,所以想了解您关于此软件是如何设置优化的呢?

附我的运算结果:
迭代数: 16
计算用时(时:分:秒:微秒): 00:00:00:470
优化算法: 麦夸特法(Levenberg-Marquardt) + 通用全局优化法
计算结束原因: 达到收敛判断标准
均方差(RMSE): 0.0492304891754422
残差平方和(SSE): 0.0363546159668
相关系数(R): 0.0732616551629511
相关系数之平方(R^2): 0.00536727011721516
决定系数(DC): 0.00536727011721524
卡方系数(Chi-Square): 0.196529389234547
F统计(F-Statistic): 0.0701510310564784

参数        最佳估算
----------        -------------
y0                 0.123377595195455
a                  0.0409168397948926
b                  1.25132017568855

论坛优秀回答者

中级

893 麦片

财富积分


5001500


0

主题

2279

帖子

194

最佳答案
  • 关注者: 127
发表于 7 天前 | 显示全部楼层
版本问题吧,1.5破解版?y0都当作参数了,明显不对。

新手

12 麦片

财富积分


050


1

主题

4

帖子

0

最佳答案
 楼主| 发表于 7 天前 | 显示全部楼层
shihe 发表于 2019-3-16 20:14
版本问题吧,1.5破解版?y0都当作参数了,明显不对。

已解决,不胜感激:handshake
您需要登录后才可以回帖 登录 | 注册

本版积分规则

关闭

站长推荐上一条 /3 下一条

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