CC BY
178
Раздел I. Системы искусственного интеллекта
А.С. Сигов, Е.С. Анцыферов, Б.И. Голубь, С.В. Ширяев
АДАПТИВНОЕ УПРАВЛЕНИЕ КАЧЕСТВОМ ИНТЕЛЛЕКТУАЛЬНОЙ ПРОДУКЦИИ
В настоящее время в диагностике и лечении сложных заболеваний все большую актуальность приобретает направление, ориентированное на выявление причинно-следственных связей в организме. Согласно этому направлению организм рассматривается как некоторая сложно организованная система, состоящая из множества взаимосвязанных подсистем [1]. При этом каждая из подсистем несет свою часть “ответственности” за состояние организма. Такому направлению отвечает разработка новых методов и средств диагностики и лечения, в частности широкого класса автоматизированных информационно-распознающих систем (АИРС) [2-11], функционирующих в условиях значительной априорной неопределенности и обеспечивающих тем самым возможность приобретения новых знаний о связях между параметрами физических полей исследуемых объектов и их состояниями. Увеличение количества методов и средств диагностики и лечения, широкое применение динамических исследований приводит к тому, что диагностика и лечение становятся необычайно сложными технологическими процессами, требующими четкого управления. При этом требуемое качество процессов диагностики и лечения определяется наборами стратегий, которыми располагает эксперт, принимающий решения. В таком случае речь идет не о качестве диагностики и лечения в общем, статистическом смысле, а о стратегической точности диагноза, когда каждому диагнозу ставится в соответствие одно единственное решение. Вместе с тем до сих пор главным управляющим звеном процесса диагностики является специалист и именно на него ложится нелегкая процедура принятия правильных решений. Поэтому для обеспечения высокого качества результатов, рациональной организации и облегчения работы специалистов целесообразно использование «интеллектуализированных» (экспертных) систем поддержки принятия решений (ЭС), что предполагает, в свою очередь, формализацию процессов управления качеством диагностики и лечения. Наборы стратегий предполагают использование огромного числа управляющих воздействий на процессы диагностики и лечения, что требует разработки методик экспертного управления качеством данных процессов с использованием автоматизированных систем копмлексного управления, построенных на основе совместного применения АИРС и ЭС, сочетающего их основные преимущества.
В этой связи важным является поиск оптимального соотношения между со-деражательным и формальным аспектами задачи. Действительно, при построении АИРС достигается предельная степень абстрагирования от содержательного аспекта знаний, которыми способен владеть лишь высококвалифицированный специалист, с одновременным эффективным решением важной задачи приобретения знаний о связях свойств анализируемых физических полей с состояниями исследуемых объектов. При построении же ЭС, необходимых для управления качеством, следует стремиться к построению таких формальных моделей, которые, обеспечивая возможность представления знаний в ЭС, минимизировали бы неизбежную потерю смыслового содержания процессов диагностики и лечения. При этом очевидно, что для более содержательного смыслового аспекта базы знаний
требуется, как правило, более простой и логически понятный механизм формализации процессов управления.
Все множество управляющих воздействий, необходимых для постановки диагноза, обозначим через и. Каждое управляющее воздействие может принимать некоторый континуум допустимых значений. В таком случае задачу управления диагностикой можно сформулировать как выбор из заданного допустимого множества значений управляющих воздействий таких, при которых достигается максимум заданной целевой функции. Вместе с тем, специфика медицинской диагностики состоит в том, что в каждом конкретном случае для получения “наилучшей” совокупности измеряемых диагностических параметров необходимо выбирать адекватно “наилучшую” совокупность управляющих воздействий. Реально добиться такой адекватности, особенно в условиях априорной неопределенности, практически невозможно. Поэтому речь должна идти об определении закона, по которому должны выбираться управляющие воздействия для каждого конкретного вида заболевания. Определение закона связано, как правило, с проведением большого числа измерительных экспериментов. В этом плане медицинская диагностика предоставляет довольно ограниченные возможности. Поэтому приходится пользоваться прежде всего такими исходными данными как число управляющих воздействий и допустимый интервал значений каждого из них.
Тип целевой функции определим как минимум средних потерь или вероятности ошибочных решений (максимум достоверности или качества диагностики). Это соответствует байесовскому принципу оптимизации, принятому при построении АИРС, где основным критерием является апостериорная плотность вероятности правильных решений или принадлежности текущих траекторий векторов измеренных диагностических параметров к установленному образу (виду заболевания). Для определения аналитического вида закона выбора управляющих воздействий Р(и) используем принцип максимума энтропии, согласно которому при неизвестной плотности распределения случайной величины следует выбирать плотность, обеспечивающую максимизацию энтропии при учете всех известных ограничений. Итак, энтропия плотности распределения
Н = -|Р(и)1пР(и)йи ,
ограничения
| Р(и )йи = 1,
и
IЬк(ЩР(и)а(и) = ак , к = 1,2,...,К
к
и
где Ьк (и) - известные непрерывно дифференцируемые функции ограничений; ак - заданные действительные числа, так называемые константы ограничений. Используя известный в математическом программировании метод множителей Лагранжа, получаем
Н =-| Р(и)
1п Р(и )-£2А (и)
к
к
к
где а0 = 1 и Ь0 (и) = 1 для всех и.
Взяв частные производные от функции Н1 по плотности распределения Р(и), будем иметь
5Н1
дР(и)
= -1
1п Р(и)-£2А (и)
к
+ 1Ш
Приравнивая правую часть данного соотношения нулю, придем к следующему соотношению:
Р(и) = ехр
(и)-1
Полученное соотношение принадлежит к целому семейству плотностей распределения, интегрируемых в пространстве управлений. Так, если Птп < и < итах , то
Р(и) = ехр(Л, -1).
Поскольку
| ехр(20 - 1)^и = 1,
то
ехр(Л) -1) = 1
и - и .
тах Ш1П
Р(и) =
и - и .
тах т1П
ит1п < и < итах
0 , ит.п > и > итах
То есть для и , отличных от нуля в некотором конечном интервале, распределение Р(и) принимает вид равномерного. К этому же семейству принадлежит бета-распределение. Путем вариации параметров этого распределения можно получить целый ряд важных с практической точки зрения распределений.
Таким образом, эксперимент по определению оптимального закона выбора значений управляющих воздействий состоит в выполнении следующих операций: установление определенных, в общем случае произвольных параметров бета- распределения; генерация для каждого диагностируемого объекта псевдослучайной последовательности (например, методом Монте-Карло) значений управляющих воздействия с последующим определением меры принадлежности к данному виду заболевания; оценка достоверности полученных результатов по обучающей выборке измеряемых диагностических параметров; коррекция параметров бета-распределения при неудовлетворительной достоверности; использование установленных экспериментально параметров бета-распределения при диагностике данного вида заболевания.
При наличии априорной информации о диапазонах изменения значений диагностических параметров в пределах исследуемого заболевания (получение такой информации вполне реально), проведение многократных обследований одних и тех же пациентов заменяется на генерацию псевдослучайных значений этих параметров по установленному закону распределения.
Предложенная методика управления качеством диагностики проходит апробацию в ОНЦ РАМН им. Н.Н. Блохина при диагностике сложных заболеваний молочной железы путем сочетания различных методов (термографии, радиоизотопы, маммография и др.). Результаты предварительных испытаний показали, что при сохранении необходимого уровня достоверности существенно снижается нагрузка на специалиста, управляющего диагностикой.
Процесс лечения будем представлять в виде чередующейся смены конечного числа состояний заболевания. При этом переходы из одного состояния в другое при лечении сложных заболеваний носят чаще всего вероятностный характер. В процессе лечения заболевание под влиянием управляющих воздействий должно пройти через несколько состояний и достигнуть одного из них, в котором вероят-
к
и
а
1
ность его наличия окажется ниже некоторого порога значимости. Представим этот процесс в следующем виде:
1) начальное распределение состояний, характеризующееся вектором вероятности Р0 и априорной вероятностью заболевания Р^1 > Р(р :
Ро
... • РА о о ... о ... ;
2) промежуточное распределение возможных состояний, достигнутое в г-й момент времени в результате лечебных воздействий и характеризующееся вектором р :
Р
„.ООО ... о ... •
3) конечное распределение состояний, характеризующееся вектором вероятности Рк и апостериорной вероятностью заболевания Р^ < Рв(о3:
Р
... о О О ... • РАЗ ... .
Аппроксимируем процесс управления лечением марковским вероятностным процессом с конечным числом состояний. Такой процесс описывается линейным векторноматричным уравнением относительно вектора вероятностей различных состояний
Р+1=РР, (1)
где Д. - матрица вероятностей переходов, зависящая от дискретного времени . , стадии заболевания, функционального состояния организма и выбранного способа лечебного воздействия в данный момент. Матрица Д в ходе управления лечением может принимать одно из т значений: В}, В2, ..., Вт. . Путем выбора последовательности матриц из возможного их множества необходимо обеспечить максимум критерия качества
3 = (А, РК), (2)
означающего вероятность одного из состояний на к-м шаге (А - целевой вектор, указывающий на конечное состояние). Согласно (1), для последовательных моментов времени будем иметь
Р1 = В0Р; Р2 = В1Р1;...; Рк = Вк-1Р-1 .
Отсюда
Рк = Вк-1Вк-2...В1В0Р0 . (3)
Последнее соотношение определяет алгоритм управления качеством. Максимум критерия качества
Р,(Р) = МАХ (А1Вк-1...Во Ро). (4)
Чтобы найти оптимальную совокупность В0,..., Вк-1, необходимо предварительно определить все значения критерия 3 и выбрать из них максимальное. Число таких вычислений равно тк. Для сокращения объема вычислений обычно используется метод динамического программирования. Согласно этому методу в рассмотрение вводится серия промежуточных критериев
Ъ (Р) = МАХ(ё, Бк- к Б р). (5)
Тогда, согласно принципу Беллмана,
Ъ (р) = МАХ^+1(р+1) = мбах Е1+1 (Б р). (6)
В В,
Производя оптимизацию последовательно в обратном времени, получим
Рк-1 р-1) = МВАХ (Бк _1Рк_1);
Вк-1
Рк-2( Рк-2)=мвах (Бк-2 Рк-2);
Вк-2
Итак, исходная к-шаговая задача сводится к последовательности одношаговых задач. Вместе с тем, чтобы выполнить одношаговую оптимизацию на (к-1)-м шаге, требуется иметь значения вектора Рк-1 .Но эти значения, как правило, неизвестны. Поэтому необходимо провести предварительное табулирование всех возможных (предполагаемых) значений этого вектора.
Таким образом, для управления качеством лечения, согласно рассмотренной методике, необходимо иметь: таблицы предполагаемых значений векторов состояний; набор матриц вероятностей переходов из одного состояния в другое; возможное число состояний. Отметим, что получение такой информации требует тщательной, кропотливой работы экспертов, обладающих высоким уровнем знаний в данной области.
Предложенная методика апробируется в ОНЦ РАМН им. Н.Н. Блохина при управлении качеством лечения злокачественных образований молочной железы. Например, при использовании интенсивной терапии (химия и облучение) были выделены три существенных состояния: исходное, промежуточное и конечное. Для управления лечением пациентки с 3-й стадией заболевания были разработаны две вероятностные матрицы, учитывающие все возможные переходы между состояниями, т.е. адекватные следующей диаграмме:
Б, Б12 Б\3 0.2 0.7 0.1 0.1 0.5 0.4
Б = Б21 Б22 Б23 Б1 = 0.1 0.4 0.5 Б2 = 0.2 0.1 0.7
Б31 Б32 Б33 0.1 0.1 0.8 0.1 0.2 0.7
В качестве начального распределения было принято Р0 = (1, 0, 0).
Для целевого вектора было принято ё = (0,0,1).
Число шагов к =2. Число возможных вариантов равно 22 = 4 :
(ё, Б1Б1Р0), (ё, Б1Б2Р0), (ё, В2В1Р0) и (ё, В2В2Р0). Результаты вычисления значений векторов промежуточных состояний Б1Р0 = (0.2; 0.7; 0.1);
В2Р0 = (0.1; 0.5; 0.4).
Компоненты векторов, соответствующие конечному состоянию, приняли следующие значения:
СВДР0) = 0.52; (BjB2 P0) = 0.63;
(B2 BjP0) = 0.56; (B2 B2 P0) = 0.67.
Согласно полученному, оптимальное воздействие состоит в двукратном использовании матрицы B2, т.е. определенным образом дозированного облучения. В результате применение предложенной методики обеспечило повышение качества лечения.
Достижение высокой эффективности управления связывается в настоящее время с использованием экспертных систем (ЭС) как неотъемлемого составного элемента практически любой системы автоматизированного управления сложными технологическими процессами. Предложенные в данной работе формализованные методики управления качеством процессов диагностики и лечения сложных заболеваний позволяют пополнять базы знаний экспертных систем о вероятностных законах выбора управляющих воздействий и их параметрах. Результаты предварительных испытаний показали, что при сохранении необходимого уровня достоверности предложенные методики позволяют существенно снизить нагрузку на специалистов, управляющих качеством процессов диагностики и лечения. Разработанные методики создают хорошие предпосылки для создания автоматизированных систем комплексного управления качеством, построенных на основе совместного использования АИРС и ЭС и сочетающих их основные преимущества.
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. Судаков К.В. Теория функциональных систем М.: Наука, 1996.
2. Куртев Н.Д., Анцыферов С.С. Структурно-стохастический метод обработки и распознавания информации тепловых изображений. // Оптический журнал, 1997, т.64, №2, с. 35-37.
3. Kurtev N.D., Antsyferov S.S. Structural-stochastic method of processing and recording information in thermal images. // J. Opt. Technol., 1997, v. 64 (2), p.102-104.
4. Анцыферов С.С. Формирование спектра тепловых изображений объектов и распознавание их образов. // Оптический журнал, 1999, т. 66, №12, с. 44-48.
5. Antsyferov S.S. Forming the spectrum of thermal images of objects and recognizing their patterns. // J. Opt. Technol., 1999, v. 66 (12), p. 1047- 1049.
6. Анцыферов С.С., Голубь Б.И. Принципы проектирования адаптивных широкоспектральных информационно-распознающих тепловизионных систем. // Сб.: «Тепловидение» / Под ред. А.С. Сигова, М.: МИРЭА, 2000, №13, с. 3-26.
7. Евтихиев Н.Н., Анцыферов С.С., Голубь Б.И. Адаптация информационно-распознающих биомедицинских систем. // Биомедицинская радиоэлектроника, 2001, № 1, с. 5-9.
8. Колесов С.Н., Анцыферов С.С., Голубь Б.И., Ширяев С.В. Построение медицинских систем распознавания тепловизионных образов. // Биомедицинская радиоэлектроника, 2001, № 1, с. 36-38.
9. Евтихиев Н.Н., Анцыферов С.С., Голубь Б.И. Технология адаптивной обработки информации тепловых широкоспектральных полей. // Наукоемкие технологии, 2002, т.3, №4, с.45-50.
10. Анцыферов С.С. Метрология виртуальных систем. // Измерительная техника, № 5, 2003, с. 17-21.
11. Евтихиев Н.Н., Анцыферов С. С. Методология оценки состояний организма по физическим полям. // Доклады 5-й Международной конференции "Радиоэлектроника в медицине", М.: ИРЭ РАН, 2003, с.3-5.
