求大佬解答一下

import numpy as np import matplotlib.pyplot as plt import math def jiankaixian(m,z,x,n): # m：模数，z:齿数，x：变位系数，n：均布多少个点 m=m z=z alpha=math.radians(25); #压力角 ha=1; c=0.25; x=x; #齿顶高系数;顶隙系数;变位系数 r= z * m / 2 #分度圆半径 ra=r+(ha+x)m #齿顶圆半径 rf=r-(ha+c-x)m #x齿根圆半径 rb=r * np.cos(alpha) #基圆半径 p=np.pim; #齿距 s=p/2+2xmnp.tan(alpha);#齿厚 theta1 = np.tan(np.arccos(rb/r)) - np.arccos(rb/r) theta2 = s/(2*r)

#第1段 渐开线
alphaK =np.linspace(0, np.arccos(rb/ra), n)    #分割点数

thetaK=np.tan(alphaK)-alphaK  - theta1 -theta2;   #极坐标系角度
rK=rb/np.cos(alphaK)   #极坐标系半径

xx = rK *np.cos(thetaK)#转为直角坐标系
yy = rK *np.sin(thetaK) 

yy_fu = -yy            #关于x轴对称


return xx,yy_fu

xx,yy = jiankaixian(1.5,44,0,10000) # 换曲线时这个函数记得换

plt.plot(xx,yy) plt.show() 大佬您能表达出来这个齿面法向的方向角吗