24小时热门版块排行榜    

Znn3bq.jpeg
查看: 1631  |  回复: 5

myb19880127

新虫 (小有名气)


[交流] 拟合非线性方程系数的精度问题

本人想拟合一个非线性方程的系数(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 ]
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

» 抢金币啦!回帖就可以得到:

查看全部散金贴

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖置顶 ( 共有1个 )

myb19880127

新虫 (小有名气)


引用回帖:
2楼: Originally posted by dingd at 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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

dingd

铁杆木虫 (职业作家)


★ ★
小木虫: 金币+0.5, 给个红包,谢谢回帖
zhangguangping: 金币+1, 谢谢参与! 2012-03-30 23:51:29
精度除了算法设置,拟合的数据质量和拟合公式的正确与否也有很大关系。把公式和数据一起放上来让大家一起看看来解决。
2楼2012-03-30 16:53:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)



小木虫: 金币+0.5, 给个红包,谢谢回帖
目标函数“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”应该是求和吧?sum(i=1,6,z(i))

数据有n=[...],公式里没出现n但出现了w(i),确认公式无误吗?
4楼2012-03-31 10:23:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

myb19880127

新虫 (小有名气)


引用回帖:
4楼: Originally posted by dingd at 2012-03-31 10:23:11:
目标函数“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”应该是求和吧?sum(i=1,6,z(i))

数据有n=,公式里没出现n但出现了w(i),确认公式无误吗?

哦,不好意思,w=n*pi/180,目的是将角度制改为弧度制。目标函数是求和的!就是使得目标函数尽量接近于0

[ Last edited by myb19880127 on 2012-3-31 at 10:40 ]
5楼2012-03-31 10:36:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)


★ ★ ★
小木虫: 金币+0.5, 给个红包,谢谢回帖
dbb627: 金币+2, 谢谢应助! 2012-04-01 21:19:08
你的目标函数z应该取绝对值吧?否则最小值远远小于0。

用1stOpt求解:
CODE:
Parameter x(4)
Constant r1=50,r=1.5;
Constant m=[41.76554,21.16515,12.83882,10.98004,14.71395,24.88945],
         n=[0,2.5,5,7.5,10,12.5];
ConstStr a=-2*x1, b=2*(x4-x2),c=2*(x3*x2-x4*x1),d=2*(x1-x3),
         e=2*(x3*x1+x2*x4-x1^2-x2^2),
         f=2*(x1^2+x2^2)+x3^2+x4^2-2*x2*x4-2*x1*x3-(r1+r)^2,
         w=n*pi/180;
MinFunction Sum(m,n)(abs(m^2+a*m+b*m*sin(w)+c*sin(w)+d*m*cos(w)+e*cos(w)+f));

目标函数值(最小): 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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 myb19880127 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 316求调剂 +4 想读研究生( ?∵ 2026-04-07 4/200 2026-04-11 11:43 by zhq0425
[考研] 调剂求助 +4 果然有我 2026-04-11 5/250 2026-04-11 11:03 by 阳光向上666
[考研] 296求调剂 +13 汪!?! 2026-04-08 14/700 2026-04-11 10:51 by 逆水乘风
[考研] 288求调剂 +15 代fish 2026-04-09 16/800 2026-04-11 10:26 by wwj2530616
[考研] 广东省 085601 329分求调剂 +14 Eddieddd 2026-04-10 14/700 2026-04-11 09:58 by bljnqdcc
[考研] 求调剂 +13 雪逢冬 2026-04-10 13/650 2026-04-11 09:58 by 猪会飞
[考研] 本9 一志愿西工大085601 324求调剂 +5 wysyjs25 2026-04-10 5/250 2026-04-10 16:57 by luoyongfeng
[考研] 285求调剂 +9 AZMK 2026-04-07 11/550 2026-04-10 15:24 by AZMK
[考研] 机械还有还有名额吗?太难了 +6 笑笑袁 2026-04-10 6/300 2026-04-10 11:54 by 高维春
[考研] 调剂 +24 不逢春 2026-04-07 25/1250 2026-04-10 09:39 by 690616278
[考研] 求调剂 +15 张zic 2026-04-05 16/800 2026-04-10 08:12 by kangsm
[考研] 311求调剂 +6 surte 2026-04-08 13/650 2026-04-09 14:00 by surte
[考研] 求调剂 +7 chenxrlkx 2026-04-05 9/450 2026-04-09 09:04 by wj165256
[考研] 388求调剂 +6 四川王涛 2026-04-07 8/400 2026-04-08 00:17 by JourneyLucky
[考研] 316求调剂 +4 15318418673 2026-04-07 4/200 2026-04-07 22:12 by hemengdong
[考研] 328求调剂 +4 ghhh88888 2026-04-06 5/250 2026-04-07 14:45 by ghhh88888
[考研] 327考研调剂推荐 +6 呜呜呜呜呢 2026-04-06 6/300 2026-04-06 21:39 by 啵啵啵0119
[考研] 362求调剂一志愿中国石油大学 +4 我要考大 2026-04-06 6/300 2026-04-06 14:11 by 无际的草原
[考研] 22408 总分320,一篇论文二作,两个国三,求调剂 +3 Leomulufu 2026-04-04 5/250 2026-04-05 19:04 by chongya
[考研] 085602调剂 初试总分335 +12 19123253302 2026-04-04 12/600 2026-04-05 08:08 by 544594351
信息提示
请填处理意见