bezier3.m

这是一次

 A=[-1,3,-3,1;3,-6,3,0;-3,3,0,0;1,0,0,0];
 %x=[1 2 3 3.6];
%y=[1 3 1 2.4];
%p=size(X);
%n=p(1,2);
%point_number=findobj(gcbf,'Tag','EditText1');
%n=str2num(get(point_number,'string'));
[x y]=ginput(4)
plot(x,y,'m+',x,y,'-.')
axis([0 10 0 10])
hold on
X=x;
Y=y;
for i=1:1001;
   t=(i-1)/1000;
   xx(i)=[t^3 t^2 t 1]*A*X;
   yy(i)=[t^3 t^2 t 1]*A*Y;
end
pause
hold on
plot(xx,yy,'r')
hold on
m=3;
while m<5
%[x1 y1]=ginput(1);
%x=[x;x1];y=[y;y1];
[x1 y1]=ginput(2);
%x=[x;x1;x];y=[y;y1;y];
%plot(x,y,'m+',x,y,'-.')
hold on
%X=[X(2:4,1);x1];
%Y=[Y(2:4,1);y1];
X=[X(1,1);x1;X(4,1)];
Y=[Y(1,1);y1;Y(4,1)];
plot(X,Y,'m+',X,Y,'-.')
for i=1:1001;
   t=(i-1)/1000;
   xx(i)=[t^3 t^2 t 1]*A*X;
   yy(i)=[t^3 t^2 t 1]*A*Y;
end
pause
hold on
plot(xx,yy,'r')
m=m+1;
end