liaochang0.m

来自「matlab讲义」· M 代码 · 共 11 行

a=[1.25 8.75 0.5  5.75 3   7.25];
b=[1.25 0.75 4.75 5    6.5 7.75];
d=[3 5 4 7 6 11];e=[20,20];w1=[5,1];w2=[2,7];
u=sqrt((w1(1)-a).^2+(w1(2)-b).^2);
v=sqrt((w2(1)-a).^2+(w2(2)-b).^2);
c=[u,v];
A=[eye(6),eye(6);ones(1,6),zeros(1,6);zeros(1,6),ones(1,6)];
B=[d,e];v=zeros(1,12);
x=lp(c,A,B,v,[],[],6);
[x(1:6)';x(7:12)'],z=c*x