Введение в цифровую обработку сигналов

Размещено на http://www.allbest.ru/

Тема. введение в цифровую обработку сигналов

цифровая обработка сигналов

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

Франсуа де Ларошфуко. Французский писатель моралист. XVII в.

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

Виль Ибрагимов. Ташкентский геофизик Уральской школы, ХХ в.

Введение

Цифровая обработка сигналов (ЦОС или DSP - digital signal processing) является одной из новейших и самых мощных технологий, которая активно внедрилась в широкий круг областей науки и техники: коммуникации, метеорология, радиолокация и гидролокация, медицинская визуализация изображений, цифровое аудио- и телевизионное вещание, разведка нефтяных и газовых месторождений, и многих других. Можно сказать, что происходит повсеместное и глубокое проникновение технологий цифровой обработки сигналов во все сферы деятельности человечества. Сегодня технология ЦОС относится к числу базовых знаний, которые необходимы ученым и инженерам всех отраслей без исключения.

Изучение принципов цифровой обработки сигналов не требует глубоких и всесторонних знаний по математике. Необходимо лишь немного опыта в элементарной алгебре, в способах декомпозиции (разложения) произвольных сигналов на составляющие и их обратной реконструкции. Почему же предмет считается сложным для понимания? Ответ в том, как материал преподносится для изучения. Одно дело писать уравнения, но другое - объяснить их значение с практической точки зрения. Именно в последнем и состоит цель данного курса.

Физические величины макромира, как основного объекта наших измерений и источника информационных сигналов, как правило, имеют непрерывную природу и отображаются непрерывными (аналоговыми) сигналами. Цифровая обработка сигналов оперирует с дискретными величинами, причем с квантованием как по координатам динамики своих изменений (во времени, в пространстве, и по любым другим изменяемым аргументам), так и по значениям физических величин. Математика дискретных преобразований зародилась в недрах аналоговой математики еще в 18 веке в рамках теории рядов и их применения для интерполяции и аппроксимации функций, однако ускоренное развитие она получила в 20 веке после появления первых вычислительных машин. В своих основных положениях математический аппарат дискретных преобразований подобен преобразованиям аналоговых сигналов и систем. Однако дискретность данных требует учета этого фактора, а его игнорирование может приводить к ошибкам. Кроме того, ряд методов дискретной математики не имеет аналогов в аналитической математике.

Стимулом развития дискретной математики является и то, что стоимость цифровой обработки данных меньше аналоговой и продолжает снижаться, а производительность вычислительных операций непрерывно возрастает. Немаловажным является и то, что системы ЦОС отличаются высокой гибкостью. Их можно дополнять новыми программами и перепрограммировать на выполнение различных операций без изменения оборудования. В последние годы ЦОС оказывает возрастающее влияние на все отрасли современной промышленности: телекоммуникации, средства информации, цифровое телевидение и пр. Интерес к научным и к прикладным вопросам цифровой обработки сигналов возрастает во всех отраслях науки и техники.

1. ПРЕДИСЛОВИЕ К ЦИФРОВОЙ ОБРАБОТКЕ СИГНАЛОВ [52]

Цифровые сигналы формируются из аналоговых операцией дискретизации - последовательным квантованием (измерением) амплитудных значений сигнала через определенные интервалы времени ?t или любой другой независимой переменной ?x. В результате равномерной дискретизации непрерывный по аргументу сигнал переводится в упорядоченную по независимой переменной последовательность чисел. В принципе разработаны методы ЦОС для неравномерной дискретизации данных, однако области их применения достаточно специфичны и ограничены. Условия, при которых возможно полное восстановление аналогового сигнала по его цифровому эквиваленту с сохранением всей исходно содержавшейся в сигнале информации, выражаются теоремами Найквиста, Котельникова, Шеннона, сущность которых практически одинакова. Для дискретизации аналогового сигнала с полным сохранением информации в его цифровом эквиваленте максимальные частоты в аналоговом сигнале должны быть не менее чем вдвое меньше, чем частота дискретизации, то есть fmax (1/2)fd, т.е. на одном периоде максимальной частоты должно быть минимум два отсчета. Если это условие нарушается, в цифровом сигнале возникает эффект маскирования (подмены) действительных частот более низкими частотами. При этом в цифровом сигнале вместо фактической регистрируется "кажущаяся" частота, а, следовательно, восстановление фактической частоты в аналоговом сигнале становится невозможным. Восстановленный сигнал будет выглядеть так, как если бы частоты, лежащие выше половины частоты дискретизации, отразились от частоты (1/2)fd в нижнюю часть спектра и наложились на частоты, уже присутствующие в этой части спектра. Этот эффект называется наложением спектров или алиасингом (aliasing). Наглядным примером алиасинга может служить иллюзия, довольно частая в кино - колесо автомобиля начинает вращаться против его движения, если между последовательными кадрами (аналог частоты дискретизации) колесо совершает более чем пол-оборота.

Преобразование сигнала в цифровую форму выполняется аналого-цифровыми преобразователями (АЦП). Как правило, они используют двоичную систему счисления с определенным числом разрядов в равномерной шкале. Увеличение числа разрядов повышает точность измерений и расширяет динамический диапазон измеряемых сигналов. Потерянная из-за недостатка разрядов АЦП информация невосстановима, и существуют лишь оценки возникающей погрешности «округления» отсчетов, например, через мощность шума, порождаемого ошибкой в последнем разряде АЦП. Для этого используется понятие отношения «сигнал/шум» - отношение мощности сигнала к мощности шума (в децибелах). Наиболее часто применяются 8-, 10-, 12-, 16-, 20- и 24-х разрядные АЦП. Каждый дополнительный разряд улучшает отношение сигнал/шум на 6 децибел. Однако увеличение количества разрядов снижает скорость дискретизации и увеличивает стоимость аппаратуры. Важным аспектом является также динамический диапазон, определяемый максимальным и минимальным значением сигнала.

Обработка цифровых сигналов выполняется либо специальными процессорами, либо на универсальных ЭВМ и компьютерах по специальным программам. Наиболее просты для рассмотрения линейные системы. Линейными называются системы, для которых имеет место суперпозиция (отклик на сумму входных сигналов равен сумме откликов на каждый сигнал в отдельности) и однородность или гомогенность (изменение амплитуды входного сигнала вызывает пропорциональное изменение выходного сигнала). Для реальных объектов свойства линейности могут выполняться приближенно и в определенном интервале входных сигналов.

Если входной сигнал x(t-t0) порождает однозначный выходной сигнал y(t-t0) при любом сдвиге t0, то систему называют инвариантной во времени. Ее свойства можно исследовать в любые произвольные моменты времени. Для описания линейной системы вводится специальный входной сигнал - единичный импульс (импульсная функция). В силу свойства суперпозиции и однородности любой входной сигнал можно представить в виде суммы таких импульсов, подаваемых в разные моменты времени и умноженных на соответствующие коэффициенты. Выходной сигнал системы в этом случае представляет собой сумму откликов на эти импульсы. Отклик на единичный импульс (импульс с единичной амплитудой) называют импульсной характеристикой системы h(n). Соответственно, отклик системы на произвольный входной сигнал s(k) можно выразить сверткой

g(k) = h(n) ? s(k-n).

Если h(n)=0 при n<0, то систему называют каузальной (причинной). В такой системе реакция на входной сигнал появляется только после поступления сигнала на ее вход. Некаузальные системы физически невозможно реализовать в реальном масштабе времени. Если требуется реализовать свертку сигналов с двусторонними операторами (при дифференцировании, преобразовании Гильберта, и т.п.), то это выполняется с задержкой (сдвигом) входного сигнала минимум на длину левосторонней части оператора свертки.

Z-преобразование. Для анализа дискретных сигналов и систем широко используется z-преобразование, которое является обобщением дискретного преобразования Фурье. Этим преобразованием произвольной непрерывной функции s(t), равномерно дискретизированной и отображенной отсчетами sk = s(k?t), ставится в соответствие степенной полином по z (или степенной полином по z-1 = 1/z), последовательными коэффициентами которого являются отсчеты функции:

sk = s(kt) TZ[s(kt)] =sk zk = S(z),

где z = +jv = rexp(-j) - произвольная комплексная переменная. Это преобразование позволяет в дискретной математике использовать всю мощь дифференциального и интегрального исчисления, алгебры и прочих хорошо развитых разделов аналитической математики.

Дискретные системы обычно описывается линейными разностными уравнениями с постоянными коэффициентами:

y(k) = ? b(n) x(k-n) - ? a(m) y(k-m), n=0, 1, … , N, m=1, 2, … , M.

Этим уравнением устанавливается, что выходной сигнал y(k) системы в определенный момент ki (например, в момент времени kit) зависит от значений входного сигнала x(k) в данный (ki) и предыдущие моменты (ki-n) и значений сигнала y(k) в предыдущие моменты (ki-m).

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

H(z) = Y(z)/X(z),

представляет собой рациональную функцию в виде отношения двух полиномов от z. Корни полинома в числителе называются нулями, а в знаменателе - полюсами функции H(z). Значения нулей и полюсов позволяют определить свойства линейной системы. Так, если все полюсы X(z) по модулю больше единицы, то система является устойчивой (не пойдет “вразнос” ни при каких входных воздействиях). Нули функции Y(z) обращают в ноль H(z) и показывают, какие колебания вовсе не будут восприниматься системой (“антирезонанс”). Систему называют минимально-фазовой, если все полюсы и нули передаточной функции лежат вне единичной окружности |z|=1 на комплексной z-плоскости. Попутно заметим, что применение z-преобразования с отрицательными степенями z-1 меняет положение полюсов и нулей относительно единичной окружности |z|=1 (область вне окружности перемещается внутрь окружности, и наоборот).

Природа сигналов. По своей природе сигналы могут быть случайными или детерминированными.

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

Случайные сигналы непредсказуемы по своим значениям во времени или в пространстве. Для каждого конкретного отсчета случайного сигнала можно знать только вероятность того, что он примет какое-либо значение в определенной области возможных значений. Закон распределения случайных значений далеко не всегда известен. Одним из самых распространенных является нормальное распределение, плотность которого имеет вид симметричного колокола. Для его описания достаточно двух первых моментов распределения случайных величин.

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

Функциональные преобразования сигналов. Одним из основных методов частотного анализа и обработки сигналов является преобразование Фурье. Различают понятия “преобразование Фурье” и “ряд Фурье”. Преобразование Фурье предполагает непрерывное распределение частот, ряд Фурье задается на дискретном наборе частот. Сигналы также могут быть заданы в наборе временных отсчетов или как непрерывная функция времени. Это дает четыре варианта преобразований - преобразование Фурье с непрерывным или с дискретным временем, и ряд Фурье с непрерывным временем или с дискретным временем. Наиболее практична с точки зрения цифровой обработки сигналов дискретизация и во временной, и в частотной области, но не следует забывать, что она является аппроксимацией непрерывного преобразования. Непрерывное преобразование Фурье позволяет точно представлять любые явления. Сигнал, представленный рядом Фурье, может быть только периодичен. Сигналы произвольной формы могут быть представлены рядом Фурье только приближенно, т.к. при этом предполагается периодическое повторение рассматриваемого интервала сигнала за пределами его задания. На стыках периодов при этом могут возникать разрывы и изломы сигнала, и возникать ошибки обработки, вызванные явлением Гиббса, для минимизации которых применяют определенные методы (весовые окна, продление интервалов задания сигналов, и т.п.).

При дискретизации и во временной, и в частотной области, обычно говорят о дискретном преобразовании Фурье (ДПФ):

S(n) = s(k) exp(-j2 kn/N),

где N- количество отсчетов сигнала. Применяется оно для вычисления спектров мощности, оценивания передаточных функций и импульсных откликов, быстрого вычисления сверток при фильтрации, расчете корреляции, расчете преобразований Гильберта, и т.п. Расчет ДПФ по приведенной формуле требует вычисления n коэффициентов, каждый из которых зависит от k элементов исходного отрезка, так что число операций не может быть меньше nk. Существует целое семейство алгоритмов, известное, как “Быстрое Преобразование Фурье” - БПФ, сокращающее число операций для вычисления коэффициентов до n log(k). “Быстрое” не следует трактовать, как “упрощенное” или “неточное”. При точной арифметике результаты расчетов ДПФ и по алгоритмам БПФ совпадают.

Известное применение находят и варианты преобразования Фурье: косинусное для четных и синусное для нечетных сигналов, а также преобразование Хартли, где базисными функциями являются суммы синусов и косинусов, что позволяет повысить производительность вычислений и избавиться от комплексной арифметики. Вместо косинусных и синусных функций используются также меандровые функции Уолша, принимающие значения только +1 и -1. И, наконец, в последнее время в задачах спектрально-временнного анализа нестационарных сигналов, изучения нестационарностей и локальных особенностей сигналов "под микроскопом", очистки от шумов и сжатия сигналов начинают получать в качестве базисов разложения вейвлеты ("короткие волны"), локализованные как во временной, так и в частотной области.

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

Необходимое условие корректного представления нелинейных и нестационарных данных заключается в том, чтобы иметь возможность формирования адаптивного базиса, функционально зависимого от содержания самих данных. Такой подход реализуется в методе преобразования Гильберта-Хуанга, хотя на данный момент без соответствующих достаточно строгих математических обоснований /54/. Хорошие результаты применения метода для решения многих практических задач позволяют надеяться, что за разработкой строгой теории метода дело не станет.

2. КЛЮЧЕВЫЕ ОПЕРАЦИИ ЦИФРОВОЙ ОБРАБОТКИ [12, 43, 53]

Существуют многочисленные алгоритмы ЦОС как общего типа для сигналов в их классической временной форме (телекоммуникации, связь, телевидение и пр.), так и специализированные в самых различных отраслях науки и техники (геоинформатике, геологии и геофизике, медицине, биологии, военном деле, и пр.). Все эти алгоритмы, как правило - блочного типа, построенные на сколь угодно сложных комбинациях достаточно небольшого набора типовых цифровых операций, к основным из которых относятся свертка (конволюция), корреляция, фильтрация, функциональные преобразования, модуляция. Эти операции уже рассматривались в "Теории сигналов и систем". Ниже приводятся только ключевые позиции по этим операциям ("повторенье - мать ученья").

Линейная свертка - основная операция ЦОС, особенно в режиме реального времени. Для двух конечных причинных последовательностей h(n) и y(k) длиной соответственно N и K свертка определяется выражением:

s(k) = h(n) ? y(k) h(n) * y(k) =h(n) y(k-n), (1.2.1)

где: ? или * - символьные обозначения операции свертки. Как правило, в системах обработки одна из последовательностей y(k) представляет собой обрабатываемые данные (сигнал на входе системы), вторая h(n) - оператор (импульсный отклик) системы, а функция s(k) - выходной сигнал системы. В компьютерных системах с памятью для входных данных оператор h(n) может быть двусторонним от -N1 до +N2, например - симметричным h(-n) = h(n), с соответствующим изменением пределов суммирования в (1.2.1), что позволяет получать выходные данные без сдвига относительно входных. При строго корректной свертке с обработкой всех отсчетов входных данных размер выходного массива равен K+N1+N2-1, и должны задаваться начальные условия по отсчетам y(k) для значений y(0-n) до n=N2, и конечные для y(K+n) до n=N1. Пример выполнения свертки приведен на рис. 1.2.1.

преобразование сигнал цифровой фильтрация

Рис. 1.2.1. Примеры дискретной свертки.

Преобразование свертки однозначно определяет выходной сигнал для установленных значений входного сигнала при известном импульсном отклике системы. Обратная задача деконволюции - определение функции y(k) по функциям s(k) и h(n), имеет решение только при определенных условиях. Это объясняется тем, что свертка может существенно изменить частотный спектр сигнала s(k) и восстановление функции y(k) становится невозможным, если определенные частоты ее спектра в сигнале s(k) полностью утрачены.

Корреляция существует в двух формах: автокорреляции и взаимной корреляции.

Взаимно-корреляционная функция (ВКФ, cross-correlation function - CCF), и ее частный случай для центрированных сигналов функция взаимной ковариации (ФВК) - это показатель степени сходства формы и свойств двух сигналов. Для двух последовательностей x(k) и y(k) длиной К с нулевыми средними значениями оценка взаимной ковариации выполняется по формулам:

Kxy(n) = (1/(K-n+1)) x(k) y(k+n), n = 0, 1, 2, … (1.2.2)

Kxy(n) = (1/(K-n+1))x(k-n) y(k), n = 0, -1, -2, … (1.2.2')

Рис. 1.2.2. Функция взаимной ковариации двух детерминированных сигналов.

Пример определения сдвига между двумя детерминированными сигналами, представленными радиоимпульсами, по максимуму ФВК приведен на рис. 1.2.2. По максимуму ФВК может определяться и сдвиг между сигналами, достаточно различными по форме.

На рис. 1.2.3 приведен аналогичный пример ФВК двух одинаковых по форме сигналов, на один из которых наложен шумовой сигнал. Мощность шума превышает мощность сигнала. Вычисление ФВК на рисунке выполнено в двух вариантах. Вариант 1 полностью соответствует формуле (1.2.2). Но в условиях присутствия в сигналах достаточно мощных шумов вычисление ФВК обычно выполняется по варианту 2 - с постоянным нормировочным множителем. Это определяется тем, что по мере увеличения сдвига n и уменьшения количества суммируемых членов в формуле (1.2.2) за счет шумовых сигналов существенно нарастает ошибка оценки ФВК, которая к тому же увеличивается за счет нелинейного увеличения значения нормировочного множителя, особенно при малом количестве отсчетов. Сохранение множителя постоянным в какой-то мере компенсирует этот эффект.

Рис. 1.2.3. ФВК двух сигналов, один из которых сильно зашумлен.

На рис. 1.2.4 приведен пример вычисления функции взаимной ковариации двух одинаковых сигналов, скрытых в шумах. ФВК позволяет не только определить величину сдвига между сигналами, но и уверенно оценить период колебаний в исследуемых радиоимпульсах.

Рис. 1.2.4. ФВК двух зашумленных радиоимпульсов.

Относительный количественный показатель степени сходства двух сигналов x(k) и y(k) - функция взаимных корреляционных коэффициентов xy(n). Она вычисляется через центрированные значения сигналов (для вычисления взаимной ковариации нецентрированных сигналов достаточно центрировать один из них), и нормируется на произведение значений стандартов (средних квадратических вариаций) функций x(k) и y(k):

xy(n) = Kxy(n)/x y). (1.2.3)

x2 = Kxx(0) = (1/(K+1))(x(k))2, y2 = Kyy(0) = (1/(K+1))(y(k))2. (1.2.4)

Интервал изменения значений корреляционных коэффициентов при сдвигах n может изменяться от -1 (полная обратная корреляция) до 1 (полное сходство или стопроцентная корреляция). При сдвигах n, на которых наблюдаются нулевые значения rxy(n), сигналы некоррелированны. Коэффициент взаимной корреляции позволяет устанавливать наличие определенной связи между сигналами вне зависимости от физических свойств сигналов и их величины.

Заметим, что в технической литературе в терминах "корреляция" и "ковариация" в настоящее время существуют накладки. Корреляционными функциями называют как функции по нецентрированным, так и по центрированным сигналам, а также и функцию взаимных корреляционных коэффициентов.

Автокорреляционная функция (АКФ, correlation function, CF) является количественной интегральной характеристикой формы сигнала, дает информацию о структуре сигнала и его динамике во времени. Она, по существу, является частным случаем ВКФ для одного сигнала и представляет собой скалярное произведение сигнала и его копии в функциональной зависимости от переменной величины значения сдвига:

Bx(n) = (1/(K-n+1))x(k) x(k+n), n = 0, 1, 2, … (1.2.5)

АКФ имеет максимальное значение при n=0 (умножение сигнала на самого себя), является четной функцией Bxy(-n)=Bxy(n), и значения АКФ для отрицательных координат обычно не вычисляются. АКФ центрированного сигнала Kx(n) представляет собой функцию автоковариации (ФАК). ФАК, нормированная на свое значение Kx(0)=?x2 в n=0:

x(n) = Kx (n)/Kx(0) (1.2.6)

называется функцией автокорреляционных коэффициентов.

Рис. 1.2.5. Автокорреляционные функции.

В качестве примера на рис. 1.2.5 приведены два сигнала - прямоугольный импульс и радиоимпульс одинаковой длительности Т, и соответствующие данным сигналам формы их АКФ. Амплитуда колебаний радиоимпульса установлена равной амплитуды прямоугольного импульса, при этом энергии сигналов будут одинаковыми, что подтверждается равными значениями максимумов АКФ. При конечной длительности импульсов длительности АКФ также конечны, и равны удвоенным значениям длительности импульсов (при сдвиге копии конечного импульса на интервал его длительности как влево, так и вправо, произведение импульса со своей копией становится равным нулю). Частота колебаний АКФ радиоимпульса равна частоте колебаний заполнения радиоимпульса (боковые минимумы и максимумы АКФ возникают каждый раз при последовательных сдвигах копии радиоимпульса на половину периода колебаний его заполнения).

Линейная цифровая фильтрация является одной из операций ЦОС, имеющих первостепенное значение, и определяется как

s(k) =h(n) y(k-n), (1.2.7)

Рис. 1.2.6. Трансверсальный цифровой фильтр.

где: h(n), n=0, 1, 2, … , N - коэффициенты фильтра, y(k) и s(k) - вход и выход фильтра. Это по сути свертка сигнала с импульсной характеристикой фильтра.

На рис. 1.2.6 показана блок-схема фильтра, который в таком виде широко известен, как трансверсальный (z - задержка на один интервал дискретизации).

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

Дискретные преобразования позволяют описывать сигналы с дискретным временем в частотных координатах или переходить от описания во временной области к описанию в частотной. Переход от временных (пространственных) координат к частотным необходим во многих приложениях обработки данных.

Самым распространенным преобразованием является дискретное преобразование Фурье. При K отсчетов функции:

S(n) =s(k) exp(-j 2 kn/K). (1.2.8)

Напомним, что дискретизация функции по времени приводит к периодизации ее спектра, а дискретизация спектра по частоте - к периодизации функции. Для дискретных преобразований s(k?t) S(n?f), и функция, и ее спектр дискретны и периодичны, а числовые массивы их представления соответствуют заданию на главных периодах Т = K?t (от 0 до Т или от -Т/2 до Т/2), и 2fN = N?f (от -fN до fN), где K, N - количество отсчетов сигнала и его спектра соответственно. При этом:

f = 1/T = 1/(Kt), t = 1/2fN = 1/(Nf), tf = 1/N, N = 2TfN = K. (1.2.9)

Соотношения (1.2.9) являются условиями информационной равноценности динамической и частотной форм представления дискретных сигналов. Другими словами: для преобразований без потерь информации число отсчетов функции и ее спектра должны быть одинаковыми.

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

Модуляция сигналов. Системы регистрации, обработки, интерпретации, хранения и использования информационных данных становятся все более распределенными, что требует коммуникации данных по высокочастотным каналам связи. Как правило, информационные сигналы являются низкочастотными и ограниченными по ширине спектра, в отличие от широкополосных высокочастотных каналов связи, рассчитанных на передачу сигналов от множества источников одновременно с частотным разделением каналов. Перенос спектра сигналов из низкочастотной области в выделенную для их передачи область высоких частот выполняется операцией модуляции. При модуляции значения информационного (модулирующего) сигнала переносятся на определенный параметр высокочастотного (несущего) сигнала.

Самые распространенные схемы модуляции для передачи цифровой информации по широкополосным каналам - это амплитудная (amplitude shift keying - ASK), фазовая (phase shift keying - PSK) и частотная (frequensy shift keying - FSK) манипуляции. При передаче данных по цифровым сетям используется также импульсно-кодовая модуляция (pulse code modulation - PCM).

3. ОБЛАСТИ ПРИМЕНЕНИЯ ЦИФРОВОЙ ОБРАБОТКИ [43]

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

Процессоры ЦОС. Обработка данных в реальном времени обычно выполняется на специальных процессорах (чипах) ЦОС. Они, как правило, имеют:

§ Встроенные умножители или умножители-накопители, работающие параллельно.

§ Отдельные шины и области памяти для программ и данных.

§ Команды организации циклов.

§ Большие скорости обработки данных и тактовые частоты.

§ Использование конвейерных методов обработки данных.

Запись, воспроизведение, использование звука.

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

Синтезаторы речи представляют собой достаточно сложные устройства генерации голосовых звуков. Микросхемы синтезаторов вместе с процессорами обычно содержат в ПЗУ словари слов и фраз в форме кадров (25 мс речи) с внешним управлением интонацией, акцентом и диалектом, что позволяет на высоком уровне имитировать человеческую речь.

Распознавание речи активно изучается и развивается, особенно для целей речевого ввода информации в компьютеры. Как правило, в режиме обучения выполняется их настройка на речь пользователя, в процессе которой система оцифровывает и создает в памяти эталоны слов. В режиме распознавания речь также оцифровывается и сравнивается с эталонами в памяти. Системы распознавания речи внедряются и в товары бытового назначения (набор телефонных номеров, включение/выключение телевизора, и пр.).

Аудиосистемы воспроизведения компакт-дисков при плотности записи выше 106 бит на мм2 обеспечивают очень высокую плотность хранения информации. Аналоговый звуковой сигнал в стереоканалах дискретизируется с частотой 44.1 кГц и оцифровывается 16-битным кодом. При записи на диск сигналы модулируются (EFM - преобразование 8-ми разрядного кода в 14-ти разрядный для надежности), при считывании сигналы демодулируются, исправляются ошибки (по возможности) и выполняется цифро-аналоговое преобразование.

Применение ЦОС в телекоммуникациях.

Цифровая сотовая телефонная сеть - двусторонняя телефонная система с мобильными телефонами через радиоканалы и связью через базовые радиостанции. Мировым стандартом цифровой мобильной связи является система GSM. Частотный диапазон связи 890-960 МГц, частотный интервал канала 200 кГц, скорость передачи информации 270 кбит/с. В мобильной связи ЦОС используется для кодирования речи, выравнивания сигналов после многолучевого распространения, измерения силы и качества сигналов, кодирования с исправлением ошибок, модуляции и демодуляции.

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

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

Литература

1. Канасевич Э.Р. Анализ временных последовательностей в геофизике. - М.: Недра, 1985.- 300 с.

2. Айфичер Э., Джервис Б. Цифровая обработка сигналов. Практический подход. / М., "Вильямс", 2004, 992 с.

3. Машеров Е. Цифровая обработка сигналов - некоторые основные понятия. http://www.nsi.ru/~EMasherow/DSP.htm

4. Давыдов А.В. Теория сигналов и систем. http://prodav.narod.ru/signals/index.html.

5. Huang, N. E., Z. Shen, S. R. Long, M. C. Wu, H. H. Shih, Q. Zheng, N.-C. Yen, С. C. Tung, and H. H. Liu, 1998: The empirical mode decomposition and the Hilbert spectrum for nonlinear and non-stationary time series analysis. Proc. R. Soc. London, Ser. A, 454, 903-995.

Размещено на Allbest.ru