Расчет шумов оконечного оборудования. Квантование по уровню в цифровых системах Дисперсия шума квантования при округлении
Основным аспектом расчета и разработки инженерных проектов является необходимость использования аналитических характеристик качества функционирования систем. Только при наличии таких характеристик система может быть объективно оценена и ее стоимость эффективно сопоставлена со стоимостью альтернативных разработок. Одной из характеристик, наобходимой для инженеров телефонной связи, является качество речи, доставляемой слушателю. Измерения качества речи усложняются из-за субъективных свойств речи, которые воспринимаются типовым слушателем. Одна из особенностей субъективного восприятия шумов или искажений в речевом сигнале связана с частотным составом, или спектром мешающих воздействий в сочетании с уровнем их мощности. Эти эффекты влияния шумов в зависимости от частоты были рассмотрены в гл.1 при введении понятий взвешивания при помощи С-контура и псофометрического взвешивания.
Последовательные ошибки квантования в ИКМ-кодере в общем случае предполагаются распределенными по случайному закону и не коррелированными друг с другом. Таким образом, совокупный эффект ошибок квантования в системах с ИКМ можно рассматривать как аддитивные шумы, имеющие субъективное воздействие, которое аналогично воздействию белого шума с ограниченной полосой. На рис. 3.9 представлена зависимость ошибок кванования от амплитуды сигнала для кодера с равномерными шагами квантования. Отметим, что если сигнал успевает измениться по амплитуде на несколько шагов квантования, ошибки квантования становятся независимыми. Если сигнал дискретизируется с частотой, намного превышающей f s , то последовательные дискреты будут часто приходиться на одни и те же шаги, что приведет к потере независимости ошибок квантования.
Ошибки, или шум квантования, возникающие при преобразовании аналогово сигнала в цифровую форму, обычно выражаются в виде средней мощности шума по отношению к средней мощности сигнала. В соответствии с этим отношение сигнал-шум квантования можно определить как
ОСШК=E{x 2 (t)}/E{ 2 }, (3.1)
где E{. } – математическое ожидание, или среднее значение,x(t) – аналоговый входной сигнал,y(t) – декодированный выходной сигнал.
При определении среднего значения шума квантования необходимо сделать три замечания.
Ошибка y(t) –x(t) ограничена по амплитуде значениемq/2, гдеq– шаг квантования. (Декодированные выходные дискреты располагаются точно посредине шага квантования.)
Можно предположить, что значения дискретов с равной вероятностью могут попадать в любую точку в пределах шага квантования (предполегается равномерная плотность вероятности, равная 1/q).
Предполагается, что амплитуды сигнала ограничены рабочим диапазоном кодера. Если значение дискрета превышает границу наивысшего шага квантования, то возникают искажения, вызванные перегрузкой.
Если для удобства предположим, что нагрузочный резистор имеет сопротивление 1 Ом, то средняя мощность шума квантования (вычислена в приложении А) определяется выражением:
Мощность шума квантования = q 2 /12. (3.2)
Если все шаги квантования имеют равные значения (равномерное квантование) и шум квантования не зависит от значений дискретов, то отношение сигнал-шум квантования (в децибелах) определяется как
ОСШК = 10lg = 10,8 + 20lg(v/q), (3.3)
где v– среднее квадратическое значение амплитуды входного сигнала. В частности, для синусоидального входного сигнала отношение сигнал-шум квантования (в децибелах) при равномерном квантовании
ОСШК = 10lg[(А 2 /2)/(q 2 /12)] = 7,78 + 20lg(А/q), (3.4)
где А – амплитуда синусоиды.
Пример 3.1 Синусоидальный сигнал с амплитудой 1 В следует преобразовать в цифровую форму таким образом, чтобы получить отношение сигнал-шум квантования не менее 30 дБ. Сколько потребуется одинаковых шагов квантования и сколько потребуется разрядов для кодирования каждого дискрета?
Решение. С помощью формулы (3.4) определяем максимальный размер шага квантованияq=10 – (30 – 7,78) / 20 = 0,078B.
Таким образом, потребуется 13 шагов квантования для каждой полярности сигнала (общее число шагов квантования 26). Число разрядов, необходимых для кодирования каждого дискрета, определяется как n=log 2 26 = 4,75 разрядов на дискрет.
При измерениях мощности шума квантования спектральные составляющие часто взвешивают тем же способом, что и шумы в аналоговых каналах. К сожалению, измерения взвешенных шумов не всегда отражают подлинное качество восприятия речи, прошедшей кодер (декодер). Если спектральное распределение шумов квантования более или менее повторяет спектральное распределение сигнала речи, эти шумы значительно менее заметны, чем шумы, некоррелированные с речью . С другой стороны, если процесс квантования создает энергию на тональных частотах, отличных от тех, которые содержатся в конкретных звуках, эти искажения становятся более заметными.
Высококачественные ИКМ-кодеры создают шум квантования, который равномерно распределен в диапазоне ТЧ и не зависит от кодируемого сигнала. В этом случае отношение сигнал-шум квантования (3.4) является хорошей мерой качества ИКМ-преобразования. В некоторых видах кодеров, рассматриваемых далее (в особенности в вокодерах), знание мощности шума квантования не приносит большой пользы. В описаны другие характеристики качества речи, прошедшей через кодер, которые лучше определяют восприятие речи слушателем.
Однако логично предположить, что промежуточные участки волны, которые оцифровщику не удалось достоверно отобразить, не могут просто так исчезнуть.
Шумы квантования
Между аналоговым сигналом и его цифровой копией в вашей системе записи существует разница, которая называется ошибками квантования , или шумами квантования .
С помощью несложных математических формул можно вычислить частоту и уровень громкости шумов квантования. Также их характер можно проследить наглядно, если проанализировать отклонения графика оцифрованной волны от оригинальной синусоиды. На рисунке справа показана разница между исходным и оцифрованным сигналом.
Шумы квантования - это неотъемлемая составляющая цифрового звука, они возникают в момент оцифровки. Для минимизации влияния этих шумов на звук в конструкциях конверторов используются специальные фильтры. Покупая оцифровщик с более дорогими характеристиками (например, 24 /192 ), многие не обращают внимания качество этих фильтров, ориентируясь лишь на красивые численные характеристики разрядности и частоты дискретизации.
Чем выше показатели конвертора , тем дороже должны быть фильтры , однако именно на них производители обычно экономят, чтобы сохранить себестоимость на низком уровне и обеспечить себе конкурентоспособность.
Алиасинг
Еще одна неприятная вещь, которая может произойти в процессе семплирования (оцифровки) звука, называется алиасингом. Алиасинг - наложение двух непрерывных сигналов разной частоты друг на друга при семплировании, в результате которого в звуке возникают искажения.
Мы можем представить алиасинг даже визуально. Вспомните вращение колес автомобилей или поездов в старых фильмах. В определенные моменты можно отчетливо заметить, что колеса как бы крутятся в обратную сторону. И это не обман зрения, этот эффект появляется в моменты, когда частота вращения колес приближается к кадровой частоте кинокамеры (обычно это 24 кадра в секунду, но когда-то это значение было на уровне 16-20). Каждая точка колеса, двигаясь по часовой стрелке, успевает пройти почти полный оборот за один кадр, оказываясь с обратной стороны исходной точки, как будто эта точка сдвинулась против часовой стрелки. И мы видим обратное вращение.
В результате алиасинга записанный сигнал отличается от ожидаемого.
В соответствии с теоремой Котельникова , для восстановления сигнала без потерь семплирование должно производиться с частотой, в два раза превышающей самую высокую частоту в записываемом спектре.
То есть, скажем, если максимальная скорость вращения колес составляет 10 оборотов в секунду, то для устранения эффекта алиасинга фиксировать этот движение нужно с частотой не менее 20 кадров в секунду. А кинокамера – этот тот же семплер, только записывающий не звук, а изображение. При указанных значениях, как бы ни крутилось колесо, камера за один его оборот успеет сделать два семпла, а значит обратного вращения мы уже не увидим.
Так что если нам надо записать звук в пределах 20 кГц (верхний порог идентифицируемых человеческим ухом частот), то семплирование должно происходить с частотой дискретизации не менее 40 кГц.
При этом половина частоты дискретизации называется числом Найквиста (Найквист и Котельников – ученые, которые независимо друг от друга занимались исследованиями данной проблемы).
Однако мы знаем, что даже если наше ухо не распознает какие-то частоты, это еще не значит, что их нет. А раз они есть, то семплер (оцифровщик) попытается их зафиксировать, работая при этом на недостаточной для записи этого спектра частоте дискретизации. И возникнет алиасинг.
Чтобы устранить негативный эффект от алиасинга, при семплировании требуется частота дискретизации с запасом более чем в два раза . Кроме того, необходимо на входе оцифровщика применять фильтры , отсекающие нежелательные частоты выше определенного значения.
Именно поэтому используемые в звукозаписи «стандартные» частоты дискретизации выше 40 кГц – 44.1 и 48 кГц: такое семплирование обеспечивает запас для устранения искажений.
В можно поочередно услышать «хорошую» и «плохую» запись пилообразной волны на частотах 440, 880 и 1760 Гц. В первом варианте были применены фильтры, а во втором отчетливо слышен алиасинг.
Сегодня уже никого не удивишь даже значениями 32 бита или 96–192 кГц. С каждым годом производители «улучшают» характеристики приборов. Но поскольку, как я уже говорил, для фильтрации более высоких частот требуются более качественные и дорогие фильтры, нередко получается, что конвертор, работающий в режиме 16/44.1, дает более качественный звук, чем конвертор 24/192. Шумы квантования, алиасинг и отсутствие хороших фильтров делают свое дело. И это мы еще опускаем возможные погрешности, связанные с повышенной нагрузкой на систему при работе с более высокими параметрами звука.
Если статья оказалась полезной, вы можете подписаться на обновления этого блога , чтобы бесплатно получать новые материалы на электронную почту. Или вступайте
При правильно выбранной частоте дискретизации, исходя из теоремы Котельникова, точность преобразования аналогового ЗС в цифровой определяется величиной шага квантования.Погрешность преобразования тем меньше, чем меньше шаг квантования. Разность между исходным и квантованным значениями сигнала в дискретные моменты времени называется шумом квантования (ошибкой квантования).
Шум квантования в отличие от флуктуационного шума, в общем случае, носит неслучайный характер. Поэтому правильнее говорить об искажениях сигнала при его аналого-цифровом преобразовании. При фиксированном максимальном уровне входного аналогового ЗС шум квантования определяется числом уровней квантования – разрядностью аналого-цифрового преобразователя (АЦП).
При кодировании двоичными числами и длине кодового слова в m разрядов количество двоичных кодовых слов r (разрешающая способность) составляет. Так при m=16, r=65536.
Поток кодовых слов на выходе АЦП характеризуется скоростью передачи данных – число бит информации переданных за 1 секунду. Скорость передачи данных есть произведение числа разрядов кодового слова на частоту дискретизации (в герцах). Объем памяти необходимый для хранения информации о реализации ЗС длительностью, определяется как произведение скорости потока данных на длительность сигнала.
При линейной импульсно-кодовой модуляции (ИКМ), т.е. при равномерном шаге квантования, мощность шума квантования определяется только шагом квантования:
где - общий динамический диапазон сигнала.
Эффективное значение ошибки квантования:
Шум квантования представляет собой, при линейной ИКМ, случайный процесс с равномерным расширением в пределах, с плотностью вероятности. Спектр шума квантования равномерный в полосе частот.
Шум квантования проявляется только при наличии сигнала. При отсутствии сигнала на входе АЦП на выходе АЦП будут иметь место квантование колебания в младшем разряде АЦП. Объясняется это наличием теплового шума входных аналоговых частей АЦП, нестабильностью питающего напряжения, дрейфом постоянной составляющей в усилителях постоянного тока и другими причинами. На выходе ЦАП (цифро-аналогового преобразователя) это квантованное колебание превращается в шум, называемый шумом паузы. Шум паузы менее равномерный, чем белый шум, характерный для аналоговых устройств, и его часто называют гранулированным. Мощность шума паузы:
на 4,7 дБ больше шума квантования.
Поскольку не зависит от уровня входного сигнала, с увеличением мощности входного отношение линейно растет до тех пор пока не возникают шумы ограничения. Уровень ограничения по входу АЦП определяется максимальным входным рабочим напряжением АЦП. Шумом ограничения называется разность между исходным и ограниченным сигналами. Система АЦП рассчитывается таким образом, чтобы ограничения не возникало т.е.
здесь R- пик-фактор сигнала; S СР – среднеквадратическое значение сигнала.
Число шагов можно определить из соотношения:
где - максимальное и минимальное значения сигнала на входе АЦП;
Шаг квантования.
С учетом выражений (9.6), (9.9),(9.10) получим выражение для мощности шума
Мощность сигнала на сопротивлении 1 Ом,тогда
или в децибелах
При m- разрядном кодировании, тогда
У гармонического сигнала пик-фактор, в этом случае
Для сигналов вещания пик-фактор зависит от жанра программы. Если в среднем считать R=13 дБ то
Если учитывать неодинаковую чувствительность слуха слушателя к составляющим шума разных частот, то отношение сигнал/шум квантования уменьшается на 8.5 дБ для сигнала в полосе частот до 15 кГц и составляет
Динамический диапазон цифрового сигнала оценивают величиной, дБ с учетом того, что получаем
Из выражения (9.15) видно, что увеличение числа разрядов на единицу приводит к улучшению отношения сигнал/шум на 6 дБ.
На рис.9.2. изображены зависимости отношения сигнал/шум для сигналов 3В при разных значениях m от уровня входного сигнала (9.17).
При 16-разрядном квантовании имеем для гармонического сигнала D=90 дБ, С-Ш=98 Б (из выражений 9.15, 9.18). Отношение С-Ш при расчете по формуле (9.17) получается равным 80дБ при кодировании сигнала максимального по уровню. При кодировании слабых сигналов отношение С-Ш меньше на величину динамического диапазона сигнала и оказывается весьма малым при D=50…60 дБ.
80 -70 -60 -50 -40 -30 -20 -10 0
Рис.9.2. Отношение сигнал/шум при ИКМ
В диапазоне оптических частот тепловой шум оказывается очень слабым. Однако в этом диапазоне при слабых сигналах существенное значение имеет "квантовый шум", вызванный дискретной природой светового излучения. Согласно квантовой теории электромагнитного поля его энергия сигнала излучается и поглошается квантами, причём энергия одного такого кванта (фотона) равна . В элементарном сигнале длительности с высокостабильной несущей частотой (когерентное одномодовое излучение) и амплитудой детерминированной может быть только средняя энергия (пропорциональная ( - среднее число фотонов на интервале Т). Конкретная же реализация элементарного сигнала имеет энергию где случайное число регистрируемых фотонов.
В современных системах оптической связи в основном используется АМ оптического несущего колебания по амплитуде или интенсивности (мощности).
Идеальная система оптической связи при изохронной передаче двоичных сообщений (1 и 0) имеет следующие характеристики:
1. Время передачи бита (тактовый интервал) постоянен и равен следовательно, скорость передачи информации
2. При передаче 1 оптическая энергия, излучаемая в виде импульсов за время передачи одного бита, где число излучённых фотонов,
Энергия одного фотона (кванта), а оптическая энергия при передаче 0 равна нулю. Оптическая энергия в месте приёма равна на тактовом интервале величине при передаче 1 и нулю при передаче 0 соответственно.
3. Вероятности передачи . В этом случае усреднённую за продолжительное время принимаемую мощность можно выразить через среднюю мощность принимаемую за время передачи бита при посылке 1. Таким образом,
Реальная система оптической связи отличается от идеальной следующим :
1. Время передачи бита информации не остаётся постоянным - этот эффект называют фазовым дрожанием цифрового сигнала.
2. Излучаемая оптическая энергия не остаётся строго одной и той же. При передаче как кодовой 1, так и кодового 0 имеет место шум передатчика, приводящий к случайным изменениям амплитуды от импульса к импульсу. Кроме того, имеет место "шум лазера", обусловленный статистической природой взаимодействия между возбуждением лазера и создаваемым потоком фотонов. Флуктуации принимаемой энергии увеличиваются ещё больше из-за изменений затухания в канале связи. Кроме того, появляются флуктуации энергии на отдельных тактовых интервалах в месте приёма, обусловленные статистической природой взаимодействия потока фотонов (оптический сигнал) и создаваемого фотодетектором (обычно это фотодиод) потока электронно-дырочных пар. Условно будем говорить в этом случае о шуме фотодетектора.
3. Весьма вероятно, что при передаче 0 излучается малый, но вполне определённый уровень энергии (шум лазера), не считая шума передатчика и канала. Отношение средней энергии, принимаемой при передаче 0, к средней энергии при передаче 1 характеризуется коэффициентом Полагают, что в идеальной системе однако обычно это не так, особенно если лазерный источник излучения смещён вблизи порога генерации.
4. Конечная длительность излучаемых импульсов и дополнительная временная дисперсия (рассеяние) при их передаче по каналу приводят к тому, что в практических системах связи происходит наложение соседних посылок, т.е. проявляется межсимвольная интерференция.
Шум лазера, о котором говорилось выше, имеет квантовую природу. Вероятность появления точно фотонов на интервале на передающей стороне определяется распределением Пуассона (см. § 2.76):
Таким образом, шум лазера - это "квантовый шум", так как проявляется во флуктуациях параметров сигнала, детерминированного по классическим представлениям. Этот шум не является аддитивным, так как зависит от самого полезного сигнала. С учётом этого в приведённой формуле следует считать, что при передаче а при передаче Как указывалось выше, при передаче 0 (отсутствие возбуждения лазера) может наблюдаться определённый, хотя и малый уровень энергии, обусловленный тем, что вероятность непоявления фотонов на этом интервале где среднее число шумовых фотонов на интервале при отсутствии возбуждения лазера. По мере увеличения средней мощности излучаемого сигнала Рпер вклад квантового шума по сравнению с другими шумами тракта передачи падает.
Шум фотодетектора имеет природу, аналогичную шуму лазера, так как падающий на фотодиод стационарный световой поток генерирует электронно-дырочные пары носителей заряда как независимые случайные события. Если за отрезок времени на фотодиод падает оптическая энергия, равная в среднем то следует ожидать, что будет создано в среднем пар носителей заряда, причём
Эффекты конечной разрядности цифровых фильтров
При анализе ЛДС предполагалось, что сигналы дискретизируются только по времени и дискретные отсчеты и коэффициенты фильтров представлялись с неограниченной точностью. Однако в реальных или цифровых системах точность вычислений ограничена и зависит от числа разрядов используемых устройств: АЦП, регистров, сумматоров, умножителей. Это обстоятельство приводит к следующим эффектам:
Шум квантования при аналого-цифровом преобразовании;
Округление результатов промежуточных вычислений;
Искажения частотных характеристик из-за квантования коэффициентов цифровых фильтров;
Переполнение разрядной сетки в процессе вычислений;
Предельные циклы малого уровня.
Шум квантования
Под шумами квантования понимают случайные ошибки между дискретными во времени отсчетами сигналов и их цифровым представлением с ограниченной разрядностью.
Соседние отсчеты шумов квантования предполагаются некоррелированными между собой. И шум квантования является «белым».
Плотность вероятности шума квантования соответствует равномерному закону распределения:
, (2.1)
где - интервал квантования по уровню.
Дисперсия шумов квантования определяется законом распределения:
. (2.2)
Если максимальное значение квантованного сигнала равно , то интервал квантования равен:
,
где - число разрядов цифрового устройства.
Спектральная плотность мощности шумов квантования определяется выражением:
,
где DТ – интервал дискретизации.
Например, коэффициент усиления приемного устройства до входа АЦП выбирается таким образом, чтобы уровень тепловых приемного устройства превышал спектральную плотность шумов квантования.
Вклад входных шумов квантования в выходной сигнал цифрового фильтра определяется выражением:
.
Соответственно, дисперсия выходного шума квантования:
Литература
Маркович И.И. Цифровая обработка сигналов в системах и устройствах: монография / И.И. Маркович; Южный федеральный университет. – Ростов н/Д: Издательство Южного федерального университета, 2012. – 236 с.
Основы цифровой обработки сигналов: учебное пособие / Ю.А. Брюханов, А.А. Приоров, В.И. Джиган, В.В. Хрящев; Яросл. гос. ун-т им. П.Г. Демидова. - Ярославль: ЯрГУ, 2013. – 344 с. (с. 152)
Карташов В.Г. Основы теории дискретных сигналов и цифровых фильтров. – М.: Высш. школа, 1982. – 109. (с. 86)
