python写一个简单的轴承故障分析代码

以下是一个简单的轴承故障分析代码示例，使用Python编写：

import numpy as np
import matplotlib.pyplot as plt

# 读取轴承振动数据
data = np.loadtxt('bearing_data.txt')

# 绘制振动信号时域图像
plt.plot(data)
plt.title('Time Domain Vibration Signal')
plt.xlabel('Time (ms)')
plt.ylabel('Vibration Amplitude')
plt.show()

# 计算频谱并绘制频谱图像
spectrum = np.fft.fft(data)  # 使用快速傅里叶变换（FFT）计算频谱
freq = np.fft.fftfreq(data.shape[-1])  # 计算对应的频率
plt.plot(freq, np.abs(spectrum))
plt.title('Frequency Domain Vibration Signal')
plt.xlabel('Frequency (Hz)')
plt.ylabel('Magnitude')
plt.show()

# 检测轴承故障
if np.max(np.abs(data)) > 0.1:  # 检测振幅是否超过阈值0.1
    print('Bearing fault detected!')
else:
    print('No bearing fault detected.')

在这个示例中，我们首先从文件中读取轴承振动数据，并绘制其时域图像和频谱图像。然后，我们检测振动信号中是否存在超过阈值的振幅，以判断是否存在轴承故障。这只是一个简单的示例，实际的轴承故障分析需要更复杂的算法和信号处理技术。