Файл: Кадыров, Х. К. Синтез математических моделей биологических и медицинских систем.pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 21.10.2024
Просмотров: 96
Скачиваний: 0
Далее, продифференцировав последнее тождество по |
и q ^ |
|
|||||||||
и подставив результат дифференцирования вместо Я в один из кор |
|
||||||||||
ней а, |
(і = 1 , 2 ) трехчлена |
|
|
|
|
|
|
|
|||
получим |
|
|
= Г + ^ |
+ |
|
|
|
і |
|||
|
|
|
|
|
|
|
|
|
|||
<ÂR( 7 1 П, q-P) + |
a £5 (qi \ |
q^) + |
aiP (i) (<7['\ q^) + |
Q (i) (<?iU)>Qz ’) = |
|
|
|||||
|
|
|
|
|
|
Q\ |
|
я1 |
|
|
|
aiR ((7i<1>' 921*) + 5 |
|
qil)) + |
щР (l) (<7j1>>qV) + |
Q (i> (<7i1)><?2 |
}) = 0» |
' |
|||||
|
|
|
|
|
|
?2 |
|
?2 |
|
|
|
|
|
|
|
i = |
1 , 2 . |
|
(4.99) |
|
|||
Учитывая, |
что |
|
„(2 ) _ |
|
.(!)„ |
n(l) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
а,- = — q\ а,- — p2 > |
|
|
|
|
|||
равенство (4.99) можно переписать в виде |
|
|
|
|
|||||||
ai [-Р «о (<7i*V ^) + |
*5 (<7і1)(72') — яѴ R |
+ |
[Q.m (q^q^) — |
|
|
||||||
|
ч1 |
|
|
|
|
|
|
|
|
|
|
, ( i ) |
Л |
= 0 , |
|
|
|
|
|
|
|
||
- q l " R(q[l\ |
|
|
|
|
|
|
|
||||
al [RА1) foiV») + |
# (^l’)^2 1 ,)1 |
+ |
[Q (1) (9 l1>9 2 1>) + |
5 (l?l1>?2 1))] — 0 . |
I |
|
|||||
|
42 |
|
|
|
|
42 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(4.100) |
|
|
Если трехчлен g2 (Я) = Я2 -f- |
^ Я |
+ q^ имеет различные корни, |
|
||||||||
т. е. |
Ф а 2, то из (4.100) следует, что каждое равенство равно нулю: |
|
|||||||||
|
/ > |
( * № |
= я\1) R № . |
9іч) - 5 ( Л 1’)- |
|
|
|
Q’?!(1) to lV 0) = 92* Я (?'!)?22)).
(4.101)
(1) fa i'V ’) = — R (<?11 )? 2 ^
42
|
|
|
Q (l) (?SU?21J) = |
— s t e i V ’)- |
|
|
; |
|
|
||||||||
|
|
|
42 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Если же известно Я-е приближение |
и |
9 21)fe |
коэффициентов |
||||||||||||||
искомого многочлена, то после |
двукратного деления на трехчлен |
||||||||||||||||
находим |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
p / „ d ) f t |
J i ) f t \ |
V w |
|
J l > k \ |
р / „ D * |
„ 0 ) f t \ |
о / „ ( О * |
. ( l ) f t \ |
), |
||||||||
а |
(q1 |
><?2 |
), |
|
, q2 |
|
), |
а (<7 ! |
, q2 |
), |
о {qi |
, <?2 |
|||||
A (i)(9 i |
|
, <?2 |
), |
а |
(I) |
/ „ < • > * |
, q2 |
), |
Ц л) (171 |
, q2 |
), |
|
|
||||
P ' |
^ ( І ) / г |
„ < 1 ) f e \ |
P ' |
|
|
y ‘ > f t \ |
n ' |
/ Д ' ) * |
„ ( D f t x |
|
|
||||||
»1 |
|
|
|
|
«2 |
|
|
|
|
|
?i' |
|
|
|
(4.102) |
ь( < ? Г .
Всоответствии с методом Ньютона находим последующие при ближения q[l) к и <7 2 ’ *.
8 ' |
115 |
Далее, |
решая систему |
|
|
|
|
|
|
|
||||
р .т (?!■’ |
- |
?ä" *) А<ЙИ + |
Р 'т (?!" *, 9І" W |
■ - - Р ( , Г ’ *, |
*), |
|||||||
|
|
- •- |
• |
|
?2 |
' |
|
|
|
|
|
|
n' /-,(*) ^ |
.92 |
^ ДлИ I п ' |
|
/л(»* ЛШ*\_ |
Qfal1’ *, ^ I)ft), |
(4.103) |
||||||
4 „(i)w i |
; л <7і |
+ Ѵ |
|
(і)(^і |
,<72 |
; — |
||||||
?2 |
|
|
|
|
?2 |
|
|
|
|
|
|
|
где |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Д /> = |
|
^ |
~ |
? Г > |
|
(4.104) |
||
|
|
|
|
Aqik) = |
|
fc+ 1 - |
qP \ |
|
||||
|
|
|
|
|
|
|
||||||
получим новое приближение |
|
|
|
(l) k |
|
|
||||||
|
|
|
|
<7iU)ft+I |
= |
? i 1)fc + |
|
|
||||
|
|
|
|
Д<7 |
|
(4.105) |
||||||
|
|
|
|
,<» k+1 |
= |
Д‘>* + |
Д ^1’ ,г |
|
||||
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
9 2 |
" " + |
Д<72 ' " |
|
|
|
Для выполнения |
описанных |
выше вычислительных процедур |
В. Э. Милн [18] предложил следующий алгоритм для осуществле ния этой последовательности.
Деление F (к) на трехчлен g2 (к) производится по схеме
Po |
P i |
- < 7 І Ч
— q™
b0 bx
P 2 |
|
T |
|
|
|
- Л |
1 |
I |
— q2l) b0 |
1 |
T ^ |
b2 |
• • • |
ЬП-2 |
P n—1 |
Pn |
— (7 )° A„_2 |
|
~ Л - 3 |
- q ? b n_ 2 |
d 0 |
d , |
|
(4.106) |
затем по схеме |
|
|
|
|
bfi— 4 |
|
bn—3 |
|||
- У |
bi |
„ü> r |
* |
|
||||||
■q\l)C0 |
|
|
„(I) n |
— |
a> г |
|||||
— <72} |
— q1W . .. |
— q1^n—6 |
6„_5 |
|||||||
|
_(» r |
|
— q2 |
<^л—6 |
— dP Ся_5 |
|||||
|
|
— ?2 |
ь 0 ... |
|||||||
|
C! |
c 2 |
- |
... |
Сл—6 |
|
Cn—5 |
|||
Для определения Aq\k) |
и Аq ^ |
|
запишем |
|
|
|||||
|
|
Ад{к) = bn - l C n— 2 ■ ■bnCn—3 |
|
|||||||
|
|
|
|
|
С 2 |
|
■ |
- 1 С л — 3 |
|
|
|
|
|
|
|
° л — 2 |
|
|
|
||
|
|
Д J b ) |
_ |
П— І1Сл— 1 — bnCn—2 |
|
|||||
|
|
~ А<^2 |
- |
с 2 |
|
- с |
— С |
• |
||
|
|
|
|
|
ь п—2 |
|
° л — 1и л — 3 |
|
bn—2
- é K .
Cn—4
-где |
|
|
|
|
|
|
ьі = Рі — я{1% -\ — qVbi-2 |
(i = |
ÖTft); |
||||
ct = bt - flf'CH - |
^ l)Ct _ 2 |
|
(/ = |
0 |
, n — 2 ); |
|
Cn—I = — ^ С л —2 |
— qCn—3 , |
|
||||
я (k + l)-e приближение |
|
|
|
|
|
|
,(D ft+i |
<7i |
+ |
Agl |
, |
|
|
<71....... = |
|
(4.107) |
||||
qi~'" ' ‘ = |
9 2 ’' " + |
a <?2 ’' |
|
|||
|
___ |
|||||
Аналогично определяются |
все остальные |
q[l) и qty (i = 1, А). |
116 .
Гл а в а 5. АНАЛИЗ ИНФОРМАЦИОННЫХ ПОТОКОВ
ИСИНТЕЗ ВЕРОЯТНОСТНЫХ МОДЕЛЕЙ
СИСТЕМ МЕДИЦИНСКОГО ОБСЛУЖИВАНИЯ
1. Метод анализа и синтеза информационных потоков
Одной из сложных и перспективных проблем современной меди цинской кибернетики является разработка методов оптимальной организации процесса медицинского обслуживания населения с целью максимального удовлетворения потребности в квалифици рованной медицинской помощи. Подобные методы можно построить на базе теории вероятности.
За последние годы появились такие новые научные направления, как теория игр, теория поиска, теория массового обслуживания и др.
Идеи и методы теории массового обслуживания за последнее время приобретают широкое применение при решении многих прикладных задач. Однако необходимо отметить, что круг практи ческих задач, решаемых методами теории массового обслуживания, непрерывно расширяется. Методы массового обслуживания с боль шим успехом могут использоваться при исследовании динамики функционирования сложных систем медицинского обслуживания, для автоматизированного управления и планирования различных медицинских мероприятий.
На начальное развитие теории массового обслуживания особое влияние оказал датский ученый А. Н. Эрланг (1878—1929). Эта теория, опираясь в основном на аппарат теории вероятности, касается изучения информационных потоков в различных систе мах.
Целью теории массового обслуживания является не изучение какого-либо конкретного обслуживания, а разработка методов ре шения, типичных для данной ситуации задач. В теории массового обслуживания разработаны математические методы определения ос новных вероятностных характеристик процессов массового обслу живания, оценки качества функционирования обслуживающей си стемы. Поэтому в современных условиях, когда заметно возрастает роль здравоохранения в жизни общества, решение вопросов организации и управления органами здравоохранения на базе тео рии массового обслуживания приобретает большую практическую ценность.
В этом аспекте важной теоретической и практической проб лемой организации и управления здравоохранения является опти мальное определение соотношения между требуемым и наличным числом обслуживающих врачей.
117
; Для исследования медицинского обслуживания'населения функ циональными органами здравоохранения эту систему детерминированно представим в виде системы уравнений Эрланга. Медицин ская система состоит из п обслуживающих медицинских работников, а поток больных, нуждающихся в обслуживании (входящий поток), удовлетворяет всем требованиям простейшего потока. При этих предположениях функционирование упомянутой системы медицин ского обслуживания можно описать при помощи марковских про цессов.
С целью применения методов теории массового обслуживания для исследования функциональных органов здравоохранения смо делируем эту систему для различных случаев.
1. Система медицинского обслуживания с потерей пациентов.
2.Система медицинского обслуживания с ожиданием.
3.Системы медицинского обслуживания с ожиданием и условной потерей пациентов.
Как увидим из последующего изложения, функционирование предполагаемых трех видов систем медицинского обслуживания будет отличаться друг от друга.
1.Пусть система медицинского обслуживания относится к системам с потерей пациентов, т. е. пациент, поступивший в си
стему обслуживания и обнаруживший занятость всех медицинских работников, покидает ее, не ожидая начала обслуживания. Система обслуживается п врачами. Если в момент поступления следующего пациента будет свободный врач, то он тут же приступит к обслужи ванию. Предположим, что время обслуживания одного пациента подчинено показательному закону с параметром ѵ, т. е. вероятность того, что время обслуживания у меньше t, можно записать следую щим образом:
Р { у < г ) = F(f) = |
(5.1) |
где ---- математическое ожидание времени обслуживания |
одного |
пациента.
В систему обслуживания поступает простейший поток пациен тов с параметром X (к — математическое ожидание числа пациентов, поступающих в систему обслуживания за единицу времени).
Основным критерием функционирования такой системы обслужи вания является вероятность отказа. Кроме того, представляет определенный практический интерес и такой критерий, как среднее число врачей, занятых в обслуживании. Первый критерий харак теризует соотношение потоков обслуженных и ушедших пациентов, т. е. полноту обслуживания входящего потока пациентов, а вто рой — степень загрузки обслуживающих врачей.
Конечной целью решения этой задачи является вычисление вероятности отказа и математическое ожидание числа занятых обслуживающих врачей. При наличии изложенных выше рассуждений функционирование таких систем обслуживания можно
118
исследовать при помощи систем однородных линейных обыкновен ных дифференциальных уравнений с постоянными коэффициентами:
Щ р - = |
Vft-1 ( 0 |
- ß + v/e) pk{t) + V (Ä + 1) pfc+1 ( 0 , !• (5.2) |
йРП}?] = |
( 0 |
— ѴЛРл (0 , |
где Po — вероятность того, что система в данный момент полностью
свободна; рк (k = 1 , га) — вероятность того, что k |
врачей |
занято |
обслуживанием пациентов; рп — вероятность того, |
что |
система |
в данный момент полностью занята. |
|
|
Итак, задача отыскания неизвестной функции вероятности при водит к решению систем га + 1 линейных однородных обыкновенных дифференциальных уравнений с постоянными коэффициентами относительно неизвестных р0 (t), рх (t), ..., рп (t) функций.
Интегрируя систему (5.2), получим все неизвестные функции
P i (t) (/ = 0, га). Для определения произвольных постоянных, вхо дящих в решения, можно использовать условие
S M 0 = i . |
(5.3) |
k = l |
|
Если нас будут интересовать предельные значения этих неиз вестных pi (t) (i = 0 , га), то, воспользовавшись предельной теоремой:
для того чтобы вероятность рк (t) при t -*■ оэ (й = 0 , га) имела пре дел рк независимо от начальных состояний, необходимо и достаточ но, чтобы к этим же пределам стремились условные вероятности Ріі при любом начальном состоянии і [19], получим систему од
нородных алгебраических уравнений (га + |
1 )-го порядка относитель |
||
но неизвестных предельных вероятностей в виде |
|
||
— Ä'Po ѵРх — Д |
|
|
|
Xpk-i — {I + vk) рк -f V (k + |
1 ) рк+1 = 0 , |
(5.4) |
|
|
1 < k < га — 1 , |
||
А,рл_, — ѵпрп = |
|
||
0 . |
|
|
|
Далее, вводя замену переменных |
|
|
|
= Pk |
(k = ТTn), |
|
систему однородных алгебраических уравнений (га + 1 )-го порядка приведем к системам неоднородных линейных алгебраических урав нений га-го порядка относительно неизвестных предельных вероят
119