| 查看: 744 | 回复: 4 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
乙凡良新虫 (小有名气)
|
[求助]
matlab拟合数据？？？？ 已有1人参与
|
||
|
目标函数如下图片所示： x，y对应值如下： x y 0.003220 5.560760 0.003420 5.561897 0.003620 5.585926 0.003820 5.667371 0.004020 5.612988 0.004220 5.670119 0.004420 5.703140 0.004620 5.747208 0.004820 5.748986 0.005020 5.780545 0.005220 5.860896 0.005420 5.855221 0.005620 5.862226 0.005820 5.935035 0.006020 5.907555 0.006220 5.961461 0.006420 5.991398 0.006620 5.991484 0.006820 5.996356 0.007020 5.989982 0.007220 6.028335 0.007420 6.020197 0.007620 6.017697 0.007820 6.134659 0.008020 6.135323 0.008220 6.110176 0.008420 6.159997 0.008620 6.135110 0.008820 6.198990 0.009020 6.178628 0.009220 6.271552 0.009420 6.257719 0.009620 6.305726 0.009820 6.305890 0.010020 6.309602 0.010220 6.286979 拟合参数，a，r, c, B的值，其中a的取值范围（0.2,0.6），r(0.3,0.5),c(-0.001,0),b(4,8)，用了matlab里的拟合工具箱，说方程太复杂拟合不出来，求matlab拟合该方程的程序？？？  12.png |
回复此楼» 猜你喜欢
欢迎发来filecode的Mz6后的代码验证其规律
已经有46人回复
-
各位道友，我要去昆明玩几天，回来见。
已经有4人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有16人回复
-
奇怪，两个人的filecode固定段从头到尾一模一样
已经有11人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有4人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有8人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有37人回复
-
好奇怪的filecode
已经有7人回复
-
filecode
已经有9人回复
-
请教兼职经验
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 使用MATLAB拟合动力学参数问题
已经有18人回复
- matlab 非线性拟合
已经有16人回复
- matlab曲线拟合及修正
已经有10人回复
- 数据关系的拟合
已经有9人回复
- matlab如何在一个图里拟合多条直线
已经有5人回复
- 关于matlab两个自变量、一个因变量曲线拟合问题
已经有3人回复
- matlab线性回归拟合数据
已经有4人回复
- matlab 拟合方法
已经有4人回复
- matlab 拟合求助
已经有8人回复
- 关于matlab拟合
已经有4人回复
- 编写matlab程序：y=aln（x）+b，求a、b及拟合的R？
已经有15人回复
- 通过实验数据拟合，求解公式中的参数
已经有14人回复
- matlab数据拟合问题求助
已经有7人回复
- 用MATLAB或者ORIGIN怎么拟合复数数据
已经有10人回复
- 坐等----MAtlab中如何有三维数据拟合曲面，并得到曲面方程？？？？？？？
已经有10人回复
- matlab拟合求值
已经有17人回复
- 【求助】用matlab最优化方法进行参数拟合
已经有17人回复
- 【求助】使用MATLAB拟合出按照公式的系数
已经有19人回复
乙凡良
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 510.5
- 散金: 55
- 帖子: 63
- 在线: 45.8小时
- 虫号: 1478502
- 注册: 2011-11-06
- 专业: 可再生与替代能源利用中的
5楼2015-05-14 15:55:34
feixiaolin
荣誉版主 (文坛精英)
-
专家经验: +518 - 应助: 942 (博后)
- 贵宾: 1.275
- 金币: 3881
- 散金: 58785
- 红花: 532
- 沙发: 11
- 帖子: 24215
- 在线: 2601.8小时
- 虫号: 2139575
- 注册: 2012-11-21
- 专业: 光学信息获取与处理
- 管辖: 数学
2楼2015-05-13 17:02:45
dingd
铁杆木虫 (职业作家)
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
【答案】应助回帖
感谢参与，应助指数 +1
|
没难度啊： 只是拟合公式是过拟合的，解不唯一。 1： 均方差(RMSE):0.0299241203402951 残差平方和(SSE):0.0322363072130566 相关系数(R): 0.991464779695399 相关系数之平方(R^2): 0.983002409376447 确定系数(DC): 0.983002409376447 卡方系数(Chi-Square): 0.00269015521761369 F统计(F-Statistic): 627.539920386496 参数 最佳估算 ---------- ------------- a 0.354996057885225 r 0.443615101337285 c -0.000999999999999999 b 5.1138108023883 2： 均方差(RMSE):0.0299241203402951 残差平方和(SSE):0.0322363072130566 相关系数(R): 0.991464779695399 相关系数之平方(R^2): 0.983002409376447 确定系数(DC): 0.983002409376447 卡方系数(Chi-Square): 0.00269015521936039 F统计(F-Statistic): 627.539920386496 参数 最佳估算 ---------- ------------- a 0.400117194038538 r 0.5 c -0.00099999999999999 b 4.29782573071645 其实原公式可以化简， 原公式：Function y=a*(1/r*ln((x-c)*a/0.0000147)+b); 可以写成： y=a*1/r*ln((x-c)*a/0.0000147)+b*a; 再写成： y=a*1/r*(ln(x-c)+ln(a/0.0000147))+b*a; 再写成： y=a*1/r*ln(x-c)+a*1/r*ln(a/0.0000147))+b*a; 最好简化为： y=r1*ln(x-c)+a1; 其中：r1=a*1/r, a1=a*1/r*ln(a/0.0000147))+b*a 解唯一稳定： 均方差(RMSE):0.0299241203402951 残差平方和(SSE):0.0322363072130566 相关系数(R): 0.9914647796954 相关系数之平方(R^2): 0.983002409376447 确定系数(DC): 0.983002409376447 卡方系数(Chi-Square): 0.0026901552191239 F统计(F-Statistic): 352.991205217404 参数 最佳估算 ---------- ------------- c -0.000999999999999999 r1 0.800234387454693 a1 9.89135968211093 |
3楼2015-05-13 22:04:49
乙凡良
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 510.5
- 散金: 55
- 帖子: 63
- 在线: 45.8小时
- 虫号: 1478502
- 注册: 2011-11-06
- 专业: 可再生与替代能源利用中的
4楼2015-05-14 15:50:44
