求包含多边形的最小圆的1stOpt解法

shamohu

附件PDF是一篇论文“求包含多边形的最小圆的遗传算法”，现示例用1stOpt求解：

代码：

DataSet;
        x,y=
        5.175635*(3+3.5*cos(0.9)) 4.741644*(2+3.5*sin(0.9))
        2.429613*(3+3.2*cos(1.75)) 5.148755*(2+3.2*sin(1.75))
        -0.297778*(3+3.5*cos(2.8)) 3.172459*(2+3.5*sin(2.8))
        1.0390690*(3+3.0*cos(4)) -0.270407*(2+3.0*sin(4))
        5.714481*(3+3.5*cos(5.6)) -0.209433*(2+3.5*sin(5.6))
EndDataSet;
Parameter x0, y0, r0;
MinFunction r0;                                                          //目标函数：半径r0最小
            For(i=1:5,x,y)(sqrt((x-x0)^2+(y-y0)^2)<=r0);  //约束：每点必须在园内

结果：
x0: 16.8774761309609
y0: 10.8490969218296
r0: 19.123088881071
约束函数
   1: sqrt(((26.7871970784569)-x0)^2+((22.4831886937579)-y0)^2)-(r0) = -3.840593646
   2: sqrt(((5.90301841118487)-x0)^2+((26.5096782044701)-y0)^2)-(r0) = 0
   3: sqrt(((0.0886717945587108)-x0)^2+((10.0644946014941)-y0)^2)-(r0) = -2.315960874
   4: sqrt(((1.0796645295386)-x0)^2+((0.0731200770461928)-y0)^2)-(r0) = -3.552713679E-15
   5: sqrt(((32.655290669483)-x0)^2+((0.0438622301932988)-y0)^2)-(r0) = 0

而原文的结果是：
x0: 3
y0: 2
r0: 3.5

不知道哪儿出问题了？

wanglei5201118

奇怪啊，可以试试其它的软件啊

wanglei5201118

版主的程序1.0版的怎么计算不了啊，说是函数表达式有错误

wanglei5201118

我计算了一下坐标并画了图，从图上看圆文算的应该有问题