@dementiaf2d

Здравствуйте, я только осваиваю matlab, нужно рассчитать КИХ фильтр.

Через канал был передан ключ для декодирования картинки. Ключ искажен при передаче через канал, помехи которого определяются КИХ фильтром 3-го порядка + белый гауссовский шум. Но ключ содержит 32 символа для определения помех, которые вносит канал. То есть известна последовательность из 32 сиволов на входе в канал и на выходе. Нужно расчитать оптимальный фильтр, который отфильтрует шумы канала.

Я пытаюсь определить коеффициенты, решая уравнение (см. картинку), Где первый множитель - матрица Тёплица из автокорреляционной функции выходного сигнала y(n)[32x32], а второй множитель - это кросс корреляция входного сигнала (x(n)) с выходным(по идее столбец из 32 элементов). С ним у меня проблемы: почему-то для двух столбцов из 32 элементов получается вектор из 63 элементов. Вы бы не могли дать совет с чего начать или что почитать. Спасибо, извините за сумбурное изложение вопроса.

Изображения

 

0

@Annza

Вообще странно, да. Хорошо бы, чтобы rxy имела размерность 32х1, тогда и ИХ будет вектором столбцом размерности 32х1. Но с другой стороны то, что взаимная корреляция будет иметь размерность 2х32-1 тоже вроде верно. Единственное, взаимная корреляция и свертка тесно связаны(http://ru.wikipedia.org/wiki/В... ая_функция), а вот в отношении свертки я знаю нюанс - свертка может быть линейной и циклической(http://www.dsplib.ru/content/conv/conv.html), и в случае циклической результат будет как раз нужный 32х1, а не 63х1 как в случае линейной.

1

@dementiaf2d

в общем, я взяла значения xcorr начиная с n=32 до 63, она симметричная относительно n=32. получила коефициенты фильтра (аж 32). Отфильтрованный ключ в итоге расшифровал картинку, но она не идеальна, теперь надо как-то оптимизировать коеффициенты.. И наверное, уменьшить порядок. Потому-что найденный в интернете код, который я не очень понимаю, расчитал фильтр 11 порядка, и декодированная картика выглядела получше с ним.

0

@vital792

это вы правильно сделали. Импульсная характеристика канала на входе фильтра Винера должна быть односторонней.
Если бы вы показали этот код, а также свой код, наверно было бы легче что нибудь сказать по этому поводу...

1

@dementiaf2d

vital792 , спасибо, кажется, мы (с товарищем) уже близки к ответу, если завтра ничего не получится, тогда скину коды.

0