matlab中sita未定义,未定义与 'double' 类型的输入参数相对应的函数 'int'。

z1 = 31;                %齿数

z2 = 65;

m = 3;                 %模数

a = deg2rad(20);                %齿轮压力角

inva = tan(a)-a;

c = 0.25;              %顶隙系数

ha = 1;                  %齿顶高系数

beta = 20;               %啮合线与x轴的夹角

E = 206e9;              %弹性模量

v = 0.25;                 %泊松比

L = 60/1000;                   %齿宽

r_1 = m*z1/2;         %分度圆半径

r_2 = m*z2/2;

h_f = 1/4;

r_b1 = m*z1*cos(a)/2;

r_b2 = m*z2*cos(a)/2;           %基圆半径

gamma = a:0.01:pi/2;

tao = deg2rad(10):0.01:deg2rad(20);

tao_C = deg2rad(10);

r_p = c*m/(1-sin(a));

b_1 = pi*m/4+ha*m*tan(a)+r_p*cos(a);

a_1 = (ha+c)*m-r_p;

phi_1 = (a_1./tan(gamma)+b_1)./r_1;

phi_2 = (a_1./tan(gamma)+b_1)./r_2;

sita_b1 = pi/(2*z1)+inva;         %基齿角的一半

sita_b2 = pi/(2*z2)+inva;

x_beta1 = r_b1*((beta+sita_b1)*cos(beta)-sin(beta));

y_beta1 = r_b1*((beta+sita_b1)*sin(beta)-cos(beta));

x_beta2 = r_b2*((beta+sita_b2)*cos(beta)-sin(beta));

y_beta2 = r_b2*((beta+sita_b2)*sin(beta)-cos(beta));

x1 = r_1.*sin(phi_1)-(a_1./sin(gamma)+r_p).*cos(gamma-phi_1);

y1 = r_1.*cos(phi_1)-(a_1./sin(gamma)+