多次运行结果不一致

zyz数模

Max  z = x12 + x22 +3x32 + 4x42 + 2x5 2 −8x1 −2x2 −3x3 − x4 −2x 5
约束条件：
0 ≤ xi ≤ 99        (i =1,,5)
x1 + x2 + x3 + x4 + x5   ≤ 400
x1 + 2x2 + 2x3 + x4 +6x5 ≤ 800
2x1 + x2 +6x3   ≤ 200
x3 + x4 +5x5 ≤ 200


function [f,g]=mengte(x); f=x(1)^2+x(2)^2+3*x(3)^2+4*x(4)^2+2*x(5)-8*x(1)-2*x(2)-3*x(3)-...
x(4)-2*x(5); g=[sum(x)-400
x(1)+2*x(2)+2*x(3)+x(4)+6*x(5)-800
2*x(1)+x(2)+6*x(3)-200 x(3)+x(4)+5*x(5)-200];

主程序：
rand('state',sum(clock));
p0=0;
tic
for i=1:10^6
    x=99*rand(5,1);
    x1=floor(x);x2=ceil(x);
    [f,g]=mengte(x1);
    if sum(g<=0)==4
        if p0<=f
            x0=x1;p0=f;
        end
    end
    [f,g]=mengte(x2);
    if sum(g<=0)==4
        if p0<=f
            x0=x2;p0=f;
        end
    end
end
x0,p0
toc
采用的是蒙特卡洛法，程序多次运行，每次运行结果都不一样，是这种算法的缺陷，还是我的程序不够好。多谢各位指导。

magic2728

zyz数模 发表于 2013-8-10 23:11
多谢提醒，因为初学matlab，没有接触lingo。这个算法运行多次，还是不能得出最优解，在自己加的一个matla ...

运行的够多程序的解肯定会越来越好的！

zyz数模

多谢各位，今天因为和姐姐跑出去玩，还看了一场电影（好开心的说！），登陆数学中国之后看到你们的建议心里真的好开心！

zyz数模

wujianjack2 发表于 2013-8-10 18:42
嗯，很好很好！我平常编程风格就是这样的。
  不过这次我没这么写有两个理由：
   1.估摸着这么写程序 ...

嘿嘿，多谢啦！初接触建模，几个师兄建议我先学习功能强大的matlab软件，其他软件就没有接触了（自己比较笨，东西多了可能掌握不了）。最近一直在看关于建模的东西，发现自己所想到的解决办法当转化为程序时真心好难

zyz数模

wujianjack2 发表于 2013-8-10 10:44
楼主你好，你的想法不错啊！
  首先，我对你的问题有一些疑惑，主要有如下两个方面：
  1.个人觉得你的 ...

多谢指出我的错误，因为初接触建模，几个师兄建议我学功能强大的matlab，其他软件就没有去研究了。我在一个matlab群里发问，一个师兄算出5100+，而我的最大值是5000+，多次运算结果在4900+处浮动，（具体数字没有列出，不好意思）多谢你的提醒

zyz数模

magic2728 发表于 2013-8-10 15:08
蒙特卡洛算法就是一种随机性的算法，所以每次运行结果不一样很正常，它的基本思想就是不断在可行域内找点试 ...

多谢提醒，因为初学matlab，没有接触lingo。这个算法运行多次，还是不能得出最优解，在自己加的一个matlab群里有一位回复了我，他的答案要优于我，算法却不一样，同时也没有给我看他的程序。不知道有没有对我程序的改进方法。至于用到floor以及ceil是因为题目要求是整数解。

wujianjack2

magic2728 发表于 2013-8-10 19:42
你的风格是只把自变量当成集里面的元素写在程序里，我的风格是把所有的参数都按照集的方式写进来，这样的 ...

嗯，大家都是为了能更好地解决楼主问题而尽力！至于风格嘛，也是可以变通学习的。

magic2728

wujianjack2 发表于 2013-8-10 18:42
嗯，很好很好！我平常编程风格就是这样的。
  不过这次我没这么写有两个理由：
   1.估摸着这么写程序 ...

你的风格是只把自变量当成集里面的元素写在程序里，我的风格是把所有的参数都按照集的方式写进来，这样的话会用到很多lingo的函数，如求和，循环等等，能体现lingo这个软件的特点和优势。而你那个却是更加明确一些。因为到底是变量还是参数，一目了然不会混淆。

wujianjack2

  4楼的回复很不错，楼主可以参考下！

wujianjack2

magic2728 发表于 2013-8-10 15:08
蒙特卡洛算法就是一种随机性的算法，所以每次运行结果不一样很正常，它的基本思想就是不断在可行域内找点试 ...

  嗯，很好很好！我平常编程风格就是这样的。
  不过这次我没这么写有两个理由：
   1.估摸着这么写程序显得更长，而且会涉及到LINGO过多的规则，不如写简短直观点，方便大家阅读；
   2.楼主既然用MATLAB解题，可能以前没有接触过LINGO，我有意要突出LINGO语言的表述可以与实际模型基本接近的优势，以此来向楼主推荐比较专业的优化软件LINGO，是别有居心的，哈哈。