Научный журнал
Фундаментальные исследования
ISSN 1812-7339
"Перечень" ВАК
ИФ РИНЦ = 1,074

ИССЛЕДОВАНИЕ БИОМЕХАНИЧЕСКИХ СИГНАЛОВ НА ОСНОВЕ СИСТЕМЫ БАЗИСНЫХ ФУНКЦИЙ ХААРА

Бакусов Л.М. 1 Насыров Р.В. 1 Старцев Ю.В. 1 Бакусова Н.С. 1
1 ФГБОУ ВПО «Уфимский государственный авиационный технический университет»
Рассматривается метод анализа результатов исследования по методу функциональных биомеханических проб. Предлагается проводить анализ на основе сокращенной выборки данных, в качестве которой рассматривается информативная часть полученного сигнала. Анализ сигнала предлагается проводить на основе функций Хаара. В качестве основы рассматривается система базисных функций ортогонального полинома Хаара. Определение коэффициентов полиномов проводится на основе матричного представления разложения полинома по векторам приближения и последующего обращения полученной матрицы. Рассматривается применение метода на примере анализа биомеханических сигналов, полученных с использованием метода функциональных биомеханических проб. Приведено разложение сигналов для разных участков полученного сигнала, соответствующих механической и функциональной составляющей ответа. Приведено разложение сигналов для разных условий съема сигналов, соответствующих наличию и отсутствию механической нагрузки на исследуемый объект.
функциональные биомеханические пробы
биомеханические сигналы
система базисных функций Хаара
1. Васин Р.А. Метод функциональных биомеханических проб (на английском языке)/ Р.А.Васин, Л.М.Бакусов, Р.В.Насыров, В.В.Мальханов// Российский журнал «Биомеханика». – 1998. – №1–2. – С.58–63.
2. Бакусов Л.М. Применение метода функциональных биомеханических проб для исследования состояния позвоночника человека/ Л.М. Бакусов, А.Е. Лепилин, Р.В. Насыров // Мехатроника. – 2001, № 6. – С.44–48.
3. Бакусов Л.М. Методы и модели причинно-структурного анализа в исследовании самоорганизующихся систем. – М.: Машиностроение, 2005. – 229 с.
4. Бакусов Л.М. Качественный анализ скрининговых функциональных биомеханических проб на основе сокращенной выборки данных.// Л.М. Бакусов, Р.В. Насыров, Т.Б. Минасов // Фундаментальные исследования. – 2013. – №12.
5. Деза Е. Энциклопедический словарь расстояний/ Е.Деза, М.М.Деза. – М.:Наука, 2008. – 444 с.
6. Математическая энциклопедия в 5 томах. Т.5. – М.: Советская энциклопедия, 1984.–1248 стб.
7. Минасов Т.Б. Структурная самоорганизация костной ткани и её механизмы как источник диагностической информации/Т.Б.Минасов, А.Е.Стрижков, Л.М.Бакусов, Р.В.Насыров. – Уфа: ПЕЧАТНЫЙ ДОМЪ, 2010.–114 с.

Исследования по методу функциональных биомеханических проб (ФБМП) приводят к получению дискретных массивов отсчетов сигналов, соответствующих ответным реакциям биообъекта на стандартизованные (дозированные) механические воздействия [1, 2, 7]. В представленной ранее статье [4] рассматривалась компонентная обработка сигналов ФБМП на основе интерполяционных полиномов Лагранжа и Гаусса с целью выбора наилучшего из этих полиномов с точки зрения информативности коэффициентов разложения. Однако применение степенных полиномов для интерполирования приводит при большом числе отсчетов к необходимости вычислять определители, члены которых различаются на много порядков, что делает вычисления неустойчивыми. Применение тригонометрического полинома Гаусса при интерполяции функций с большим числом отсчетов требует включения в число базисных функций компонентов с высокими частотами, которые фактически не наблюдаются в реальном физическом процессе (являются артефактами). Это заставляет искать более удобные в вычислительном отношении, а также мало чувствительные к числу отсчетов методы анализа и соответствующие им множества базисных функций.

Как будет показано далее, на роль такого множества может претендовать система ортогональных функций Хаара [6]. Особенностью базиса Хаара является то, что каждая функция является причинно–обусловленной. Более того, каждая компонента одного порядка учитывает момент «включения» очередного причинно–обусловленного преобразования [3]. Кроме того, базис Хаара, как известно, является полным и ортогональным, что весьма существенно облегчает вычисления, а также фактически воспроизводит сигнал в виде ступенчатой функции независимо от числа отсчетов. Это представляется важным в свете построения общей методики компонентного анализа сигнала, содержащего данные реакций ФБМП.

Поскольку компонентный анализ данных ФБМП имеет конечной целью создание методик массового скринингового обследования, то очевидно, что способы обработки и содержательной интерпретации таких данных должны быть, согласно сложившейся практике, простыми, допускать удобное графическое представление. Все процедуры анализа признаков, объединения (интеграции) данных и их представления врачу–специалисту для принятия решений должны быть в максимальной степени автоматизированы. Ниже будет построена дискретная схема анализа на основе системы дискретных ортогональных функций Хаара.

Разложение вектора отсчетов наблюдаемого сигнала bak1.tif по дискретным базисным функциям Хаара связано с решением системы линейных алгебраических уравнений вида

bak2.tif, (1)

где H – квадратная матрица порядка nхn, составленная из столбцов, каждый из которых представляет собой набор последовательных значений базисных функций bak3.tif, bak4.tif – вектор коэффициентов разложения. Решение системы (1) имеет вид

bak5.tif. (2)

Обращение [H] можно упростить, используя свойство полноты и ортогональности базиса Хаара.

Домножая обе части (1) на [H–1], получим

bak6.tif. (3)

Учитывая ортогональность, заключаем, что [HTH] – диагональная матрица с одинаковыми диагональными элементами hii=n, i=1,…,n.

С учетом этого получаем

bak7.tif, где bak8.tif. (4)

Тогда окончательно получаем

bak9.tif. (5)

Легко оказать, что (5) эквивалентно (2). Действительно,

bak10.tif (6)

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

Рассмотрим пример. Матрицу H базисных функций Хаара легко задать, используя рекуррентное определение

bak11.tif (7)

где 0 ≤ log2n, 1 ≤ n ≤ 2l; n – количество функций в базисе. Сразу отметим, что для полноты системы функций необходимо n брать равным 2, 4, 8, …, т.е. 2k.

При n=8 [6] приводит к следующему виду матрицы H (8):

bak12.tif,

bak13.tif. (8)

Матрица D, согласно (4), примет вид (9)

bak14.tif, bak15.tif. (9)

Тогда

bak16.tif. (10)

Рассмотрение выражений (7)–(10) с точки зрения вычислительной сложности показывает, что с ростом n (числа отсчетов) эта сложность растёт линейно. Другой аспект интерполирования сигналов в базис Хаара связан при размерности матрицы равной, числу отсчётов, с тем, что интерполированный сигнал в точности совпадает с исходной системой отсчетов. Однако при этом мы всё же получаем ценные данные о компонентном составе сигнала в виде набора коэффициентов разложения по базисным функциям. Тем самым открывается возможность использовать качественные меры различия сигналов, такие, как χ2 – расстояние, расстояние Махалонобиса, спектральное расстояние [5], которые в перспективе позволяют построить качественные шкалы «тяжести» патологических изменений состояний костной ткани.

Применение системы ортогональных функций Хаара для анализа сигналов ФБМП

Обработка сигналов, полученных в результате ФБМП, проводилась по n = 16 точкам с соответствующими размерами матриц в выражениях (7), (8), (9), (10).

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

Таблица 1

Оценка различия коэффициентов разложения по базису Хаара для двух последовательных участков сигнала

bak17.tif

Различия между значениями коэффициентов разложения находятся в интервале от 50 % до 433 %, что является весьма существенным. Выраженный «пороговый» характер позволяет использовать эти различия в качестве информативного признака разделения механической и функциональной части ответа на тестовое воздействие.

Также представлял интерес анализ чувствительности коэффициентов разложения к некоторым характерным изменениям в механическом и функциональном состоянии кости. Для исследования этого факта реализовали измерение с механической нагрузкой на кость по типу «сдавливания». Результаты вычислений коэффициентов и их сравнения приведены в табл. 2.

Таблица 2

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

bak18.tif

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

Выводы

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

Рецензенты:

Васильев В.И., д.т.н., профессор, заведующий кафедрой вычислительной техники и защиты информации, ФГБОУ ВПО «Уфимский государственный авиационный технический университет» Министерства образования и науки РФ, г. Уфа;

Куликов Г.Г., д.т.н., профессор, заведующий кафедрой «Автоматизированные системы, управления» ФГБОУ ВПО «Уфимский государственный авиационный технический университет» Министерства образования РФ, г. Уфа.

Работа поступила в редакцию 24.06.2014.


Библиографическая ссылка

Бакусов Л.М., Насыров Р.В., Старцев Ю.В., Бакусова Н.С. ИССЛЕДОВАНИЕ БИОМЕХАНИЧЕСКИХ СИГНАЛОВ НА ОСНОВЕ СИСТЕМЫ БАЗИСНЫХ ФУНКЦИЙ ХААРА // Фундаментальные исследования. – 2014. – № 9-2. – С. 247-251;
URL: http://fundamental-research.ru/ru/article/view?id=34833 (дата обращения: 14.12.2019).

Предлагаем вашему вниманию журналы, издающиеся в издательстве «Академия Естествознания»
(Высокий импакт-фактор РИНЦ, тематика журналов охватывает все научные направления)

«Фундаментальные исследования» список ВАК ИФ РИНЦ = 1.074