Файл: Кадыров, Х. К. Синтез математических моделей биологических и медицинских систем.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