- •Часть 2
- •1.1. Клиническая оценка эффективности
- •1.2 Инженерные методы оценок качества приборов.
- •1.3 Информативность методик обследования
- •2 Электроды
- •Характеристики электродов
- •3 Усилители электробиосигналов (убс)
- •3.1. Построение усилителей биосигналов (убс)
- •3.2 Согласование усилителя с электродами
- •3.4. Методы подавления сетевой помехи
- •3.5. Подавление сетевой помехи вычитанием в усилителе
- •3.6. Типовое построение схем убс (Рис 2.7-2.9)
- •3.7. Шумовые характеристики убс
- •3.8. Защита убп от статических разрядов и импульса дефибриллятора.
- •3. Современное построение эк
- •4. Выбор частоты квантования ацп
- •4.1. Влияние неограниченности спектра сигнала
- •4.3. Искажение сигнала во входном фильтре. Метод парных эхо
- •4 .4. Погрешности амплитудного квантования ацп
- •4.5. Погрешности за счет выходного интерполирующего фильтра
- •4.6. Пример выбора частоты квантования
- •5. Простейшие элементы цифровой фильтрации
- •5.1 Общие свойства цифровой фильтрация
- •5.2 Примеры простых, легко реализуемых фильтров
- •5.3 Использованием пропуска и добавления отсчетов
- •5.4 Комбинаторика структур простейших фильтров.
- •5.6. Дифференцирующие фильтры
- •7.16. Распознавание образов
- •6. Морфологический, контурный анализ, измерение
- •6.1 Выделение комплексов
- •6.2 Измерение параметров
- •6.3 Обнаружители r пика
- •6.3.2. Технологии "Фильтр - ограничитель - фильтр (фоф)
- •6.4 Измерение амплитуды и длительности
- •Измерение моментов точек начала и конца
- •7. Математические методы обработки биосигналов
- •7.1. Удобными считаются аналитические функции
- •7.4. Преобразование Фурье
- •7.5 Равенство Парсеваля
- •7.6. Системы базисных функций
- •7.7. Разложение функций по ортогональным полиномам
- •7.8. Векторное представление сигнала.
- •7.9. Интегральное преобразование Фурье
- •7.10 Использование интегральных преобразований
- •7.12. Теорема отсчетов как преобразование Фурье
- •7.13. Фильтрация
- •7.14 . Интерполяция
- •7.15. Экстраполяция
- •8. Статистические методы
- •8.2 Оптимальная фильтрация сигнала в шуме
- •8.3 Оценка параметров. Потенциальная точность
- •8.3.1. Байесовский метод
- •8.3.2 Разделение сигнальной и шумовой функции
- •8.3.3 Ошибки случайного смещения
- •8.3.4 Ошибки неоднозначности
- •8.3.5 Метод функции правдоподобия. Проблема априорной недостаточности
- •8.4. Фильтрация Калмана - Бьюиси
- •8.5 Проверка гипотез
- •8.6 Статистические методы распознавания образов
- •8.8 Персептрон Розенблата
- •Приложение 1
- •Доказательство принципа взаимности
- •Приложение 2.
- •Приложение 3.
- •Приложение 4
- •Приложение 5
- •5.1 Проверка гипотез
- •5.2 Последовательный анализ
- •Приложение 6 Технический уход за электродами и их подготовка.
- •1) Удаление с электродов старого слоя хлористого серебра
- •2) Хлорирование серебряных электродов
- •Приложение 7
- •Приложение 8
- •8.1. Работа процессора в реальном времени
- •Приложение 9
- •Приложение 10
- •9.1. Требования к защите медицинских приборов от питающей сети 220в 50Гц.
- •6.4 Обеспечение безопасности
- •9.2. Безопасность при комплексировании приборов
- •9 .3 Безопасность при работе с электродефибрилятором
- •9.4. Измерение токов утечки
8.3.2 Разделение сигнальной и шумовой функции
Однозначное выделение пика апостериорной вероятности достигается при достаточно большом отношении сигнал/шум. Для оценки этого желаемого отношения наш сигнал у(t) и его преобразованную величину q(τ) мысленно разделим на две составляющие:
Функции g и h называются сигнальной и шумовой. Сигнальная функция является автокореляционной функцией эталона (для нашего случая она представлена на рис 8.1б). Ее значение в максимальной точке gмакс =2Е/N. Шумовая составляющая h(τ) является случайной функцией от τ, точнее, шумом n(t), отфильтрованным фильтром, согласованным с сигналом (Рис 8.2г). Дисперсия h(τ) так же равна значению 2Е/N. Величине 2Е/N0 целесообразно дать отдельное обозначение: R=2Е/N0. Этот параметр равен пиковому значению сигнальной функции g(τ), среднему квадрату шумовой функции h(τ) и отношению мощности сигнала к мощности шума в q(τ).
8.3.3 Ошибки случайного смещения
Мы можем сформировать оценку точечную или интервальную. С интервальной оценкой все просто, если нас удовлетворяет величина доверительной вероятности нахождения оценки. В качестве точечной оценки принимаем максимум q(τ). Этот максимум не совпадает с истинным значением τ0 т.к. существуют малые шумовые сдвиги вокруг истинного значения. По этому для точечной оценки обязательно надо указывать дисперсию оценки и вероятность аномальных ошибок.
Для
нахождения τ^ (как и для всякого определения
точки максимума) берется производная
q(τ) и приравнивается нулю.
.
Процедура завершена, τ^ найдена. Необходимо еще найти дисперсию этой оценки. Для этого сигнальную функцию g(τ) разложим в ряд Тейлора вокруг точки τ0:
g(τ-τ0)=(2/N0)(1+ (1/2)g11(τ0)х(τ-τ0)2).
После взятия производной q(τ) в целом и приравнивания нулю можно выделить ошибку оценки (τ-τ0):
Дисперсия ошибки (τ-τ0) находится как ее средний квадрат:
Как мы видим, дисперсия оценки τ пропорциональна уровню шума N0 и пропорциональна ширине пика сигнальной функции g(τ), точнее обратно пропорционально квадрату ее второй производной в точке τ0. Сопутствующие преобразования показаны на рис 8.3. Теперь остается найти вероятность аномальных ошибок (неоднозначности).
8.3.4 Ошибки неоднозначности
Если значение R недостаточно велико, например меньше 6, то возникают ложные пики ру(τ). Явно видна возможность двух родов ошибок измерения: 1 - ошибки неоднозначности (аномальные) и 2 - случайные погрешности в области точного значения параметра (нормальные). На рис 8.4 представлено изменение ру(τ) для разных значений R. Вероятность появления аномальных ошибок определяется значением отношения сигнал/шум R и числом интервалов размещения сигнальной функции g(τ) на априорном диапазоне существования τ (а также выбранным порогом обнаружения). Если интервал существования τ содержит М позиций и они равновероятны, то на каждой позиции установленный порог может быть пересечен шумовым пиком R с конкретной вероятностью Рлт. Нам необходимо определить вероятность ложного пересечения хотя бы на одном их М наблюдаемых интервалов. Вероятность не пересечения на М интервалах равна (1-Рлт)М. Вероятность пересечения хотя бы на одном их М интервалов равна: 1-(1-Рлт)М или примерно МРлт.. Уровень ложных тревог растет линейно с увеличением числа возможных позиций сигнала.
