chenyangbu

金蟲(chóng) (著名寫(xiě)手)

應(yīng)助: 1 (幼兒園)
金幣: 5320.8
散金: 2666
紅花: 10
帖子: 1048
在線: 206.9小時(shí)
蟲(chóng)號(hào): 2831799
注冊(cè): 2013-11-27
性別: GG
專業(yè): 計(jì)算機(jī)應(yīng)用技術(shù)

[求助] 求助大神，為何不同機(jī)器運(yùn)行MATLAB結(jié)果不同已有3人參與

求助：不同機(jī)器運(yùn)行MATLAB結(jié)果不同

我調(diào)用MATLAB優(yōu)化工具箱的庫(kù)函數(shù)fmincon，使用相同的初始解（可行解），對(duì)同一個(gè)問(wèn)題進(jìn)行局部搜索（算法為序列二次規(guī)劃，即SQP），但在不同機(jī)器上得到的結(jié)果不同。

一共有五臺(tái)機(jī)器 (為了方便，簡(jiǎn)稱為A1，A2，A3,  B1，B2),  B1-B2配置相同，均為新買(mǎi)的機(jī)器，A1-A3為以前的機(jī)器。A1-A3上運(yùn)行結(jié)果相同, B1-B2上結(jié)果相同。

在A1-A3以及B1-B2上均安裝了MATLAB2014b和MATLAB2015a，說(shuō)明MATLAB版本對(duì)結(jié)果沒(méi)有影響。

A1-A3的系統(tǒng)分別為Windows 7, Windows Server 2008, Windows Server 2003;
B1-B2的系統(tǒng)分別為Windows 7, Windows 10;
都有win7,說(shuō)明系統(tǒng)對(duì)結(jié)果應(yīng)該沒(méi)有影響。

但A1-A3與B1-B2的結(jié)果不同，說(shuō)明機(jī)器對(duì)結(jié)果有影響？

調(diào)用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為所求解函數(shù)的目標(biāo)函數(shù)，@mycon為所求解函數(shù)的約束函數(shù)；lb，ub分別為自變量的下界和上屆（lb=[0,0] ,ub=[3,4]).

初始解均為：[2.32952 3.17849]

兩種結(jié)果如下
1）A1-A3上的運(yùn)行結(jié)果如下：
                                                                                 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
目標(biāo)函數(shù)值是：-3.508013271595291

2）B1-B2的運(yùn)行結(jié)果如下：

                                                                        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
目標(biāo)函數(shù)值是：-3.508013271595249。

求教大神，是否可能因?yàn)椴煌瑱C(jī)器的浮點(diǎn)數(shù)運(yùn)算部件不同導(dǎo)致會(huì)有最終結(jié)果有精度的差異？

回復(fù)此樓

» 猜你喜歡

309求調(diào)劑已經(jīng)有3人回復(fù)
070300化學(xué)354求調(diào)劑已經(jīng)有4人回復(fù)
330分求調(diào)劑已經(jīng)有3人回復(fù)
299求調(diào)劑已經(jīng)有8人回復(fù)
一志愿北京理工大學(xué)本科211材料工程294求調(diào)劑已經(jīng)有6人回復(fù)
300求調(diào)劑，材料科學(xué)英一數(shù)二已經(jīng)有8人回復(fù)
招收生物學(xué)/細(xì)胞生物學(xué)調(diào)劑已經(jīng)有5人回復(fù)
070305高分子化學(xué)與物理 304分求調(diào)劑已經(jīng)有7人回復(fù)
289求調(diào)劑已經(jīng)有13人回復(fù)
一志愿哈爾濱工業(yè)大學(xué)材料與化工方向336分已經(jīng)有9人回復(fù)

1樓 2016-03-04 21:37:25

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

回帖支持 ( 顯示支持度最高的前 50 名 )

cmit

捐助貴賓 (著名寫(xiě)手)

應(yīng)助: 6 (幼兒園)
金幣: 3050.9
散金: 5350
紅花: 13
沙發(fā): 1
帖子: 2314
在線: 291.6小時(shí)
蟲(chóng)號(hào): 2966476
注冊(cè): 2014-02-10
性別: GG
專業(yè): 計(jì)算機(jī)應(yīng)用技術(shù)

cpu是不是不一樣，指令集有區(qū)別?

發(fā)自小木蟲(chóng)Android客戶端

贊一下(2人)

回復(fù)此樓

5樓2016-03-05 07:00:03

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

ll550

木蟲(chóng) (職業(yè)作家)

應(yīng)助: 259 (大學(xué)生)
金幣: 4586.3
散金: 230
紅花: 34
帖子: 3331
在線: 767小時(shí)
蟲(chóng)號(hào): 3180771
注冊(cè): 2014-05-05
專業(yè): 馬、恩、列、斯思想研究

這個(gè)精度是可以接受的吧

發(fā)自小木蟲(chóng)Android客戶端

贊一下(1人)

回復(fù)此樓

livelong

3樓2016-03-04 22:49:17

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

chenyangbu

金蟲(chóng) (著名寫(xiě)手)

應(yīng)助: 1 (幼兒園)
金幣: 5320.8
散金: 2666
紅花: 10
帖子: 1048
在線: 206.9小時(shí)
蟲(chóng)號(hào): 2831799
注冊(cè): 2013-11-27
性別: GG
專業(yè): 計(jì)算機(jī)應(yīng)用技術(shù)

引用回帖:

3樓: Originally posted by ll550 at 2016-03-04 22:49:17
這個(gè)精度是可以接受的吧

因?yàn)槲业脑创a以后會(huì)公開(kāi)，想讓其他人可以完全重現(xiàn)我的論文~

贊一下(1人)

回復(fù)此樓

10樓2016-03-05 12:53:29

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

普通回帖

chenyangbu

金蟲(chóng) (著名寫(xiě)手)

應(yīng)助: 1 (幼兒園)
金幣: 5320.8
散金: 2666
紅花: 10
帖子: 1048
在線: 206.9小時(shí)
蟲(chóng)號(hào): 2831799
注冊(cè): 2013-11-27
性別: GG
專業(yè): 計(jì)算機(jī)應(yīng)用技術(shù)

另外，五臺(tái)機(jī)器均為64位機(jī)器

贊一下

回復(fù)此樓

2樓2016-03-04 21:39:46

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

人民海軍

木蟲(chóng) (職業(yè)作家)

應(yīng)助: 101 (高中生)
金幣: 904.2
散金: 10685
紅花: 35
帖子: 4761
在線: 418.8小時(shí)
蟲(chóng)號(hào): 1829469
注冊(cè): 2012-05-22
性別: GG
專業(yè): 概率論與隨機(jī)分析

操作系統(tǒng)是不是都64位的？

發(fā)自小木蟲(chóng)Android客戶端

贊一下

回復(fù)此樓

Letbygonesbebygones.

4樓2016-03-05 05:51:01

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

IEEEYLQ

�。。。�！

發(fā)自小木蟲(chóng)Android客戶端

6樓2016-03-05 07:09:10

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

cooooldog

鐵桿木蟲(chóng) (著名寫(xiě)手)

ส็็็

應(yīng)助: 237 (大學(xué)生)
金幣: 6101.9
散金: 1114
紅花: 39
帖子: 1380
在線: 553.8小時(shí)
蟲(chóng)號(hào): 506699
注冊(cè): 2008-02-18
專業(yè): 模式識(shí)別

【答案】應(yīng)助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
chenyangbu: 金幣+14, ★★★★★最佳答案 2016-03-08 17:38:18

引用回帖:

5樓: Originally posted by cmit at 2016-03-05 07:00:03
cpu是不是不一樣，指令集有區(qū)別?

的確；這種浮點(diǎn)數(shù)上的微小差別，只能是CPU硬件級(jí)別的差異。

從代差時(shí)間看，主要還是CPU升級(jí)的影響，樓主把CPU的型號(hào)報(bào)出來(lái)看看

贊一下

回復(fù)此樓

ส็็็็็็็็็็็็็็็็็็็็

7樓2016-03-05 09:01:06

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

chenyangbu

金蟲(chóng) (著名寫(xiě)手)

應(yīng)助: 1 (幼兒園)
金幣: 5320.8
散金: 2666
紅花: 10
帖子: 1048
在線: 206.9小時(shí)
蟲(chóng)號(hào): 2831799
注冊(cè): 2013-11-27
性別: GG
專業(yè): 計(jì)算機(jī)應(yīng)用技術(shù)

引用回帖:

7樓: Originally posted by cooooldog at 2016-03-05 09:01:06
的確；這種浮點(diǎn)數(shù)上的微小差別，只能是CPU硬件級(jí)別的差異。

從代差時(shí)間看，主要還是CPU升級(jí)的影響，樓主把CPU的型號(hào)報(bào)出來(lái)看看...

各CPU型號(hào)為：
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

贊一下

回復(fù)此樓

8樓2016-03-05 12:51:51

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

chenyangbu

金蟲(chóng) (著名寫(xiě)手)

應(yīng)助: 1 (幼兒園)
金幣: 5320.8
散金: 2666
紅花: 10
帖子: 1048
在線: 206.9小時(shí)
蟲(chóng)號(hào): 2831799
注冊(cè): 2013-11-27
性別: GG
專業(yè): 計(jì)算機(jī)應(yīng)用技術(shù)

引用回帖:

4樓: Originally posted by 人民海軍 at 2016-03-05 05:51:01
操作系統(tǒng)是不是都64位的？

操作系統(tǒng)都是64位的~

贊一下

回復(fù)此樓

9樓2016-03-05 12:52:16

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

最具人氣熱帖推薦 [查看全部]	作者	回/看	最后發(fā)表

[考研] 300求調(diào)劑，材料科學(xué)英一數(shù)二 +8	leaflight 2026-03-24	8/400	2026-03-29 01:31 by fmesaito
[考研] 085602 化工專碩 338分求調(diào)劑 +12	路癡小琪 2026-03-27	12/600	2026-03-28 15:41 by L135790
[考研] 化學(xué)0703 調(diào)劑 306分一志愿211 +4	26要上岸 2026-03-28	4/200	2026-03-28 15:30 by 1018329917
[考研] 346求調(diào)劑一志愿070303有機(jī)化學(xué) +3	蘿卜燉青菜 2026-03-28	3/150	2026-03-28 14:11 by 唐沐兒
[考研] 304求調(diào)劑 +6	曼殊2266 2026-03-27	6/300	2026-03-28 14:10 by 唐沐兒
[考研] 0703化學(xué) +10	妮妮ninicgb 2026-03-27	10/500	2026-03-28 14:07 by 噠噠噠呱呱呱
[考研] 0703化學(xué)求調(diào)劑 +9	奶油草莓. 2026-03-22	10/500	2026-03-28 13:30 by 唐沐兒
[考研] 一志愿南京航空航天大學(xué)材料學(xué)碩求調(diào)劑 +3	@taotao 2026-03-28	3/150	2026-03-28 10:26 by JourneyLucky
[考研] 一志愿南昌大學(xué)324求調(diào)劑 +7	hanamiko 2026-03-27	7/350	2026-03-28 09:56 by 李上岸0921
[考研] 308求調(diào)劑 +7	墨墨漠 2026-03-27	7/350	2026-03-28 07:43 by 熱情沙漠
[考研] 330一志愿中國(guó)海洋大學(xué) 化學(xué)工程 085602 有讀博意愿求調(diào)劑 +3	wywy.. 2026-03-27	4/200	2026-03-28 03:32 by fmesaito
[考研] 265求調(diào)劑11408 +3	劉小鹿lu 2026-03-27	3/150	2026-03-27 20:53 by nihaoar
[考研] 287求調(diào)劑 +10	land xuxu 2026-03-26	10/500	2026-03-27 15:33 by 帕爾馬拉特
[考研] 考研化學(xué)308分求調(diào)劑 +10	你好明天你好 2026-03-23	12/600	2026-03-27 14:43 by shangxh
[考研] 298調(diào)劑 +3	jiyingjie123 2026-03-27	3/150	2026-03-27 11:57 by wxiongid
[考研] 0854電子信息求調(diào)劑 +7	α____ 2026-03-22	9/450	2026-03-25 13:37 by α____
[考研] 285求調(diào)劑 +3	AZMK 2026-03-24	3/150	2026-03-25 12:23 by userper
[考研] 一志愿北化315 求調(diào)劑 +3	akrrain 2026-03-24	3/150	2026-03-24 19:35 by 了了了了。。
[考研] 292求調(diào)劑 +4	鵝鵝鵝額額額額?/a> 2026-03-24	4/200	2026-03-24 16:41 by peike
[考研] 石河子大學(xué)（211、雙一流）碩博研究生長(zhǎng)期招生公告 +3	李子目 2026-03-22	3/150	2026-03-22 21:01 by 怎么釋?xiě)?/a>