24小时热门版块排行榜    

Znn3bq.jpeg
查看: 1679  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[有机交流] 如何实现卤原子转化 10+3 BT20230424 2026-05-15 5/250 2026-05-16 16:20 by czyzsu
[考研] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +5 cjf4bx70cj 2026-05-14 6/300 2026-05-16 16:17 by 0i5p09z61n
[硕博家园] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 x0mp7owy2b 2026-05-15 3/150 2026-05-16 12:49 by vcdazktkjx
[论文投稿] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 xx7gd5zq4e 2026-05-15 5/250 2026-05-16 12:24 by h3oerqvkv9
[考博] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 l7k6xnh0yc 2026-05-14 4/200 2026-05-16 11:36 by h3oerqvkv9
[考研] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 l7k6xnh0yc 2026-05-14 6/300 2026-05-16 11:29 by h3oerqvkv9
[硕博家园] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 cjf4bx70cj 2026-05-14 6/300 2026-05-16 11:16 by h3oerqvkv9
[硕博家园] 申请博士 +3 呃?呃 2026-05-15 3/150 2026-05-16 11:01 by a4742549
[公派出国] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 l7k6xnh0yc 2026-05-14 4/200 2026-05-16 04:44 by k37jurhrau
[硕博家园] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 l7k6xnh0yc 2026-05-14 5/250 2026-05-16 04:29 by k37jurhrau
[找工作] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 l7k6xnh0yc 2026-05-14 3/150 2026-05-16 01:49 by k37jurhrau
[文学芳草园] 窗边初夏的小雨 +8 阿美_Lml888 2026-05-09 11/550 2026-05-15 23:54 by WASM
[考博] 西南大学考核制博士 +4 lijunjie84 2026-05-11 7/350 2026-05-15 23:20 by 同仁堂教主
[基金申请] 精华III评审感受-评审感受-评审感受 +14 ferrarichen 2026-05-11 18/900 2026-05-15 11:12 by cmhchen
[教师之家] 教学课件你会给同学吗 +8 硕士研究生吗 2026-05-13 8/400 2026-05-14 22:23 by 常规沥青
[有机交流] 求助2,4-二氯-5-嘧啶甲醛的合成方法 20+3 光吃不拉 2026-05-14 5/250 2026-05-14 20:15 by 一切都是空工
[高分子] 本人最近太闲了,谁有问题可以提,每天会统一回复 +8 一切都是空工 2026-05-12 19/950 2026-05-14 20:03 by 一切都是空工
[考博] 材料类只有一篇综述能申博么 +4 乐逍遥谷 2026-05-13 4/200 2026-05-14 12:05 by zhyzzh
[基金申请] 请问大佬b0816评完了吗 +3 市民华南虎 2026-05-12 7/350 2026-05-14 07:41 by 市民华南虎
[论文投稿] 求助大佬sci投稿哪个好中 +3 江沅188 2026-05-12 4/200 2026-05-13 14:35 by 江沅188
信息提示
请填处理意见