sketcher.m

来自「薛定宇高等应用数学问题的MATLAB求解例子源码」· M 代码 · 共 14 行

   function sketcher(vis)
   x=[]; y=[]; i=1; h=[];
   axis([0,1 0 1])
   while 1
      [x0,y0,but]=ginput(1);
      if but==1,
         x=[x,x0]; y=[y,y0];
         h(i)=line(x0,y0); set(h(i),'Marker','o'); i=i+1;
      else, break; end
   end
   if nargin==1, delete(h); end  % 若需要，可以删除样本点标识
   xx=[x(1):(x(end)-x(1))/100: x(end)];
   yy=interp1(x,y,xx,'spline'); line(xx,yy)