- 大家好。我已经检测了数百个这种形式的信号，其中的峰值超过了某个阈值。
我定义峰宽为FWHM（半高全宽）。然而，我已经对信号的谷部拟合了一个三次多项式，因此我定义峰值为从基线到峰值的距离，基线与峰值在同一索引处。
我计算峰宽为信号与最大值一半的线相交的最大距离。看起来是这样的：
roots = indeces_of_intersectionsintersection_lengths=[abs(y - x) for x, y in it.combinations(roots, 2)]calculated_width = max(intersection_lengths)我遇到的问题是无法一致地计算峰宽，这是因为有时线会与不同峰上的点相交。
我已经限制了这条相交线定义的域：
域 = [峰值左侧一点，峰值右侧一点]
但这个域限制对所有峰值都是相同的。
我考虑过如何让这个域根据不同的峰值进行变化，但不确定如何实现。我的代码几乎完全自动化，我必须保持这种状态。
回答：
将我的问题发布在这里帮助我意识到一个简单的解决方案：
more_than_peak=[x for x in it.ifilter(lambda x: x if x>peaks[i] else 0, roots)]less_than_peak=[x for x in it.ifilter(lambda x: x if x<peaks[i] else 0, roots)]if len(more_than_peak)>0 or len(less_than_peak)>0: width = min(more_than_peak)-max(less_than_peak)在这里，我找到峰值索引左侧和右侧的交点。然后，我找到峰值右侧最小的交点和左侧最大的交点（x轴向右增加）。如此简单且快速！
