1.分析结果
包络图是第三行。第一行水平时域信号,第二行垂直时域信号,第三行对第1行信号在故障频点包络计算后的结果。
上方是最终的视图。右下角我把光标对准了低频位置的那个故障峰,与理论值是一致的.
2.源码
2.1 测试例程源码
import gp_xjtu_data
import gp_calc_envelope
import gp_xjtu_tools
def test_gp_calc_envelope():
(x, sig_h, sig_v, fullpathOfFile) = gp_xjtu_data.get_demo_data();
descToDraw = []
descToDraw.append(x);
descToDraw.append(sig_h);
descToDraw.append("Time Domain Signal (XJTU-SY Bearing_Datasets 35Hz12KN 1-1)")
descToDraw.append("Time(s) 256000sa/s, 32768pt.")
descToDraw.append("g(acc)");
descToDraw.append(fullpathOfFile + " horizon");
signals = []
signals.append(descToDraw.copy());
descToDraw = []
descToDraw.appe