| 查看: 4567 | 回复: 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。 求教大神，是否可能因为不同机器的浮点数运算部件不同导致会有最终结果有精度的差异？ |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有5人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有5人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有3人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有3人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有4人回复
-
怎么看青基中了没有啊
已经有6人回复
-
国社科又开始会评了，不知道这次命运如何
已经有11人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有3人回复
-
有多少人是今天查系统知道结果的？
已经有19人回复
-
基金未中，这种答复是模板吗？
已经有7人回复
ll550
木虫 (职业作家)
- 应助: 259 (大学生)
- 金币: 4586.3
- 散金: 230
- 红花: 34
- 帖子: 3331
- 在线: 767小时
- 虫号: 3180771
- 注册: 2014-05-05
- 专业: 马、恩、列、斯思想研究
【答案】应助回帖
★
感谢参与，应助指数 +1
chenyangbu: 金币+1 2016-03-08 17:38:46
感谢参与，应助指数 +1
chenyangbu: 金币+1 2016-03-08 17:38:46
|
哥们儿，小数点后12位的差异 MATLAB采用的double来进行运算的，所以一般如果输入数据精度够的话，输出基本差别1e-12-1e-16.你这个完全是可以接受的。 这个是计算机本身精度的问题。 当然也有可能是你的输入数据精度不够导致的。 |
chenyangbu
12楼2016-03-05 13:22:27
chenyangbu
金虫 (著名写手)
- 应助: 1 (幼儿园)
- 金币: 5326.6
- 散金: 2666
- 红花: 10
- 帖子: 1050
- 在线: 206.9小时
- 虫号: 2831799
- 注册: 2013-11-27
- 性别: GG
- 专业: 计算机应用技术
2楼2016-03-04 21:39:46
ll550
木虫 (职业作家)
- 应助: 259 (大学生)
- 金币: 4586.3
- 散金: 230
- 红花: 34
- 帖子: 3331
- 在线: 767小时
- 虫号: 3180771
- 注册: 2014-05-05
- 专业: 马、恩、列、斯思想研究
3楼2016-03-04 22:49:17
4楼2016-03-05 05:51:01
10