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

Формула изобретения

Способ определения физической работоспособности человека, включающий определение величины порога анаэробного обмена (ПАНО), отличающийся тем, что величину ПАНО определяют через показатели электропроводности биологически активных точек (БАТ), для чего исследуют репрезентативные точки-пособники правых ветвей меридианов легких Р(п) и тонкого кишечника IG(п) и левых ветвей меридианов почек R(л) и печени F(л), а количественное значение ПАНО определяют в единицах измерения л/мин по формуле:

ПАНО=2,64-4,56·Р(п)+4,61·IG(п)-1,15·R(л)+9,15·F(л),

где 2,64 - свободный член (регуляторная константа);

Р(п), IG(п), R(л), F(л) - показатели электропроводности, измеренные в репрезентативных точках-пособниках ветвей меридианов легких, тонкого кишечника, почек и печени соответственно;

п - правая ветвь меридиана;

л - левая ветвь меридиана;

-4,56; 4,61; -1,15 и 9,15 - весовые коэффициенты указанных ветвей меридианов.

Описание изобретения к патенту

Изобретение относится к области спортивной медицины и медицинской техники, а именно к рефлексодиагностике, и может быть использовано для оценки физической работоспособности (ФР) человека по электрическим параметрам биологически активных точек (БАТ).

В настоящее время в спорте для определения ФР применяют различные виды физической нагрузки, в качестве которой используют восхождение на ступеньку определенной высоты с определенной частотой (Гарвардский степ-тест), велоэргометр, с помощью которого изменяют сопротивление вращению педалей, и тредбан («бегущую дорожку») [1-4].

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

Однако при выполнении восхождения на ступеньку в работе участвуют, в основном, мышцы нижних конечностей, что нарушает стандартное требование ВОЗ к нагрузочным тестам, согласно которому в работу должно вовлекаться не менее 2/3 мышечной массы. При этом выполнение шага назад после восхождения на ступеньку является неестественным движением, достаточно сложно по координации и требует предварительного навыка. Кроме того, тестирование требует много времени из-за смены высоты ступеньки и темпа восхождения (не менее 20 мин), что неудобно при массовых обследованиях.

Следующая группа тестов связана с выполнением испытуемыми максимальных нагрузок и оценкой максимальной аэробной мощности спортсмена путем определения величины порога анаэробного обмена (ПАНО). Под ПАНО понимают «интенсивность нагрузки, выше которой у обследуемого развивается метаболический ацидоз» [1-3].

Известен способ определения ПАНО в процессе однократной нагрузки со ступенчатым повышением мощности педалирования на велоэргометре или скорости бега на тредбане, при котором используют протокол с трехминутными (по длительности работы) ступенями педалирования или бега при числе таких ступеней в нагрузке не менее 4. При этом прирост интенсивности работы в единицах потребления кислорода составляет 15-20% от максимального потребления кислорода (МПК) на каждую ступень.

ПАНО определяют по уровню ступенчато увеличивающейся мощности педалирования (или скорости бега), при которой изменения выбранного показателя достигают критического значения. Чем выше мощность работы (потребление кислорода), тем выше ФР исследуемого. При регистрации потребления кислорода на каждой ступени нагрузки по результатам обследования строят график зависимости потребления кислорода от мощности мышечной работы, на котором находят критическую точку перегиба кривой, соответствующую уровню МПК в %. Начало метаболического ацидоза определяют по резкому изменению (излому) физиологических кривых на графиках зависимости содержания лактатов, буферных оснований в крови (рН) и легочной вентиляции от мощности мышечной работы. Поэтому для обозначения данного показателя ряд авторов использует другие термины: «анаэробный порог», «аэробно-анаэробный переход», «лактатный порог» и «вентиляционный порог» и разные способы его определения, в результате чего у одних и тех же испытуемых результаты исследований не совпадают [1-3].

Таким образом, при определении ПАНО у одного и того же человека разными методами (по перелому кривой потребления кислорода, по перелому кривой накопления лактатов или по перелому кривой легочной вентиляции) в результате получаются разные уровни ФР по отношению к МПК, что снижает объективность данного способа [2]. Кроме того, недостатком известного способа является изнуряющий характер процедуры. 6% спортсменов прекращают работу, не достигнув уровня критической мощности.

Известен также способ определения ФР посредством теста PWC-170 (от первых букв английского термина «Physical Working Capaciti» - физическая работоспособность) [1, 2], основанный на определении мощности мышечной нагрузки при повышении ЧСС до уровня 170 уд/мин. Выбор ЧСС, равной 170 уд/мин, обусловлен тем, что физиологически характеризует начало оптимальной зоны функционирования кардио-респираторной системы при нагрузке. Существует два варианта определения ФР по реакции пульса на нагрузку: посредством оценки ЧСС при выполнении испытуемым стандартной мышечной нагрузки и посредством нахождения величины мощности той нагрузки, при которой ЧСС увеличивается до некоторого стандартного уровня. В основе способа лежит принцип линейной зависимости между ЧСС и мощностью нагрузки (до уровня 170 уд/мин), что позволяет по двум значениям нагрузки умеренной интенсивности путем экстраполяции линейной зависимости определить ФР. Чем больше PWC-170, тем выше ФР. Способ реализован в приборе «Кардиолидер», который позволяет, изменяя мощность мышечной работы, повысить ЧСС до нужной величины 170 уд/мин) [2].

Недостатками данного способа являются его относительная сложность, большие затраты времени (20-30 мин) и графический способ расчета величины PWC-170, который не вполне точен. Кроме того, он требует дорогостоящей аппаратуры.

Известен также способ определения ФР посредством теста Новакки [1-4], основанный на определении времени, в течение которого испытуемый способен выполнить ФН определенной, зависящей от его массы тела, продолжительности. Процедура тестирования следующая. Исходную нагрузку, равную 1 Вт/кг, увеличивают через каждые 2 мин на 1 Вт/кг до тех пор, пока испытуемый откажется выполнять работу. В момент отказа потребление кислорода близко или равно МПК, а ПАНО определяют в процентах от МПК.

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

Наиболее близким аналогом является способ определения физической работоспособности человека, включающий определение порога анаэробного обмена (ПАНО) во время выполнения ступенчато-возрастающей нагрузки «до отказа» по заявке на изобретение № 930536648 [5], который принят автором в качестве прототипа. Сущность прототипа заключается в том, что у спортсменов и лиц, не занимающихся спортом, в возрасте 17-30 лет во время ступенчато-возрастающей нагрузки на велоэргометре (тредмиле, гребном тренажере или степ-тесте) при длительности ступени 1-2 мин начальной мощности 50-100 Вт и приросте мощности на 25-50 Вт непрерывно через 10, 20, 30 или 40 сек синхронно регистрируют объем легочной вентиляции VE и дополнительно объем выделяемого углекислого газа VCO₂ , объем потребляемого кислорода VO₂ и эксцесс СО ₂ ЕхсСО₂ и затем по величине максимального потребления кислорода выделяют зоны наиболее вероятного проявления ПАНО 30-70% от МПК О₂ для лиц, не занимающихся спортом, и 40-90% для спортсменов и максимумы перегибов VO₂, VE, VCO ₂ и ЕхсСО₂, при этом если 4 максимальных перегиба укладываются в ±30 сек, то их среднее время принимается за момент проявления ПАНО, если один из максимумов не укладывается в заданный интервал, то его значения не учитываются и усредняется время трех максимумов, а если в ±30 сек укладываются только два максимума, то за основу для определения момента проявления ПАНО берется максимум, определяемый по VCO₂.

Общим существенным признаком прототипа и заявляемого объекта является определение порога анаэробного обмена (ПАНО) - основного показателя физической работоспособности человека.

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

Целью предлагаемого изобретения является расширение функциональных возможностей известного способа, повышение его объективности и информативности посредством количественного определения значения ПАНО через показатели электропроводности репрезентативных биологически активных точек (БAT).

Поставленная цель достигается тем, что ПАНО определяют через показатели электропроводности БАТ, для чего исследуют репрезентативные точки-пособники правых ветвей меридианов легких (Рп) и тонкого кишечника (IGп) и левых ветвей меридианов почек (Rл) и печени (Fл), а количественное значение ПАНО определяют в единицах измерения (л/мин) по формуле (1):

где 2,64 - свободный член (регуляторная константа);

Р(п), IG(п), R(л), F(л) - показатели электропроводности, измеренные в репрезентативных точках-пособниках ветвей меридианов легких, тонкого кишечника, почек и печени, соответственно;

п - правая ветвь меридиана;

л - левая ветвь меридиана;

-4,56; 4,61; -1,15 и 9,15 - весовые коэффициенты указанных ветвей меридианов.

Для достижения поставленной цели были проведены экспериментальные исследования, в которых приняло участие 42 юных спортсмена обоего пола в возрасте от 15 до 20 лет Училища олимпийского резерва № 1, специализирующихся в циклических видах спорта: плавании и гребле.

Сущность исследования заключалась в одновременном обследовании испытуемых двумя независимыми способами: нагрузочное тестирование на велоэргометре [1-4] и электропунктурная диагностика (ЭПД) по Накатани [6, 7]. ЭПД по методу Накатани, основанная на определении электропроводности репрезентативных БАТ 12 парных меридианов, соответствующих текущему состоянию внутренних органов и систем человека, широко применяется в спортивной и восстановительной медицине.

Для измерения электропроводности БАТ использован аппаратно-программный комплекс «АКРО», разработанный в Физико-техническом институте РАН им. А.Ф. Иоффе, Санкт-Петербург. В качестве репрезентативных БАТ, в проекциях которых производят электропунктурные измерения, используют точки-пособники. В качестве исследуемого электрического показателя, определяющего уровень активности меридианов, выбрано электрокожное сопротивление (ЭКС) БАТ.

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

Массив данных для расчета величины ПАНО формируют из 25 переменных, из которых 24 независимые переменные (результаты электропунктурных измерений БАТ в реальном масштабе времени), а 25-я зависимая переменная (функция 24 переменных) - ПАНО.

На основании полученных экспериментальных данных рассчитывают линейное уравнение множественной регрессии вида:

где а[0] - свободный член (регуляторная константа);

а[1]÷a[24] - расчетные значения весовых коэффициентов по каждому меридиану;

X[1]÷X[24] - измеренные значения электропроводности в 24 - х БАТ.

На основании полученных экспериментальных данных по пакету стандартных, программ «Statistica - 6.0» и специальной компьютерной программе «SSPS-10.0», проводят регрессионный анализ, рассчитывают регуляторную константу и значения весовых коэффициентов при всех меридианах, после чего выбирают те из них, которые в наибольшей степени влияют на значение ПАНО. В результате регрессионного анализа полученных экспериментальных данных нагрузочного тестирования (велоэргометр) и ЭПД получают математические модели для определения ПАНО, каждая из которых представляет собой количественную взаимосвязь ПАНО с электропроводностью БАТ значимых меридианов через регуляторную константу и весовые коэффициенты. В зависимости от числа задействованных меридианов математические модели различаются значениями коэффициента множественной регрессии R и коэффициента достоверности модели R sg. Результаты регрессионного анализа представлены в таблице 1, где Р, IG, R, F - меридианы легких, тонкого кишечника, почек и печени;

п - правосторонние меридианы; л - левосторонние меридианы.

Таблица 1
Результаты регрессионного анализа
№ модели	Константа	Меридианы	Весовой коэффициент	R	R sg
1	3,457	Рп	-5,858	0,479	0,229
2	2,887	Рп	-6,586	0,657	0,432
		IGп	3,535
3	2,997	Рп	-4,183	0,739	0,546
		IGп	4,224
		Rл	-5,510
4	2,635	Рп	-4,559	0,846	0,715
		IGп	4,606
		Rл	-1,154
		Fл	9,154

Из таблицы 1 видно, что в результате множественного регрессионного анализа данных получают 4 и более математических моделей, каждая из которых содержит регуляторную константу и меридианы, которые через весовые коэффициенты связаны в наибольшей степени с ПАНО, а также расчетные значения коэффициентов R и R sg. Высокое значение коэффициента множественной регрессии R (0,85) позволяет заключить, что электропроводность БАТ очень тесно связана с величиной ПАНО. По критерию достоверности R sg достаточно точной является математическая модель № 4, которая принята нами для определения ПАНО и более полно представлена в таблице 2.

Таблица 2
Результаты регрессионного анализа для построения математической модели порога анаэробного обмена (ПАНО)
Регуляторная константа	Меридиан	Весовой коэффициент	Стандартная ошибка	t	р
2,64	-	-	0,291	9,07	0,000
Рп	-4,56	0,02	-2,61	0,017
IGп	4,61	0,01	4,61	0,000
Rл	-1,15	0,03	-4,30	0,000
Fл	9Д5	0,03	3,36	0,003
Примечание: t - критерий Стьюдента, а р - уровень значимости.

Из таблицы 2 видно, что наибольшее влияние на ПАНО оказывают меридианы печени F, тонкого кишечника IG, легких Р и почек R (весовые коэффициенты 9,15; 4,61; -4,56 и -1,15 соответственно).

На основании данных таблицы 2 получают математическую модель (1), описывающую ПАНО через значения электропроводности БАТ: ПАНО=2,64-4,56·Pп+4,61·IGп-1,15·Rл+9,15·Fл (R sq=0,71),

где R sq коэффициент достоверности математической модели.

Весовые коэффициенты при меридианах со знаком плюс или минус указывают на количественную взаимосвязь меридианов с ПАНО.

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

Вторые сомножители в данной формуле - это текущие значения электропроводности ВАТ в кОм.

Определение ПАНО человека в заявляемом способе осуществляют следующим образом. Предварительно вводят в память компьютера весовые коэффициенты наиболее значимых меридианов Р, IG, R и F выявленной количественной взаимосвязи между ПАНО и значениями электропроводности БАТ. Посредством измерительных электродов (активного и пассивного) измеряют электропроводность 4-х БАТ значимых меридианов на кисти и стопе обследуемого. При этом пассивный электрод помещают в руку обследуемого, а активный электрод перемещают последовательно по кожным проекциям 4-х БАТ на его кисти и стопе. Условия проведения ЭПД стандартные: сидя, натощак, сняв с себя металлические предметы, без постороннего шума [5-10]. Время измерения электропроводности 4-х БАТ составляет не более 5-10 мин. Измеренные значения электропроводности согласно формуле 1 (4 измерения по одному человеку) в реальном масштабе времени вводят в компьютер, куда предварительно также вводят значения весовых коэффициентов при выбранных меридианах и регуляторную константу.

В компьютере их подвергают статистической обработке:

1) преобразуют измеренные значения электропроводности БАТ по каждому значимому меридиану в сигнал, воспринимаемый компьютером;

2) вычисляют произведения значений электропроводности БАТ на весовые коэффициенты по каждому значимому меридиану;

3) формируют сигнал количественной оценки ПАНО в л/мин посредством суммирования регуляторной константы и полученных произведений.

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

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

Вышеперечисленные блоки конструктивно оформлены в виде портативного прибора, сопряженного с компьютером 1 через порт 12.

Устройство работает следующим образом.

Тестирующее постоянное напряжение заданной величины с блока 3 тестового напряжения подают на измерительные электроды 6 и 7 блока 5 через блок А измерения, выход которого подключен к входу блока 8, обеспечивающего преобразование измеряемых величин токов через БАТ в сигнал, воспринимаемый компьютером 1. Исследуемые сигналы, содержащие информацию об активности значимых меридианов Рп, IGп, Rл и Fл обследуемого в данный момент времени, с выхода блока 8 преобразования посылают на первый вход блока 10. На второй вход блока 10 из блока 9 памяти поступают значения весовых коэффициентов при выбранных меридианах и регуляторная константа для определяемого ПАНО. В блоке 10 вычисляют значение ПАНО по формуле (1) и через блок 11 направляют в порт 12 компьютера 1 и далее на принтер 2.

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

13 - ввод программы и начальных данных (дата и время обследования);

14 - ввод в память компьютера формулы (1) для определения показателя ПАНО через значения электропроводности БАТ, весовые коэффициенты при них и регуляторную константу;

15 - ввод в память компьютера исходных данных обследуемого (ФИО, пол и возраст);

16 - измерение значений электропроводности 4-х репрезентативных БАТ обследуемого (ЭКС в кОм) и ввод их в компьютер;

17 - преобразование измеренных значений силы тока в мкА в ЭКС в кОм (при необходимости);

18 - аналого-цифровое преобразование данных;

19 - вычисление количественного значения ПАНО в л/мин по формуле (1), введенной в память компьютера;

20 - запись на жесткий диск компьютера и распечатка результатов.

Описание работы алгоритма.

После ввода в компьютер программы и начальных данных оператором 13 и формулы для вычисления показателя ПАНО оператором 14 устройство, реализующее заявленный способ, готово к работе. Далее оператором 15 регистрируют основные данные обследуемого и вводят их в компьютер. Оператором 16 производят измерение электропроводности 4-х репрезентативных БАТ значимых меридианов Рп, IGп, Rл, Fл и вводят их в компьютер. При использовании прибора электропунктуры, измеряющего силу тока в мкА (например, «Диакомс-7» и «Риодораку»), оператором 17 производят преобразование измеренных значений силы тока в ЭКС. При использовании комплекса «АКРО» необходимость в операторе 17 отпадает.

Затем оператором 18 производят аналого-цифровое преобразование, после чего все операции выполняют в цифровом коде. Оператором 19 производят вычисление произведений электропроводности 4-х БАТ на весовые коэффициенты по каждому значимому меридиану и формируют сигнал количественной оценки исследуемого показателя ПАНО посредством суммирования полученных произведений и регуляторной константы. Оператором 20 производят запись полученных результатов на жесткий диск компьютера и их распечатку.

Предлагаемый способ прошел апробацию и показал достаточно высокую сходимость результатов. В таблице 3 для сравнения приведены значения ПАНО, полученные путем нагрузочного тестирования (велоэргометр) и ЭПД для 10 человек: 5 - мужчины и 5 - женщины (из выборки 42).

Таблица 3
Сравнительная оценка значений ПАНО в л/мин по результатам нагрузочного тестирования и ЭПД (n=10)
Показатель	Способ определения	Количественные значения ПАНО, л/мин
ПАНО	Велоэргометр	3,78	3,485	3,08	2,52	1,996
ЭПД	3,5	3,3	2,8	2,3	1,7
Велоэргометр	2,3	2,23	2,15	1,7	1,47
ЭПД	2,4	2,1	2,1	1,5	1,5

Анализ табл.3 показал, что максимальное расхождение данных, полученных двумя независимыми способами, составляет порядка 5%, что свидетельствует о достаточной точности предлагаемого способа.

Для подтверждения осуществления предлагаемого способа приводим расчеты определения ПАНО по формуле (1) в л/мин на примере двух конкретных спортсменов (мужчины и женщины), специализирующихся в плавании. Подставляя в предлагаемую формулу (1) рассчитанные значения весовых коэффициентов и регуляторной константы и полученные на основании данных ЭПД значения ЭКС БАТ, рассчитывают количественное значение ПАНО.

Спортсмен Т-ев Ярослав, 15 лет, 1 разряд:

ПАНО=2,64-4,56·Рп+4,61·IGп-1,15·Rл+9,15·Fл=2,64-4,56·100+4,61·200-1,15·50+9,15·250=2,8 л/мин;

Спортсменка Ж-ва Анна, 15 лет, кандидат в мастера спорта: ПАНО=2,64-4,56·Pп+4,61·IGn-1,15·Rл+9,15·Fл=2,64-4,56·372+4,61·532-1,15·160+9,15·204=2,4 л/мин.

Таким образом, в результате проведенных исследований с использованием современных технических средств и математического аппарата достоверно (р<0,05) установлена и экспериментально проверена количественная взаимосвязь между показателем ПАНО и электропроводностью 4-х репрезентативных БАТ. На основании установленной взаимосвязи разработана математическая модель, позволяющая с высокой степенью достоверности определить количественное значение ПАНО только через электропроводность 4-х репрезентативных БАТ значимых меридианов и получить информацию о их влиянии на значение ПАНО.

Использование предлагаемого способа расширяет функциональные возможности известного способа и повышает его объективность и информативность. Кроме того, предлагаемый способ существенно сокращает время определения ПАНО до 5-10 мин вместо 20-30 мин.

Структурный раздел «Краткое описание чертежей»

Перечень чертежей:

фиг.1 - блок-схема устройства электропунктурного определения физической работоспособности человека;

фиг.2 - блок-схема алгоритма компьютерной программы.

Фиг.1 содержит следующие блоки:

1 - персональный компьютер;

2 - принтер;

3 - блок тестового напряжения;

4 - блок измерения;

5 - блок измерительных электродов, включающий активный электрод 6 и пассивный электрод 7;

8 - блок преобразования «аналог-код»;

9 - блок памяти весовых коэффициентов по каждому значимому меридиану для определения ПАНО;

10 - блок вычисления;

11 - блок гальванической связи;

12 - порт компьютера,

соединенные электрически последовательно таким образом, что выход блока 3 через блок 4 подключен к входам блоков 5 и 8, выход блока 8 соединен с входом блока 10, второй вход блока 10 соединен с выходом блока 9, а выход блока 10 подключен к входу блока 11, выходы которого соединены со вторым входом блока 5 и компьютером 1 через его порт 12, при этом компьютер 1 подключен к принтеру 2 и блоку 5 измерительных электродов.

Вышеперечисленные блоки конструктивно оформлены в виде портативного прибора, сопряженного с компьютером 1 через порт 12.

Фиг.2 включает в себя следующие операторы:

13 - ввод программы и начальных данных (дата и время обследования);

14 - ввод в память компьютера математической модели для определения показателя ПАНО через значения электропроводности БАТ, весовые коэффициенты при них и регудяторную константу;

15 - ввод в память компьютера исходных данных обследуемого (ФИО, пол и возраст);

16 - измерение и ввод значений электропроводности 4-х репрезентативных БАТ значимых меридианов обследуемого;

17 - преобразование измеренных значений силы тока в мкА в ЭКС в кОм (при необходимости);

18 - аналого-цифровое преобразование данных;

19 - вычисление количественного значения показателя ПАНО в принятых единицах измерения по введенной в память компьютера формуле (1), включающее вычисление произведений электропроводности БАТ на весовые коэффициенты по каждому значимому меридиану по показателю ПАНО и суммирование полученных произведений и регуляторной константы;

20 - запись на жесткий диск компьютера и распечатка результатов, соединенные электрически последовательно.

Литература

1. Аулик И.В. Определение физической работоспособности в клинике и спорте. - М.: Медицина, 1979.

2. Карпман В.Л., Белоцерковский З.Б., Гудков И.А. Тестирование в спортивной медицине. - М.: ФиС, 1988.

3. Зациорский В.М. Основы спортивной метрологии. - М.: ФиС, 1970.

4. Дубровский В.И. Спортивная медицина. М.: Владос, 1999.

5. Фалалеев А.Г., Дрычкин А.В., Черенина С.В. Способ определения физической работоспособности по заявке на изобретение № 93053648, А61В 5/091, А61В 5/08, 10.12.1996.

6. Портнов Ф.Г. Электропунктурная рефлексотерапия. Рига: Зинатне, 1988.

7. Бойцов И.В. Электропунктурная диагностика по "риодораку". - Витебск, 1996.