function myfun = pbpCircleintpol(x0, y0, x1, y1)
%x0,y0,x1,y1分别为起始点横、纵坐标和终止点的横、纵坐标
%以下程序是将图形限制在area区域内，并且将网格间距调整为1个单位长度 area = max(max(x0,y0),max(x1,y1)); plot([0 0],[area area]); grid on
set(gca,'XTick',[0:1:area]) set(gca,'YTick',[0:1:area]) hold on;
%以下部分程序是画出需要加工的工件轮廓 r = sqrt(y0^2 + x0^2); %求圆的半径 theta0 = atan(y0/x0); theta1 = atan(y1/x1);
theta = theta0:0.01:theta1 dx = r * cos(theta); dy = r * sin(theta); plot (dx, dy); hold on;
%以下部分程序是模拟走刀路径 num = abs(x0-x1) + abs(y0-y1); px = x0; py = y0;
rxy = r; %当前位置离圆心距离，初始化为r
for i=1:num lastX = px;
lastY = py; %lastX，lastY为走刀之前的位置坐标 if (rxy