24小时热门版块排行榜    

查看: 4555  |  回复: 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的回帖

cooooldog

铁杆木虫 (著名写手)

ส็็็

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
chenyangbu: 金币+14, ★★★★★最佳答案 2016-03-08 17:38:18
引用回帖:
5楼: Originally posted by cmit at 2016-03-05 07:00:03
cpu是不是不一样,指令集有区别?

的确;这种浮点数上的微小差别,只能是CPU硬件级别的差异。

从代差时间看,主要还是CPU升级的影响,楼主把CPU的型号报出来看看
ส็็็็็็็็็็็็็็็็็็็็
7楼2016-03-05 09:01:06
已阅   回复此楼   关注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的回帖

人民海军

木虫 (职业作家)

操作系统是不是都64位的?

发自小木虫Android客户端
Letbygonesbebygones.
4楼2016-03-05 05:51:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 基金不中,共勉 +8 eulota 2026-08-26 8/400 2026-08-27 12:34 by lice_1987
[基金申请] 为什么资助数各大高校都创新高,自己申请怎么就这么难 +5 Kittylucky 2026-08-27 5/250 2026-08-27 12:10 by 姜太公烤鱼
[教师之家] 导师吐槽:我怎么摊上了这么个极品研究生! +5 苏东坡二世 2026-08-23 5/250 2026-08-27 10:39 by 炜雨
[基金申请] 怎么看青基中了没有啊 +5 叶九微 2026-08-26 5/250 2026-08-27 10:35 by l_zh2008
[基金申请] 看板上这么多中的,有点像50人群里49个人都是骗子的那种感觉…… +4 a089 2026-08-26 5/250 2026-08-27 09:24 by Howard28
[基金申请] 让我中一个面上吧! +13 大萍1987 2026-08-20 17/850 2026-08-27 08:20 by 大萍1987
[基金申请] 国自然面上复盘~欢迎讨论 +7 晴天加油 2026-08-26 8/400 2026-08-27 07:53 by longfie172
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +7 process2012 2026-08-23 10/500 2026-08-26 19:23 by hmhminy
[基金申请] 系统查不到 +10 董八千 2026-08-26 10/500 2026-08-26 16:30 by Equinoxhua
[基金申请] 怎么查啊 +5 huang1991js 2026-08-26 5/250 2026-08-26 16:09 by Equinoxhua
[基金申请] 范进中举一文的中心思想 +9 炎黄贵胄 2026-08-22 10/500 2026-08-26 15:40 by semaglutide
[基金申请] 国际合作可查了,中了面上 +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 国合现在查不到了吗? +10 chengyan1220 2026-08-24 20/1000 2026-08-26 08:57 by peasantsprig
[基金申请] 牛来!米来!面来! +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
[基金申请] 人气不行了 +11 fansofjerry 2026-08-21 11/550 2026-08-25 11:04 by 孤独的英雄6
[教师之家] 跳槽后在研项目怎么办? +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[基金申请] 时间戳今天,20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
信息提示
请填处理意见