clear allclcN = 1000;T = 02piN2pi;r = 50; 基圆半径为50h = 5; 齿高z1 = 6; 齿数z2 = z1;d = 100; 中心距ri_1 = dz1z1+z2; 节圆半径h1 = 08; 小齿齿高m = 6; 小齿齿数z_1 = z1m;z_2 = z_1; 极坐标系下的齿廓l = h1sinz1mT;u = hsinz1T + r;c

注：本程序为MATLAB代码，需要在MATLAB环境下运行。

本程序演示了如何通过极坐标系下的齿廓，计算出齿廓和节圆的交点，并旋转得到节点时的啮合点。然后利用共轭原理，通过计算共轭齿廓方程，生成两个齿轮的共轭齿廓并绘制出来。

程序中的参数可以根据实际情况进行修改，具体含义如下：

• N：极坐标系下绘制齿廓的点数。
• T：极坐标系下的角度序列。
• r：基圆半径。
• h：齿高。
• z1：齿数。
• z2：齿数，与z1相等。
• d：中心距。
• ri_1：节圆半径。
• h1：小齿齿高。
• m：小齿齿数。
• z_1：小齿轮齿数，等于z1*m。
• z_2：大齿轮齿数，等于z_1。
• y1：原始齿轮的共轭齿廓方程。
• y2：共轭齿轮的共轭齿廓方程。

程序中采用了MATLAB自带的函数fzero对齿廓和节圆的交点进行数值求解，请注意选择合适的初始值。