| 查看: 1517 | 回复: 5 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[交流]
拟合非线性方程系数的精度问题
|
|||
|
本人想拟合一个非线性方程的系数（4个未知数），精度要求很高，期望要到0.01mm。用MATLAB中优化函数lsqnonlin、fmincon、fminunc等函数和遗传算法计算时还得不到满意的精度.不知道上述提到的计算方法大家在实际工程应用中都能达到多少精度啊！并且为了能达到我想要的精度还有什么好的方法不？ 目标函数为求和sum(i=1,6,z(i)),z(i)=m(i)^2+a*m(i)+b*m(i)*sin(w(i))+c*sin(w(i))+d*m(i)*cos(w(i))+e*cos(w(i))+f,求它的最小值，接近与0.其中a=-2*x(1); b=2*(x(4)-x(2));c=2*(x(3)*x(2)-x(4)*x(1));d=2*(x(1)-x(3)); e=2*(x(3)*x(1)+x(2)*x(4)-x(1)^2-x(2)^2); f=2*(x(1)^2+x(2)^2)+x(3)^2+x(4)^2-2*x(2)*x(4)-2*x(1)*x(3)-(R+r)^2;R=50;r=1.5 w=n*pi/180 m=[41.76554 21.16515 12.83882 10.98004 14.71395 24.88945] n=[0 2.5 5 7.5 10 12.5] x(1)~x(4)为所求未知数，理想解为x(1)=422.443 ;x(2)=20 ;x(3)=62.918 ;x(4)=46.95,结果精度要求0.01 麻烦各位大侠了！ [ Last edited by myb19880127 on 2012-3-31 at 10:43 ] |
回复此楼» 猜你喜欢
求国际会议网站
已经有1人回复
-
求取一些关于纳米材料和纳米技术相关的英文PPT。
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有56人回复
-
【复旦大学】二维材料方向招收2026年博士研究生1名
已经有0人回复
-
北京纳米能源与系统研究所 王中林院士/曹南颖研究员课题组2026级硕/博/博后招生
已经有10人回复
-
荷兰Utrecht University超快太赫兹光谱王海教授课题招收2026 CSC博士生
已经有19人回复
-
反铁磁体中的磁性切换：两种不同的机制已成功可视化
已经有0人回复
-
26申博推荐：南京航空航天大学国际前沿科学研究院光学方向招收博士生！
已经有0人回复
-
求标准粉末衍射卡号 ICDD 01-076-1802
已经有0人回复
» 本主题相关价值贴推荐，对您同样有帮助:
在matlab中拟合曲面的表达式
已经有29人回复- 求助：非线性方程系数回归问题？
已经有10人回复
- 求非线性方程中的参数问题
已经有13人回复
- 求助关于线性拟合的
已经有1人回复
» 抢金币啦！回帖就可以得到:
-
Ei期刊青年编委招募（工程设计方向）
+1/380
-
天津大学精密测试技术及仪器全国重点实验室柔性电子技术实验室专任副研究员招聘启事
+1/76
-
以色列理工-生物质塑料等催化转化及流体力学方向---全奖博士研究生和科研助理
+2/64
-
南京大学自旋全国重团队陆显扬课题组招聘博士后
+2/44
-
中国科学院理化所微纳材料与技术前沿交叉研究中心诚聘英才加盟
+1/42
-
2026年博士招生--北京理工大学交叉学科（航空，力学，能动，计算机等方向）
+1/35
-
北京科技大学鲁启鹏招收2026年博士生1名
+1/31
-
能够检测核磁、LCMS的机构或个人请跟我联系
+1/15
-
长春理工大学和西安工业大学主动光电探测成像技术重点实验室招收博士生
+1/12
-
荷兰奈梅亨大学招收2026 CSC博士: 非线性控制与神经调控
+1/7
-
SCI辅助 3 纯 4 自 7 己 0 写 9 非 ③ 中 1 介 ⑦ 优0 惠
+1/7
-
海南大学国家高层次人才团队2026年博士招生
+1/6
-
哈尔滨工业大学（深圳）赵怡潞课题组诚招博士后
+1/6
-
南京航空航天大学航天学院黄护林教授课题组博士研究生招生（工程热物理专业）
+1/6
-
诚招博士后及研究人员
+1/5
-
招收26年秋季入学博士生（北科大高精尖学院 力学超材料/机器学习/增材制造相关方向）
+1/5
-
国家纳米科学中心鄢勇课题组26年博士招生
+1/5
-
SCI辅 助 3 纯 4 自 7 己 0 写 9 非 ③ 中 1 介 ⑦ 优0 惠
+1/5
-
澳大利亚南昆士兰大学（UniSQ）量子点课题组 招收CSC全奖博士生
+1/3
-
海南大学化学院 招聘 材料与电化学方向——研究助理，博士（2026年入学）
+1/1
★ ★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
dbb627: 金币+2, 谢谢应助！ 2012-04-01 21:19:08
小木虫: 金币+0.5, 给个红包，谢谢回帖
dbb627: 金币+2, 谢谢应助！ 2012-04-01 21:19:08
|
你的目标函数z应该取绝对值吧？否则最小值远远小于0。 用1stOpt求解： 目标函数值(最小): 0.288024770736229 x1: 422.733761186019 x2: 19.3667780059412 x3: 62.9959407161319 x4: 46.9203589652975 将你的理想解(x(1)=422.443 ;x(2)=20 ;x(3)=62.918 ;x(4)=46.95)代入验证，目标函数值(最小)=5.506768407 从理论上看，结果已经好于你的理想值了。 |
6楼2012-03-31 10:49:49
2楼2012-03-30 16:53:41
|
目标函数为z=m(i)^2+a*m(i)+b*m(i)*sin(w(i))+c*sin(w(i))+d*m(i)*cos(w(i))+e*cos(w(i))+f,求它的最小值，接近与0.其中a=-2*x(1); b=2*(x(4)-x(2));c=2*(x(3)*x(2)-x(4)*x(1));d=2*(x(1)-x(3)); e=2*(x(3)*x(1)+x(2)*x(4)-x(1)^2-x(2)^2); f=2*(x(1)^2+x(2)^2)+x(3)^2+x(4)^2-2*x(2)*x(4)-2*x(1)*x(3)-(R+r)^2;R=50;r=1.5 w=n*pi/180 m=[41.76554 21.16515 12.83882 10.98004 14.71395 24.88945] n=[0 2.5 5 7.5 10 12.5] x(1)~x(4)为所求未知数，理想解为x(1)=422.443 ;x(2)=20 ;x(3)=62.918 ;x(4)=46.95,结果精度要求0.01 麻烦各位大侠了！ [ Last edited by myb19880127 on 2012-3-31 at 10:38 ] |
3楼2012-03-31 09:57:53
4楼2012-03-31 10:23:11