Формирование акустического импульса, без эха Опции

[phantom]

Допустим есть у нас трубопровод. На некотором расстоянии от одного конца установлен динамик излучающий импульсы. Около динамика, в трубе, установлен микрофон, который эти имульсы записывает. При попытке посылки импульса получаем резкое искажение формы сигнала записанного микрофоном от излучаемого динамиком (появляются хвосты). Происходит это от того, что во время формирования импульса на него накладываются отражения от ближайшего конца трубопровода. Вопрос: существуют ли адаптивные алгоритмы, способные сформировать заданную форму импульса вне зависимости от характера начальных отражений (поскольку динамик может устанавливаться в трубы разной начальной конфигурации)? Основной критерий формирователя - малая длина полученого импульса.

[phantom]

fontp
Пробовали запихнуть в Вашу программу наши сигналы - тот, что подается на динамик, и записанный в начале трубы, причем в качестве входного пробовали шум, кусок синусоиды и еще несколько вариантов. Проблема в том, что отклик достаточно длинный - пару тысяч отсчетов, так что при небольшом количестве отсчетов ИХ фильтра естественно ничего не получается, а побольше его сделать проблематично - даже для тестовых данных (случайный сигнал на входе и свернутый со случайной ИХ) длиной 10200 и порядком фильтра 1000 fir_covar возвращает код ошибки (-1,-2), а коэффициенты считаются только до 544-го (остальные нули). Чтоб считались все коэффициенты и не было ошибки, нужно, чтобы количество отсчетов сигнала и отклика было раз в 100 больше порядка фильтра (для тестового сигнала, для наших так и не получилось ничего). В то же время из статьи можно понять, что вроде достаточно данных всего в 2 раза больше порядка. То ли есть какие-то нюансы в реализации алгоритма, то ли я идею неправильно понимаю?

[fontp]

fontp
Пробовали запихнуть в Вашу программу наши сигналы - тот, что подается на динамик, и записанный в начале трубы, причем в качестве входного пробовали шум, кусок синусоиды и еще несколько вариантов. Проблема в том, что отклик достаточно длинный - пару тысяч отсчетов, так что при небольшом количестве отсчетов ИХ фильтра естественно ничего не получается, а побольше его сделать проблематично - даже для тестовых данных (случайный сигнал на входе и свернутый со случайной ИХ) длиной 10200 и порядком фильтра 1000 fir_covar возвращает код ошибки (-1,-2), а коэффициенты считаются только до 544-го (остальные нули).

Куда Вам столько много? Это подгонка. Порядок нужно ограничивать в вызове функции.

Обычно порядок увеличивают только до тех пор пока ошибка уменьшается резко. Обычно ошибка снижается сначала более резко, а потом менее резко. На перегибе останавливаются :-)
Дальше - это всё подгонка под шум. Так моделируются FIR-системы, есть там специальная даже стат-наука по выбору порядка.

"Моя программа" это всего лишь быстрый алгоритм решения ковариационных уравнений. Можете выписать нормальные уравнения и решать их любой другой программой для решения линейных систем уравнений. При обычной точности вычислений возможно решать системы порядка несколько десятков, особенно если они недостаточно хорошо обусловлены


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