Из каких значений индексного массива находится определенный диапазон ошибок?

Question

Есть ли встроенная функция numpy для проверки, из какого индекса сигнал (массив) не оставляет определенного диапазона ошибок?

---

Работа с цифровыми фильтрами, мне нужно, чтобы определить длину импульсной характеристики для использования в scipy.signal.filtfilt. Довольно легко с Финитный импульсный отклик (FIR) фильтры, но вроде невозможного с Бесконечный импульсный отклик (IIR) фильтры.
Однако было бы сделать расчет точки, из которой импульсная характеристика не покидает определенную группа ошибок:

На данный момент я использую быстрый и грязный обходной путь, проверяя обращенным массив вручную для первого значения вне диапазона ошибки:

def ringing_time(sig, th): 
    return len(sig) - np.argmax(np.abs(sig[::-1]) > th) 

есть ли быстрый встроенный numpy подход для этого?

Answer 1

В общем, нет. Вы используете некоторые специфические для сигнала знания, которые не универсальны (факт, что огибающая распадается со временем). Ваше решение хорошее, я думаю, если вы хотите сделать это численно.

Вы можете сделать что-то вроде принятия уравнения импульсного отклика, символически дифференцировать его с помощью sympy и применить метод Ньютона, но вам придется иметь дело с его запуском в нужном месте, чтобы он не попадал в один из локальные минимумы. В целом я думаю, что вы в хорошем месте.