Файл: Кадыров, Х. К. Синтез математических моделей биологических и медицинских систем.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 21.10.2024

Просмотров: 80

Скачиваний: 0

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

При анализе эмпирических распределений важное значение име­ ет вычисление показателей симметрии и эксцесса, для чего необ­ ходимо определить моменты разного порядка. Если мы будем рас­ сматривать распределения

X = Х

* 1. * 2.

 

Рі> Рг<

где х{ — конкретные значения

случайной величины, полученные

в результате эксперимента;

pt — соответствующие вероятности,

то моменты разного порядка определяются по формуле

=

2

(*; —a)kPh

 

t= 1

где а — произвольно выбранная константа.

Моменты по отношению к константе а разделяются на начальные, когда а = 0; центральные, когда а = М\ условные, когда а ф М,

а Ф 0.

Центральные моменты применяются в анализе эмпирических распределений для характеристики дисперсии, асимметрии и экс­ цесса. Центральные моменты разного порядка вычисляются по сле­ дующим формулам:

рх = а г — М — первого порядка, р2 = а 2 — (М)г — второго порядка,

р3 = а3ЗМа2-j- 2М3— третьего порядка,

р4 = а 4 — 4Мх3 + 6М2а2+ 3Мі — четвертого порядка,

где

і= І

Далее воспользуемся конкретными моментами третьего и четвер­ того порядков. В случае симметрического распределения нечетные центральные моменты равны нулю. Поэтому, если третий момент не равен нулю, то это свидетельствует об асимметрии распределе­ ния. Исходя из центрального момента третьего порядка, определим асимметрию

А =

і = і

 

А— характеризует асимметрию распределения.

Спомощью параметра А мы получаем некоторые предположения

оформе исследуемого статистического ряда. Распределение счита­ ется почти симметричным, если |Л| < 0,1, и сильно асимметричным,

если \А\ >■ 0,5. Если А > 0, то асимметрия положительна,

при А <

•< 0 — отрицательна. Далее, с помощью центрального

момента

59


четвертого порядка р,4 определяем эксцесс Е по формуле

Р а 4 — Ш а 3 + 6М2а 3 + ЗМ4

К - (М)3]4

На основании параметра Е мы получаем сведения о крутизне кривой. Если Е <. О, то кривая вогнута по сравнению с нормальной кривой, если Е > 0 — она заострена сверху. Распределение счита­ ется близким к нормальному, если |Д| < 0,1, и значительно откло­ няющимся от нормального, если |Д| > 0,5.

Далее выдвигается гипотеза о том, что закон распределеная'рассматриваемых случайных процессов близок к закону нормального распределения. В математической статистике разработано мно­ жество критериев для оценки достоверности принятых допущений о распределениях. Некоторые из них справедливы лишь для определе­ ния моделей, другие применимы для широкого класса законов рас­ пределений. Мы будем использовать для проверки любого распреде­ ления критерий хи-квадрат Пирсона, который определяется по

формуле

(2.18),

где

pt — значения вероятности теоретического

распределения;

р\ — значения вероятности эмпирического распре­

деления;

п — число

наблюдений.

Для

определения

значений теоретического распределения ста­

тистический ряд необходимо разделить на классы, а затем опреде­ лять вероятности совпадения случайных величин на участках между двумя классами. Эта вероятность находится по формуле

Pt {ßi-i < М < ß,} = Ф (ß,) - Ф (ß,-i),

где ф (ßt) — значение функции распределения, соответствующее высшей границе класса і; Ф (ß(—і) — значение функции распреде­

ления, соответствующее высшей границе

предшествующего

класса

і — 1.

численных значений функций

распределений ф

(ß,) и

Для

Ф (ßi—і)

разработаны специальные таблицы. Численные значения

эмпирического распределения определяются следующим отноше­ нием:

где с,- — число наблюдений, появляющихся в і классе.

По формуле (2.30) определяется численное значение X2, а затем по числу степеней свободы из специальной таблицы находится веро­ ятность близости эмпирической кривой распределения к нормаль­ ному закону распределения. Итак, при помощи критерия Пирсона оценивается степень согласованности эмпирического закона распре­ деления с нормальным законом распределения.

Из изложенного выше очевидно, что исследование случайных величин, описывающих некоторые случайные процессы и определе­ ние законов их распределения, требует выполнения большого коли­ чества вычислительных работ, чего нельзя осуществить без при­ влечения современных электронных вычислительных машин.

60



Г л а в а 3. МЕТОДЫ СИНТЕЗА МАТЕМАТИЧЕСКИХ МОДЕЛЕЙ СЛОЖНЫХ СИСТЕМ

1.Синтез математических моделей одномерных систем

Алгоритм определения веса влияющих факторов. Исследование мно­ гомерной системы с помощью методов математического моделирова­ ния приводит к решению сложных математических задач. В связи с этим перед исследователем естественно возникает мысль о поиске методов упрощения математической модели без ущерба для основ­ ного содержания рассматриваемого явления.

Как известно, степень сложности математической модели в ка­ кой-то мере обусловлена степенями свободы рассматриваемой много­ мерной системы. Следовательно, уменьшение степеней свободы исследуемой системы является одним из источников упрощения мате­ матической задачи определения математической модели рассматрива­ емой системы. Вместе с тем, всякое снижение степеней свободы у ис­ следуемой системы приводит к ухудшению качества математической модели, т. е. модель становится неточной и не полностью отража­ ет реальность. Однако роль каждого фактора в процессе формирова­ ния исследуемого процесса различна, т. е. каждый фактор, который является отдельным элементом большой системы, по-разному воз­ действует на выходную функцию большой системы. Следовательно, критерий оценки веса каждого фактора должен строиться так, чтобы можно было получить количественный показатель оценки веса каж­ дого фактора в многомерной системе.

В природе большинство рассматриваемых многомерных систем обусловлено случайными процессами. Поэтому выдвигаемый крите­ рий оценки веса влияния каждого фактора должен быть построен на основе законов случайности, что сильно усложняет процесс выработ­ ки критериев.

Одной из распространенных в практике теорией случайных про­ цессов является теория корреляции. Сущность этой теории заключа­ ется в следующем. Пусть задан случайный процесс, описываемый случайной величиной Y, значения которой обусловлены влиянием случайных величин

X » Хя, х а, х п.

Будем считать распределения этих случайных величин извест­ ными и подчиненными нормальному закону. Взаимосвязь случай­ ных величин исследуем по принципу «черного ящика». Иными слова­ ми, на интегральное воздействие факторов, описываемых случайными величинами

— Х 1{Хц',

. . . ; Х\m),

61


Х2

%22’

• • • і

%2т) t

Xn Xn {xn\y Xn2,

. . . \

(3.1)

^dmjt

«черный ящик» отвечает определенной реакцией, описываемой слу­ чайной величиной

Ѵ = Ѵ І У ѵ У і ...........

У m V

Механизм организации реакции «черного ящика» на внешнее воз­ действие остается неизвестным. Задача состоит в раскрытии функ­ ционального механизма организации реакции «черного ящика» на внешние раздражения.

Пусть мы располагаем достаточной информацией, характеризую­ щей влияние факторов, а также значением выхода, который легко установить экспериментально. В этом случае по критерию оценки веса влияющих факторов можно выбрать коэффициенты корреля­ ции, характеризующие взаимоотношения каждого соответствующе­ го влияющего фактора с реакцией системы. Для описания методов корреляционного анализа предположим, что для изучения механиз­ ма организации реакции «черного ящика» на различные воздействия произведено т экспериментов при одинаковых условиях. Результа­ ты экспериментов сведены в табл. 10.

Т а б л и ц а 10. Результаты экспериментов

В первом столбце указаны номера опытов, во втором приведена соответствующая реакция «черного ящика», фиксированная в соответ­ ствующем эксперименте. В последующих столбцах помещены при­ даваемые каждому фактору значения раздражений соответственно в каждом эксперименте. Например, значение у является количест­ венным изображением реакции «черного ящика» на интегральное раз­ дражение влияющих факторов в /-м эксперименте, а х^ — количе­ ственным значением раздражения, вносимым і-м фактором в у'-м эксперименте.

Далее предположим, что все результаты эксперимента описыва­ ются случайными величинами и законы распределения их подчинены нормальному закону распределения. Из теории вероятности извест­ но, что если основные численные характеристики случайных вели-

62

чин подчинены нормальному закону распределения, то они опреде­ ляются проще:

 

 

і=1

1

m

----

Дисперсию находим в виде

i=i

 

 

 

1 т

Г

Аналогично определяется среднеквадратическое отклонение

au = V D (у)

= V D ( х / ) =

(/ = 1, п).

Приведенные показатели являются основными характеристика­ ми соответствующих случайных величин, однако они не могут харак­ теризовать взаимоотношения случайных величин.

Для изучения взаимоотношения исследуемых случайных вели­ чин необходимо рассмотреть другие характеристики, отличающиеся от указанных выше. Одной из таких численных характеристик яв­ ляется соотношение

тт

kyxi — ДГ £ Уі ~ ДГ S Уі Хі‘ ---- т £ ХП и ~ 1>я).

называемое коэффициентом корреляции. Эта величина в определен­

ной мере

дает

представление^ взаимосвязи случайных величин

*/ и у (І =

1

«).

Однако с помощью коэффициента корреляции нельзя количест­ венно оценить силу или слабость корреляции этих случайных вели­ чин. Он указывает лишь на тесноту их связи. Коэффициент корреля­ ции — это показатель степени сложности законов зависимости меж­ ду двумя случайными величинами. Следовательно, он не может быть хорошим показателем для оценки веса влияющих факторов. Поэто­ му критерием оценки веса влияющих факторов необходимо вы­ брать другую численную характеристику. Для этого рассмотрим закон, устанавливающий взаимосвязь между двумя случайными

63