ЦИФРОВАЯ ОБРАБОТКА СИГНАЛОВ. Аппроксимация сигналов. Цифровая обработка сигналов
 Скачать 360.5 Kb.
|
|
14.5. ДЕЦИМАЦИЯ И ИНТЕРПОЛЯЦИЯ ЦИФРОВЫХ СИГНАЛОВ [43] Применительно к цифровым сигналам децимация – уменьшение частоты дискретизации данных с сохранением в новом сигнале всей полезной информации. Интерполяция обратна децимации – увеличение частоты дискретизации без изменения информации. Цифровая децимация и интерполяция широко используется в современных системах обработки данных для сжатия и восстановления данных, для уменьшения объемов памяти хранения данных, для увеличения скорости передачи данных, и т.п. Простой, но мало производительный подход – восстановить сигнал в аналоговой форме (ЦАП) и заново оцифровать его (АЦП) с новой частотой дискретизации. Цифровые методы позволяют выполнить эту операцию в более эффективной форме. Децимация с целым шагом. Кратная компрессия частоты дискретизации снижает частоту дискретизации входного сигнала x(k) с fd до fd/M путем отбрасывания М-1 отсчетов в каждой последовательной серии из М-отсчетов, т.е. из М-отсчетов оставляется только 1: y(m) =  x(k) (k-mM), m = 0, 1, 2, … , K/M. (14.5.1)Естественно, что частота Найквиста fN входного сигнала x(k) компрессора для выходного сигнала y(k) также уменьшается в М раз и становится равной fN' = fN/M для выходного сигнала. Для полного сохранения после компрессии полезной информации, содержащейся в сигнале x(k), максимальная частота полезной информации во входном сигнале не должна превышать значения fmax fN/2M. В противном случае децимация будет некорректной и в новом главном частотном диапазоне выходного сигнала произойдет искажение спектра полезной информации за счет сложения со спектрами боковых диапазонов. Пример корректной децимации сигнала с М=2 и спектры входного и децимированного сигнала приведены на рис. 14.5.1.  Рис. 14.5.1. Входные сигналы кроме полезной информации могут содержать статистические шумы и помехи, распределенные по всему частотному диапазону. При децимации шумы и помехи в частотном диапазоне от fN/M до fN входного сигнала зеркально отражаются от fN' нового главного частотного диапазона и их суммирование со спектром нового главного диапазона и полезного сигнала может приводить к увеличению уровня шумов и искажению информации. Для исключения этого эффекта перед конверсией сигнала необходимо выполнять его низкочастотную фильтрацию со срезом на частоте fN/2M.  Рис. 14.5.2. На рис. 14.5.2 приведен спектр Z2n сигнала zk с рисунка 14.5.1 (спектр Zn), в который для наглядности эффекта зеркального отражения условно введен только высокочастотный шум на интервале fs - fN, и соответствующий данному спектру сигнал z2k. При децимации сигнала z2 с М=2 сначала была выполнена его низкочастотная фильтрация с частотой среза fs, что полностью сняло модельный шум, и при восстановлении сигнала (интерполяции) получен сигнал y1k, полностью соответствующий сигналу zk. При выполнении децимации без предварительной фильтрации восстанавливается сигнал y2k, который по своей форме отличается как от сигнала y1k = zk, так и от входного сигнала z2k. Интерполяция с целым шагом. Экспандер частоты дискретизации увеличивает частоту дискретизации входного сигнала xk с fd до Lfd путем введения (L-1) нулевых отсчетов после каждого отсчета входного сигнала. При этом форма спектра выходного сигнала yk остается без изменения, но частотная шкала спектра сжимается в L раз и в границы главного диапазона спектра входного сигнала ±fN заходят боковые диапазоны спектра выходного сигнала (зеркальные частоты). Это наглядно можно видеть на рис. 14.5.3 сравнением спектров Xn для входного сигнала xk, и Yn для экспандированного сигнала xk с L=2. Следовательно, фактическая частота Найквиста выходного сигнала становится равной fN/L. Для исключения зеркальных частот и распределения энергии отсчетов xk по L выходным интервалам экспандированный сигнал пропускается через фильтр нижних частот со срезом на частоте fN/L и с коэффициентом L для компенсации распределения энергии отсчетов по интервалам L. Результат операции можно видеть на сигнале yk по сравнению с исходным сигналом zk (рис. 14.5.1), децимацией которого с М=2 был получен сигнал xk.  Рис. 14.5.3. |