[转帖]数据排序函数

linda

[转帖]数据排序函数function shujupaixu(x,y)
%x=[-150.15 -190.68 -60.696 -106.53 -228.15 -262.03 -292.56 -319.19 -340.96 -357.69 -358.76 -373.12 -371.73 -364.5 -286.34 -293.94 -299.68 -307.31 -308.56 -307.84 -299.54 -294.38 -285.93 -276.33 -265.93 -254.85 -242.6 -232.73 -223.34 -215.98 -208.78 -196.96 -190.87 -188.9 -188.11 -191.77 -172.35 -178.78 -136.85 -51.993 0.5328 47.901 89.581 122.44 157.47 195.22 234.39 241.67 268.07 293.04 311.74 331.06 322.74 326.05 324.15 317.29 304.38 285.53 260.99 231.26 197.36 159.99 118.56 73.696 28.82 -15.352 -60.696 -60.696 -60.696 -60.696 -60.696 -60.696 -60.696];
%y=[-355.61 -336.36 -377.16 -369.43 -313.26 -287.48 -259.95 -230.29 -198.79 -162.69 -133.88 -69.54 -20.059 22.862 55.472 64.554 75.369 87.294 101.24 115.55 131.07 141.16 149.5 158.62 164.99 168.41 168.67 168.25 162.08 159.19 154.19 134.8 122.71 114.06 105.71 82.116 74.341 61.697 34.991 33.975 30.926 28.617 28.735 31.622 12.644 -8.3258 -31.478 -80.975 -114.25 -123.29 -145.16 -163.84 -168.29 -190.42 -208.23 -223.45 -242.59 -265.13 -287.6 -308.57 -325.97 -339.57 -355.18 -372.31 -381.45 -381.36 -377.16 -377.16 -377.16 -377.16 -377.16 -377.16 -377.16];
X(1)=x(1); Y(1)=y(1); e=length(x);
for i=1:(e-1)
  x=x(2:e-i+1); y=y(2:e-i+1); c=zeros(size(x));  c(:)=x(1); d=zeros(size(y)); d(:)=y(1);
  for j1=1:(e-i), T(j1)=sqrt((x(j1)-c(j1))^2+(y(j1)-d(j1))^2); end
  [a,index]=sort(T);
  if index==1, X(1+i)=x(index(1)); Y(1+i)=y(index(1));
  else
     X(1+i)=x(index(2)); Y(1+i)=y(index(2)); T=T(2:(e-i));
     U=x(2); x(2)=x(index(2)); x(index(2))=U;
     U=y(2); y(2)=y(index(2)); y(index(2))=U;
  end
end
X(e+1)=X(1); Y(e+1)=Y(1);
plot(X,Y); xlabel('X'),ylabel('Y'); title('按照距离排序后的二维封闭图形');

[ 本帖最后由 ChaChing 于 2009-12-14 13:34 编辑 ]