斜齿轮啮合刚度matlab程序，势能法，有偿

求一个如下类似的斜齿轮啮合刚度matlab程序，势能法，有偿
参数自设，求一个求刚度的程序，标注清晰

参考文献：齿轮磨损故障动态响应特征与诊断指标研究