Discrete value approximation of the signals by a cubic spline - filters

Abstract


The cubic spline - filter for the signal discrete values approximation is considered. The filter coefficients are determined on the basis of a frequency criterion. The frequency characteristics of a cubic spline - filter and a series of parabolic spline - filters are presented. Errors of approximation of the discrete values of harmonic signal are determinated. The error of approximation of discrete values of the Gaussian signal by considered filters are analyzed. Areas of approximating spline - filters application for on - line signal processing are described.

Full Text

Алгоритмы аппроксимации сигнала часто используются при сжатии измерительных сигналов во временной области. Они позволяют на предварительном этапе его преобразования определять коэффициенты функций, аппроксимирующих его дискретные значения. Эти коэффициенты часто являются информационными параметрами, например коэффициенты параболической или кубической аппроксимации сигналов определяют значения первой и второй производных таких сигналов. Кроме того аппроксимация дискретных значений сигнала позволяет снизить погрешности, связанные с наложением копий частотных спектров сигнала (aliasing), возникающим при его дискретизации. При аппроксимации дискретных данных обычно используется алгоритм аппроксимации параболой дискретных данных x[n], заданных в моменты времени t[n] с использованием метода наименьших квадратов по пяти дискретным точкам [1]. Однако такого рода алгоритмы при работе в режиме реального времени (online) требуют значительных вычислительных ресурсов, и с этой точки зрения более перспективными являются алгоритмы, интерполирующие или аппроксимирующие дискретные значения сигнала простыми функциями на ограниченном отрезке его наблюдения. К таким функциям относятся кусочно-линейная интерполяция, а также сплайн-функции, обеспечивающие удовлетворительную погрешность аппроксимации дискретных значений сигнала. Однако большинство известных алгоритмов определения коэффициентов сплайн-функций также требует значительных вычислительных ресурсов и поэтому используется в основном при обработке всей реализации сигнала в режиме раздельного времени (offline). В связи с этим актуальной является задача разработки простых алгоритмов определения коэффициентов сплайн-аппроксимаций дискретных значений сигналов в режиме реального времени. Одним из методов аппроксимации сигнала является метод аппроксимации сплайн-функциями, представляющими собой гладкие кривые, «сшитые» на границах участков аппроксимации вместе со своими несколькими производными. Параболические сплайны «сшиты» на границах дискретных участков по 0-й и 1-й производным, а кубические - по 0-й, 1-й и 2-й производным. Вторая производная у параболической сплайн-функции и третья - у кубической сплайн-функции на границах участков претерпевают разрыв с конечным скачком. Алгоритмы параболической сплайн-аппроксимации были рассмотрены в [2]. Параболическая функция, аппроксимирующая сигнал Ux(t), на его n-м дискретном участке описывается выражением , (1) где a2[n], a1[n], a0[n] - постоянные коэффициенты. Рассмотрим задачу определения коэффициентов кубического сплайна. Он описывается на n-м участке аппроксимации уравнением, аналогичным (1): . (2) Кубическая сплайн-аппроксимация дискретных значений сигнала характеризуется меньшей погрешностью аппроксимации по сравнению с параболической [2], однако при ее использовании на каждом участке аппроксимации необходимо определять на один коэффициент больше. Рассмотрим метод определения коэффициентов кубической сплайн-аппроксимации. Производные сигнала Uy(t) определяются выражениями (3) (4) Из (3) и (4) следует: (5) Из (2)-(5) определяются граничные значения (при t = td): где td - интервал дискретизации сигнала Uy(t). Следовательно, коэффициенты ai[n] кубического сплайна связаны соотношениями (6) Z-преобразование этих соотношений имеет вид (7) Пусть передаточная функция фильтра, определяющая младший коэффициент a0[n] кубической функции (2), имеет вид [2] Используя Z-преобразование этой функции, получаем . Подставляя полученное соотношение в (7), определяем систему (8) Определитель (8) равен Определители коэффициентов ai равны Отсюда определяются выражения для коэффициентов ai[z]: (9) Из рассмотрения (9) можно сделать вывод, что для минимальной фазовой погрешности аппроксимации функция F[z] кубического сплайн-фильтра должна делиться без остатка на трехчлен [1]. Методика определения коэффициентов функции F[z] кубического сплайн-фильтра соответствует методике [2] определения аналогичных коэффициентов функции для параболического сплайн-фильтра. В качестве примера определим коэффициенты функции для пятиточечного кубического сплайн-фильтра с симметричной весовой характеристикой: (10) Z-преобразование этого выражения имеет вид . (11) В связи с тем, что это выражение аналогично соответствующему выражению для пятиточечного параболического сплайн-фильтра [2], из уравнений, определяющих частотную характеристику фильтра и ее производную, получаем: (12) Еще одно соотношение, связывающее коэффициенты bi функции фильтра, определяем исходя из (9) для коэффициента a1[n]: (13) Последний член в (13) определяет остаток при делении двух многочленов, и он должен быть равен нулю. Это выполняется при условии . (14) Уравнения (12) и (14) определяют систему: (15) Решение (15) определяет значения коэффициентов сплайн-фильтра: (16) Таким образом, младший коэффициент a0[n] 5-точечного кубического сплайна определяется уравнением (17) Выражения для остальных трех коэффициентов кубического сплайна легко определяются исходя из (9): (18) Выражения (17) - (18) с целью более простой реализации [3] целесообразно записать в виде (19) где Z-преобразование амплитудно-частотной характеристики (АЧХ) кубического сплайна может быть определено из (11) выражением (20) АЧХ кубического сплайна определяется из (20) при использовании подстановки : где - относительная частота входного сигнала; N - число интервалов дискретизации на периоде входного сигнала Ux . Амплитудно-частотная характеристика кубического сплайна, а также аналогичные характеристики параболических сплайнов [2] изображены на рис. 1. Рис. 1. Амплитудно-частотные характеристики цифровых фильтров, реализующих сплайн-аппроксимацию: 1 - четырехточечным параболическим сплайном, 2 - пятиточечным параболическим сплайном, 3 - шеститочечным параболическим сплайном, 4 - пятиточечным кубическим сплайном Сравнив частотные характеристики параболических и кубического сплайн-фильтров, можно сделать вывод, что кубический сплайн-фильтр при малом числе дискретных значений на периоде сигнала (N<5) обеспечивает меньшую погрешность аппроксимации этого сигнала, чем параболические сплайн-фильтры. Например, пятиточечный параболический сплайн-фильтр при 2,5 дискретных отсчетах на периоде сигнала () имеет погрешность по амплитуде около 50 %, а пятиточечный кубический сплайн-фильтр - около 30 %. Представляет интерес анализ погрешности аппроксимации кубическим сплайном дискретных значений конкретного сигнала. В качестве примера рассмотрена аппроксимация сигнала гауссовой формы, характерного для сигналов приборов для анализа состава и свойств веществ (хроматографа, спектроанализатора, термического анализатора). Аппроксимация такого сигнала кубическим сплайном представлена на рис. 2. Как можно заметить, погрешность аппроксимации сигнала гауссовой формы кубическим сплайном достаточно мала при небольшом числе дискретных значений (в данном примере - 8 значений на гауссов пик). Рис. 2. Аппроксимация дискретных значений сигнала UX [n] пятиточечным кубическим сплайном U5k (штриховая линия) График погрешности аппроксимации такого сигнала разными сплайнами представлен на рис. 3. Рис. 3. Погрешность аппроксимации гауссова сигнала сплайнами: 1 - четырехточечным параболическим сплайном, 2 - пятиточечным параболическим сплайном, 3 - шеститочечным параболическим сплайном, 4 - пятиточечным кубическим сплайном При рассмотрении этих графиков можно сделать вывод, что кубический пятиточечный сплайн имеет более чем в два раза меньшую погрешность (≈2 %) аппроксимации гауссова сигнала по сравнению с пятиточечным параболическим сплайном (5 %). Используемый подход может быть применен и при определении коэффициентов кубического сплайна с числом точек, большим пяти, например коэффициентов семи- или девятиточечного сплайна. Многоточечные сплайны целесообразно использовать для сглаживания сигнала в случаях его аппроксимации при наличии аддитивных помех. Применительно к обработке измерительных сигналов достаточным является использование пятиточечных кубических сплайн-фильтров. Необходимо заметить, что используемый в данной работе подход позволяет определить коэффициенты кубической сплайн-аппроксимации сигнала при произвольном числе его дискретных значений. Достоинством кубических сплайнов по сравнению с параболическими является более качественная аппроксимация сигнала, а также возможность оценки его 1-й, 2-й и 3-й производных.

About the authors

Petr K Lange

Samara State Technical University

244, Molodogvardeyskaya st., Samara, 443100, Russian Federation
(Dr. Sci. Techn.), Professor

Marat B Ungarov

Samara State Technical University

244, Molodogvardeyskaya st., Samara, 443100, Russian Federation
Postgraduate Student

Oleg V Matveev

Samara State Technical University

244, Molodogvardeyskaya st., Samara, 443100, Russian Federation
Postgraduate Student

References

  1. Лайонс Р. Цифровая обработка сигналов. - М.: Бином, 2011. - 656 с.
  2. Ланге П.К. Сплайн-аппроксимация дискретных значений сигналов с применением методов цифровой фильтрации // Сб. труд. Самарского гос. тех. ун-та. Сер. Физ.-матем. науки. - Самара: Сам. гос. техн. ун-т, 2003. - Вып. 19. - С. 134-138.
  3. Методы сплайн - функций / Завьялов Ю.С., Квасов Б.И., Мирошниченко В.Л. - М.: Наука, 1980. - 353 с.

Statistics

Views

Abstract - 58

PDF (Russian) - 13

Cited-By


Article Metrics

Metrics Loading ...

PlumX

Dimensions

Refbacks

  • There are currently no refbacks.

Copyright (c) 2016 Samara State Technical University

Creative Commons License
This work is licensed under a Creative Commons Attribution 4.0 International License.

This website uses cookies

You consent to our cookies if you continue to use our website.

About Cookies