| 查看: 736 | 回复: 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 |
回复此楼» 猜你喜欢
应该是93bebmhtak前后十一个字符比较关键
已经有32人回复
-
2019年青年基金涵评意见，大家看看几个A，几个B?
已经有5人回复
-
职称评审，求友友推荐见刊最快的期刊
已经有4人回复
-
奇怪，两个人的filecode固定段从头到尾一模一样
已经有10人回复
-
综述论文作为代表作会不会影响评审专家的印象分？
已经有13人回复
-
是这周出结果还是下周出结果？
已经有3人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有9人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有4人回复
-
帮忙看看fileCode
已经有13人回复
-
为什么网上很多人说本周 12号出结果
已经有7人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 使用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人回复
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
乙凡良
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 510.5
- 散金: 55
- 帖子: 63
- 在线: 45.8小时
- 虫号: 1478502
- 注册: 2011-11-06
- 专业: 可再生与替代能源利用中的
5楼2015-05-14 15:55:34
