bezier2.m

来自「这是一次」· M 代码 · 共 46 行

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