Какова роль удержания нулевого порядка в гибридной аналоговой / цифровой системе выборочных данных?


14

Я признаю, я задаю этот вопрос риторически. Мне любопытно, какие ответы вернутся из этого.

Если вы решите ответить на этот вопрос, убедитесь, что вы хорошо понимаете теорему Шеннона-Найквиста. Особенно реконструкция. Также будьте осторожны с "гочами" в учебниках. Инженерное представление о функции импульса Дирака является достаточным. Вам не нужно беспокоиться обо всем «распределении», импульс Дирака как зарождающаяся дельта-функция достаточно хорош:

δ(t)=limτ01τrect(tτ)

где

rect(t){0if |t|>121if |t|<12

Проблемы, связанные с точностью, шириной слова в образцах и квантованием при преобразовании, не имеют отношения к этому вопросу. Но масштабирование от входа к выходу имеет значение.

Я напишу свой собственный ответ на него в конце концов, если кто-то еще не представит точный и педагогически полезный ответ. Я мог бы даже вознаградить за это (может с тем же успехом потратить то маленькое количество повторений, которое у меня есть).

Имейте в этом.


Вы заинтересованы услышать о псевдонимах в первую очередь?
Deadude

нет. Я предполагаю, что все правила теоремы выборки соблюдаются. то есть, ни содержимое, ни энергия во входных данных непрерывного времени не отбираются, которые находятся на уровне или выше . теперь, помните, есть разница между «псевдонимами» и «изображениями». fs2
Роберт Бристоу-Джонсон

Насколько я помню, удержание нулевого порядка - это просто задержка между выборками в цифровой системе, и, очевидно, может повлиять на аналоговую сторону вещей между одной выборкой и следующей
KyranF

@KyranF, это немного больше, чем это.
Роберт Бристоу-Джонсон

@ robertbristow-johnson из ответов, данных Тимо, действительно выглядит более запутанным, чем я думал. Удачи с этим!
KyranF

Ответы:


6

Настроить

Мы рассматриваем систему с входным сигналом , и для ясности мы называем значения как напряжения, где это необходимо. Наш образец период , а соответствующий образец ставка .x ( t ) T f s1 / Tx(t)x(t)Tfs1/T

Для преобразования Фурье мы выбираем соглашения давая обратное преобразование Фурье Обратите внимание, что с этими соглашениями является функцией переменной Лапласа .

X(i2πf)=F(x(t))x(t)ei2πftdt,
x(t)=F1(X(i2πf))X(i2πf)ei2πftdf.
Xs=iω=i2πf

Идеальный отбор проб и реконструкция

Начнем с идеальной выборки: в соответствии с теоремой выборки Найквиста-Шеннона для заданного сигнала который ограничен полосой частот , то есть то исходный сигнал может быть полностью восстановлен по выборкам , где . Другими словами, учитывая условие ширины полосы сигнала (называемое критерием Найквиста ), достаточно знать его мгновенные значения в эквидистантных дискретных моментах времени.е < 1x(t)X(i2πf)=0,f<12fs x[n]x(nT)nZ

X(i2πf)=0,when|f|12fs,
x[n]x(nT)nZ

Теорема выборки также дает явный метод для проведения реконструкции. Остановимся на этом так, как это будет полезно в дальнейшем: оценим преобразование Фурье сигнала по его сумме Римана с шагом : где . Давайте перепишем это как интеграл, чтобы определить ошибку, которую мы делаем: х ( т ) Т Х ( я 2 π е ) ~ Е п = - х ( п Δ т ) е - я 2 π е п Δ т Δ т , Δ т = Т Е n = - x ( n T ) e -X(i2πf)x(t)T

X(i2πf)n=x(nΔt)ei2πfnΔtΔt,
Δt=T x(t)n = - Tδ(t-nT)n = - δ(f-k/T)
n=x(nT)ei2πfnTT=n=x(t)ei2πftTδ(tnT)dt=X(i2πf)F(Tn=δ(tnT))(1)=k=X(fk/T),
где мы использовали теорему свертки о произведении и функцию выборки , тот факт, что преобразование Фурье функции выборки является , и выполнил интеграл по дельта-функциям.x(t) n=Tδ(tnT)n=δ(fk/T)

Обратите внимание, что левая часть в точности равна , где - преобразование Фурье с дискретным временем соответствующего дискретизированного сигнала , где - безразмерная частота дискретного времени.X 1 / T ( i 2 π f T ) x [ n ] x ( n T ) f TTИкс1/T(я2πеT)Икс1/T(я2πеT)Икс[N]Икс(NT)еT

Здесь мы видим существенную причину критерия Найквиста: это именно то, что требуется, чтобы гарантировать, что условия суммы не перекрываются. С помощью критерия Найквиста сумма сводится к периодическому расширению спектра от интервала до всей действительной линии.[-еs/2,еs/2]

Поскольку DTFT в имеет такое же преобразование Фурье в интервале что и наш исходный сигнал, мы можем просто умножить его на прямоугольную функцию и верните исходный сигнал. С помощью теоремы свертки это сводится к свертыванию гребня Дирака с преобразованием Фурье прямоугольной функции, которая в наших соглашениях имеет вид где нормализованная функция sinc имеет вид [ - е с / 2 , е с / 2 ] т е с т ( е / е ы ) Р ( г е с т ( е / е ы ) ) = 1 / T сек я н с ( т / Т ) , s i n c ( x ) sin ( π(1)[-еs/2,еs/2]ресT(е/еs)

F(rect(f/fs))=1/Tsinc(t/T),
x ( t ) = n = - x [ n ] s i n c ( t / T - n ) .
sinc(x)sin(πx)πx.
Затем свертка просто заменяет каждую дельту Дирака в гребне Дирака sinc-функцией, сдвинутой в положение дельты, давая Это интерполяционная формула Уиттекера-Шеннона .
(2)x(t)=n=x[n]sinc(t/Tn).

Неидеальная выборка

Для перевода вышеприведенной теории в реальный мир, наиболее трудной частью является обеспечение ограничения полосы пропускания, что необходимо сделать перед выборкой. Для целей этого ответа мы предполагаем, что это было сделано. Оставшаяся задача состоит в том, чтобы затем взять образцы мгновенных значений сигнала. Поскольку для реального АЦП потребуется конечное количество времени для формирования аппроксимации к выборке, обычная реализация будет сохранять значение сигнала в схеме выборки и удержания, из которой формируется цифровая аппроксимация.

Несмотря на то, что это очень похоже на удержание нулевого порядка, это особый процесс: значение, полученное из выборки и удержания, действительно является точно мгновенным значением сигнала, вплоть до приближения, что сигнал остается постоянным для продолжительность зарядки конденсатора, в котором хранится значение выборки. Это обычно хорошо достигается системами реального мира.

Следовательно, мы можем сказать, что реальный АЦП, игнорирующий проблему ограничения полосы пропускания, является очень хорошим приближением к случаю идеальной выборки, и, в частности, «лестница», возникающая из выборки и удержания, не вызывает никакой ошибки в Выборка сама по себе.

Неидеальная реконструкция

Для реконструкции цель состоит в том, чтобы найти электронную схему, которая выполняет сумму синусов, появляющуюся в . Поскольку sinc имеет бесконечную протяженность во времени, совершенно очевидно, что это не может быть точно реализовано. Кроме того, формирование такой суммы сигналов даже в разумном приближении потребует нескольких подсхем и быстро станет очень сложным. Поэтому обычно используется гораздо более простое приближение: в каждый момент дискретизации выводится напряжение, соответствующее значению выборки, и поддерживается постоянным до следующего момента дискретизации (хотя см. Дельта-сигма-модуляция для примера альтернативного метода). Это удержание нулевого порядка , и оно соответствует замене sinc, которую мы использовали выше, функцией прямоугольника. 1 / Т т е с т ( т / Т - 1 / 2 ) ( 1 / Т т е с т ( т / Т - 1 / 2 ) ) * ( Е п = - Т х [ п ] δ ( t - n T ) ) , 1 / T T(2)1/Trect(t/T1/2) . Оценка свертки используя свойство определения дельта-функции, мы видим, что это действительно приводит к классической лестничной форме с непрерывным временем. Множитель вводит для отмены введенного в . Необходимость такого фактора также очевидна из того факта, что единицы импульсного отклика равны 1 / время.

(1/Trect(t/T1/2))(n=Tx[n]δ(tnT)),
1/TT(1)

Сдвиг на просто гарантирует причинность . Это сводится только к сдвигу вывода на 1/2 выборки относительно использования (что может иметь последствия в системах реального времени или когда требуется очень точная синхронизация с внешними событиями ), который мы будем игнорировать в дальнейшем.1 / Т т е с т ( 1 / Т )1/2T1/Trect(1/T)

Сравнивая с , мы заменили прямоугольную функцию в частотной области, которая оставила базовую полосу полностью нетронутой, и удалили все высокочастотные копии спектра, называемые изображениями , с помощью преобразования Фурье функции . Это, конечно, 1 / T r e c t ( t / T ) s i n c ( f / f s ) .(1)1/Trect(t/T)

sinc(f/fs).

Обратите внимание, что логика несколько отличается от идеального случая: там мы определили нашу цель, которая заключалась в удалении изображений в частотной области, и вывели последствия во временной области. Здесь мы определили, как восстанавливать во временной области (поскольку это то, что мы знаем, как делать), и вывели последствия в частотной области.

Таким образом, результатом удержания нулевого порядка является то, что вместо прямоугольного оконного управления в частотной области мы в итоге получаем sinc как оконную функцию. Следовательно:

  • Частотная характеристика больше не ограничена полосой пропускания. Скорее он затухает как , причем верхние частоты являются изображениями исходного сигнала1/f
  • в основной полосе отклик уже значительно затухает, достигая около -4 дБ при1/2fs

В целом, удержание нулевого порядка используется для аппроксимации функции sinc во временной области, появляющейся в формуле интерполяции Уиттекера-Шеннона . При выборке аналогично выглядящий выборка и удержание является техническим решением проблемы оценки мгновенного значения сигнала и не дает никаких ошибок сам по себе.

Обратите внимание, что при реконструкции также не теряется информация, поскольку мы всегда можем отфильтровать высокочастотные изображения после первоначального удержания нулевого порядка. Потери усиления также могут быть скомпенсированы обратным sinc-фильтром до или после ЦАП. Таким образом, с более практической точки зрения удержание нулевого порядка используется для построения начального реализуемого приближения к идеальной реконструкции, которое затем может быть дополнительно улучшено, если это необходимо.


Тимо это интересно. Вы столкнулись с последствиями политики Википедии. Проверьте эту старую версию статьи Википедии о теореме выборки . вместо того, чтобы прятаться за формулой суммирования Пуассона, он просто показывает, как выборка генерирует изображения, и явно то, что требуется для восстановления исходного непрерывного сигнала. и вы можете понять, почему в функции выборки присутствует этот фактор. T
Роберт Бристоу-Джонсон

Интересно, что старая версия статьи в Википедии, на мой взгляд, более понятна. Вычисления почти такие же, как я пишу выше, за исключением того, что они дают немного больше деталей.
Тимо

Во всяком случае, я не совсем уверен, зачем это нужно, чтобы понять, зачем нужен фактор : я думаю, что то, что я пишу в ответе, является достаточным условием для необходимости (технически, это условие согласованности, но мы уже предположим, что реконструкция возможна). Теперь, конечно, понимание всегда субъективно. Например, здесь можно считать более глубокой причиной появления фактора то, что существу становится мерой интегрирования когда принимается предел . Т Т Т д т Т 0TTTTdtT0
Тимо

Я полагаю, что вы называете «почему» появление 1 / T в представлении гребня Дирака в виде суммы комплексных экспонент, в en.wikipedia.org/w/… ? Это, конечно, один из способов выразить это, и оно напрямую связано с ролью как меры. T
Тимо

1
Я не могу помочь, но думаю, что вы должны просто добавить ответ, который вы ищете. Комментарии не для расширенного обсуждения.
Дэвид

4

Удержание нулевого порядка играет роль аппроксимации дельта- и -функций, появляющихся в теореме выборки, в зависимости от того, что подходит.sinc

Для ясности рассмотрим систему АЦП / ЦАП с сигналом напряжения. Все следующее относится, однако, к любой системе отбора проб с соответствующей сменой единиц измерения. Я также предполагаю, что входной сигнал уже магически ограничен по полосе, чтобы соответствовать критерию Найквиста.

Начните с выборки: в идеале можно выбрать значение входного сигнала в один момент. Поскольку реальным АЦП требуется конечное количество времени для формирования их приближения, мгновенное напряжение аппроксимируется выборкой и удержанием (мгновенное значение аппроксимируется временем переключения, используемым для зарядки конденсатора). По сути, удержание преобразует проблему применения дельта-функционала к сигналу в задачу измерения постоянного напряжения.

Отметим здесь, что разница между умноженным входным сигналом на последовательность импульсов или удержанием нулевого порядка, применяемым в те же моменты времени, является просто вопросом интерпретации, поскольку АЦП, тем не менее, будет хранить только удерживаемые мгновенные напряжения. Одно можно реконструировать из другого. Для целей этого ответа я приму толкование того, что дискретизированный сигнал является непрерывным сигналом в виде где является опорное напряжение АЦП / ЦАП, есть число битов, являются образцы представлены обычным образом как целые числа, аVrefnxkΔtxk

x(t)=ΔtVref2nkxkδ(tkΔt),
VrefnxkΔtпериод выборки. Эта несколько нетрадиционная интерпретация имеет то преимущество, что я всегда рассматриваю непрерывный сигнал, а выборка здесь просто означает представление его в терминах чисел , которые действительно являются выборками в обычном смысле.xk

В этой интерпретации спектр сигнала в основной полосе точно такой же, как и у исходного сигнала, а эффективная свертка с помощью последовательности импульсов дает эффект репликации этого сигнала, что делает спектр периодическим. Реплики называются изображениями спектра. То, что нормировочный коэффициент необходим, можно увидеть, например, с учетом смещения постоянного тока импульса длительностью 1 В длительностью : его DC-смещение, определенное как -компонент преобразования Фурье, равно Чтобы получить тот же результат из нашей выборочной версии, мы действительно должны включить коэффициент . *Д т е = 0 х ( 0 ) = Д т 0 1 В д т = 1 В Д т . Δ тΔtΔtf=0

x^(0)=0Δt1Vdt=1VΔt.
Δt

Тогда идеальная реконструкция означает построение электрического сигнала, который имеет тот же спектр основной полосы частот, что и этот сигнал, и не содержит компонентов на частотах вне этого диапазона. Это то же самое, что свертывание импульсного поезда с соответствующей . Это довольно сложно сделать с помощью электроники, поэтому часто аппроксимируется прямоугольной функцией AKA, удерживающей нулевой порядок. По сути, для каждой дельта-функции значение выборки сохраняется в течение периода выборки.с я н сsincsinc

Чтобы увидеть, какие последствия это имеет для восстановленного сигнала, я заметил, что удержание в точности эквивалентно свертыванию последовательности импульсов с прямоугольной функцией Нормализация этой прямоугольной функции определяется требованием правильного воспроизведения постоянного напряжения, или, другими словами, если при выборке измерялось напряжение , то же напряжение выводится при реконструкции.V1

rectΔt(t)=1Δtrect(tΔt).
V1

В частотной области это равносильно умножению частотной характеристики на преобразование Фурье прямоугольной функции, которое равно Обратите внимание, что усиление на постоянном токе равно . На высоких частотах распадается как , и, следовательно, ослабляет изображения спектра.1sinc1/f

rect^Δt(f)=sinc(πΔtf).
1sinc1/f

В конце концов, возникающая в результате удержания нулевого порядка, ведет себя как фильтр нижних частот сигнала. Обратите внимание, что никакая информация не теряется в фазе дискретизации (при условии, что критерий Найквиста), и в принципе, ничего не теряется при реконструкции: фильтрация в основной полосе с помощью может быть компенсирована обратным фильтром (и это действительно иногда делается, см., например, https://www.maximintegrated.com/en/app-notes/index.mvp/id/3853 ). Скромное затухание обычно требует некоторой формы фильтрации для дальнейшего ослабления изображений.s i n c - 6 d B / o c t a v e s i n csincsinc6dB/octavesinc

Отметим также, что мнимый генератор импульсов, который может физически воспроизводить последовательность импульсов, использованную в анализе, будет выводить бесконечное количество энергии при восстановлении изображений. Это также может вызвать некоторые волосатые эффекты, такие как то, что АЦП при повторной выборке на выходе ничего не увидит, если он не будет полностью синхронизирован с исходной системой (он будет в основном производить выборку между импульсами). Это ясно показывает, что даже если мы не можем точно ограничить полосу выходного сигнала, некоторое приблизительное ограничение полосы частот всегда необходимо для регуляции полной энергии сигнала, прежде чем он может быть преобразован в физическое представление.

Подвести итоги:

  • в обоих направлениях удержание нулевого порядка действует как приближение к дельта-функции, или ее ограниченная полосой форма, функции .sinc
  • с точки зрения частотной области, это приближение к кирпичному фильтру, который удаляет изображения и, следовательно, регулирует бесконечное количество энергии, присутствующей в идеализированной последовательности импульсов.

* Это также ясно из размерного анализа: единицами преобразования Фурье сигнала напряжения являются тогда как Дельта-функция имеет единицы , которые отменяют единицу времени, полученную из интеграла в преобразовании.1/сVs=VHz,1/s


когда таймер позволит мне, я вознагражу за это, Тимо. Есть некоторые вещи, которые мне нравятся: например, иметь усиление постоянного тока = 1, что соответствует уравнению. 1 на вашей сентенции цитате, но путь слишком много учебников облажаться с усилением , что они не знают , что делать с. и кажется, что вы понимаете, что ZOH не имеет никакого отношения к любому возможному S / H на входе АЦП. это хорошо. я все еще буду ждать немного более точного ответа. и не беспокойтесь о . Я предполагаю, что это то же самое для АЦП и ЦАП. V refTVref
Роберт Бристоу-Джонсон

@ Робертбристоу-Джонсон: спасибо за добрые слова! Можете ли вы указать немного, в каком направлении вы ищете более строгий? Больше деталей, больше математического ответа или что-то совсем другое?
Тимо

я предполагаю математическую обработку с чистыми и последовательными математическими обозначениями. я бы предложил быть согласным с Оппенгеймом и Вильским или что-то в этом роде. возможно, так что преобразования Лапласа и Фурье имеют согласованную и совместимую запись . обсудите, что говорит теорема о сэмплировании и как она отличается в реальности и где ZOH вступает в это.
T1fs
x[n]x(nT)
F{x(t)}=X(j2πf)+x(t)ej2πft dt
Роберт Бристоу-Джонсон

Хорошо, позвольте мне на самом деле попробовать написать другой ответ, так как редактирование этого, чтобы изменить нотацию на то, что вы предпочитаете и т. Д., Вероятно, оставит немного беспорядка. Сначала я исправлю небольшую ошибку, потому что это беспокоит меня ...
Тимо

Я был немного смущен и медлителен при розыгрыше и не нажал значок щедрости, чтобы наградить вашу награду. согласно правилам: если вы не начисляете награду в течение 7 дней (плюс льготный период), ответ с наибольшим количеством голосов, полученный после того, как награда началась с минимальной оценкой 2, будет присуждена в размере половины суммы вознаграждения. Если два или более отвечающих критериям ответа имеют одинаковый балл (т.е. их баллы связаны), самый старый ответ присуждается за вознаграждение. Если нет ответа, отвечающего этим критериям, награда не присуждается никому. - согласно этим правилам вы должны получить его в течение недели.
Роберт Бристоу-Джонсон

3

Преобразование Фурье :

X(j2πf)=F{x(t)}+x(t) ej2πft dt

Обратное преобразование Фурье:

x(t)=F1{X(j2πf)}=+X(j2πf) ej2πft df

Функция прямоугольного импульса :

rect(u){0if |u|>121if |u|<12

Функция "Sinc" ("sinus cardinalis") :

sinc(v){1if v=0sin(πv)πvif v0

Определить частоту дискретизации , , как величина , обратная периоду дискретизации .fs1TT

Обратите внимание:

F{rect(tT)}=T sinc(fT)=1fs sinc(ffs)

Гребень Дирака (он же «функция выборки», он же «функция Ша») :

IIIT(t)n=+δ(tnT)

Гребень Дирака является периодическим с периодом . Ряд Фурье :T

IIIT(t)=k=+1Tej2πkfst

Дискретный непрерывный сигнал :

идеально подобранный сигнал с помощью гребня Дирака

xs(t)=x(t)(TIIIT(t))=x(t)(Tn=+δ(tnT))=T n=+x(t) δ(tnT)=T n=+x(nT) δ(tnT)=T n=+x[n] δ(tnT)

где .x[n]x(nT)

Это означает, что определяется исключительно выборками и периодом выборки и полностью теряет любую информацию о значениях за промежутки времени между экземплярами выборки. представляет собой дискретную последовательность чисел и является сокращенной записью для DSP сорта для . Хотя верно, что для , значение для любого не являющегося целым числом, не определено.xs(t)x[n]Tx(t)x[n]ИксNИксs(T)знак равно0NT<T<(N+1)TИкс[N]N

NB: дискретный сигнал и все операции с дискретным временем над ним, такие как -преобразование , дискретное преобразование Фурье (DTFT) , дискретное преобразование Фурье (DFT) , являются «независимыми» относительно частоты дискретизации или периода дискретизации . После того, как вы в дискретном времени домен, вы не знаете (или уход) о . Это только с Найквиста-Шеннона сэмплирования и реконструкции теоремы , что и вместе взятые.Икс[N]ZTИкс[N]TИкс[N]T

Преобразование Фурье являетсяИксs(T)

Иксs(J2πе)F{Иксs(T)}знак равноF{Икс(T)(TIIIT(T))}знак равноF{Икс(T)(TΣКзнак равно-+1TеJ2πКеsT)}знак равноF{ΣКзнак равно-+Икс(T) еJ2πКеsT}знак равноΣКзнак равно-+F{Икс(T) еJ2πКеsT}знак равноΣКзнак равно-+Икс(J2π(е-Кеs))

Важное примечание о масштабировании: функция выборки и дискретизированный сигнал имеют коэффициент который вы не увидите почти во всех учебниках. Это педагогическая ошибка авторов этих учебников по нескольким (связанным) причинам: TIIIT(T)Иксs(T)T

  1. Во-первых, пропуская изменяет размер выборочного сигнала по сравнению с выборки сигнала .TИксs(T)Икс(T)
  2. Этот фактор понадобится где-то в цепочке сигналов. Эти учебники, которые исключают его из функции выборки, в конечном итоге помещают его в часть восстановления теоремы выборки, обычно в качестве усиления полосы пропускания фильтра восстановления. Это сбивает с толку. Кто-то может разумно спросить: «Как мне спроектировать LPF с кирпичной стеной с усилением полосы пропускания ?»TT
  3. Как будет показано ниже, оставление здесь приводит к аналогичной ошибке масштабирования для чистой передаточной функции и чистой частотной характеристики удержания нулевого порядка (ZOH). Все учебники по цифровым (и гибридным) системам управления, которые я видел, допускают эту ошибку, и это серьезная педагогическая ошибка.T

Обратите внимание, что DTFT для и преобразование Фурье для дискретизированного сигнала при правильном масштабировании практически идентичны:Икс[N]Иксs(T)

DTFT:

ИксDTFT(ω)Z{Икс[N]}|Zзнак равноеJωзнак равноИксZ(еJω)знак равноΣNзнак равно-+Икс[N] е-JωN

Можно показать, что

ИксDTFT(ω)знак равноИксZ(еJω)знак равно1TИксs(J2πе)|езнак равноω2πT


Приведенная выше математика верна независимо от того, правильно ли выбрана или нет. "правильно выбран", если может быть полностью восстановлен из выборок и знания частоты дискретизации или периода выборки. Sampling теорема говорит нам , что необходимо восстановить или реконструировать от и .Икс(T)Икс(T)Икс(T)Икс[N]Икс(T)Икс[N]T

Если является ограниченной полосой до некоторой bandlimit , это означает ,Икс(T)В

Икс(J2πе)знак равно0для всех|е|>В

полосовой спектр

Рассмотрим спектр дискретизированного сигнала, составленного из смещенных изображений оригинала:

Иксs(J2πе)знак равноΣКзнак равно-+Икс(J2π(е-Кеs))

Исходный спектр может быть восстановлен из выборочного спектра если ни одно из сдвинутых изображений , перекрывают соседних соседей. Это означает, что правый край изображения (то есть ) должен находиться полностью слева от левого края ( ) -ое изображение (то есть ). Пересчитано математически,Икс(J2πе)Иксs(J2πе)Икс(J2π(е-Кеs))КИкс(J2π(е-Кеs))К+1Икс(J2π(е-(К+1)еs))

Кеs+В<(К+1)еs-В

что эквивалентно

еs>2В

Если мы производим выборку с частотой дискретизации, которая в два раза превышает пропускную способность, ни одно из изображений не перекрывается, исходный спектр , то есть изображение, где может быть извлечено из с фильтром нижних частот кирпичной стены, который сохраняет исходное изображение (где ) немасштабированным и удаляет все другие изображения. Это означает, что он умножает исходное изображение на 1 и умножает все остальные изображения на 0.Икс(J2πе)Кзнак равно0Иксs(J2πе)Кзнак равно0

Икс(J2πе)знак равнопрямоугольник(ееs)Иксs(J2πе)знак равноЧАС(J2πе) Иксs(J2πе)

фильтр реконструкции

Фильтр реконструкции является

ЧАС(J2πе)знак равнопрямоугольник(ееs)

и имеет акаузальный импульсный отклик :

час(T)знак равноF-1{ЧАС(J2πе)}знак равноеsсинк(еsT)

Эта операция фильтрации, выраженная как умножение в частотной области, эквивалентна свертке во временной области:

Икс(T)знак равночас(T)Иксs(T)знак равночас(T)T ΣNзнак равно-+Икс[N] δ(T-NT)знак равноT ΣNзнак равно-+Икс[N] (час(T)δ(T-NT))знак равноT ΣNзнак равно-+Икс[N] час(T-NT))знак равноT ΣNзнак равно-+Икс[N] (еsсинк(еs(T-NT)))знак равноΣNзнак равно-+Икс[N] синк(еs(T-NT))знак равноΣNзнак равно-+Икс[N] синк(T-NTT)

Это ясно объясняет, как исходный восстанавливается из выборок и знания частоты дискретизации или периода выборки.Икс(T)Икс[N]


Так что то, что выводится из практического цифро-аналогового преобразователя (ЦАП), не является ни

ΣNзнак равно-+Икс[N] синк(T-NTT)

который не нуждается в дополнительном лечении для восстановления , ниИкс(T)

Иксs(T)знак равноΣNзнак равно-+Икс[N] Tδ(T-NT)

который с помощью идеальной кирпичной стены LPF восстанавливает , изолируя и сохраняя изображение основной полосы частот и отбрасывая все другие изображения.Икс(T)

Выход ЦАП

Что получается из обычного ЦАП, если не выполняется обработка или масштабирование для оцифрованного сигнала, так это значение сохраняемое при постоянном значении до тех пор, пока не будет выведен следующий образец. Это приводит к кусочно-постоянной функции :Икс[N]

ИксDAC(T)знак равноΣNзнак равно-+Икс[N] прямоугольник(T-NT-T2T)

Обратите внимание на задержку периода выборки примененную к функции . Это делает это причинно-следственной. Это просто означает, что12прямоугольник()

xDAC(t)=x[n]=x(nT)whennTt<(n+1)T

Заявлено по-другому

xDAC(t)=x[n]=x(nT)forn=floor(tT)

где - это функция пола , определенная как наибольшее целое число, не превосходящее .floor(u)=uu

Этот выход ЦАП напрямую моделируется как линейная не зависящая от времени система (LTI) или фильтр, который принимает идеально дискретизированный сигнал и для каждого импульса в идеально дискретизированном сигнале выводит эту импульсную характеристику:xs(t)

hZOH(t)=1Trect(tT2T)

Подключить, чтобы проверить это ...

xDAC(t)=hZOH(t)xs(t)=hZOH(t)T n=+x[n] δ(tnT)=T n=+x[n] (hZOH(t)δ(tnT))=T n=+x[n] hZOH(tnT))=T n=+x[n] 1Trect(tnTT2T)=n=+x[n] rect(tnTT2T)

Выход ЦАП , как выход системы LTI с импульсной характеристикой согласуется с кусочно-постоянной конструкцией выше. И входом в эту систему LTI является выборочный сигнал разумно масштабированный таким образом, чтобы изображение основной полосы частот точно таким же, как спектр исходного дискретизированного сигнала . То естьxDAC(t)hZOH(t)xs(t)xs(t)x(t)

X(j2πf)=Xs(j2πf)forfs2<f<+fs2

Спектр исходного сигнала такой же, как у дискретизированного спектра, но все изображения, появившиеся в результате дискретизации, отбрасываются.

Передаточная функция этой системы LTI, которую мы называем удержанием нулевого порядка (ZOH) , является преобразованием Лапласа импульсного отклика:

HZOH(s)=L{hZOH(t)}+hZOH(t) est dt=+1Trect(tT2T) est dt=0T1T est dt=1T1sest|0T=1esTsT

Частотная характеристика получается путем подстановкиj2πfs

HZOH(j2πf)=1ej2πfTj2πfT=ejπfTejπfTejπfTj2πfT=ejπfTsin(πfT)πfT=ejπfTsinc(fT)=ejπfTsinc(ffs)

Это указывает на линейный фазовый фильтр с постоянной задержкой, равной половине периода выборки, , и с коэффициентом усиления, который уменьшается с увеличением частоты . Это мягкий эффект фильтра нижних частот. При постоянном токе коэффициент усиления равен 0 дБ, а при Найквисте коэффициент усиления равен -3,9224 дБ. Таким образом, в изображении основной полосы частот некоторые высокочастотные компоненты немного уменьшены.T2ff=0f=fs2

Как и в случае дискретизированного сигнала , в дискретизированном сигнале имеются изображения с целыми кратными частоты дискретизации, но эти изображения значительно уменьшены по амплитуде (по сравнению с изображение основной полосы), потому чтопроходит через ноль, когда для целого числа , которое не равно 0, что находится прямо в середине этих изображений.xs(t)xDAC(t)|HZOH(j2πf)|f=kfsk

Заключительное:

  1. Удержание нулевого порядка (ZOH) - это линейная не зависящая от времени модель восстановления сигнала, выполняемая практическим цифро-аналоговым преобразователем (ЦАП), который поддерживает выходную константу при значении выборки, , до тех пор, пока не будет обновлено следующий образец .x[n]x[n+1]

  2. Вопреки распространенному заблуждению, ZOH не имеет ничего общего со схемой выборки и удержания (S / H), которую можно обнаружить перед аналого-цифровым преобразователем (АЦП) . Пока ЦАП поддерживает выходное значение постоянным в течение каждого периода выборки, не имеет значения, имеет ли АЦП сигнал S / H или нет, эффект ZOH сохраняется. Если ЦАП выводит что-то отличное от кусочно-постоянного вывода (например, последовательность узких импульсов, предназначенных для аппроксимации импульсов Дирака), изображенных выше как , то эффект ZOH отсутствует (что-то еще вместо этого) есть ли цепь S / H, предшествующая АЦП или нет.xDAC(t)

  3. Чистая передаточная функция ZOH - а чистая частотная характеристика ZOH - Во многих учебниках пропущен фактор в знаменателе передаточной функции, и это ошибка.

    HZOH(s)=1esTsT
    HZOH(j2πf)=ejπfTsinc(fT)
    T

  4. ZOH значительно уменьшает изображения дискретизированного сигнала , но не устраняет их. Для устранения изображений нужен хороший фильтр нижних частот, как и раньше. Кирпичные стены LPF являются идеализацией. Практический ФНЧ также может ослаблять изображение в основной полосе (которое мы хотим сохранить) на высоких частотах, и это затухание должно учитываться как затухание, возникающее в результате ZOH (которое составляет менее 3,9224 дБ затухания). ZOH также задерживает сигнал на половину периода выборки, что может потребоваться принять во внимание (наряду с задержкой фильтра LPF против формирования изображения), особенно если ZOH находится в контуре обратной связи.xs(t)


Я признаю, что ваш ответ чище и немного тщательнее, чем мой. Я все еще оставался, задаваясь вопросом, что было большим показом? Может быть, вы хотели бы подчеркнуть трюм нулевого порядка в качестве модели в DAC-выход?
Тимо

в вашем ответе есть некоторые ошибки. например, он не показывает 1/2 выборки задержки в частотной характеристике. извините, что так получилось, что наша награда (которая была моей и теперь должна быть вашей ) ушла в туалет.
Роберт Бристоу-Джонсон

eiπfT

@ Тимо, теперь у тебя в два раза больше, чем у меня. когда ты собираешься опубликовать награду, на которую я могу нанести удар? :-)
Роберт Бристоу-Джонсон

Справедливо, я должен попытаться что-то придумать: D
Тимо
Используя наш сайт, вы подтверждаете, что прочитали и поняли нашу Политику в отношении файлов cookie и Политику конфиденциальности.
Licensed under cc by-sa 3.0 with attribution required.