Page 141 - 数学建模算法与应用

第五章 整数规划的探讨

度究竟怎样呢 ?
6
下面就随机取样采集 10 个点计算，应用概率理论来估计一下可信度。不失
一般性，假定一个整数规划的最优点不是孤立的奇点。假设目标函数落在高值区
6
的概率分别为 0.01、0.00001，则当计算 10 个点后，至少有一个点能落在高值区
的概率分别为

解1.首先编写 M 文件 mente. m 定义目标函数f和约束向量函数 g，程序如下：

function [f，g] =mengte(x)；
f=x(1)~2+x(2)2+3*x(3)2+4*x(4)~2+2*x(5)^2-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]；

2. 编写如下 Matlab 程序求问题的解。
rand(‘state’，sum(clock))；名初始化随机数发生器
po =0；
tic

for i =1：10
X=randint(1，5，[1，99］)；
[f，g] =mengte(x)；
if all(g< =0)
if p0<f

xO=x；pO=f；
end
end

end
x0 ，po
Toc %

131