打算用或正在用isight或sigiht-fd的朋友，别急，先看看这道简单的整数优化问题

TBE_Legend

打算用或正在用isight的朋友，别急，先看看这道简单的整数优化问题：

题目见如下帖子中，题目虽简单，但很考验优化算法。帖子中1stopt这个很小的优化软件表现很彪悍！ 我用ModeFrontier和Optimus都没有找到解。

http://forum.simwe.com/thread-913182-1-1.html

想看看isight能否求解，求解时间和结果又如何？ 会sight的朋友给算算吧，程序应该很简单。

多谢~！

jacicka1984

原文：
“用1到9这九个数组成等式( )/( )( )+( )/( )( )+( )/( )( )=1
用1至9这9个数字怎么组合相加后等于1,不可以重复使用数字”
请问楼主，这个问题有什么意义吗？做理工科的可能和专门学数学的想法不一样。
  这个问题最难实现的是9个框框里面的数字 不能重复呗，也就是约束如何定义的过程。这种约束在任何软件里面都没法定义，因为它不是“大于”“等于”“小于”的描述方式。
   你要实现这个约束，为什么不编写一小段程序（Matlab\VB.....），把约束关系在程序里面体现出来，然后用MDO软件去运行。
   须知所有的MDO软件的目的是实现外部程序及软件的嵌入，然后搭建流程，再优化。
   另：不管是做复杂的算法研究还是详细的产品优化设计，一定要知道你问题定义及优化的目的是什么，脱离了实际问题，一切工具都是空架子。

TBE_Legend

原文：
“用1到9这九个数组成等式( )/( )( )+( )/( )( )+( )/( )( )=1
用1至9这9个数字怎么组合相加后等于1,不可以重复使用数字”
请问楼主，这个问题有什么意义吗？做理工科的可能和专门学数学的想法不一样。
...
jacicka1984 发表于 2010-1-21 10:16

多谢关注~！


我是工科的，搞振动控制。

这个问题和我现在做的一个问题类似，都是整数优化问题，但是这个问题比较简单。我的问题比较复杂（用matlab+cfx+ansys（emag））用mF和optimus得到结果都不是很好（减震效果达到了，但是设备需要的功率太大了），而且内容涉及专利，不方便贴图。但是1stopt怎么整合cae软件还不清楚。

不管背景如何，优化问题到头来就是个数学问题。这个问题难就难在约束太强。

就算是空架子，这么简单的一个小问题，闲时做下又何妨？呵呵

多谢~！

jacicka1984

期待有人编10行以内的程序把这个约束定义好。

TBE_Legend

期待有人编10行以内的程序把这个约束定义好。
jacicka1984 发表于 2010-1-21 10:51

按最笨的办法来就行。或驱动matlab编个循环。


abs( x1 - x2 ) * abs( x1 - x3 ) * abs( x1 - x4 ) * abs( x1 - x5 ) * abs( x1 - x6 ) * abs( x1 - x7 ) * abs( x1 - x8 ) * abs( x1 - x9 ) * abs( x2 - x3 ) * abs( x2 - x4 ) * abs( x2 - x5 ) * abs( x2 - x6 ) * abs( x2 - x7 ) * abs( x2 - x8 ) * abs( x2 - x9 ) * abs( x3 - x4 ) * abs( x3 - x5 ) * abs( x3 - x6 ) * abs( x3 - x7 ) * abs( x3 - x8 ) * abs( x3 - x9 ) * abs( x4 - x5 ) * abs( x4 - x6 ) * abs( x4 - x7 ) * abs( x4 - x8 ) * abs( x4 - x9 ) * abs( x5 - x6 ) * abs( x5 - x7 ) * abs( x5 - x8 ) * abs( x5 - x9 ) * abs( x6 - x7 ) * abs( x6 - x8 ) * abs( x6 - x9 ) * abs( x7 - x8 ) * abs( x7 - x9 ) * abs( x8 - x9 )

zxkchina

TBE_Legend,诚如斯言,强约束是个麻烦的地方。不过那个1stopt用零点几秒，5个迭代就找到了，也太夸张了。倒想研究一下那个算法。我做DOE时候，使用cs算法，生成20个样本，迭代1000步最后得到了19个可行解，光是DOE找可行样本就用了零点几秒呢。

另外请大家看一下偶个性签名的最后一行，或许和此问题亦有些关联呢。

TBE_Legend

会isight的朋友帮忙看看吧，多谢~！

TBE_Legend

TBE_Legend,诚如斯言,强约束是个麻烦的地方。不过那个1stopt用零点几秒，5个迭代就找到了，也太夸张了。倒想研究一下那个算法。我做DOE时候，使用cs算法，生成20个样本，迭代1000步最后得到了19个可行解，光是DOE找可 ...
zxkchina 发表于 2010-1-21 16:47

1stopt的算法很神秘，只是说“其计算核心是基于七维高科有限公司科研人员十数年的革命性研究成果【通用全局优化算法】(Universal Global Optimization - UGO)”

再找不到其他内容了，这也是我刚开始看不起这个软件的一个原因，觉得他可能是骗人的，呵呵。但现在不服不行，的确很彪悍，Mathematica的优化也左不过他，Matlab版也有一些题目，他的官方网站上也有一些例题。跑题了跑题了。。。

声明一下： 此次开贴的目的在于如何设置isight或mF或optimus等软件求解我提出的这个问题，绝无意惹起“谁比谁强”的无意义的争论。

我也不是任何软件公司的代理。