频谱分析_滚动轴承的频谱分析代码

滚动轴承的频谱分析代码需要用到信号处理的知识，以下是一个简单的示例代码：

import numpy as np
import matplotlib.pyplot as plt

# 读取数据
data = np.loadtxt('data.txt')

# 定义采样频率
fs = 10000

# 傅里叶变换
freqs = np.fft.fftfreq(len(data), 1/fs)
fft = np.fft.fft(data)

# 取振动信号频率范围内的振动信号分量
freq_range = (0, 2000)
mask = (freqs >= freq_range[0]) & (freqs <= freq_range[1])
freqs = freqs[mask]
fft = fft[mask]

# 绘制频谱图
plt.plot(freqs, np.abs(fft))
plt.xlabel('Frequency (Hz)')
plt.ylabel('Amplitude')
plt.show()

其中，`data.txt` 是存储振动数据的文本文件，每行为一个采样点的振动信号值。`fs` 是采样频率。通过傅里叶变换将时域信号转换为频域信号，取出振动信号频率范围内的分量，并绘制频谱图。可以根据实际情况调整频率范围和绘图参数。