我试图计算最大值处的宽度,但使用peak_widths时,它返回信号内的所有峰值宽度。
from scipy.signal import chirp, find_peaks, peak_widths
import matplotlib.pyplot as plt
x = np.linspace(0, 6 * np.pi, 1000)
y = np.sin(x) + 0.6 * np.sin(2.6 * x)
peaks, _ = find_peaks(y)
results_half = peak_widths(y, peaks, rel_height=0.5)
results_half[0]
plt.plot(y)
plt.plot(peaks, y[peaks], "x")
plt.hlines(*results_half[1:], color="C2")
plt.show()
这是来自 scipy.signal.peak_widths 的示例。
我只想得到最大值的宽度。
我想用最大峰位置的索引来找到最大峰的FWHM,如何找到最大值对应的宽度?
i_max_peak = peaks[np.argmax(y[peaks])]
y_max = y[i_max_peak]
我想得到最大峰值的宽度值。
1条答案
按热度按时间wxclj1h51#
我想你已经成功了一半。
您已经通过
[np.argmax(y[peaks])]
得到了最大峰值的索引,剩下的是通过results_half[0]
得到宽度。所以,我相信,这会给予你你想要的: