自编的绘制任意多边形程序

function dbx(x,y)
a=min(x);b=max(x);c=min(y);d=max(y);
[X,Y]=meshgrid(linspace(a,b,100),linspace(c,d,100));k=size(X)
Z=X+Y*i;
z=x+y*i;n=length(x);P=zeros*Z;Q=P;R=P;T=P;
P=abs(Z-z(1));Q=abs(Z-z(2));R=abs(Z-z(3));T=abs(Z-z(4));S=P+Q+R+T;S=reshape(S,1,k(1)*k(2));
m=find(S==min(S));
Z=reshape(Z,1,k(1)*k(2));Zmin=Z(m);
a1=angle(z-Zmin);
b1=sort(a1);
c1=zeros(size(a1));
for j=1:n
c1(j)=find(a1==b1(j));
end
c1
for j=1:(n-1)
hold on,plot([x(c1(j)),x(c1(j+1))],[y(c1(j)),y(c1(j+1))]);hold on,plot([x(c1(j)),x(c1(j+1))],[y(c1(j)),y(c1(j+1))]);
if j==n-1 plot([x(c1(j+1)),x(c1(1))],[y(c1(j+1)),y(c1(1))]);hold off
end
end