Способ построения спектра n-мерных неразделимых цифровых сигналов

Изобретение относится к области обработки n-мерных цифровых сигналов и может быть использовано для анализа и синтеза фильтров, n-мерных неразделимых цифровых сигналов, для преобразования и обработки статических и динамических изображений, стеганографии, спектрального анализа многомерных дискретных полевых структур. Техническим результатом является унифицированное получение спектра и его локализация в пространстве и (или) времени. Способ основан на одномерном вейвлет преобразовании, при котором в ходе прямого преобразования производят последовательные, итеративные вычисления в соответствии с выражением: fk-1=fk-Sk, где: k - номер спектральной составляющей по масштабному признаку протяженности, fk - остаточный сигнал, а исходный сигнал fk, при k=K, где К - максимальное, исходное значение индекса k, Sk - элементы исходного сигнала выбирают в соответствии с индексом k, сумма всех полученных в ходе прямого преобразования спектральных составляющих Sk является исходным сигналом: f k = k = 1 K S k , каждая спектральная составляющая Sk является суперпозицией одномерного вейвлет преобразования векторов, образующих исходный сигнал по каждому из n-измерений пространства, в котором определен цифровой сигнал, на базе функции типа меандр, с определенным масштабным коэффициентом по протяженности и дискретно варьирующейся амплитудой, а разложение исходного n-мерного сигнала на спектральные составляющие производят по взаимно параллельным и перпендикулярным элементам исходного сигнала. 18 ил.

 

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

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

Известен способ получения спектра цифровых сигналов с помощью преобразования Фурье (Лайонс Р. Цифровая обработка сигналов: Второе издание. Пер. с англ. - М.: ООО «Бином-пресс», 2009, стр.63). При преобразовании Фурье происходит разложение сигнала, заданного на плоскости, образованной осями времени и амплитуды, на спектр, отображаемый на плоскости, заданной осями частоты и амплитуды, на основе базиса гармонических функций. К недостаткам преобразования Фурье можно отнести то, что в стандартном виде оно применимо к сигналам, расположенным в пространстве, образованном двумя ортогональными осями: амплитудой и временем, такие сигналы обычно называют одномерными.

Известен способ получения спектра в случае многомерного сигнала на основе ортогональных преобразований, в частности дискретное косинусное преобразование (ДКП) и его модификации, применяемые для обработки изображений (Стивен Смит. Цифровая обработка сигналов. Практическое руководство для инженеров и научных работников. Пер. с англ. Линовича А.Ю., Витязева С.В., Гусинского И.С. - М.: Додека - XXI, 2011, стр.461). Ортогональные преобразования являются линейным преобразованием евклидова векторного пространства, сохраняющем неизменным скалярное произведение векторов. В ортогональном и нормированном базисе ортогональному преобразованию соответствует ортогональная матрица. Ортогональное преобразование образует группу вращений данного евклидова пространства вокруг начала координат. В трехмерном пространстве ортогональное преобразование сводится к повороту на некоторый угол вокруг некоторой оси, проходящей через начало координат, если определитель соответствующей ортогональной матрицы равен +1. Если же этот определитель равен -1, то поворот дополняется зеркальным отражением относительно плоскости, проходящей через начало координат и перпендикулярной оси поворота. В двумерном пространстве ортогональное преобразование определяет поворот на некоторый угол вокруг начала координат или зеркальное отражение относительно некоторой прямой, проходящей через начало координат. Недостатком ДКП является то, что при его применении происходит разделение многомерного сигнала на компоненты, в результате чего утрачивается неразделимость многомерного сигнала, а полученное отображение сигнала в частотную область не во всех случаях соответствует физическому спектру сигнала.

Известен способ получения спектра как одномерного, так и многомерного цифрового сигнала путем вейвлет преобразования (Чобаноу М. Многомерные многоскоростные системы обработки сигналов. - М.: Техносфера, 2009, стр.210; Шарк Г.-Г. Применение вейвлетов для ЦОС. М., Техносфера, 2007. - 192 с., стр.14), который является наиболее близким по совокупности существенных признаков к заявляемому способу и принят в качестве прототипа. Вейвлет преобразование является способом, разбивающим сигнал на составляющие с различными частотами, каждая из которых исследуется с разрешением в подходящем масштабе. Вейвлет преобразование позволяет получить спектр n-мерного сигнала, как цифрового, так и аналогового, на основе заранее выбранных базисных n-мерных функций. В случае действительной или комплексной n-мерной функции для вейвлет преобразования выбирается анализирующая (базисная) функция g, осциллирующая относительно поверхности (в случае одномерного сигнала, около оси), образованной бесконечным множеством точек с нулевой амплитудой, что можно выразить выражением: + g ( u ) d u = 0 , где u может быть многомерным пространством. При этом функция g(u) является быстро затухающей при u→±∞. Такую функцию g в литературе называют анализирующим, базисным или материнским вейвлетом. Взяв за основу функцию вида g ' ( u , u ' , a ) = 1 a g ( u u ' a ) , полученную посредством сдвигов функции g(u), на величины u' и сопровождаемую сжатием (a<1) или растяжением (a>1) этой функции, получают вейвлет преобразование вида: L g f ( a , u ' ) = + f ( u ) g ' ( u , u ' , a ) d u . Для получения спектра неразделимых сигналов при вейвлет преобразовании размерность базисного вейвлета, по которому производится разложение, должна совпадать с размерностью разделяемого на спектр сигнала. Вейвлет преобразование предоставляет свободу в выборе базисного вейвлета, с одной стороны, это является преимуществом, т.к. есть возможность выбрать наиболее оптимальную базисную функцию разложения данного сигнала на спектр, с другой стороны, недостатком: формирование библиотек вейвлетов является трудоемкой процедурой, осложненной принципиальной невозможностью выбора универсального базиса даже для сходных цифровых сигналов, например изображений.

Задачей, на решение которой направлено заявляемое изобретение является унифицированное получение спектра n-мерного неразделимого цифрового сигнала и его локализации в пространстве и (или) времени.

Указанная задача решается за счет того, что в способе построения спектра n-мерных неразделимых цифровых сигналов, основанном на одномерном вейвлет преобразовании, в ходе прямого преобразования производят последовательные, итеративные вычисления в соответствии с выражением: fk-1=fk-Sk; где k - номер спектральной составляющей по масштабному признаку, значение индекса k соответствует длине (периоду) взаимно перпендикулярных и параллельных элементов спектральной составляющей, fk - остаточный сигнал (исходный сигнал fk, где k=K - максимальное (исходное) значение индекса k), Sk - элементы исходного сигнала, отобранные по масштабному признаку, сумма всех полученных в ходе прямого преобразования спектральных составляющих Sk является исходным сигналом: f K = k = 1 K S k , каждая спектральная составляющая Sk является суперпозицией одномерного вейвлет преобразования векторов, образующих исходный сигнал по каждому из n-измерений, на базе функции типа меандр, с определенным масштабным коэффициентом по протяженности и дискретно варьирующейся амплитудой, таким образом, разложение исходного n-мерного сигнала на спектральные составляющие происходит не по выбранному заранее n-мерному базису (базисной функции), а по взаимно параллельным и перпендикулярным элементам исходного сигнала.

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

f k 1 = f k S k ( 1 )

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

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

Для раскрытия способа необходимо уточнить обозначения измерений пространств, в которых применим рассматриваемый способ. Будем использовать понятия векторного поля и зададим функции координат положения: f ( x 1 ) - для одномерного; f ( x 1 , x 2 ) - для двухмерного; f ( x 1 , x 2 , x 3 ) - для трехмерного; … f ( x 1 , x 2 , x 3 , , x n ) - для n-мерного. Координаты вектора поля в указанной точке: { ν 1 } - для одномерного; { ν 1 , ν 2 } - для двухмерного; { ν 1 , ν 2 , ν 3 } - для трехмерного; … { ν 1 , ν 2 , ν 3 , ν p } - для p-мерного. В частном случае координаты положения могут соответствовать пространственно-временным координатам. При этом в практике обработки цифровых сигналов координаты положения и направлений вектора поля часто имеют различную размерность. Например, если в качестве сигнала рассматривать сигнал, отображаемый на дисплее осциллографа, то координата времени будет соответствовать координате положения (f(t) - одномерное), координате вектора поля будет соответствовать амплитуда ( { ν 1 } - одномерное). В данном случае координаты положения и координаты вектора поля имеют одинаковую размерность, но различные единицы измерения. Следующие примеры описывают случаи, когда размерности положения и вектора поля различны. Положение каждого пикселя цифрового изображения, сформированного на двумерной плоскости, характеризуется двумя координатами f ( x 1 , x 2 ) , в случае представления изображения в виде оттенков серого цвета, мы будем иметь одномерное пространство для вектора поля { ν 1 } , в случае цветного изображения мы можем разделить спектр по базису RGB, получив трехмерное пространство для вектора поля {R, G, B}.

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

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

Разделение исходного сигнала fK на спектральные составляющие происходит, начиная с самых протяженных (долгопериодических) элементарных всплесков, последовательно по всем направлениям положения. В каждой последующей итерации в соответствии с выражением (1) протяженность элементарного всплеска имеет значение на единицу меньше, чем в предыдущей. При выделении из исходного сигнала взаимно перпендикулярных элементарных всплесков учитывается принцип суперпозиции волн, т.е. сложения их амплитуд в месте пересечения. По этой причине последовательность обхода направлений положения, для выполнения условия однозначности разложения, должна быть одинаковой для каждой спектральной составляющей. Каждая спектральная составляющая Sk является совокупностью (суммой) элементарных всплесков протяженностью k, выделенных в остаточном сигнале fK.

Заявляемый способ иллюстрируется фиг.1-18.

Фиг.1 - последовательное разложение на спектр сигнала с одномерным пространством положения и одним вектором поля предлагаемым способом.

Фиг.2 - последовательное разложение на спектр сигнала с двумерным положением и одномерным вектором поля предлагаемым способом. Вектор поля направлен вверх от плоскости, образуемой осями x1, x2. В соответствии с выражением (1) исходный сигнал f10 и остаточные составляющие сигнала f4, f1, f0, причем f0 - нулевой остаток. В данном случае выделение спектральных составляющих S5, S2, S1 осуществляется по протяженности всплесков, причем отдельно выделяются всплески сначала в направлении x1, затем в направлении x2. Следует обратить внимание, что максимальная протяженность всплеска соответствует максимальному размеру (числу элементов дискретизации) разделяемой на спектр области из направлений положения. Для случая, приведенного на фиг.2, максимальный размер соответствует направлению X) с числом дискретов 10. Это означает, что максимальная протяженность всплесков возможна только вдоль оси x1, по оси x2, максимальная длина всплесков, ограничена размером шесть. Для случая, рассмотренного ранее, фиг.1, максимальная протяженность всплесков ограничена числом элементов дискретизации по оси x, т.е. двенадцатью.

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

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

Фиг.5 - сигнал (а) и диаграммы на фазовой плоскости, полученные методом преобразования Фурье в движущемся окне (с) и предлагаемым способом (b).

Фиг.6 - спектрограмма сигнала (а) фиг.5, полученная при помощи предлагаемого способа.

Фиг.7 - диаграмма плотности сигнала (а) фиг.5 на фазовой плоскости, полученная при помощи предлагаемого способа.

Фиг.8 - спектрограмма сигнала (а) фиг.5, полученная при помощи преобразования Фурье по всему ряду данных.

Фиг.9 - сигнал (а) и диаграммы на фазовой плоскости, полученные методом преобразования Фурье в движущемся окне (с) и предлагаемым способом (b).

Фиг.10 - спектрограмма сигнала (а) фиг.9, полученная при помощи предлагаемого способа.

Фиг.11 - диаграмма плотности сигнала (b) фиг.9 на фазовой плоскости, полученная при помощи предлагаемого способа.

Фиг.12 - спектрограмма сигнала (а) фиг.9, полученная при помощи преобразования Фурье по всему ряду данных.

Фиг.13 - исходное двумерное поле, заданное функцией

f ( x 1 , x 2 ) = S i n ( x 1 3 ) + S i n ( x 2 3 ) , x 1 , x 2 [ 0 , 1 2 π 5 ] .

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

Фиг.15 - исходное двумерное поле, заданное функцией

f ( x 1 , x 2 ) = S i n [ x 1 2 + x 2 2 ] , x 1 , x 2 [ 7 7 π 5 , 7 7 π 5 ] .

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

Фиг.17 - распределение плотности ненулевых элементов в спектральной составляющей, в зависимости от протяженности всплеска для всплесков вдоль оси x1. Плотность ненулевых элементов в спектральной составляющей рассчитывалась как число ненулевых элементов в поле спектральной составляющей, отнесенное к полному числу (площади) всей спектральной составляющей. Очевидно значительное повышение плотности, соответствующее выделению всплесками радиальных волн фиг.15 и фиг.16.

Фиг.18 - распределение плотности ненулевых элементов в спектральной составляющей, в зависимости от протяженности всплеска для всплесков вдоль оси x2. Плотность ненулевых элементов в спектральной составляющей рассчитывалась как число ненулевых элементов в поле спектральной составляющей, отнесенное к полному числу (площади) всей спектральной составляющей. Очевидно значительное повышение плотности, соответствующее выделению всплесками радиальных волн фиг.15 и фиг.16.

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

Порядок последовательности действий способа строго определен и детерминирован.

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

Осуществление изобретения в самом простом случае: одномерного пространства положения и одного вектора поля для дискретного сигнала фиг.1. Частным случаем такого сигнала является измерение некоторой одномерной величины (например, тока) во времени.

В соответствии с выражением (1) на диаграммах обозначены остаточные составляющие сигнала f3, f2, f1, f0 и спектральные составляющие S4, S3, S2, S1. Сигнал f12 - является исходным сигналом, а f0 - нулевым остатком. По причине отсутствия в спектре составляющих с номерами (длиной) 12…5 их диаграммы не приведены, т.к. спектральные составляющие с соответствующими номерами равны нулю, и соответственно неизменны значения остаточного сигнала. Данный пример очень показателен с точки зрения практических результатов предлагаемого способа:

- если положить, что ось x соответствует времени, то сумма S4, S3, S2, S1 в каждый отдельный момент времени позволяет получить спектр сигнала, при этом не требовалось производить вычисления спектра в сдвигаемом окне, как это требуется при Фурье преобразовании для получения отображения сигнала на фазовую плоскость;

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

Результаты данного спектрального преобразования практически соответствуют результатам вейвлет преобразования на базе единичного всплеска типа меандр. При переходе к большим мерностям принципиальное различие подходов становится очевидным. Рассмотрим случай двумерного положения ( f ( x 1 , x 2 ) ) и одномерного вектора поля { ν 1 } фиг.2.

Вектор поля направлен вверх от плоскости, образуемой осями x1, x2. В соответствии с выражением (1) исходный сигнал f10 и остаточные составляющие сигнала f4, f1, f0, причем f0 нулевой остаток. В данном случае выделение спектральных составляющих S5, S2, S1 осуществляется по протяженности всплесков, причем отдельно выделяются всплески сначала в направлении x1, затем в направлении x2.

Следует обратить внимание, что максимальная протяженность всплеска соответствует максимальному размеру (числу элементов дискретизации) разделяемой на спектр области из направлений положения. Для рассматриваемого случая фиг.2 максимальный размер соответствует направлению x1, с числом дискретов 10. Это означает, что максимальная протяженность всплесков возможна только вдоль оси x1, по оси x2 максимальная длинна всплесков ограничена размером шесть. Для случая, рассмотренного ранее, фиг.1, максимальная протяженность всплесков ограничена числом элементов дискретизации по оси x, т.е. двенадцатью.

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

Обобщенный алгоритм формирования спектральных составляющих, позиционируемый как реализация предлагаемого способа (прямое преобразование) для N-мерного пространства положения x1, x2…xN и одномерного вектора поля (алгоритм 1).

1. На поле fk выбирают начальное значение K = max ( x 1 , x 2 x N ) как максимальное значение определенности (т.е. наличия значений) сигнала в направлениях положения. В следующих пунктах алгоритма последовательно в порядке х1, х2…xN выделяют перпендикулярные и параллельные элементарные всплески.

2. Задают начальное значение координаты выделения элементарных всплесков n=1.

3. Проверяют, если протяженность поля в выбранном направлении xn меньше k, то переходим к пункту 4, иначе переходят к пункту 6.

4. Выделяют вдоль направления xn элементарные всплески протяженностью k, всплеск должен удовлетворять следующим условиям:

- границами всплеска считаются элементы, имеющие либо нулевое, либо отличное по знаку (в случае знакопеременного сигнала) значение от значений всплеска;

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

- у элементарного всплеска единичная протяженность по всем направлениям положения кроме xn.

5. Если в пункте 4 были локализованы элементарные всплески, то вычитают их из остаточного сигнала fk и добавляют к спектральной составляющей Sk в соответствии с выражением fk-1=fk-Sk.

6. Если значение n меньше N, то n=n+1, и переходят к пункту 3, иначе переходят к пункту 7.

7. Если значение k больше 1, то k=k-1, и переходят к пункту 2, иначе окончание алгоритма.

Результатом вычисления по алгоритму 1 является совокупность спектральных составляющих, образованных суммами перпендикулярных и параллельных элементарных всплесков в направлениях x1, x2…xN. Таким образом, разложение исходного сигнала f m a x ( x 1 , x 2 x N ) осуществляется по функциям «суммы элементарных всплесков», содержащихся в самом исходном сигнале. Иллюстрация к алгоритму 1 для двумерного сигнала на фиг.1.

Спектральные составляющие, полученные предлагаемым способом, сами имеют размерность исходного сигнала: положения p ( x 1 , x 2 , , x n ) и одномерного вектора поля { ν 1 } . При разложении в спектр предлагаемым способом для многомерных полевых структур с размерностью поля { ν 1 , ν 2 , ν 3 , ν n } расчет спектральных составляющих необходимо производить отдельно для каждого направления ν1, ν2, ν3, …νp, получая отдельные наборы спектральных составляющих для каждого направления поля. Полученные наборы спектральных составляющих можно рассматривать как векторные составляющие спектра поля и объединять путем сложения векторных компонент, перенося в пространство, имеющее размерность исходного поля.

Однозначность разложения определяется соблюдением условий последовательности действий. При сложении взаимно перпендикулярных всплесков в силу суперпозиции, в случае пересечения, область пересечения является суммой пересекающихся всплесков фиг.3. Поэтому при последовательном вычитании перпендикулярных всплесков без соблюдения выбранной заранее покоординатной последовательности может возникнуть неоднозначность разложения. Таким образом, необходимо строго выполнять покоординатную последовательность вычитания перпендикулярных всплесков, т.е. придерживаться заранее выбранной последовательности (алгоритм 1, пункт 2 и пункт 6).

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

Обобщенный алгоритм формирования исходного сигнала из спектральных составляющих, позиционируемый как реализация способа (обратное преобразование) для N-мерного пространства положения х1, х2…xN и одномерного вектора поля (алгоритм 2).

Исходный сигнал f m a x ( x 1 , x 2 x N ) является суммой полученных спектральных составляющих Sk по алгоритму 1, поэтому его можно определить как:

f m a x ( x 1 , x 2 x N ) = k = 1 m a x ( x 1 , x 2 x N ) S k . ( 2 )

В соответствии с выражением (2) получение исходного сигнала осуществляется суммированием спектральных составляющих Sk. Каждая спектральная составляющая является массивом данных размерности N и размера (по направлениям) исходного сигнала, суммирование осуществляется по каждому элементу спектральных составляющих Sk, имеющих одинаковые координаты. Таким образом, значение каждого элемента исходного сигнала f max ( x 1 , x 2 x N ) можно определить как:

f max ( x 1 , x 2 x N ) [ x 1 ' , x 2 ' , x N ' ] = k = 1 max ( x 1 , x 2 x N ) S k [ x 1 ' , x 2 ' , x N ' ] , ( 3 )

где: x 1 ' , x 2 ' , x N ' - соответствующие координаты элемента дискретного сигнала и спектральной составляющей.

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

{ O ( X ) = N ( N 1 ) [ k x i 2 + k x i 2 ] , N > 1 ( 4 )

где: N - размерность положения x1, x2, …xN, k x i 2 - размер сигнала в направлении xi.

Для одномерного сигнала N=1 выражение имеет вид:

O ( X ) = [ k 2 + k 2 ] . ( 5 )

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

O ( Д П Ф ) = 2 k 2 , ( 6 )

где: k - размер массива данных.

Для быстрого преобразования Фурье (БПФ):

O ( Д П Ф ) = k L o g 2 ( k ) . ( 7 )

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

Возможность реализации предлагаемого способа на кольце целых чисел позволяет позиционировать способ как достаточно легко реализуемый только аппаратными средствами. Предлагаемый способ может быть практически реализован программными средствами на базе ЭВМ или полностью аппаратными средствами на базе логических элементов, в частности ПЛИС.

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

Способ построения спектра n-мерных неразделимых цифровых сигналов, основанный на одномерном вейвлет преобразовании, отличающийся тем, что в ходе прямого преобразования производят последовательные, итеративные вычисления в соответствии с выражением: fk-1=fk-Sk, где k - номер спектральной составляющей по масштабному признаку протяженности, значение индекса k соответствует длине и (или) периоду взаимно перпендикулярных и параллельных элементов спектральной составляющей, fk - остаточный сигнал, а исходный сигнал fk, при k=K, где К - максимальное исходное значение индекса k, Sk - элементы исходного сигнала выбирают в соответствии с индексом k, сумма всех полученных в ходе прямого преобразования спектральных составляющих Sk является исходным сигналом: f k = k = 1 K S k , каждая спектральная составляющая Sk является суперпозицией одномерного вейвлет преобразования векторов, образующих исходный сигнал по каждому из n-измерений пространства, в котором определен цифровой сигнал, на базе функции типа меандр, с определенным масштабным коэффициентом по протяженности и дискретно варьирующейся амплитудой, а разложение исходного n-мерного сигнала на спектральные составляющие производят по взаимно параллельным и перпендикулярным элементам исходного сигнала.



 

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

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

Изобретение относится к области медицины и может быть использовано для автоматизированного исследования электроэнцефалограмм (ЭЭГ) человека. .

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

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

Изобретение относится к средствам сжатия графики, изображений и видеоинформации. .

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

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

Группа изобретений относится к способу и устройствам для сжатия и отображения изображения. Техническим результатом является исключение ошибок аппроксимации вычисления обратного дискретного косинусного преобразования с использованием вычислений с фиксированной запятой. Способ содержит этапы, на которых масштабируют каждый коэффициент в матрице закодированных коэффициентов с размерами 8×8 одним из множителей для формирования матрицы масштабированных коэффициентов, повторно применяют масштабированные одномерные преобразования с фиксированной точкой для преобразования матрицы масштабированных коэффициентов в матрицу преобразованных коэффициентов, сдвигают вправо преобразованные коэффициенты для формирования матрицы скорректированных коэффициентов, причем каждый скорректированный коэффициент в матрице скорректированных коэффициентов аппроксимирует соответствующее значение в матрице значений, которые можно сформировать с применением идеального двумерного обратного дискретного косинусного преобразования к матрице закодированных коэффициентов, отображают блок пикселей с размерами 8х8, причем каждый пиксель в блоке пикселей включает в себя значение компонента пикселя с учетом скорректированного коэффициента. 4 н. и 23 з.п. ф-лы, 12 ил., 3 табл.

Звуковое кодирующее устройство (10), приспособленное для кодирования фреймов квантованного звукового сигнала для получения кодированных фреймов, где фрейм включает ряд звуковых образцов временной области. Звуковое кодирующее устройство (10) включает этап анализа предиктивного кодирования (12) для определения информации о коэффициентах синтезирующего фильтра и фрейма области предсказания, основанного на фрейме звуковых образцов. Звуковое кодирующее устройство (10) далее включает преобразователь, вводящий временное совмещение имен (14), для преобразования перекрывающихся фреймов области предсказания в частотную область для получения спектров фрейма области предсказания, где преобразователь, вводящий временное совмещение имен (14), приспособлен для преобразования перекрывающихся фреймов области предсказания способом критической выборки. Кроме того, звуковое кодирующее устройство (10) включает кодирующее устройство, уменьшающее избыточность (16) для кодирования спектров фрейма области предсказания для получения кодированных фреймов, основанных на коэффициентах, и кодированных спектров фрейма области предсказания. 6 н. и 15 з.п. ф-лы, 20 ил.

Изобретение относится к области цифровой обработки сигналов и может быть использовано при обработке видео- и аудиосигналов в реальном масштабе времени. Техническим результатом является обеспечение выполнения различных подклассов дискретных полиномиальных преобразований (ДПП) и для реализации КИХ-фильтров с использованием заданной системы функций на конечном интервале длины N=2c. Устройство содержит систему из n N-разрядных регистров сдвига, Т блоков вычисления системы функций ДПП и Т блоков комбинационных сумматоров для N ƒ-разрядных двоичных чисел, где Т - количество элементов образа заданного подкласса ДПП, представленных (ƒ-1+log2N)-разрядными двоичными числами, ƒ=n+r, n - количество двоичных разрядов входа устройства, r - максимальное количество двоичных разрядов постоянных коэффициентов для заданного подкласса ДПП. 4 ил.

Изобретение относится к области вычислительной техники, в частности к устройствам, выполняющим операцию быстрого преобразования Хартли массива действительных чисел, и может быть использовано в системах и устройствах цифровой обработки сигналов. Технический результат заключается в повышении помехоустойчивости. Бортовой спецвычислитель содержит процессор, четыре АЦП, ОЗУ, двухпортовое ОЗУ, контроллер мультиплексного канала, два приемопередатчика мультиплексного канала, ПЗУ, системный контроллер, четыре приемопередатчика RS-485 и узел сброса. 6 з.п. ф-лы, 35 ил.

Изобретение относится к средствам для измерения времени прихода сигналов с двухпозиционной угловой манипуляцией на приемной позиции. Техническим результатом изобретения является повышение вычислительной эффективности и повышение точности измерения. Данный технический результат достигается за счет исключения неопределенности частоты приема, позволяющего исключить необходимость двумерного поиска аргументов максимума двумерной дискретной кросс-корреляционной функции (ДККФ) для сигнала с двухпозиционной угловой манипуляцией и обойтись, таким образом, поиском аргумента максимума одномерной ДККФ. Формируют двухуровневый модулирующий сигнал и дифференциально декодируют. Отображают каждый символ дифференциально декодированного сигнала на соответствующее количество отсчетов преобразованного в цифровую форму принятого сигнала относительно шкалы времени приемной позиции путем повтора значений дифференциально декодированного сигнала с повышением его частоты выборки до частоты выборки преобразованного в цифровую форму принятого сигнала. Формируют второй цифровой поток данных, который преобразуют с использованием БПФ в значения второго спектра. 2 н.з. и 3 з.п. ф-лы, 3 ил.

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

Изобретение относится к области цифровой обработки сигналов. Техническим результатом изобретения является создание ядра сопроцессора быстрого преобразования Фурье реального времени для автономного, параллельного с работой DSP - процессора (процессора цифровой обработки сигнала), выполнения быстрых преобразований Фурье комплексных массивов и некоторых сопутствующих операций, которое имеет следующие преимущества: ввод/вывод выполняются в реальном времени, параллельно с обработкой; входные/выходные данные для пользователя входных/выходных данных располагаются в прямом порядке; для расчетов и хранения данных в прямом порядке дополнительная память не требуется; форматы действительных/мнимых компонент входных и выходных данных: 32-разрядная плавающая точка (стандарт IEEE-754), 32-разрядная фиксированная точка (целое число, дополнительный код), 16-разрядная фиксированная точка (целое число, дополнительный код); формат вычислений: 32-разрядная плавающая точка; максимальный размер непосредственно выполняемого преобразования - 8192 К, минимальный - 16 К, предельный размер наращиваемого преобразования - 256 К. 5 з.п. ф-лы, 4 ил.

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

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

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