Децимация (обработка сигналов)
Децима́ция (от лат. decimatio, от decem — «десять») — уменьшение частоты дискретизации дискретного во времени сигнала путём прореживания его отсчётов.
Отсчёт — численное значение напряжения сигнала в определённый момент времени.
Термин децимация в обработке сигналов происходит от изначального значения этого слова. Однако имеется существенное различие: если в Древнем Риме при децимации каждого десятого казнили, то при децимации сигналов каждый N-й отсчёт, напротив, остаётся.
Прореживание отсчётов может осуществляться путём их выборочного удаления из исходного массива или на основе частичного суммирования[1] (накопления со сбросом, фильтрации) в фиксированных временных интервалах — стробах (отсюда происходит другое название этой операции — дополнительное стробирование отсчётов АЦП[2]).
Децимация удалением отсчётов сигналов
При таком виде децимации из исходной последовательности отсчётов
- a0, a1, a2, …
берётся каждый N-й отсчет (N — целое число):
- a0, aN, a2N, … ; N > 1
остальные отсчёты отбрасываются. Преобразование спектра при децимации существенно зависит от спектра исходного сигнала:
- Если исходный сигнал не содержит частот, превышающих частоту Найквиста децимированного сигнала, то форма спектра полученного (децимированного) сигнала совпадает с низкочастотной частью спектра исходного сигнала. Частота дискретизации, соответствующая новой последовательности отсчётов, в N раз ниже, чем частота дискретизации исходного сигнала, и спектр полученного сигнала масштабирован по оси абсцисс относительно спектра исходного сигнала.
- Если исходный сигнал содержит частоты, превышающие частоту Найквиста децимированного сигнала, то при децимации будет иметь место алиасинг (наложение спектров).
Таким образом, для сохранения спектра необходимо до децимации удалить из исходного сигнала частоты, превышающие частоту Найквиста децимированного сигнала. Эта операция производится цифровыми фильтрами.
Децимация на основе дополнительного стробирования отсчетов АЦП
Суть этого метода децимации сводится к тому, что по серии отсчётов АЦП формируется один суммарный отсчёт, жёстко привязанный к сетке импульсов такта АЦП[2].
В случае децимации отсчётов напряжений видеосигнала при этом результат прореживания описывается выражением:
- [math]\displaystyle{ y[n] = \sum_{k=0}^{M-1} x[nM+k]\, }[/math],
где x[•] — входные отсчеты напряжений сигнала до децимации, M — длительность строба.
Для гармонических сигналов[2]
- [math]\displaystyle{ y[n] = \sum_{k=0}^{M-1} x[nM+k]\ e^{-i 2\pi f kT}, n=0,1,.., N }[/math],
где T — период дискретизации АЦП (интервал между отсчётами).
Если [math]\displaystyle{ 2\pi f kT = k\pi/2 }[/math], то имеет место [math]\displaystyle{ e^{-i 2\pi f kT} = 1, e^{-i \pi/2}, e^{-i \pi},... }[/math] и следовательно[2]
- [math]\displaystyle{ Re (y[n]) = \sum_{k=0}^{M-1}( Re (x[nM+k]) cos { k \pi/2 } + Im (x[nM+k]) sin { k \pi/2 })\, }[/math],
- [math]\displaystyle{ Im (y[n]) = \sum_{k=0}^{M-1}( Im (x[nM+k]) cos { k \pi/2 } - Re (x[nM+k]) sin { k \pi/2 })\, }[/math].
При [math]\displaystyle{ Im (x[nM+k])=0 }[/math] получим
- [math]\displaystyle{ Re (y[n]) = x[nM]-x[nM+2]+x[nM+4]-..., }[/math],
- [math]\displaystyle{ Im (y[n]) = -x[nM+1]+x[nM+3]-x[nM+5]-... }[/math].
Такая обработка позволяет осуществить прореживание отсчётов сигнала без энергетических потерь, декоррелировать сигналы за счёт перехода к их укрупнённому представлению[3], осуществить цифровую фильтрацию и формирование квадратурных составляющих напряжений сигналов (I/Q-демодуляция), реализовать сверхрэлеевское разрешение импульсных сигналов по времени прихода[2].
В случае, если аналоговый сегмент не позволяет достаточно эффективно обеспечить предварительную (anti-aliasing) фильтрацию, указанный метод децимации может быть модифицирован в виде:
- [math]\displaystyle{ y[n] = \sum_{k=0}^{M-1} x[nM+k]h[k]\ e^{-i 2\pi f kT}, n=0,1,.., N }[/math],
где [math]\displaystyle{ h[k] }[/math] — вектор весовых множителей.[4]
В качестве примера при [math]\displaystyle{ 2\pi f kT = k\pi/2 }[/math] следует указать процедуру децимации с нечётной длительностью строба:[5]
- [math]\displaystyle{ Re (y[n]) = -6[nM+1]+32x[nM+3]-52x[nM+5]+32x[nM+7]-6x[nM+9], }[/math]
- [math]\displaystyle{ Im (y[n]) = x[nM]-17x[nM+2]+46x[nM+4]-46x[nM+6]+17x[nM+8]-x[nM+10]. }[/math]
Децимация с использованием КИХ фильтров
Альтернативным вариантом дополнительного строобирования отсчётов АЦП является их низкочастотная фильтрация с использованием фильтров с конечной импульсной характеристикой (КИХ или FIR). При этом по массиву входных отсчётов также формируется только каждый M-й выходной отсчёт как взвешенная сумма напряжений исходных отсчётов с весами в виде дискретной импульсной характеристики КИХ -фильтра :
- [math]\displaystyle{ y[n] = \sum_{k=0}^{K-1} x[nM-k]\cdot h[k], }[/math]
где h[•] — импульсная характеристика, K — её длительность; x[•] — входные отсчеты напряжений сигнала до децимации.
Децимация с дробным коэффициентом прореживания отсчётов
Данный вид децимации необходим, например, в тех случаях, когда частота дискретизации сигналов некогерентна с частотой обрабатываемого радиосигнала.
В этом случае для децимации с коэффициентом M/L, где M, L ∈ ℤ; M > L, сначала необходимо провести интерполяцию отсчётов c помощью интерполирующего фильтра порядка L, а затем выполнить их прореживание с коэффициентом M, например, с помощью описанной процедуры дополнительного стробирования отсчетов АЦП[2]. Как правило, обе операции совмещают в одном фильтре.
Возможна также децимация с иррациональными коэффициентами пересчёта.[6]
См. также
Примечания
- ↑ Антипов В.Н., Горяинов В.Т., Кулин А.Н. и др. Радиолокационные станции с цифровым синтезированием апертуры антенны. — М.:: Радио и связь, 1988. — С. 42 — 43. — 304 с.
- ↑ 2,0 2,1 2,2 2,3 2,4 2,5 Слюсар В.И. Синтез алгоритмов измерения дальности М источников при дополнительном стробировании отсчетов АЦП.// Известия вузов. Сер. Радиоэлектроника.- Том 39, № 5. — 1996. — С. 55 — 62.
- ↑ Харкевич А.А. Теория информации. Опознание образов. Избранные труды в трёх томах. Т. 3. -. — М.:: Наука, 1973. — С. 85 — 89. — 524 с.
- ↑ Слюсар В. И. Развитие схемотехники ЦАР: некоторые итоги. Часть 2.// Первая миля. Last mile (Приложение к журналу "Электроника: наука, технология, бизнес"). – N2. - 2018. - C. 76 - 80.[1] Архивная копия от 20 июня 2018 на Wayback Machine
- ↑ Слюсар В. И., Живило Е. А. Цифровая фильтрация, эквивалентная тандемному квадратурному дециматору. //VI Международный научно-технический симпозиум «Новые технологии в телекоммуникациях» (ГУИКТ-Карпаты '2013), 21 — 25 января 2013. — Карпаты, Вышков.- C. 41 — 43. [http://web.archive.org/web/20160406103605/http://slyusar.kiev.ua/VYSHKIV_2013_2.pdf Архивная копия от 6 апреля 2016 на Wayback Machine]
- ↑ Milić, Ljiljana. Multirate Filtering for Digital Signal Processing (англ.). — New York: Hershey, 2009. — P. 192. — ISBN 978-1-60566-178-0.. — «Generally, this approach is applicable when the ratio Fy/Fx is a rational, or an irrational number, and is suitable for the sampling rate increase and for the sampling rate decrease.».
Литература
- Oppenheim, Alan V.[англ.]; Schafer, Ronald W.; Buck, John R. Discrete-Time Signal Processing. — 2nd. — Prentice Hall, 1999. — ISBN 0-13-754920-2.
- Proakis, John G. Digital Signal Processing: Principles, Algorithms and Applications (англ.). — 3rd. — India: Prentice-Hall, 2000. — ISBN 8120311299.