Способ и устройство для цифрового сжатия и восстановления сигналов

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

 

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

Известны и широко применяются способы и реализующие их устройства для цифрового сжатия и восстановления сигналов, основанные на полиноминальной, сплайновой и Фурье-аппроксимации сигналов [1, 2, 3, 4], и последующей интерполяции данных. Их основным недостатком является низкая точность, особенно для знакопеременных сигналов, и значительное влияние на получаемые результаты случайных шумов, неизбежно наложенных на сигнал.

Из известных наиболее близким по технической сущности является способ сжатия и восстановления сигналов, основанный на представлении выборочных данных в виде конечной линейной комбинации экспонент (колебаний, характеризуемых амплитудами, начальными фазами, частотами и коэффициентами затухания), параметры которых рассчитываются по методу наименьших квадратов Прони (С.Л.Марпл - мл. Цифровой спектральный анализ и его приложения. Пер. с англ.- М.: Мир, 1990. - 584 с., стр.365-391 [5]). Данный способ реализован в системе для измерения параметров и спектрального анализа слабых сигналов (Ионов С.В., Мясникова М.Г., Цыпин Б.В., Шаповал И.О. Модульная система для измерения параметров и спектрального анализа слабых сигналов // Измерительная техника. - 2011. - №4 [5]).

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

В способе-прототипе для представления в сжатом виде непрерывного сигнала x(t), заданного своими дискретными значениями х[1],…,x[N], полученными с постоянным шагом дискретизации Т, используется функция дискретного времени (n), представляющая собой сумму р экспонент (колебаний):

где 1≤n≤N, N>2р;

hk=Akexp(jθk) - комплексная амплитуда, представляющая независимый от времени параметр экспоненты;

zk=exp[(αk+j2πfk)T] - комплексная экспонента, описывающая параметр экспоненты, зависимый от времени.

Так как исходные данные x(n) являются действительными числами, то экспоненты hk и zk появляются, в общем случае, комплексно-сопряженными парами равной амплитуды.

Известный способ реализует следующую последовательность действий.

1. Исходя из априори известного вида сигнала x(t) или его предварительного анализа (см. [4] гл.8) и требуемой точности восстановления, задают число р экспонент (колебаний), используемых для представления сигнала.

2. Непрерывный сигнал x(t) дискретизируют с постоянным шагом T, то есть представляют его дискретными значениями х[1],…, x[N], причем число дискретных значений N>>2p.

3. Рассчитывают параметры экспонент hk и zk в следующей последовательности:

3.1. Составляют по методу наименьших квадратов систему уравнений

из решения которой находят коэффициенты авторегрессии аk.

3.2.Подставляют значения коэффициентов авторегрессии аk в характеристическое уравнение

3.3. zp+a1zp-l+…+ap=0

и находят его корни zk, представляющие параметры экспонент, зависимые от времени.

3.4. Минимизируя сумму квадратов ошибок по каждому параметру hk, составляют и решают вторую систему уравнений, вычисляя ее корни hk, представляющие параметры экспонент, независимые от времени:

3.5.Рассчитывают параметры колебаний: частоты

,

коэффициенты затухания амплитуды амплитуды Ak=|hk| [В] и начальные фазы θk=arghk [рад].

4. Восстанавливают по полученным значениям частот fk, коэффициентов затухания амплитуды αk, амплитуд Аk и начальных фаз θk исходный сигнал.

Известный способ при значительном уровне аддитивного шума, наложенного на сигнал, часто не обеспечивает удовлетворительных результатов. Значения коэффициентов затухания и частот часто намного превышают их истинные значения. Для решения этой проблемы предлагается ([4, стр.376, 384-390]) использовать значения р, превышающие число необходимых, исходя из вида сигнала. При практическом применении способа такое решение приводит к существенному возрастанию трудоемкости и продолжительности расчетов, что препятствует реализации способа на микроконтроллерах и использованию в составе мобильных систем реального времени. Кроме того, недостатком известного способа является то, что одновременно с восстановлением исходного сигнала происходит и восстановление сопровождающих его шумов.

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

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

В заявляемом способе для представления в сжатом виде непрерывного сигнала x(t), заданного своими дискретными значениями х[1],…,x[N], полученными с постоянным шагом дискретизации Т, также используется функция дискретного времени , представляющая собой сумму р экспонент:

где 1≤n≤N, N>2р;

hk=Akехр(jθk) - комплексная амплитуда, представляющая независимый от времени параметр колебания;

zk=exp[(αk+j2πfk)T] - комплексная экспонента, описывающая параметр колебания, зависимый от времени.

Так как исходные данные х(n) являются действительными числами, то экспоненты hk и zk появляются, в общем случае, комплексно-сопряженными парами равной амплитуды.

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

1. Исходя из априори известного вида сигнала x(t) или его предварительного анализа и требуемой точности восстановления, задают число р экспонент (колебаний), используемых для представления сигнала.

2. Непрерывный сигнал частот fk, коэффициентов затухания амплитуды αk, амплитуд Аk и начальных фаз θk дискретизируют с постоянным шагом Т, то есть представляют его дискретными значениями х[1],…, x[N], причем число дискретных значений N>>2p.

3. Разбивают последовательность дискретных значений сигнала на равные блоки длительностью m, число которых М несколько больше 1р.

4. Проводят численное интегрирование сигнала в каждом блоке.

5. Используют полученные интегральные значения для расчета параметров экспонент (колебаний) hk и zk в следующей последовательности:

5.1. Составляют по методу наименьших квадратов систему уравнений

из решения которой находят коэффициенты авторегрессии .

5.2.Подставляют значения коэффициентов авторегрессии в характеристическое уравнение

и находят параметры , представляющие параметры колебаний, зависимые от времени.

5.3.Минимизируя сумму квадратов ошибок по каждому параметру , составляют и решают вторую систему уравнений, вычисляя ее корни , представляющие параметры колебаний, независимые от времени:

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

, .

5.5. Исключают из дальнейшего рассмотрения колебания с малой относительно остальных амплитудой и колебания с большим относительно остальных коэффициентом затухания амплитуды.

6. Восстанавливают по оставшимся значениям частот fk, коэффициентов затухания амплитуды αk, амплитуд Аk и начальных фаз θk исходный сигнал.

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

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

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

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

Устройство для цифрового сжатия и восстановления сигналов содержит аналого-цифровой преобразователь 1, интегратор 2, запоминающее устройство 3, блок 4 вычисления коэффициентов авторегрессии, блок 5 вычисления параметров экспонент, зависимых от времени, блок 6 вычисления параметров экспонент, независимых от времени, блок 7 расчета параметров колебаний, блок 8 восстановления сигнала и таймер 9. Аналого-цифровой преобразователь 1, интегратор 2 и запоминающее устройство 3 соединены последовательно. Управляющие входы интегратора 2 и запоминающего устройства 3 подключены к выходу таймера 9. Выход запоминающего устройства 3 соединен с входом блока 4 вычисления коэффициентов авторегрессии и первым входом блока 6 вычисления параметров экспонент, независимых от времени. Выход блока 4 вычисления коэффициентов авторегрессии соединен с входом блока 5 вычисления параметров экспонент, зависимых от времени, и входом блока 7 расчета параметров колебаний, к выходу которого подключен блок 8 восстановления сигнала.

Все элементы, входящие в состав устройства, могут быть реализованы в виде отдельных функциональных узлов, например, блоки 2-7, 9 на программируемых логических интегральных схемах (ПЛИС), блок 8 - ЦАП или все - программным способом при использовании микроконтроллеров, оснащенных аналого-цифровыми преобразователями.

Работает устройство для цифрового сжатия и восстановления сигналов следующим образом. Аналого-цифровой преобразователь 1 осуществляет дискретизацию с шагом Т и квантование непрерывного входного сигнала x(t), преобразуя его в последовательность N цифровых отсчетов х[1],…,x[N]. Частота дискретизации высокая и определяется возможностями аналого-цифрового преобразователя 1. В интеграторе 2 выполняется численное интегрирование дискретизированного сигнала на интервалах времени m, задаваемых таймером 9. Число интервалов М устанавливается до начала работы не меньшим, чем удвоенное число экспонент, требуемых для представления сигнала с заданной точностью, и зависит от ожидаемого вида сигнала x(t). Последовательность интегрированных значений сигнала накапливается в запоминающем устройстве 3 за время существования сжимаемого кадра (фрагмента) сигнала. Затем накопленные значения передаются в блок 4 вычисления коэффициентов авторегрессии , а в запоминающем устройстве 3 накапливаются последовательность интегрированных значений следующего фрагмента сигнала. В блоке 4 вычисления коэффициентов авторегрессии используют интегральные значения сигнала для нахождения корней системы уравнений, составленных по методу наименьших квадратов.

Вычисленные значения коэффициентов авторегрессии передаются в блок 5 вычисления параметров экспонент , зависимых от времени. Результаты вычислений являются корнями характеристического уравнения и характеризуют частоты fk и коэффициенты затухания амплитуд αk колебаний, представляющих сигнал. Значения параметров экспонент , зависимых от времени, совместно с интегральными значениями сигнала поступают в блок 6 вычисления параметров экспонент , независимых от времени. В блоке 6 вычисления параметров экспонент , независимых от времени, происходит определение корней системы уравнений, составленных по методу наименьших квадратов для подгонки независимых от времени параметров экспонент, представляющих проинтегрированный сигнал, при условии, что зависимые от времени параметры этих экспонент известны.

Блок 7 расчета параметров колебаний на основе полученных значений и находит параметры колебаний - частоты fk, коэффициенты затухания амплитуды αk, амплитуды Аk и начальные фазы θk. При необходимости в блоке 7 расчета параметров колебаний исключают из дальнейшего рассмотрения колебания с малой относительно остальных амплитудой и колебания с большим относительно остальных коэффициентом затухания амплитуды.

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

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

а) проинтегрированные значения сигнала после интегратора 2;

б) накопленная за время сжимаемого фрагмента сигнала последовательность проинтегрированных значений после запоминающего устройства 3;

в) значения коэффициентов авторегрессии и значения параметров экспонент , зависимых от времени;

г) значения зависимых от времени параметров экспонент и независимых от времени параметров экспонент;

д) значения параметров колебаний - частоты fk, коэффициенты затухания амплитуды αk, амплитуды Аk и начальные фазы θk.

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

Использованные источники

1. Айфичер, Эммануил С., Джервис, Барри У. Цифровая обработка сигналов. - М.: Вильямс, 2004. - 992 с.

2. Попов Б.А., Теслер Г.С. Вычисление функций на ЭВМ. Справочник. - Киев: Наукова думка, 1984. - 600 с.

3. Патент RU №94028881 от 20.06.1996.

4. Патент RU №2099720 от 20.12.1997.

5. Марпл - мл. С.Л. Цифровой спектральный анализ и его приложения. Пер. с англ. - М.: Мир, 1990. - 584 с.

6. Ионов С.В., Мясникова М.Г., Цыпин Б.В., Шаповал И.О. Модульная система для измерения параметров и спектрального анализа слабых сигналов // Измерительная техника. - 2011. - №4.

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

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

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



 

Похожие патенты:

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

Изобретение относится к устройствам и способам кодирования и декодирования, используемым в системе связи, в которой сообщение кодируется и передается и принимается и декодируется.

Изобретение относится к кодирующему и декодирующему устройствам и их способам, используемым в системе связи с масштабируемой схемой кодирования. .

Изобретение относится к устройству и способу кодирования речи, в частности к устройству и способу кодирования речи для выполнения поиска фиксированной кодовой книги.

Изобретение относится к устройству обработки аудиосигналов для эффективного кодирования и декодирования аудиосигналов разных видов и способу для этого. .

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

Изобретение относится к способу и устройству для обработки аудиосигнала, в частности аудиосигнала, принятого в виде цифрового носителя, широковещательного сигнала или т.п.

Изобретение относится к способам кодирования и декодирования массивов цифровых данных, основанных на принципах цифровой голографии, предназначенных для их хранения на электронных носителях и передачи по каналам связи.

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

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

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

Изобретение относится к вычислительной технике

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

Изобретение относится к способам кодирования и декодирования аудиосигнала без потерь. Техническим результатом является повышение эффективности кодирования и снижение требуемого для кодирования объема памяти. Раскрыты устройство и способ основанного на контексте арифметического кодирования, а также устройство и способ основанного на контексте арифметического декодирования. Устройство, основанное на контексте арифметического декодирования может определять контекст текущего кортежа из N элементов, который должен быть декодирован, определять контекст старшего бита (MSB), соответствующий MSB-символу текущего кортежа из N элементов, и определять вероятностную модель с использованием контекста кортежа из N элементов и MSB-контекста. Затем устройство основанного на контексте арифметического декодирования может выполнять декодирование для MSB на основе определенной вероятностной модели и выполнять декодирование для младшего бита (LSB) на основе битовой глубины LSB, извлеченного из процесса декодирования для управляющего кода. 28 н. и 51 з.п. ф-лы, 29 ил.

Объектом изобретения является кодирование/декодирование цифровых сигналов, использующее, в частности, перестановочные коды, сопровождающиеся вычислением комбинаторных выражений. Согласно изобретению, эти комбинаторные выражения представлены разложениями на степени простых множителей и определяются путем считывания в памяти заранее записанных представлений разложений выбранных целых чисел. Технический результат - обеспечение эффективного декодирования перестановочных кодов. 3 н. и 5 з.п. ф-лы, 9 ил., 18 табл., 7 прил.

Объектом изобретения является кодирование/декодирование цифровых сигналов, использующее, в частности, перестановочные коды, сопровождающиеся вычислением комбинаторных выражений. Согласно изобретению, эти комбинаторные выражения представлены разложениями на степени простых множителей и определяются путем считывания в памяти заранее записанных представлений разложений выбранных целых чисел. Технический результат - обеспечение эффективного декодирования перестановочных кодов. 3 н. и 10 з.п. ф-лы, 9 ил., 11 табл., 7 прил.

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

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

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

Наверх