func7p4-try.m

来自「数值分析最常用的四十种算法」· M 代码 · 共 18 行

%使用穷举法解决四村问题
a=[0,0];b=[0,3],c=[8,1];d=[10,5]
x1=[0,0];x2=[0,0];h=0.2;xo1=[0,0];xo2=[0,0]
s=norm(a-x1)+norm(b-x1)+norm(x2-x1)+norm(c-x2)+norm(d-x2)
start=now
x1=[0.8,1.4];x2=[7.8,1.4]
        
                    s=norm(a-x1)+norm(b-x1)+norm(x2-x1)+norm(c-x2)+norm(d-x2)

%使用穷举法，时间消耗是惊人的。本题循环25*25*50*50次，且其中包含开方运算，在C4-1.7G/256M的PC上需要20多秒；通过最优化，计算
%几何等方法简化算法,是今后的努力方向。
%Answer:距离 S=16.7555   最佳点 xo1=[0.8000,1.8000]  xo2=[4.8000,2.2000]