Znn3bq.jpeg
汕头大学海洋科学接受调剂
查看: 4377  |  回复: 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的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

cmit

捐助贵宾 (著名写手)

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

发自小木虫Android客户端
5楼2016-03-05 07:00:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ll550

木虫 (职业作家)

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

chenyangbu

金虫 (著名写手)

引用回帖:
3楼: Originally posted by ll550 at 2016-03-04 22:49:17
这个精度是可以接受的吧

因为我的源码以后会公开,想让其他人可以完全重现我的论文~
10楼2016-03-05 12:53:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

chenyangbu

金虫 (著名写手)

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

人民海军

木虫 (职业作家)

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

发自小木虫Android客户端
Letbygonesbebygones.
4楼2016-03-05 05:51:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
6楼2016-03-05 07:09:10
已阅   回复此楼   关注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的回帖

chenyangbu

金虫 (著名写手)

引用回帖:
7楼: Originally posted by cooooldog at 2016-03-05 09:01:06
的确;这种浮点数上的微小差别,只能是CPU硬件级别的差异。

从代差时间看,主要还是CPU升级的影响,楼主把CPU的型号报出来看看...

各CPU型号为:
A1: Intel(R) Core(TM) i7-2600 3.40GHz
A2: Intel(R) Xeon(R) E5-2640 v2 2.00GHz
A3: Intel(R) Xeon(R) X5670 2.93GHz

B1-B2:
Intel(R) Core(TM) i7-4790 3.60GHz
8楼2016-03-05 12:51:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chenyangbu

金虫 (著名写手)

引用回帖:
4楼: Originally posted by 人民海军 at 2016-03-05 05:51:01
操作系统是不是都64位的?

操作系统都是64位的~
9楼2016-03-05 12:52:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 chenyangbu 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 化学070300 求调剂 +23 哈哈哈^_^ 2026-04-12 23/1150 2026-04-14 16:30 by zhouxiaoyu
[考研] 335求调剂 +19 想上岸呀!! 2026-04-12 21/1050 2026-04-14 16:23 by Art1977
[考研] 297,工科调剂? +6 河南农业大学-能 2026-04-14 6/300 2026-04-14 16:04 by Art1977
[考研] 考研求调剂 +12 子木呐 2026-04-12 13/650 2026-04-14 01:19 by 王珺璞
[考研] 一志愿华南理工大学331分材料求调剂 +10 天下ww 2026-04-09 11/550 2026-04-13 23:25 by pies112
[考研] 0856专硕求调剂 希望是a区院校 +24 好好休息好不好 2026-04-09 27/1350 2026-04-13 22:22 by pies112
[考研] B区0809 ,数一英一,290 求调剂 +3 泠潍1111 2026-04-12 4/200 2026-04-13 20:35 by 学员JpLReM
[考研] 一志愿211 0703化学 346分求调剂 +26 土豆er? 2026-04-09 29/1450 2026-04-13 15:15 by 独醉梦孤城
[考研] 生物学调剂,一志愿西南大学348,Top期刊一区二作、二区三作,三等奖学金三次 +5 candyyyi 2026-04-09 5/250 2026-04-13 09:02 by 可淡不可忘
[考研] 307求调剂 +10 tzq94092 2026-04-10 10/500 2026-04-12 08:18 by wise999
[考研] 求调剂,262机械专硕 +8 嗯yyl 2026-04-08 8/400 2026-04-12 02:31 by 秋豆菜芽
[考研] 22408调剂315分 +3 zhuangyan123 2026-04-09 3/150 2026-04-12 00:25 by 蓝云思雨
[考研] 331求调剂 +5 王国帅 2026-04-11 5/250 2026-04-11 22:56 by 溪涧流水
[考研] 343求调剂 +9 王国帅 2026-04-10 9/450 2026-04-11 20:31 by dongdian1
[考研] 求调剂,一志愿大连理工大学354分 +5 雨声余生 2026-04-11 6/300 2026-04-11 16:12 by 雨声余生
[考研] 289 分105500药学专硕求调剂(找B区学校) +6 白云123456789 2026-04-09 8/400 2026-04-10 21:13 by zhouxiaoyu
[考研] 求调剂 材料与工程 324分 专硕 +19 翩翩一书生 2026-04-10 21/1050 2026-04-10 11:41 by wp06
[考研] 292求调剂 +9 笑笑袁 2026-04-09 9/450 2026-04-10 10:05 by LHGeng
[考博] 博士自荐 +7 可可小胖 2026-04-08 7/350 2026-04-10 08:28 by kimhero
[考研] 349学科化学045106求调剂,化学类都可以 +8 保好懂懂 2026-04-08 8/400 2026-04-09 14:03 by xulei3024
信息提示
请填处理意见