| 查看: 3211 | 回复: 8 | |||
[求助]
1stopt三维空间曲线拟合 已有1人参与
|
|
请教一下，已知三组数据，如何拟合出三维空间曲线，并得到他的参数方程，目前我用的1stopt，得到的是一个曲面，想问下到底如何能得到曲线拟合。如果有别的方法，如matlab，也请不吝赐教，谢谢！ 0 2.125 414897.577994617 1.5915 2.165 455102.831620648 3.1831 2.215 525696.470359034 4.7746 2.325 591448.219202163 6.3662 2.475 648298.441256297 7.9577 2.705 697003.786876068 9.5493 3.015 737917.784853322 11.1408 3.45 772047.838155198 12.7324 4.14 800063.260857717 14.3239 5.265 822044.806251477 15.9155 7.25 838474.029280583 17.507 11.035 848481.557824979 19.0986 14.79 852087.532478039 Function: z = p1+p2*x+p3*x^2+p4*x^3+p5*x^4+p6*x^5+p7/y+p8/y^2+p9/y^3+p10/y^4+p11/y^5 Algorithms: 准牛顿法(BFGS) + 通用全局优化法 Root of Mean Square Error (RMSE): 101.518407556763 Sum of Square Error (SSE): 133977.831947195 Correlation Coef. (R): 0.999999763438562 R-Square: 0.99999952687718 Determination Coef. (DC): 0.999999523682752 Parameters Name Parameter Value =============== =============== p1 126891.285393683 p2 12605.6079267492 p3 16813.6042796314 p4 -2990.12528858069 p5 190.825152062815 p6 -4.05893118643963 p7 1662713.16080946 p8 -1920920.79169137 p9 15965140.4020001 p10 -90888904.3994398 p11 118054087.577405  1.jpg |
回复此楼» 猜你喜欢
PRB投稿，Acknowledgment sent to author状态十天了
已经有3人回复
-
PVA溶解
已经有2人回复
-
物理学I论文润色/翻译怎么收费?
已经有148人回复
-
基金申请
已经有44人回复
-
CSC与新西兰维多利亚大学PhD奖学金项目
已经有0人回复
-
新西兰Robinson研究所 招聘CSC公派访问人员
已经有0人回复
-
帮我的英语口语老师找学生
已经有0人回复
-
什么时候开奖？
已经有13人回复
-
散金币祈福
已经有94人回复
-
青基已中
已经有1人回复
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
2楼2016-01-07 09:38:11
lipenggg
专家顾问 (正式写手)
-
专家经验: +65 - 应助: 122 (高中生)
- 金币: 4841.3
- 红花: 43
- 帖子: 722
- 在线: 435.7小时
- 虫号: 3929201
- 注册: 2015-06-17
- 性别: GG
- 专业: 流体力学
- 管辖: 计算模拟
3楼2016-01-07 09:49:24
4楼2016-01-07 18:56:47
5楼2016-01-07 18:58:56
|
上传的网页图片看不了，可以看看这个网页，http://www.ilovematlab.cn/thread-340717-1-1.html 里面的一个回复给出了三维空间曲线和参数方程，可是没有给方法 |
6楼2016-01-07 19:01:59
瀚海行舟: 回帖置顶 2016-01-16 19:50:17
|
感谢你在QQ上的热心帮忙，把你的回答贴出来，方便大家以后参考  。TiTle"三位函数参数方程拟合-编程"; Parameter a0,a1,a2,a3,b0,b1,b2,b3,c0,c1,c2,c3; Variable t,x[output],y[output],z[output]; StartProgram; var i: integer; begin for i := 0 to DataLength - 1 do begin x:=a0+a1*t+a2*t^2+a3*t^3; y:=b0+b1*t+b2*t^2+b3*t^3; z:=c0+c1*t+c2*t^2+c3*t^3; end; end; EndProgram; data; 1 0 2.125 414897.578 2 1.5915 2.165 455102.8316 3 3.1831 2.215 525696.4704 4 4.7746 2.325 591448.2192 5 6.3662 2.475 648298.4413 6 7.9577 2.705 697003.7869 7 9.5493 3.015 737917.7849 8 11.1408 3.45 772047.8382 9 12.7324 4.14 800063.2609 10 14.3239 5.265 822044.8063 11 15.9155 7.25 838474.0293 12 17.507 11.035 848481.5578 13 19.0986 14.79 852087.5325 迭代数: 16 计算用时(时:分:秒:微秒): 00:00:00:632 优化算法: 麦夸特法(Levenberg-Marquardt) + 通用全局优化法 计算结束原因: 达到收敛判定标准 均方差(RMSE): 3465.37346903185 残差平方和(SSE): 468343717.914924 相关系数(R): 0.999981637656814 相关系数之平方(R^2): 0.999963275650804 决定系数(DC): 0.999963275650804 F统计(F-Statistic): 299517.793313736 参数 最佳估算 -------------------- ------------- a0 -1.59079509167053 a1 1.59114273175162 a2 5.58198121310261E-5 a3 -2.21755578554155E-6 b0 0.888656538960089 b1 1.19767463959602 b2 -0.279347661876009 b3 0.0206738835030876 c0 334648.005832586 c1 70428.8906870194 c2 -1334.11539353067 c3 -80.0928008306951 |
7楼2016-01-16 19:49:17
8楼2016-10-09 09:01:04
9楼2016-10-09 09:23:57
投票:
10