24小时热门版块排行榜    

查看: 1608  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 0856求调剂285 +3 吕仔龙 2026-02-28 3/150 2026-03-01 07:59 by ms629
[考研] 材料类求调剂 +7 wana_kiko 2026-02-28 7/350 2026-03-01 07:55 by ms629
[考研] 291分工科求调剂 +3 science饿饿 2026-03-01 3/150 2026-03-01 07:49 by ms629
[考研] 298求调剂 +5 axyz3 2026-02-28 5/250 2026-03-01 06:45 by 刘兵
[考研] 285求调剂 +6 满头大汗的学生 2026-02-28 6/300 2026-03-01 06:29 by Trying]
[考研] 材料调剂 +4 爱擦汗的可乐冰 2026-02-28 4/200 2026-03-01 00:38 by 猫猫球alter
[考研] 307求调剂 +4 73372112 2026-02-28 6/300 2026-03-01 00:04 by ll247
[考研] 304求调剂 +3 52hz~~ 2026-02-28 5/250 2026-03-01 00:00 by 52hz~~
[考研] 295求调剂 +5 19171856320 2026-02-28 5/250 2026-02-28 21:39 by gaoxiaoniuma
[考研] 材料学调剂 +5 提神豆沙包 2026-02-28 5/250 2026-02-28 21:34 by gaoxiaoniuma
[考研] 求调剂 +4 repeatt?t 2026-02-28 4/200 2026-02-28 21:16 by gaoxiaoniuma
[考研] 085600材料工程一志愿中科大总分312求调剂 +8 吃宵夜1 2026-02-28 10/500 2026-02-28 20:27 by L135790
[考研] 298求调剂 +8 人间唯你是清欢 2026-02-28 11/550 2026-02-28 20:26 by L135790
[考研] 276求调剂 +3 路lyh123 2026-02-28 4/200 2026-02-28 19:45 by 路lyh123
[高分子] 求环氧树脂研发1名 +3 孙xc 2026-02-25 11/550 2026-02-28 16:57 by ichall
[考研] 265分求调剂不调专业和学校有行学上就 +4 礼堂丁真258 2026-02-28 6/300 2026-02-28 16:18 by 求调剂zz
[考研] 0856调剂 +3 刘梦微 2026-02-28 3/150 2026-02-28 13:22 by houyaoxu
[考研] 寻找调剂 +3 LYidhsjabdj 2026-02-28 3/150 2026-02-28 12:59 by miniwendy
[硕博家园] 博士自荐 +6 科研狗111 2026-02-26 9/450 2026-02-28 12:32 by seaskyy
[硕博家园] 【博士招生】太原理工大学2026化工博士 +4 N1ce_try 2026-02-24 8/400 2026-02-26 08:40 by N1ce_try
信息提示
请填处理意见