Отправлено SМ 11 ноября 2003 г. 23:12 В ответ на: А можно чуть подробнее+ отправлено AndrI 11 ноября 2003 г. 21:07

Под квадратурной дискретизацией я понимаю следующее. Тоже, что и все обычно понимают.

Два АЦП работают на частоте 100 герц. Цифруют одно и то-же. Но со сдвигом по фазе друг относительно друга на 90 градусов. Их сигналы назовем a и b. Вместо двух АЦП естественно применить один, но выборки с него идут так: 0ms->a, 2.5ms->b, 10ms->a, 12.5ms->b и так далее. Раз в 10 ms собираем из а и b комплексный сигнал x = a+j*b. Аналоговая предфильтрация необходима потому как в оцифровываемом сигнале не должно быть составляющих вне диапазона 50..150 герц. Они вызовут спектральные наложения. Поэтому я предложил для примера полосу 75..125. А вообще-то, как заметил SVI - чем уже, тем лучше. Именно в аналоге. Такая хитрая дискретизация осуществляет гетеродинирование входного сигнала с "несущей" частотой 100 герц. Т.е. перенос спектра. Затем необходимо произвести децимацию (прореживание) для перехода к новой частоте дискретизации - к 20 герцам (Ваши раз в пять периодов). Для этого сам бог велел применить CIC-фильтр подходящего порядка. На выходе дециматора получится опять-же комплексный сигнал x1=a1+j*b1=decimate(x). Этот, прореженный, с выборками с частотой уже в 20 герц, комплексный сигнал надо пропустить через ФНЧ с частотой среза 0.33 герца. С рассчетом этого ФНЧ надеюсь Вы справитесь. С выхода ФНЧ получится окончательный сигнал. Опять комплексный. x2 = a2+j*b2 = filter(decimate(x)). Представляющий комплексную огибающую сигнала с несущей частотой 100 герц и добротностью 150 (полосой 0.66 герца). Теперь осталось разделить полученную комплексную огибающую на запчасти. Посчитав A=abs(x2) вы получите амплитудную составляющую комплексной огибающей (видимо это то, что вам и надо). А посчитав phi=arg(x2) - фазовую (на всяк случай - вдруг пригодится).

Эту схему можно заменить и обычным оцифровыванием сигнала (с частотой дискретизации >200 герц) с дальнейшим его домножением на комплексную 100-герцовую экспоненту exp(j*2*pi*100*t) для сдвига спектра (гетеродинирования). Дальнейшая обработка - прореживание, фильтрация, модуль/аргумент - аналогично предыдущему случаю, но требования к дециматору пожестче будут в части подавления наложений. Но это ресурсов больше отъест, зато меньше аналоговых предфильтров.

• У меня такое чувство, что природу не обмануть+ — AndrI   (12.11.2003 20:10, 1374 байт)
  • Вы противоречите сами себе (+) — SМ   (12.11.2003 21:43, 1561 байт)
    • Простите за дремучесть, но все же...+ — AndrI   (13.11.2003 09:20, 1844 байт)
      • Ответ: (+) — SМ   (13.11.2003 23:05, 2125 байт)
        • Спасибо, я думая толк будет+ — AndrI   (14.11.2003 21:17, 936 байт)
          • Ответ: (+) — SМ   (14.11.2003 23:47, 1460 байт)
            • Ответ+ — AndrI   (15.11.2003 20:25, 1199 байт)
              • Ответ: (+) — SМ   (16.11.2003 11:08, 742 байт)
                • Ответ+ — AndrI   (16.11.2003 20:03, 543 байт)
        • Хотя, чувствуется, советуй, не советуй, толку от этого... Абсолютный ноль. — SМ   (13.11.2003 23:07, пустое)
• И я добавлю :) — SVI   (12.11.2003 10:26, 311 байт)
  • Кстати, по поводу 50 герц - ее лучше вообще задавить отдельно, хотя-бы Т-мостом. Применительно к конкретно этому случаю. — SМ   (12.11.2003 10:44, пустое)
    • По моему, не стоит — SVI   (12.11.2003 11:58, 213 байт)
      • Во первых, вопрос сколько этих 50-ти герц в сигнале было изначально. Во вторых - оно прям на частоту Найквиста попало, что не есть хорошо при квадратурной дискретизации. Однако если вариант с домножением на комплексную exp, то без проблем. — SМ   (12.11.2003 19:34, пустое)
• И, на всякий случай, еще раз повторю (+) — SМ   (12.11.2003 00:05, 508 байт)
  • Кажется дошло - спасибо большой+ — AndrI   (12.11.2003 20:22, 625 байт)