Отправлено bp 06 апреля 2006 г. 15:15

Задача вот такая вырисовывается
Fs = 10000 Hz
У сигнала полезная полоса примерно от 30 до 2000 гц
Необходимо произвести сдвиг на 90 градусов, но чтобы в этой полосе амплитуда не пострадала (неравномерность АЧХ не более 0,1%)
Попробовал Гильберта в fdatool, но вот требование по амплитуде не удовлетворить.
Есть еще какие-нибудь другие способы ?

• Вот тут поискал вручную (движок форума в корзину небытия) и нашёл то, что мне подсказвали. (+) — уни   (06.04.2006 19:37 81.91.36.3, 773 байт, ссылка)
  • спасибо, сейчас посмотрю — bp   (07.04.2006 10:11 81.3.189.234, пустое)
• небольшое дополнение (+) — bp   (06.04.2006 18:10 81.3.189.234, 398 байт)
  • Может тогда так?+ — andy_P   (07.04.2006 14:11 80.82.63.185, 300 байт)
    • Да, это интересно. Сейчас попробую промоделировать. — bp   (07.04.2006 14:29 81.3.189.234, пустое)
  • Так может, сигнал u2 у Вас узкополосный, и соответственно, интересует именно его корреляция с сигналом U2? Тогда постоянный сдвиг фазы по всему спектру и не нужен, достаточно сдвинуть в полосе сигнала u2... — Sokol   (06.04.2006 18:27 62.118.172.222, пустое)
    • Там где "U2" читать "U1" :) — Sokol   (06.04.2006 18:29 62.118.172.222, пустое)
      • у них у обоих полоса одинаковая — bp   (06.04.2006 21:49 83.143.195.14, пустое)
• Возьмите ИХ длиной ~6000 отсчетов и будет Вам счастье (только что проверил на матлабе - метод синтеза - remez, полоса от 0 до 2.2 , Fs=10, вес 1). Или несчастье :-) — andy_P   (06.04.2006 16:39 80.82.63.185, пустое)
  • да, я это пробовал, примерно столько же получилось — bp   (06.04.2006 18:02 81.3.189.234, пустое)
• Ответ: — TAP_   (06.04.2006 16:32 83.149.201.129, 81 байт)
  • Вы имеете ввиду что я упираюсь в физический предел ? — bp   (06.04.2006 18:01 81.3.189.234, пустое)
    • нет — TAP_   (07.04.2006 09:26 83.149.201.129, 84 байт)
• Ничего не скажу про этот самый tool, который fda, за его неимением, но в таких случаях принято играть с порядком, в самых отчаянных случаях начинают играть с четностью нечетностью длины, наряду с порядком, и извращаться с методами рассчета(окна, Ремез,...) — -=ВН=-   (06.04.2006 16:11 193.125.71.140, пустое)
• Если бы это надо было сделать 1 раз :) , то можно вычислить коэффициенты спектра в алгебраической форме и поменять их местами, при мнимой части знак поменять вроде и ОБПФ - вот без среза НЧ, а по делу - ... — уни   (06.04.2006 16:09 87.224.180.187, 403 байт)
  • -j на положительных, +j на отрицательных - вся ЧХ ПГ:-) ) — -=ВН=-   (06.04.2006 16:25 193.125.71.140, пустое)
    • Я как-то делал нечто похожее ... по ссылке — уни   (06.04.2006 16:32 87.224.180.187, пустое, ссылка)
      • Ой, форум... Да в этом деле и обсуждать нечего:-) — -=ВН=-   (06.04.2006 16:37 193.125.71.140, пустое)
        • Да это я для красоты... почитал учебник... про Гильберта, там формулы.. ну, думаю, хочу в живую посмотреть, вот и посмотрел... а форум - картинку больше некуда выложить. А фильтры Гильберта - это конечно отдельная тема. — уни   (06.04.2006 16:40 87.224.180.187, пустое)
          • А-а-а, да не в ней ничего отдельного, кроме ИХ-ЧХ:-). Фильтры, как фильтры. Подчиняются всем фильтровым законам:-) — -=ВН=-   (07.04.2006 10:21 193.125.71.140, пустое)
            • Я про методы проектрования ФГ, вон выше привёл ссылку. Как только нужно нечто специальное, то общие утилиты расчёта коэффициентов, как тут оказывается, не удовлетворяют потребности. — уни   (07.04.2006 10:34 87.224.180.187, пустое)
              • Да ладно. Все абсолютно удовлетворяет:-) — -=ВН=-   (07.04.2006 10:42 193.125.71.140, пустое)
                • :-) А-а-а, ну тогда забыли. Пусть автор решает. — уни   (07.04.2006 10:45 87.224.180.187, пустое)
                  • А чего забыли-то? Борьба с узкими перех. полосами только одна - увеличение порядка. И все сведется к этому. Или автору надо переформулировать задачу. — -=ВН=-   (07.04.2006 10:52 193.125.71.140, пустое)
                    • Где я спорил то? :) Я просто сам хотел синтезировать фильтр по Frequency-Response Masking Technique и посмотреть, чем это отличается от типовых приёмов с увелич. частоты дискр. или изм. длины или уже замены ТЗ. Не получилось у меня... пока. Я знаю, что теория нерушима, но ведь сколь угодно близко к идел. ФГ можно разыми способами приближаться. Я так думаю. — уни   (07.04.2006 10:57 87.224.180.187, пустое)
                      • Ну не спорили и славно. А про типовые приемы с увеличением частоты дискретизации Вы сами додумались?:-) Интересна ведь не полоса в герцах, а относительная полоса, отнесенная к частоте дискретиз. Так что ее увеличением Вы добровольно отстригаете себе гениталии. А методы раасчета - Ваш, с импортным назвнием, алгоритм Ремеза, просто оконный способ - не так сильно между собой отличаются. Это уже от безысходности, когда и их начинают перебирать. — -=ВН=-   (07.04.2006 11:10 193.125.71.140, 160 байт)
                        • В моём случае (по ссылке, а потом наверх) случай был для меня безысходный почти, пока SM посоветовал комплексн. фильтры. — уни   (07.04.2006 11:20 87.224.180.187, пустое)
                          • Так елки-палки, ну это же абсолютно очевидно, что коплексный полосовой фильтр (при должных требованиях к полосе) служит оч. хорошим заменителем ПГ, выравнивателю ГЗ, при том условии, что преобразование Гильберта нужно в полосе. — -=ВН=-   (07.04.2006 11:30 193.125.71.140, пустое)
                            • Дык Вам очевидно, а я только на свет можно сказать народился. Мне это в голову не приходило и не учили :( навёрстываю вот. — уни   (07.04.2006 12:30 87.224.180.187, пустое)
                        • Ну не надо так кипятиться. Я всегда имею в виду относительность, т.к. частотный план задаю относительным при расчёте коэф-тов. А что я отрезаю (но не себе) я итак знаю, т.к. это приводит к плохому подавл. боковой полосы, например, при дальнейшем переносе на ВЧ методом ОБП. Я всю эту фигню моделировал и промучился вполне достаточно, думаю, чтобы знать о чём говорю и что при этом Вам можно было и предположить в хорошем смысле. Вот. — уни   (07.04.2006 11:17 87.224.180.187, пустое)
                          • Да разве я кипячусь? Я Вас за язык, что ли, тянул, заставлял говорить о типовости повышения частоты дискр?:-) Вы только расскажите, что я должен предполагать в хорошем смысле, и я обязательно буду в этом смысле предполагать:-) — -=ВН=-   (07.04.2006 11:26 193.125.71.140, пустое)
                            • Во-первых, меня на гениталии не разведёшь, если бы Вы сходили бы по ссылке, то не постили бы сюда про "безысходность" - её вообще нет, есть просто отсутсвие знания, либо физическая (либо программная) нереализуемость. Во-вторых, о догадливости и предположениях. Если Вы сделали предположение ... — уни   (07.04.2006 12:52 87.224.180.187, 1628 байт)
                              • Ну не хватало мне еще по ссылкам лазить:-) Нет никаких предположений о безисходности. В не знаю какой раз повторю - все решают полосы, они определяют порядок, метод рассчета в этом деле вещь десятая. Если он не совсем идиотский. Ну даст один метод 3088 порядок, а второй 4012 - много легче будет? По комплексным фильтрам - о них, во-первых, додуматься не велик труд. Во-вторых, пишут на каждом углу, в этой конфе тоже, мной в том числе и неоднократно:-). Захотите - найдете. А игнорировать я могу, — -=ВН=-   (07.04.2006 13:27 193.125.71.140, пустое)
                                • 3088~=3988:-)) — -=ВН=-   (07.04.2006 13:51 193.125.71.140, пустое)
                                • Конец отстригло, ну да бог с ним:-) — -=ВН=-   (07.04.2006 13:28 193.125.71.140, пустое)
                                  • Вот говорят, что случаный сигнал сжать не возможно, т.к. в пространствах, в которых работают сейчас алгоритмы сжатия с потерями (в т.ч. вейвлеты) эти самые случайные сигналы такими же и остаются. От того и бытует такое мнение. (+) — уни   (07.04.2006 14:33 87.224.180.187, 2367 байт)
                                    • По порядку. Готов не поскупиться и дать Вам даже 3800-4200, и сомневайтесь на здоровье:-) — -=ВН=-   (07.04.2006 15:20 193.125.71.140, 1002 байт)
  • А...Вы про неравномерность, тогда да... я только "играл формулой частотного плана" и порядком. — уни   (06.04.2006 16:17 87.224.180.187, пустое)
• Продифференцировать. — geekkoo   (06.04.2006 16:06 194.67.79.1, пустое)
  • А если не поможет, то проинтегрировать :-) — AntZ   (06.04.2006 16:46 212.5.70.188, пустое)
  • Муддро. Причем на каждой частоте индивидуально. — -=ВН=-   (06.04.2006 16:07 193.125.71.140, пустое)
• Учение Гильберта всесильно, потому что оно верно — AntZ   (06.04.2006 15:53 212.5.70.188, пустое)
  • с этим никто не спорит, просто не знаю как его адаптировать под свою задачу — bp   (06.04.2006 16:01 81.3.189.234, пустое)

E-mail: info@telesys.ru