已知圆弧槽形，反求砂轮/铣刀廓形

clc ;clear all ;close all;
r1=15;
r2=30;
b=30;%%螺旋角
c=75;%%刀具轴向与工件轴线最短距离
d=60;
%u=10;
%u1=0:5:75;
u=0:10:70;
n=length(u);


for u1=1:n
syms a
eqn=r2*(a/180*pi)*(cotd(b))^2*cosd(u(u1)-a)-c*sind(u(u1)-a)+sind(u(u1))*(c+r1*cosd(u(u1)-a)+r2*(1-cosd(a)))==0;
a=vpasolve(eqn);

xt=r2-r1*cosd(u(u1));
yt=r1*sind(u(u1));

x=xt*cosd(a)-yt*sind(a);
y=xt*sind(a)+yt*cosd(a);
z=r2*(a/180*pi)*cotd(b);

x0=x-c;
z0=(-y)*sind(d)+z*cosd(d);
y0=y*cosd(d)+z*sind(d);

R=sqrt(x0^2+y0^2);
xx(u1)=R;
yy(u1)=z0;
end 
xx=double(xx)
yy=double(yy)
plot(yy,xx);

%fprintf('%6.2f is u ; %6.8f is a; %6.8f is z; %6.8f is R;%6.8f ',u,a,z0,R,y);
fprintf('%6.2f is u1 ; %6.8f is n',u1,n);