Перечень ссылок. 
Разработка алгоритма обработки электромиограммы

• 1. Гехт Б. М. Теоретическая и клиническая электромиография. — Л.: Наука, 1990. — 229 с.
• 2. Васильева-Линецкая Л.Я., Роханский А. О., Галацан А. В., Черепащук Г. А., Степанов А. М., Шабалдас Д. А. Автоматизированная система исследований электромиографических сигналов человека // Открытые информационные и компьютерные информационные технологии. — Харьков, 1998. — Вып. 2 — с.215−220.
• 3. Иванова-Смоленская И.А., Кандель Э. И., Андреева Е. А. и др. Спектральный электромиографический анализ эссенциального тремора// Журн. невропатол. и психиатр. — 1986. —Т. 86. Вып. 7. — С. 975—980.
• 4. Бабкин Л. С., Гехт Б. М., Полуказаков С. Я., Федотов В. Л. Автоматический анализ игольчатой ЭМГ в дифференциальной диагностике нервно-мышечных заболеваний// Журн. невропатол. и психиатр.—1988.—Т. 86, Вып. II.—С. 1623—1628.
• 5. Коуэн X. Л., Брумлик Дж. Руководство по электромиографии и электродиагностике: Пер. с англ.—М.: Медицина, 1975.— 192 с.
• 6. Методичні вказівки до лабораторних робіт з дисципліни «Автоматизація обробки і аналізу біомедичної інформації» Упоряд.: Жемчужкіна Т.В., Козіна О.А. — Харків: ХНУРЕ, 2007. — 92 с.

ПРИЛОЖЕНИЯ

Приложение 1.

Спектральный анализ ЭМГ. Код программы спектрального анализа ЭМГ.

close all; %закрываем все открытые ранее графические окна.

%Сигнал 1.

%вызов диалогового окна.

fname=uigetfile ({'*.txt','Teкстовый файл'});

fid=fopen (fname,'rt');

y=fscanf (fid,'%f',[1,Inf]);

figure (1).

fd = 1600;

T=1/fd;

t= 0: T:(length (y)-1)*T;

plot (t, y).

title ('Исходный сигнал').

xlabel ('Время (t)');

Y = fft (y); %Быстрое преобразование Фурье по 1600 точкам.

PSDPyy = abs (Y).^2; %Спектральная плотность мощности сигнала (квадрат модуля).

f = -fd/2:fd/(length (PSDPyy)-1):fd/2; % шкала частоты для построения — достаточно половины точек.

f=fftshift (f);

figure (2).

plot (f (1:fix (end/2)), PSDPyy (1:fix (end/2))) % Вывод спектральной плотности мощности сигнала.

title ('Спектральная плотность мощности сигнала').

xlabel ('Частота Гц');

fres=f (2)-f (1); %вычисляем разрешающую способность по частоте.

n5=fix (5/fres+1);%преобразование частоты в номер элемента масива.

n150=fix (150/fres+1);

n300=fix (300/fres+1);

n800=fix (800/fres+1);

SVLF=trapz (f (n5:n150), PSDPyy (n5:n150)).

SVLF=trapz (f (n150:n300), PSDPyy (n150:n300)).

SVLF=trapz (f (n300:n800), PSDPyy (n300:n800)).

Приложение 2.

Рисунок В1 — График исходного сигнала 4pat.txt.

Рисунок В2 — Спектральная плотность мощности сигнала 4pat. txt

Рисунок В3 — График исходного сигнала 5hf. txt

Рисунок В4 — Спектральная плотность мощности сигнала 5hf. txt

Рисунок В5 — График исходного сигнала 5lf. txt

Рисунок В6 — Спектральная плотность мощности сигнала 5lf. txt

Рисунок В7 — График исходного сигнала 6hf. txt

Рисунок В8 — Спектральная плотность мощности сигнала 6hf. txt

Рисунок В9 — График исходного сигнала 6lf. txt

Рисунок В10 — Спектральная плотность мощности сигнала 6lf. txt