MATLAB 的渐开线齿轮的范成法仿真

clear;clc
hd=pi/180;m=2;%模数
z=input('渐开线齿轮齿数z =');
phi0=20;%压力角
x=0;
x0=0;y0=0;r=m*z/2;%分度圆半径
p=pi*m;%齿距
s=2.5*m*tan(phi0* hd);h=(2*s+p)/4;
c=x0-2*p-h;
x1(1,1)=c;y1(1,1)=-(r+(1.25 +x)* m)+y0;
x1(2,1)= x1(1,1)+ s;
y1(2,1)= y1(1,1)+ 2.5* m;
x1(3,1)= x1(2,1)+ (p /2-s);
y1(3,1)= y1 (2,1);x1 (4,1)= x1 (3,1)+ s;
y1(4,1)= y1(3,1)-2.5* m;
for i =5:20
x1(i,1)= x1(i-4,1)+ p;
y1(i,1)= y1(i-4,1);
end
j = 0;
for d_phi = 0:(6* hd):4.398226 /2;j =j + 1;
for i = 1:20;  
x1(i,j)= x1(i,1)-r* d_phi;
y1(i,j)= y1(i,1);
s2 = y1(i,j)- y0;s1 = x1(i,j)-x0;
r1(i,j)= sqrt((s1)^2 + (s2)^2);
phi(i,j)= atan(s1 /s2);
x2(i,j)= r1(i,j)* sin(phi(i,j)-d_phi)+ x0;