考虑一个由装有振动传感器的电机驱动的六齿轮传动系统，如下图所示。电机轴上的齿轮1与齿轮2啮合，传动比为17:1。最终的齿轮比，即齿轮1和齿轮2和齿轮3和齿轮4之间的比率是51:1。齿轮5，也在电机轴上，以10:1的传动比与齿轮6啮合。电机转速为180rpm，振动传感器采样率为50khz。为了获得只包含齿轮5和齿轮6的啮合组件的信号，通过指定它们的齿轮比为17和51英寸，滤除齿轮1和2和，3和4的信号组件订单．与轴旋转(order = 1)相对应的信号分量总是隐式地包含在计算中。

RPM = 180;fs = 50e3;t =（0：1 / fs：（1/3）-1 / fs）';%样品时间OrderList = [17 51];f = rpm / 60 * [1 orderlist 10];

在实践中，您将使用诸如从加速度计获得的振动信号的测量数据。对于此示例，生成TSA信号X，这是安装在电机上的振动传感器的模拟数据。

X = sin2 *f(1)*t + sin2 *f(1)*t +．..%电机轴转动及谐波3*sin(2* f(2) t) + 3*sin(2* f(2) t) +．..%齿轮啮合振动和谐波齿轮1和24*sin(2* f(3)*t) + 4*sin(2* f(3)*t) +．..%齿轮啮合振动和谐波齿轮3和42 * sin(2 *π* 10 * f (1) * t);齿轮5和齿轮6的啮合振动%

使用采样时间、转速和要过滤掉的网格顺序计算TSA信号的残差。

Y = tsaresidual (X, t, rpm, orderList);

输出Y是包含齿轮网信号和齿轮5和6的谐波的矢量。

在图中可视化残差信号，原始TSA信号及其幅度谱。

Tsaresidual（X，FS，RPM，OrderList）

从振幅谱图中观察到以下分量:

在这个例子中,sineWavePhaseMod.mat包含相位调制正弦波的数据。XT是否有正弦波数据和时间表rpm.使用的是60 rpm。正弦波具有32 Hz的频率，并滤除未经调制的正弦波，使用32作为订单．

加载数据和所需的变量。

加载（“sineWavePhaseMod.mat”，“XT”，“转”，“订单”)头(XT, 4)

ans =4×1时间表时间数据______________ _______ 0 sec 0 0.00097656 sec 0.2011 0.0019531 sec 0.39399 0.0029297 sec 0.57078

请注意时间值XT是严格递增的、等距的、有限的。

计算残差信号及其幅度谱。设置价值'领域'来“频率”由于订单处于Hz。

[Y S] = tsaresidual (XT, rpm命令,'领域'，“频率”）

Y =1024×1时间表时间数据______________ _________ 0秒2.552e-15 0.00097656秒0.051822 0.0019531秒0.10116 0.0029297秒0.14566 0.0039062秒0.18317 0.0048828秒0.21188 0.0058594秒0.23039 0.0068359秒0.23776 0.0078125秒0.2336 0.0087891秒0.21803 0.0097656秒0.19174 0.010742秒0.1559 0.011719秒0.11215 0.012695秒0.062503 0.013672SEC 0.0092782 0.014648秒-0.045032⋮

S =1024×1复杂-0.0000 + 0.00000 i -0.0000 + 0.00000 i -0.0000 + 0.00000 i -0.0000 + 0.00000 i -0.0000 + 0.00000 i⋮

输出Y是一个包含剩余信号的时间表，即相位调制信号，而年代是一个包含残差信号振幅谱的向量吗Y．

在这个例子中,sineWaveAmpMod.mat包含调幅正弦波的数据。X是一个矢量，其振幅调制的正弦波数据在轴转速为60rpm时获得。未调制的正弦波的频率为32赫兹，振幅为1.0单位。

加载数据，绘制调幅TSA信号的残余信号X．为了获得剩余信号，指定频率为32hz in，滤除未调制的正弦波订单．设置价值'领域'来“频率”．

加载（“sineWaveAmpMod.mat”，'X'，“t”，“转”，“orderList”) tsaresidual (X, t, rpm, orderList,'领域'，“频率”）;

从图中观察残差信号和原始信号的波形和幅值谱。