2.4.2. Дискретизация сигнала по времени и явление наложения частот
Импульсная модуляция. Для численного анализа отсчеты большинства непрерывных сигналов
будут производиться через некоторый фиксированный интервал
и полученные таким образом дискретизованные сигналы будут затем использоваться для цифровых вычислений. Дискретизованный сигнал можно рассматривать как результат умножения первоначального непрерывного сигнала на сигнал
состоящий из бесконечного ряда единичных импульсов, или дельта-функций
Это дает дискретизованный, или импульсно-модулированный, сигнал
Следовательно, воспользовавшись теоремой о свертке (П2.1.8), находим
где
является преобразованием Фурье от
Используя для
выражение (2.2.16), преобразуем (2.4.8):
Равенство (2.4.9) показывает, что дискретизованный, или импульсно-модулированный, сигнал
имеет периодическое преобразование Фурье с периодом
и если
обращается в нуль при
то
является просто периодически повторяемой функцией
как показано на рис.
Это означает, что можно восстановить
по
умножив
на
где
Так как умножение в частотной области соответствует свертке во временной области, то отсюда следует, что
Функция
является идеальным фильтром для восстановления непрерывного сигнала
из дискретизованного сигнала
Иначе говоря, функция
является идеальной интерполирующей функцией для равноотстоящих ординат, и формулу (2.4.11) иногда называют интерполяционной формулой Уиттекера.
Наложение частот. Если интервал отсчета таков, что
убывает до нуля, не доходя до
как в случаях
или в на рис. 2.11, то можно восстановить
по
. С другой стороны, если
не равна нулю за частотой
то частотные компоненты от частот выше
присутствуют в
в диапазоне частот —
как, например, в случае
на рис. 2.11. Настота
называется частотой Найквиста и является
наивысшей частотой, которую можно обнаружить на данных, полученных с интервалом отсчета А.
Рис. 2.11. (см. скан) Преобразования Фурье входного сигнала и дискретизованных сигналов для различных интервалов отсчета.
Если, например,
сек, то частота Найквиста равна 5 гц. Преобразование Фурье
дискретизованного сигнала на 4 гц будет состоять из вкладов преобразования
на 4 гц, на
гц, на
гц, на
гц, на
гц и т.
Все эти частоты, кроме первой, называются обычно двойниками (aliases) частоты 4 гц, а их влияние на преобразование
Фурье — явлением наложения частот (aliasing). Следовательно, при дискретизации по времени непрерывных временных рядов нужно надлежащим образом позаботиться о выборе достаточно высокой частоты отсчетов
чтобы избежать искажающего влияния наложения частот на
Явление наложения частот возникает в ряде практических ситуации, например при использовании стробоскопа или в кинофильмах. Так, если в фильме колеса телеги приходят в движение, то впачале видно, что они вращаются в направлении движения, затем при возрастании скорости кажется, что направление вращения меняете”
обратное и скорость колес уменьшается до полной остановки, затем они начинают вращаться с возрастающей скоростью в направлении движения и т. д.
Пример. Чтобы проиллюстрировать обсуждаемые в этом разделе вопросы, предположим, что желательно вычислить длину записи Т и интервал отсчета
необходимые для достижения некоторых целей. Предположим, известно, что изучаемый сигнал содержит две синусоидальные компоненты на частотах 100 и 99 гц. Тогда, если мы хотим различить эти пики в преобразовании Фурье, взятом от конечной записи, нам нужно, как показывает (2.4.5), взять
порядка
гц, т. е. Т должно быть порядка 1 сек. Чтобы оценивать частоты порядка 100 гц, величина
должна быть по меньшей мере 100 гц и, следовательно,
мсек. Таким образом, нужно взять по крайней мере 200 точек.
Если бы нам захотелось различить две частотные компоненты на 999 и 1000 гц, необходимая длительность записи была бы все еще 1 сек, однако интервал отсчета в этом случае нужно было бы взять 0,5 мсек, так что потребовалось бы 2000 точек.
Следовательно, длина записи Т определяет степень различимости пиков в преобразовании Фурье, а интервал отсчета
определяет максимальную частоту, которую можно различать.
ЛИТЕРАТУРА
(см. скан)