| 查看: 2150 | 回复: 18 | ||||
[求助]
小女子求助求助，有数据，需要求方程，已崩溃，求大神指点 已有1人参与
|
||||
|
各位大神，求助！急求！急求！急求！急求！最近要进行数据拟合，得到其方程，我最近看了可以使用1STOPT,就自己尝试了一下，完全是新手啊，我的数据还都是离散的三维数据，压根没有方程直接可以拿过来用然后求系数，也没有可以参考的方程。别人说可以用1STOPT的自动搜索功能，就用了，但问题有三个： 1，得出的方程很复杂，参数有10个左右，我就两个自变量，一个时间，一个温度，因变量是在不同温度和时间的组合下测得的指标含量（举个例子：蛋白质含量等），是不是自动搜索功能得到的方程都是这么复杂？能不能精简方程呢？ 2，自动搜索得到方程的R2为0.852048，但看三维图时候有太多的点都不在曲面上，怎么回事？为什么R2还会这么大？ 3，像我这种只有数据，没有方程，也没有可以供参考的方程，用什么办法可以得到方程，能不能用MATLAB做，怎么做呢？麻烦大神帮我指点指点，不胜感激，拜托拜托啦！ 小女子求助大神了，已经查了好几天资料了，没有头绪，快崩溃了，谢谢啦！！！ 数据如下：第一列x是时间，第二列y是温度，第三列z是指标 10 75 91.78 20 75 90.28 30 75 92.94 40 75 85.88 50 75 85.43 60 75 84.39 90 75 78.60 120 75 77.04 10 80 87.39 20 80 92.60 30 80 79.82 40 80 82.68 50 80 78.95 60 80 77.04 90 80 74.89 120 80 75.85 10 85 86.13 20 85 86.61 30 85 80.97 40 85 67.10 50 85 79.18 60 85 72.05 90 85 73.52 120 85 71.17 10 90 89.16 20 90 80.05 30 90 78.31 40 90 73.13 50 90 80.74 60 90 74.80 90 90 76.87 120 90 68.13 10 95 88.93 20 95 81.60 30 95 75.31 40 95 71.52 50 95 58.31 60 95 65.60 90 95 66.79 120 95 66.52 10 100 77.47 20 100 74.34 30 100 65.94 40 100 64.97 50 100 62.98 60 100 61.58 90 100 59.09 做出的一张图如下：附件中 问题是这么多的点都不在曲面上，为什么R2还这么高？  bw166h1411689_1444911705_636.jpg |
回复此楼
dingd
铁杆木虫 (职业作家)
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
5楼2015-10-21 10:59:20
lipenggg
专家顾问 (正式写手)
-
专家经验: +65 - 应助: 122 (高中生)
- 金币: 4841.3
- 红花: 43
- 帖子: 722
- 在线: 435.7小时
- 虫号: 3929201
- 注册: 2015-06-17
- 性别: GG
- 专业: 流体力学
- 管辖: 计算模拟
2楼2015-10-20 13:09:40
lipenggg
专家顾问 (正式写手)
-
专家经验: +65 - 应助: 122 (高中生)
- 金币: 4841.3
- 红花: 43
- 帖子: 722
- 在线: 435.7小时
- 虫号: 3929201
- 注册: 2015-06-17
- 性别: GG
- 专业: 流体力学
- 管辖: 计算模拟
3楼2015-10-20 13:26:45
lipenggg
专家顾问 (正式写手)
-
专家经验: +65 - 应助: 122 (高中生)
- 金币: 4841.3
- 红花: 43
- 帖子: 722
- 在线: 435.7小时
- 虫号: 3929201
- 注册: 2015-06-17
- 性别: GG
- 专业: 流体力学
- 管辖: 计算模拟
|
试试TableCurve 4.0，有公式库： z=a+b/x+c*ln(y)+d/(x^2)+e*(ln(y))^2+f*(ln(y))/x+g/(x^3)+h*(ln(y))^3+i*(ln(y))^2/x+j*(ln(y))/(x^2); z=a+b/x+c/y+d/(x^2)+e/(y^2)+f/(x*y)+g/(x^3)+h/(y^3)+i/(x*y^2)+j/(x^2*y); z=a+b/x+c*y+d/(x^2)+e*(y^2)+f*y/x+g/(x^3)+h*(y^3)+i*(y^2)/x+j*y/(x^2); z=a+b/x+c/(x^2)+d/(x^3)+e*y+f*(y^2)+g*(y^3); |
6楼2015-10-21 16:19:28
4楼2015-10-21 10:47:29
7楼2015-10-21 19:23:28
8楼2015-10-21 19:37:44
dingd
铁杆木虫 (职业作家)
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
9楼2015-10-21 20:30:57
10楼2015-10-22 15:59:49
