Способ и устройство для аппроксимации сигналов



Способ и устройство для аппроксимации сигналов
Способ и устройство для аппроксимации сигналов
Способ и устройство для аппроксимации сигналов
Способ и устройство для аппроксимации сигналов
Способ и устройство для аппроксимации сигналов
Способ и устройство для аппроксимации сигналов
Способ и устройство для аппроксимации сигналов
Способ и устройство для аппроксимации сигналов
Способ и устройство для аппроксимации сигналов
Способ и устройство для аппроксимации сигналов
Способ и устройство для аппроксимации сигналов
Способ и устройство для аппроксимации сигналов

 


Владельцы патента RU 2467385:

Открытое акционерное общество "Научно-исследовательский институт физических измерений" (RU)

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

 

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

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

Из известных способов аппроксимации наиболее близким по технической сущности является способ интерполяции экспоненциальными суммами, основанный на алгоритме частных и разностей (Рутисхаузер Г. Алгоритм частных и разностей. - М.: Изд-во иностр. литер. 1960. - 93 с. [5], стр.54-60).

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

В способе-прототипе для аппроксимации непрерывной функции F(t) используется сумма экспонент (конечная сумма Дирихле) вида

совпадающая с функцией F(t) в 2n точках tν=t0+νh, ν=0, 1,…, 2n-1, заданных с постоянным шагом дискретизации h. Задача сводится к представлению в виде суммы простейших дробей рациональной функции f(z), представляющей z-преобразование исходной функции F(t), для которой

ν=0, 1,…, 2n-1.

Тогда [5, с.54]

где , .

На фиг.1, 2, 3 и 4 представлены формы отдельных компонент аппроксимирующей функции F(tν) при различных видах αk. Так как исходные данные для аппроксимации - значения функции F(t) являются действительными числами, то a k и αk могут быть как действительными числами, так и встречаться комплексно-сопряженными парами a k, и αk, . Если αk=0, то компонента представляет собой постоянное напряжение, изображенное на фиг.1. Если αk - действительное число, то компонента представляет собой экспоненту, изображенную на фиг.2. Если действительные части комплексно-сопряженной пары αk, равны нулю, то компонента представляет собой гармоническое колебание, изображенное на фиг.3. Если действительные части комплексно-сопряженной пары αk, не равны нулю, то соответствующая им компонента сигнала представляет собой затухающее гармоническое колебание, изображенное на фиг.4. Вид a k влияет только на фазу гармонического колебания.

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

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

2. Непрерывную функцию F(t) дискретизируют с постоянным шагом дискретизации h, то есть представляют ее последовательностью , содержащей 2n дискретных цифровых значений F(tν)=sν.

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

3.1. По значениям sν вычисляют частные и разности , σ=1, 2,…, c помощью правил ромба [5, с.13]

при условии, что , ν=0, 1,…,2n, а все числа равны нулю.

3.2. Находят последовательности и , а все полагают равными нулю.

3.3.Определяют полюсы функции f(z) по найденным в п.3.2 последовательностям с помощью прогрессивной формы алгоритма частных и разностей [5, стр.37-41,43-53], считая все равными нулю,

3.4. Находят по результатам п.п.3.2 и 3.3 коэффициенты в разложении функции f(z) на простейшие дроби с помощью модификаций алгоритма частных и разностей [5, стр.28-31, стр.57-60].

4. Вычисляют по полученным значениям и искомые параметры аппроксимирующего агрегата αk=h-1ln(λk) и , k=1,…, n.

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

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

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

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

В заявляемом способе для аппроксимации непрерывного сигнала F(t) также используется сумма экспонент (конечная сумма Дирихле) вида

совпадающая с исходным сигналом F(t) в 2n точках tν=t0+νh, ν=0, 1,…, 2n-1, заданных с постоянным шагом дискретизации h. Для этого сигнал F(t) дискретизируют и оцифровывают, представляя его N отсчетами сигнала , сигнала F(t) в моменты времени t0, t1,…, tN-1, tν=t0+νh, ν=0, 1,…, N-1, причем N>>2n. Отсчеты сигнала разбивают на М одинаковых блоков , , число которых не меньше удвоенного числа экспонент, а N-1=Mm. Интегрирование P(t) по интервалу [t0+µmh, t0+(µ+1)mh], где µ=0,…, M-1, дает

µ=0,…, M-1.

То есть исходный интервал дискретизации h заменяется большим интервалом дискретизации mh.

По найденным отсчетам интегрированного сигнала находим z-преобразование этого сигнала

где , .

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

и , k=1,…, r.

Число r экспонент, используемых для аппроксимации, не превышает целой части числа М/2 и определяется в процессе выполнения алгоритма частных и разностей из условия: меньше наперед заданной погрешности.

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

1. Непрерывный сигнал F(t) дискретизируют и оцифровывают с постоянным шагом дискретизации h, то есть представляют его N отсчетами в моменты времени t0, t1,…, tN-1, tν=t0+νh, ν=0, 1,…, N-1, причем число отсчетов N много больше удвоенного числа требуемых для аппроксимации экспонент.

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

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

, , а N-1=Mm.

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

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

5.1. По значениям вычисляют частные и разности , k=1, 2,…, с помощью правил ромба

при условии, что , µ=0, 1…, M-1, а все числа равны нулю.

5.2. Находят последовательности и , где r не превышает целой части числа М/2 и определяется в процессе выполнения алгоритма частных и разностей из условия: меньше наперед заданной погрешности. Число r ограничивает количество экспонент, необходимых для аппроксимации сигнала с заданной погрешностью.

5.3. Полагают равным нулю и по найденным в п.5.2 последовательностям и определяют полюсы функции .

5.4. Находят по результатам п.п.5.2 и 5.3 коэффициенты в разложении функции на простейшие дроби.

6. Вычисляют по полученным значениям и искомые параметры аппроксимирующего агрегата и k=1,…, r.

6.1. При измерении параметры колебательных компонент сигнала - частоты fk [Гц], коэффициенты затухания βk-1], амплитуды Ak [В] и начальные фазы θk [рад] - вычисляют по полученным значениям αk и a k, соответственно:

fk=(2π)-1Im(αk), βk=Re(αk),

,

или непосредственно по значениям и :

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

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

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

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

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

Работает устройство для аппроксимации сигналов следующим образом. Аналого-цифровой преобразователь 1 осуществляет дискретизацию и оцифровку аппроксимируемого непрерывного сигнала F(t), преобразуя его в последовательность N цифровых отсчетов в моменты времени t0, t1,…, tN-1, tν=t0+νh, ν=0, 1,…, N-1. Частота дискретизации высокая и определяется возможностями аналого-цифрового преобразователя 1. В интеграторе 2 выполняется численное интегрирование дискретизированного сигнала на интервалах времени М, задаваемых блоком 5 задания интервалов интегрирования. Число интервалов М устанавливается до начала работы не меньшим, чем удвоенное число экспонент, требуемых для аппроксимации сигнала с заданной точностью, и зависит от ожидаемого характера сигнала F(t). Последовательность интегрированных значений сигнала , где µ=0,…, М-1, поступает в блок 3 вычисления последовательностей частных и разностей, где с использованием алгоритма частных и разностей вычисляют частные и разности , k=1, 2,…, по которым находят последовательности и . В процессе вычисления этих последовательностей каждое очередное значение сравнивается пороговой схемой 6 с наперед заданной погрешностью ε. Если очередное значение становится меньше ε, вычисление прекращается. Если это условие не достигается, то вычисление прекращается, когда число членов последовательности становится равным целой части M/2. В блоке 4 вычисления параметров аппроксимирующего агрегата по значениям , определяют полюсы функции , находят коэффициенты в разложении функции на простейшие дроби и вычисляют искомые параметры аппроксимирующего агрегата по следующим формулам

и , k=1,…, r.

Моделирование известного и предлагаемого способов аппроксимации в среде «MAPLE» показало, что при использовании точных значений сигнала среднее квадратичное отклонение аппроксимирующей функции от сигнала для обоих способов не превышает сотых долей процента. При использовании квантованных значений сигнала с четырьмя десятичными разрядами в качестве исходных данных для аппроксимации в известном способе проявляется неустойчивость алгоритма и в ряде случаев аппроксимации не возможна, в то время как в предлагаемом способе среднее квадратичное отклонение аппроксимирующей функции от сигнала не превышает десятых долей процента. При наложении на аппроксимируемый сигнал нормального шума с отношением сигнал/шум до 500 и использовании оцифрованных значений сигнала с четырьмя десятичными разрядами в качестве исходных данных для аппроксимации известный способ практически перестает работать, в то время как предлагаемый способ сохраняет то же значение среднего квадратичного отклонения аппроксимирующей функции от сигнала. Известный способ аппроксимации не работает, если число отсчетов сигнала превышает число точек аппроксимации.

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

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

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

2. Патент RU 2010324, G06F 15/31. Устройство для аппроксимации функций. Опубл. 30.03.1994 г.

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

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

5. Рутисхаузер Г. Алгоритм частных и разностей. - М.: Изд-во иностр. литер., 1960. - 93 с.

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

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

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



 

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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