Отправлено R2D2 23 ноября 2005 г. 10:25 В ответ на: Нет, не безобразие. Есть ф-ция polystab, которая полюса на место загоняет, не меняя АЧХ, но корежа ФЧХ. А то ведь (+) отправлено SM 23 ноября 2005 г. 01:34

Можно ж было вставить в алгоритм условие устойчивости. Ну нет подходящего фильтра, значит нет. Зачем чо попало выдавать?

• Как это "чо попало" - эта ф-ция и не предназначена для синтеза устойчивых фильтров. Это просто преобразование, обратное к freqz. И более ничего. — SM   (23.11.2005 10:53 213.141.159.26, пустое)
  • А почему все-таки не хотите синезировать all-pass фильтр с заданной ФЧХ прямо в частотной области? Просто и со вкусом. Потом и обработку в частотной области можно сделать будет с использованием overlap-add или overlap-save приема... Или задержка не устраивает? — homekvn   (23.11.2005 11:44 212.185.161.237, пустое)
    • Не понял про что Вы. Так же я не понял как с помощью iirgrpdelay синтезировать аллпасс с заданной ФЧХ. — R2D2   (23.11.2005 14:42 195.46.182.4, пустое)
      • Попробую объяснить. (+) — homekvn   (23.11.2005 15:30 212.185.161.237, 4207 байт)
        • А про то, как при помощи iirgrpdelay желаемый БИХ-фильтр получить... Ну неужели не догадались? (+) — homekvn   (23.11.2005 15:39 212.185.161.237, 362 байт)
          • Не все! — R2D2   (23.11.2005 15:47 195.46.182.248, 192 байт)
            • И я не знаю. Спросите у MathWorks - они знают. Но вобщем, согласен, данная проблема останется. Но для Гильбертовского фильтра сделана другая функция hilbiir, которую Вы почему-то отшили. — homekvn   (23.11.2005 16:12 212.185.161.237, пустое)
              • А не получается нормальнйы гильберт там для моей задачи и я решил отказаться от условия нулевой АЧХ на нулевой частоте. — R2D2   (23.11.2005 16:52 195.46.191.12, пустое)
                • А Гильберт всегда нулевой на нулевой частоте. Как ни крутитесь. Их у него антисимметричная. Сделайте (+) — SM   (23.11.2005 17:01 213.141.159.26, 582 байт)
                  • Еще, как вариант, разбейте сигнал на полосы, в каждой перейдите на более низкую Fs, и считайте мгновенные мощности отдельно. — SM   (23.11.2005 17:02 213.141.159.26, пустое)
                    • Так это тогда банки фильтров получатся... Тогда проще сразу в частотной области лепить, потому как (+) — homekvn   (23.11.2005 17:24 212.185.161.237, 227 байт)
                      • А гадость Вы в любом случае получите. Все дело в ее размере. — -=ВН=-   (23.11.2005 17:48 194.190.181.231, пустое)
                      • Человек хотел мгновенную мощность мерять. Я думаю, что склеить мощности по всем полосам не проблема, как и дефекты Гильбертования поправить, сделав полосы достаточно пересекающимися. — SM   (23.11.2005 17:30 213.141.159.26, пустое)
                        • Ну, может, для мгновенной мощности и сойдет, но так вообще-то вопрос открытым остается: если я сделаю гильберта на пониженной в два раза частоте, в ветке, отвечающей за ВЧ, то что я там увижу, когда потом все склею? — homekvn   (23.11.2005 17:36 212.185.161.237, пустое)
                          • Ну, я же в контексте задачи, описанной товарищем (внизу страницы, после флейма :) )... Про склейку и не думал, если склеивать, по-Вашему наверное лучше выйдет, через БПФ. — SM   (23.11.2005 17:47 213.141.159.26, пустое)

E-mail: info@telesys.ru