python能量谱代码
下面是一个计算信号的能量谱的Python代码示例:
```python
import numpy as np
import matplotlib.pyplot as plt
# 生成示例信号
fs = 1000  # 采样频率
t = np.arange(0, 1, 1/fs)  # 时间范围
f1 = 10  # 第一个频率成分
f2 = 20  # 第二个频率成分
x = np.sin(2*np.pi*f1*t) + np.sin(2*np.pi*f2*t)  # 信号
# 计算能量谱
freqs, Pxx = plt.psd(x, NFFT=1024, Fs=fs, pad_to=1024, detrend='linear')
# 绘制能量谱图
frequency函数计算频数plt.plot(freqs, 10 * np.log10(Pxx))
plt.xlabel('Frequency (Hz)')
plt.ylabel('Power Spectral Density (dB/Hz)')
id(True)
plt.show()
```
上述代码首先生成了一个包含两个频率成分的示例信号。然后使用`plt.psd`函数计算信号的能量谱,其中`NFFT`表示FFT的窗口大小,`Fs`表示采样频率,`pad_to`表示FFT的长度,`d
etrend`表示去趋势项的方法。最后使用`plt.plot`函数绘制能量谱图,并调整坐标和显示网格。

版权声明:本站内容均来自互联网,仅供演示用,请勿用于商业和其他非法用途。如果侵犯了您的权益请与我们联系QQ:729038198,我们将在24小时内删除。