精选优质文档-倾情为你奉上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 onset(gca,XTick,0:1:area) set(gca,YTick,0:1:area)hold on;%以下部分程序是画出需要加工的工件轮廓r = sqrt(y02 + x02); %求圆的半径theta0 = atan(y0/x0);theta1 = atan(y1/x1);theta = theta0:0.01:theta1dx = r * cos(theta);dy = r * sin(theta);plot (dx, dy);hold on;%以下部分程序
