| 查看: 1679 | 回复: 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 ] |
» 猜你喜欢
【修正版】单环用磁约束低速电子实验方案(简化版)
已经有0人回复
桂林理工大学物理学专业招收调剂,还有三个名额!!!
已经有22人回复
物理学I论文润色/翻译怎么收费?
已经有289人回复
考博自荐
已经有4人回复
山东大学第二批博士研究生招生
已经有0人回复
中国科学院东莞材料科学与技术研究所-2026年博士招生-吴昊研究员-磁学与自旋电子学
已经有0人回复
《电磁学》教材推荐
已经有1人回复
【急招】合肥工大核聚变材料计算方向2026级工程博士生
已经有4人回复
大豆异黄酮分离
已经有0人回复
湖南大学材料学院急招2026年博士生,临时增加一名博士联培指标
已经有10人回复
天津理工大学晶体材料全国重点实验室刘红军教授课题组招收博士生1-2名
已经有0人回复
» 本主题相关价值贴推荐,对您同样有帮助:
在matlab中拟合曲面的表达式
已经有29人回复
求助:非线性方程系数回归问题?
已经有10人回复
求非线性方程中的参数问题
已经有13人回复
求助关于线性拟合的
已经有1人回复
» 抢金币啦!回帖就可以得到:
河南大学化学与分子科学学院博士招生
+2/418
北京科技大学能源与环境工程学院王庆功教授团队,招收全日制博士生1-2名,5月21日前
+1/91
南京农业大学人工智能学院 农业智能传感器与检测技术实验室招收2026级博士研究生
+1/34
邀请函 | 华南理工大学集成电路技术前沿论坛(含电子信息)诚邀您参加!
+1/14
求siesta、Transiesta教程资料
+1/10
杨亚研究员课题组(广西大学-纳米能源研究中心)招收2026年物理学专业学术型博士
+1/8
邀请函 | 华南理工大学集成电路技术前沿论坛(含电子信息)诚邀您参加!
+1/8
211高校2026年入学博士名额
+1/6
南京邮电大学邓超教授课题组招收2026博士生(5月18日前有效)
+1/6
知名外资仪器厂家急招Application Scientist(2名)
+1/6
华南师范大学(211)- 生物光子学研究院招收生物学、物理学、神经科学等专业博士
+1/5
福建师范大学优青团队招收申请-考核制博士1名
+1/5
哈尔滨工业大学(深圳)——急招能动/材料/机械——26年9月份博士/博后/RA
+1/4
中国科学院半导体研究所半导体量子材料与器件课题组多岗位招聘
+1/4
天津大学张小岩教授团队招收教授/副教授/副研究员岗位(有编制和名额)
+1/4
天津理工大学国家杰青鲁统部团队第二批博士生招生简章-名额有限
+1/4
广东工业大学聚合物与胶体国家优青团队青年教师招聘启事(高分子合成及自组装方向)
+1/4
华侨大学(厦门)——急招化工/材料/化学/能动——26年9月份博士
+1/3
请教:国际会议报账的问题
+1/2
清华大学化学系赵亮课题组诚聘博士后
+1/1
|
目标函数为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
2楼2012-03-30 16:53:41
4楼2012-03-31 10:23:11
5楼2012-03-31 10:36:37
★ ★ ★
小木虫: 金币+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












回复此楼
精华III: