Отправлено Newka 21 апреля 2003 г. 10:37

Я в данный момент делаю интерполяцию (повышение Фдискр2=Фдискр1*N в целое число) путем вставления нулей между отсчетами смгнала взятыми с Фдискр1 и последующей фильтрацией ФНЧ с частотой среза немного больше Фмах в исходном сигнале. Это чтобы порезать спектр на частотах вокруг Фдискр1*К. Так вот я заметил что при увеличение N мне требуется больший порядок ФНЧ для уменьшения искажения в интерполированном сигнале. Вопрос1: из-за чего такое? предположительный ответ - больше спектров при ув. N начинают залазить на мой сигнал. Вопрос2: что применить ? Полифазные фильтры? или попробовать кубическую (полиномами)интерполяцию. Объясните что нить про это :) или линки плс.

• ________***** Во первых всем Спасибо! Во вторых исщо вопрос. ****________ — Newka   (22.04.2003 10:19, 620 байт)
  • Пример хорошо реализуемой интерполяции кубическим сплайном по ссылке. Применяется например при передискретизации изображений. Да, надо хранить четыре отсчета, и получать результаты меж двух средних точек. — SM    (22.04.2003 10:29, пустое, ссылка)
    • Кстати, КИХ тоже имеет задержку. И не на 1-2-3 сэмпла, а на пол-ИХ. — SM    (22.04.2003 10:31, пустое)
      • Кстати о КИХах, какой максимальный порядок фильтра можно делать на 16 битных коэффициентах — Newka   (22.04.2003 11:10, 248 байт)
        • Вообще квантование коэффициентов считается грубо так. — st256   (22.04.2003 11:49, 437 байт)
        • Такой порядок - как правильно сказал TAP каскадами. Или блочно через FFT. Хотя c FFT все сильно зависит от требуемой ИХ фильтра. — SM   (22.04.2003 11:46, пустое)
        • Ответ: а зачем такого порядка Не проще ли делать многокаскадные фильтры — TAP   (22.04.2003 11:13, пустое)
• Ответ: Polifaznii mnogokaskadnii filter (+) — TAP   (21.04.2003 22:03, 269 байт)
• А не проще ли считать, что сигнал представляет самый нижний уровень вэйвлет-преобразования и затем его интерполировать с помощью него обратного? На каждом этапе длина фильтра сохраняется и спектр сигнала вполне съедобно сохраняется тоже. — Вадим Кудрявцев   (21.04.2003 19:20, пустое)
• Смотрите ссылку. — [забывший пароль]GroundCtrl   (21.04.2003 15:14, 45 байт, ссылка)
• Ответ1: Растет добротность фильтра и пропорционально растет длина КИХ фильтра. Ответ2. Можно попробывать интерполяцию по Котельникову (sinx/x). А если забивать, то хотя бы предыдущим значением, а не нулем. Порядок фильтра НЧ должен упасть. Возможны комбинации забития каким-нибудь многочленом и слабеньким ФНЧ. — st256   (21.04.2003 11:37, пустое)
  • Ответ: — [забывший пароль]GroundCtrl   (21.04.2003 15:50, 559 байт)
    • Ответ: — AntZ   (21.04.2003 17:16, 568 байт)
  • если забивать не нулём то z преобразование не будут равны. — Newka   (21.04.2003 12:52, пустое)
    • Кроме того спектры будут убывать к Fd/2. Ну и что с того? — st256   (21.04.2003 13:05, пустое)
      • Сигнал будет искажен, то есть в спектре сигнала будут изменения. Можно конечно повторять одно и то-же значение, несколько раз и фильтровать. Но это будет другой сигнал. Простой пример в маолабе - 16 самплов синуса от 0 до 2pi. После upsampilng по методу st256 (повторяем каждый сампл два раза) в спектре сигнала от 0 до pi появляются новые частоты. При забивании нулями, идет только алиасинг на частоты >pi. — AntZ   (21.04.2003 13:18, пустое)
        • про забивание нулей — [забывший пароль]GroundCtrl   (21.04.2003 16:01, 371 байт)
          • Вы попробуйте сделать это, например в матлабе. Вы увидите, что в спектре появляются высокие частоты, которых не было. Жаль не могу привести картинку из Matlab. Полоса сигнала меняется. Думаю лучше всех это может разяснить BH. Как я понимаю, Вы хотите сказать, что фазовая характеристика меняется, а частотная нет. Несколько строчек на матлаб могут опровергнуть Вашу гипотезу. — AntZ   (21.04.2003 16:13, 176 байт)
            • Ответ: Опоздал. — ВН   (21.04.2003 18:29, 71 байт)
              • Пока выяснилось следующее. — st256   (21.04.2003 18:46, 309 байт)
                • Ответ: Вот ведь.... — ВН   (21.04.2003 19:58, 2206 байт)
                  • Как все просто. Спасибо за объяснение. Мне становится стыдно, что я не распознал достаточно простую операцию фильтрации, это же очевидно, что IR фильтра растянутый единичный импульс. — AntZ   (22.04.2003 10:13, пустое)
                    • Блин, Антоша! Это - ЦПШ..... — st256   (22.04.2003 10:31, пустое)
                  • Я попробую блеснуть кучей знаний :) — Newka   (22.04.2003 10:08, 224 байт)
                  • Я не злопамятный. Просто я злой и память у меня хорошая. — st256   (22.04.2003 05:10, 2376 байт)
                    • Ответ: — ВН   (22.04.2003 12:01, 930 байт)
                      • Не въехали? Повторяю — st256   (22.04.2003 12:38, 771 байт)
                • Меня Ваше предложение подвигло на эксперименты, в результате которых я узнал кое-что новое. Чего греха таить, я тоже первое что подумал - зачем добивать нулями, а не линейно интерполировать, когда впервые столкнулся с интерполяцией. То что иногда интуитивно кажется хорошим, часто таковым не является — AntZ   (21.04.2003 18:57, пустое)
                  • Антоша, ну зачем цепляться к человеку, если человек всячески Вас игнорирует? Придется показать Вам зубы. — st256   (22.04.2003 05:40, 1322 байт)
                    • Безусловно, я ошибся что появляются новые составляющие, но я был прав, что спектр искажается. Ваше письмо, появилось после обстоятельного ответа BH, который справедливо поставил на место нас обоих. Меряться ... по поводу образования - не будем. — AntZ   (22.04.2003 10:09, пустое)
                      • Правильно, Антоша. Вам мериться нечем. Поэтому Вы и пишите, что ВН поставил на место нас обоих. В выше означеном случае он лишь завизировал Ваши "дополнения" к Котельникову. И показал определенный навык дисскуссии. — st256   (22.04.2003 10:16, пустое)
                    • Мда... У меня вот нету В/О - ну и что? :-)))) — SM    (22.04.2003 09:54, пустое)
                      • Да ничего. Мой учитель, написавший книгу "Цифровые РПУ" Фима Побережский, тоже закончил что-то не профильное. — st256   (22.04.2003 10:27, 425 байт)
              • Интересно узнать у BH, есть ли какие либо преимущества интерполяции "повторением" самплов? — AntZ   (21.04.2003 18:32, пустое)
            • Ещё раз внимательно прочитайте мой ответ. — [забывший пароль]GroundCtrl   (21.04.2003 16:51, 268 байт)
              • Блин, не туда нажал... — [забывший пароль]GroundCtrl   (21.04.2003 17:10, 369 байт)
                • Понял. Просто Вы криво выразились, а я криво понял :-))) — AntZ   (21.04.2003 17:19, пустое)
                  • В интернет-баталиях такое часто случается :))) — [забывший пароль]GroundCtrl   (21.04.2003 17:30, пустое)
              • Тогда поясните фразу "Т.е полоса сигнала остаётся прежней, но сам сигнал, после фильтрации, получается сильно искажённым" и что Вы имете в виду под полосой сигнала. Если сигнал является "сильно искаженным", то следовательно у него меняется либо амплитуда от частоты, либо фаза от частоты, либо обе характеристики. При этом "полоса сигнала" остается неизменной?. — AntZ   (21.04.2003 17:05, 51 байт)
                • Ответ: — [забывший пароль]GroundCtrl   (21.04.2003 17:23, 788 байт)
                  • Из Вашего ответа, следовало, что Вы говорите о методе st256. При upsampling'e методом st256 сигнал портится так, что восстановить его становится очень проблематично... Я не вижу никих выгод, и наоборот массу недостатков интерполяции путем "повторения" самплов. — AntZ   (21.04.2003 17:29, пустое)
        • No comment. — st256   (21.04.2003 13:52, пустое)
  • Забивание предыдушим значением - это новое слово в науке. В качестве эффективного решения передлагается линейная или нелинейная интерполяция. Это все сводится к переходу цифра->аналог->цифра. Если устраивает решение "digital-only", то фильтрацию после upsample часто выгодно заменить на полифазный фильтр — AntZ   (21.04.2003 12:03, пустое)
• Я делал недавно ресамплинг с 9600 на 8000. (upsampling 5, downsampling 6). Есть модель в симулинке, могу бросить. Полифаза, как мне и рекомендовали здесь пару недель назад, получилась в разы эффективней. — AntZ   (21.04.2003 10:51, 253 байт)
  • Залейте сюда плс — Newka   (21.04.2003 11:16, пустое)
    • Ловите, послал — AntZ   (21.04.2003 11:32, пустое)

E-mail: info@telesys.ru