| 查看: 1713 | 回复: 6 | |||
[求助]
数据拟合多项式
|
|
文献上有一色散参数：x=10^(-3)*[365 376 387 400 413 428 443 459 477 496 517 534 563 590 620 653 689 729 775 827 886 954 1033 1127 1240 1378 1550 1771]; y=[4.43 4.56 4.71 4.86 5.02 5.2 5.38 5.28 5.8 6.03 6.28 6.55 6.85 7.15 7.48 7.82 8.21 8.57 8.6 8.31 8.21 8.95 10.2 11.2 12.5 14 16 18.3];我用matlab拟合成10次多项式，结果为：1.0e+006 * 0.0072 -0.0687 0.2844 -0.6746 1.0139 -1.0070 0.6687 -0.2930 0.0812 -0.0128 0.0009 拟合数据和原始数据符合较好！但用所给系数作图却和原数据相差极大。不知错在哪里！请高手指教！谢谢！ |
回复此楼» 猜你喜欢
基金申请
已经有44人回复
-
CSC与新西兰维多利亚大学PhD奖学金项目
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有230人回复
-
新西兰Robinson研究所 招聘CSC公派访问人员
已经有0人回复
-
帮我的英语口语老师找学生
已经有0人回复
-
什么时候开奖？
已经有13人回复
-
散金币祈福
已经有94人回复
-
青基已中
已经有1人回复
-
散金币祈福
已经有113人回复
-
求助，如何提取ELK的rt-TDDFT在某一时刻的自旋密度分布和电子密度分布
已经有1人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 从excel读取数据多项式拟合出错
已经有3人回复
- 求助：数据拟合问题
已经有4人回复
- R软件如何做局部多项式拟合的预测
已经有10人回复
- R软件如何做多项式拟合线
已经有3人回复
- 【致H兄】数列题的万能多项式拟合求法
已经有28人回复
- 帮帮忙看看我的难题
已经有15人回复
- 【软件】一个科研数据处理软件和一个化学分子量，元素含量，配平等计算软件
已经有83人回复
- 【求助】怎样用chebyshev多项式拟合XRPD 背景?
已经有6人回复
- 【原创】多项式和非多项式曲线拟合的MATLAB实现视频教程
已经有10人回复
- 【一个数据拟合的问题，比较麻烦】
已经有17人回复
change0618
铁杆木虫 (著名写手)
方丈大师
- 计算强帖: 1
- 应助: 44 (小学生)
- 金币: 17724.5
- 红花: 17
- 帖子: 2413
- 在线: 546.7小时
- 虫号: 496517
- 注册: 2008-01-19
- 专业: 化学反应工程
2楼2012-03-24 22:01:38
change0618
铁杆木虫 (著名写手)
方丈大师
- 计算强帖: 1
- 应助: 44 (小学生)
- 金币: 17724.5
- 红花: 17
- 帖子: 2413
- 在线: 546.7小时
- 虫号: 496517
- 注册: 2008-01-19
- 专业: 化学反应工程
【答案】应助回帖
★ ★
dbb627: 金币+2, 感谢应助 2012-03-26 09:15:35
dbb627: 金币+2, 感谢应助 2012-03-26 09:15:35
|
x=10^(-3)*[365 376 387 400 413 428 443 459 477 496 517 534 563 590 620 653 689 729 775 827 886 954 1033 1127 1240 1378 1550 1771]; y=[4.43 4.56 4.71 4.86 5.02 5.2 5.38 5.28 5.8 6.03 6.28 6.55 6.85 7.15 7.48 7.82 8.21 8.57 8.6 8.31 8.21 8.95 10.2 11.2 12.5 14 16 18.3]; pp=spline(x,y) xi=10^(-3)*linspace(365,1771,1001); yi=ppval(pp,xi); plot(x,y,'o',xi,yi,'-') |
3楼2012-03-24 22:02:21
715211229
木虫 (正式写手)
- 应助: 52 (初中生)
- 金币: 3652.8
- 散金: 20
- 红花: 6
- 帖子: 339
- 在线: 249.2小时
- 虫号: 1357002
- 注册: 2011-07-27
- 性别: GG
- 专业: 传动机械学
4楼2012-03-25 12:10:07
5楼2012-10-23 17:32:20
715211229
木虫 (正式写手)
- 应助: 52 (初中生)
- 金币: 3652.8
- 散金: 20
- 红花: 6
- 帖子: 339
- 在线: 249.2小时
- 虫号: 1357002
- 注册: 2011-07-27
- 性别: GG
- 专业: 传动机械学
6楼2012-10-23 20:25:58
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
【答案】应助回帖
★ ★ ★
fegg7502: 金币+3, 应助指数+1, 3ks 2012-10-24 10:39:54
fegg7502: 金币+3, 应助指数+1, 3ks 2012-10-24 10:39:54
|
试下这个模型公式： y = p1*x^(p2+p3/x)+p4*x^p5+p6*exp(-0.5*((x-p7)/p8)^2) 均方差(RMSE): 0.0849994832064847 残差平方和(SSE): 0.202297540070345 相关系数(R): 0.999698462012221 相关系数之平方(R^2): 0.9993970149496 决定系数(DC): 0.999396761182227 卡方系数(Chi-Square): 0.0151137306236624 F统计(F-Statistic): 4736.33786978735 参数 最佳估算 ---------- ------------- p1 -6.60984611539542 p2 -19.31271541339 p3 17.948245176272 p4 11.0125462832269 p5 0.889759125054543 p6 5.22077463093514 p7 0.983832288151778 p8 0.158992885145476  v3.JPG |
7楼2012-10-24 09:25:35
