24小时热门版块排行榜    

查看: 1889  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +55 医学老男孩 2026-08-13 123/6150 2026-08-19 23:44 by stan_et
[基金申请] 今天基金会出结果吗?20260819 +13 kkkl_v 2026-08-19 13/650 2026-08-19 22:37 by 宝贝虫子
[基金申请] 时间戳变了,能看出什么问题? +12 基诺咪客 2026-08-17 16/800 2026-08-19 22:04 by kk 的小木虫
[基金申请] filecode,4个jtjc了 +11 ziyangfang 2026-08-19 13/650 2026-08-19 21:28 by aasahr
[教师之家] 为什么余额宝的年化利率越来越低?主要原因有哪些? +5 瞬息宇宙 2026-08-15 5/250 2026-08-19 20:23 by super2002521
[基金申请] 今天放榜没戏了吧 +4 yuleib84 2026-08-19 5/250 2026-08-19 19:52 by hhs666
[基金申请] 今天放榜吗? +14 布布和一二 2026-08-19 15/750 2026-08-19 18:07 by gltch
[基金申请] 2027广东省杰青 +4 奶牛小黑 2026-08-15 10/500 2026-08-19 11:02 by wanfengnew
[基金申请] 快农历七夕节了,轻松一下,男人悄悄话,女施主请不要进来。 +6 Tide man 2026-08-14 7/350 2026-08-19 09:56 by ZJTJZ
[基金申请] 什么时候开奖? +6 CrisMessi 2026-08-18 6/300 2026-08-19 08:06 by Equinoxhua
[基金申请] 重要消息,中午系统在维护 +10 yuleib84 2026-08-18 11/550 2026-08-19 07:30 by 羽毛枫f
[基金申请] 明天放榜? +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[论文投稿] 投稿咨询 +4 wwm09 2026-08-17 6/300 2026-08-18 15:36 by wwm09
[基金申请] 时间戳又变了8-15 +14 archvillain 2026-08-15 26/1300 2026-08-18 13:37 by phantomgost
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +13 Ldrop2023 2026-08-13 16/800 2026-08-18 12:25 by 淀粉搬运工
[基金申请] 93BebMhtakh前后11位开头都是大写 +4 且听虎啸 2026-08-17 5/250 2026-08-18 00:49 by 蔡棒棒菂
[基金申请] 感觉是下周放榜了 +6 angus9576 2026-08-17 11/550 2026-08-17 23:57 by angus9576
[基金申请] 今天系统多次维护,明天很可能放榜! +8 zju2000 2026-08-16 9/450 2026-08-17 12:20 by lmz0216
[精细化工] 招聘 金属平磨液,抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
信息提示
请填处理意见