写一个将凯西西储大学轴承数据输入绘出其时域图和频域图的简单python程序
时间: 2024-02-20 11:01:15 浏览: 26
好的,以下是一个使用Python绘制凯西西储大学轴承数据的时域图和频域图的简单程序:
```python
import numpy as np
import matplotlib.pyplot as plt
from scipy.fftpack import fft
# 读取数据
data = np.loadtxt("kcxy_data.txt")
# 绘制时域图
plt.subplot(2, 1, 1)
plt.plot(data[:, 0], data[:, 1])
plt.xlabel("Time (s)")
plt.ylabel("Amplitude")
plt.title("Time Domain")
# 计算频域
n = len(data)
k = np.arange(n)
T = n / 2000
frq = k / T
frq = frq[range(n//2)]
Y = fft(data[:, 1]) / n
Y = Y[range(n//2)]
# 绘制频域图
plt.subplot(2, 1, 2)
plt.plot(frq, abs(Y))
plt.xlabel("Frequency (Hz)")
plt.ylabel("Amplitude")
plt.title("Frequency Domain")
plt.show()
```
其中,我们假设数据文件名为“kcxy_data.txt”,数据文件应该是一个两列的文本文件,第一列是时间,第二列是振幅。
在程序中,首先使用`numpy`库读取数据文件。然后,我们使用`matplotlib`库绘制时域图,使用`scipy.fftpack`库计算数据的FFT,最后绘制频域图。
希望这个程序能够帮到你!