24小时热门版块排行榜    

查看: 4599  |  回复: 17
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

chenyangbu

金虫 (著名写手)

[求助] 求助大神,为何不同机器运行MATLAB结果不同 已有3人参与

求助:不同机器运行MATLAB结果不同

我调用MATLAB优化工具箱的库函数fmincon,使用相同的初始解(可行解),对同一个问题进行局部搜索(算法为序列二次规划,即SQP),但在不同机器上得到的结果不同。

一共有五台机器 (为了方便,简称为A1,A2,A3,  B1,B2),  B1-B2配置相同,均为新买的机器,A1-A3为以前的机器。A1-A3上运行结果相同, B1-B2上结果相同。

在A1-A3以及B1-B2上均安装了MATLAB2014b和MATLAB2015a,说明MATLAB版本对结果没有影响。

A1-A3的系统分别为Windows 7, Windows Server 2008, Windows Server 2003;
B1-B2的系统分别为Windows 7, Windows 10;
都有win7,说明系统对结果应该没有影响。

但A1-A3与B1-B2的结果不同,说明机器对结果有影响?

调用fmincon的核心代码如下:
options=optimset('Algorithm','sqp','GradObj','off','GradConstr','off','TolCon',0,'Display','iter','TolFun',1e-4,'MaxFunEvals',50);
[X,FVAL,EXITFLAG,OUTPUT,LAMBDA,GRAD,HESSIAN]=fmincon(@myFun,x0,[],[],[],[],lb,ub,@mycon,options);

其中,@myFun为所求解函数的目标函数,@mycon为所求解函数的约束函数;lb,ub分别为自变量的下界和上届(lb=[0,0] ,ub=[3,4]).

初始解均为:[2.32952 3.17849]

两种结果如下
1)A1-A3上的运行结果如下:
                                                                                     Norm of First-order
   Iter      F-count     f(x)         Feasibility  Steplength       step            optimality
    0       3   -5.508010e+00   2.000e+00                                          1.000e+00
    1       6   -3.508006e+00   2.685e-08   1.000e+00   2.000e+00   2.000e+00
    2       9   -3.508013e+00   2.927e-12   1.000e+00   6.319e-06   1.000e+00
    3      12   -3.508013e+00   0.000e+00   1.000e+00   1.946e-12   1.279e-08

得到的解为:2.329520197477608   1.178493074117684
目标函数值是:-3.508013271595291


2)B1-B2的运行结果如下:
                                          
                                                                            Norm of First-order
Iter   F-count     f(x)        Feasibility  Steplength     step  optimality
    0       3   -5.508010e+00   2.000e+00                          1.000e+00
    1       6   -3.508006e+00   2.685e-08   1.000e+00   2.000e+00   2.000e+00
    2       9   -3.508013e+00   2.984e-12   1.000e+00   6.319e-06   1.694e+00
    3      12   -3.508013e+00   0.000e+00   1.000e+00   1.972e-12   1.179e-06

得到的解为:2.329520197477604   1.178493074117646
目标函数值是:-3.508013271595249。


求教大神,是否可能因为不同机器的浮点数运算部件不同导致会有最终结果有精度的差异?
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

人民海军

木虫 (职业作家)

Letbygonesbebygones.
4楼2016-03-05 05:51:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 18 个回答

chenyangbu

金虫 (著名写手)

另外,五台机器均为64位机器
2楼2016-03-04 21:39:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ll550

木虫 (职业作家)

livelong
3楼2016-03-04 22:49:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cmit

捐助贵宾 (著名写手)

cpu是不是不一样,指令集有区别?

发自小木虫Android客户端
5楼2016-03-05 07:00:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +6 LwdutQ8HoqWP 2026-09-13 6/300 2026-09-17 01:50 by Tql5LQhh5rLK
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +5 Aj1rhIDL5ixY 2026-09-14 5/250 2026-09-16 21:39 by MHf40yA0Uhyb
[博后之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +5 6F5UbRU2I5hL 2026-09-14 5/250 2026-09-16 20:38 by MHf40yA0Uhyb
[考研] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +7 QUjhNVAcOSff 2026-09-13 8/400 2026-09-16 18:37 by rRvxEj3zorCy
[考博] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 6F5UbRU2I5hL 2026-09-14 4/200 2026-09-16 17:00 by i3KBRD19ERkG
[找工作] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 s3fFTmArrBt6 2026-09-14 4/200 2026-09-16 16:51 by i3KBRD19ERkG
[公派出国] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +5 s3fFTmArrBt6 2026-09-14 5/250 2026-09-16 16:24 by i3KBRD19ERkG
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 s3fFTmArrBt6 2026-09-13 5/250 2026-09-16 15:51 by i3KBRD19ERkG
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 s3fFTmArrBt6 2026-09-13 5/250 2026-09-16 15:48 by i3KBRD19ERkG
[考博] 上海工程技术大学激光智能制造课题组|2027级博士研究生招生公告 +7 水士口 2026-09-16 8/400 2026-09-16 13:47 by kkklens
[考博] 上海工程技术大学激光智能制造课题组招收博士研究生 +8 两三岁ss 2026-09-14 8/400 2026-09-16 09:46 by 水士口
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 Aj1rhIDL5ixY 2026-09-14 4/200 2026-09-16 07:48 by i0IGWqQuchto
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 Aj1rhIDL5ixY 2026-09-14 3/150 2026-09-16 07:12 by LH5NkK5Ud4bh
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 23jxep3nCNZb 2026-09-14 3/150 2026-09-16 06:24 by LH5NkK5Ud4bh
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +5 s3fFTmArrBt6 2026-09-13 7/350 2026-09-16 03:26 by DgNGHc3h5tPl
[考博] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +5 3n8v2C8RimXI 2026-09-13 5/250 2026-09-16 02:50 by DgNGHc3h5tPl
[考博] 申博发邮件 +4 Lmengk 2026-09-15 4/200 2026-09-16 00:05 by 木心想继续深造
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 s3fFTmArrBt6 2026-09-14 5/250 2026-09-15 16:59 by CgyNCDVNhGVg
[教师之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 6F5UbRU2I5hL 2026-09-14 4/200 2026-09-15 07:50 by C79jjtjAKjEn
[教师之家] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 LwdutQ8HoqWP 2026-09-13 4/200 2026-09-14 16:13 by Aj1rhIDL5ixY
信息提示
请填处理意见