Реализация

Простейшая реализация - дифференцирующая цепочка, разность между текущим и предыдущим отсчетом.

Более сложная реализация заключается в имитации отклика цепочки.

В работе "Siddharth Mathur, Vocal-Tract Modeling: Fractional Elongation of Segment Lengths in a Waveguide Model With Half-Sample Delays" приводится следующее разностное уравнение:

где:

Fm - значение давления прямой волны.

Fprev - значение давления прямой волны на предыдущем шаге.

Pprev - значение выходного давления на предыдущем шаге.

r - радиус излучающего поршня, м

Код см. Бегущие волны - Излучение.

Замечание:

Есть некоторое сомнение в правильности формулы для L'. Мне кажется, что надо домножить на π.