Файл: Кадыров, Х. К. Синтез математических моделей биологических и медицинских систем.pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 21.10.2024
Просмотров: 74
Скачиваний: 0
следовательно, неравенство (2.17) справедливо. Таким образом, сочетание метода кодирования с методом упаковки дает то преиму щество, что в ячейку можно записать хотя бы два кодовых числа.
Для объяснения сути комбинированного метода приведем конк ретный пример. Проанализируем вариационный ряд предыдущего примера с точки зрения применения комбинированного метода.
Перепишем вариационный ряд: 120, 116, ПО, 112, 160, ПО, ...
.., 115 (я = 5000). Дано
max [хъ х2, . . . , хп) = 160,
min {xlt х2, . . . , ха) = 110,
N = max [хи х2, ... , хп) — min {хи х2, ... , хп) = 501о = 628.
Присвоим код каждому члену ряда диапазона по восьмеричной системе счисления, начиная от 018 до 638. Итак, kg = £код = 2.
На ЭЦВМ типа «Минск-22» для данного ряда определим количе ство ячеек по формуле (2.6):
Nз |
5 0 + 1 |
5000 |
842. |
|
37 |
I |
37 I |
||
|
3 • 2 |
[ 3 |
- 2 j |
|
Следовательно, вместо 5000 используется всего лишь 842 ячейки. Если при вводе использовать только кодирование, то, как следует из формулы (2.15), будет использовано
5000
N. = 885
37
4 ■3
m i n
ячеек.
Для метода упаковки это количество вычисляется по формуле
(2. 12):
|
n ; = |
5000 |
5000 |
1734. |
|
||
|
37 |
3 |
|
||||
|
|
|
|
|
|||
|
|
4 |
• 3 |
|
|
|
|
Таким образом, простой метод |
упаковки |
позволяет |
умень |
||||
шить |
число ячеек меньше, |
чем в три раза, а |
метод кодирования |
||||
с упаковкой кодированных |
чисел |
и комбинированный |
метод — |
||||
почти |
в шесть раз. |
|
|
|
|
|
|
3.Алгоритмы поиска законов распределения
ипараметров вариационной статистики
Алгоритм поиска законов распределения для случайных величин. При исследовании различных физических, технических, биологи ческих и медицинских задач сталкиваются со случайными явления ми — случайными вариациями. Эти вариации всегда связаны с наличием каких-либо вероятностных факторов, влияющих на исход наблюдения.
50
Особенно сильно проявляются случайности в различных меди цинских исследованиях. Изучение факта случайности в медицин ских исследованиях с помощью математического моделирования на базе математического аппарата случайных процессов является весь ма перспективным и многообещающим. Математическое моделирова ние, описывающее случайные процессы, дает возможность опреде лить все основные численные характеристики сложных медико-био логических процессов [13].
На практике приходится оперировать с ограниченным количе ством экспериментальных данных. Следовательно, к методам обоб щения результатов эксперимента следует предъявлять требования, которые могли бы выявить характерные черты наблюдаемого явле ния, а все несущественные элементы, связанные со случайными от клонениями, отстранить. В связи с этим возникает задача сглажива ния или выравнивания статистических данных. При обработке ста тистических материалов, полученных в результате эксперимента, необходимо подбирать теоретическую кривую распределения. Задача нахождения теоретической кривой для статистических рядов со стоит в .подборе плавной кривой распределения, наилучшим образом описывающей данное распределение.
Задача о наилучшем выравнивании статистических рядов, как и вообще задача о наилучшем аналитическом представлении эм пирических данных, в значительной мере неопределенная и решение ее зависит от того, что считать наилучшим. Одним из универсальных математических методов для выбора наилучшей кривой является метод наименьших квадратов. При этом вопрос о том, в каком именно классе функций следует искать наилучшее приближение, решается уже не математически, а на основании соображений, связанных с природой исследуемого объекта, с учетом характера полученной эмпирической кривой. Принципиальный вид выбранной теоретичес кой кривой зависит от некоторых параметров. Задача выравнивания статических рядов переходит в задачу выбора параметров. В каче стве критерия выбора неизвестных параметров, наилучшим образом обеспечивающих описание исследуемого процесса, можно восполь
зоваться |
следующими |
характеристиками случайных величин [7]: |
|
а) |
математическим ожиданием исследуемых случайных величин, |
||
определяемым в случае, если случайные величины дискретны: |
|||
|
|
/И (х) = |
V Хірі. |
|
|
|
/=і |
Для |
непрерывных |
случайных |
величин |
|
|
|
оо |
|
|
М (х) = |
j xf (X) dx, |
— оо
где xt — конкретные значения случайной величины; pt — вероят ность принятия случайными величинами конкретного численного значения; f (х) — плотность распределения случайных величин; / (х) dx — элемент вероятности;
4* |
51 |
■б) центральными моментами случайных величин разного по
рядка: |
|
|
|
|
|
|
|
Ml = |
м [X ] = 2 |
[*, — М (*)] Рі = |
О, |
|
|||
|
|
;=і |
|
|
|
|
|
м2 = |
М [ > ] = |
V |
[xt - м (х)Г Рі = |
а , - M l, |
|
||
|
|
i=i |
|
|
|
|
|
Мз = |
м [X3] = |
2 |
Iхі — Мх?Рі = “ з — з а2Мх + 2 М\, |
|
|||
|
|
і=1 |
|
|
|
|
|
м4 = М [ * 4] = S |
|
— М |
Рі = а 4 — 4 0 ^ + 7 а 2/И* — 3 М Х4, |
||||
|
|
і=і |
|
|
|
|
|
где |
|
|
|
|
п |
I |
|
|
|
|
|
|
|
||
|
|
|
|
а / = 2 |
х^ ‘'; |
|
|
|
|
|
|
|
*=1 |
|
|
|
|
|
|
м х = 2 |
|
|
|
|
|
|
|
|
«=1 |
|
|
в) коэффициентом |
асимметрии, |
характеризующим |
меру асим |
||||
|
|
|
9 |
_ _рз_ _ |
«з — Шхаг+ 2М-1 |
|
|
метрии законов распределений случайных величин: |
|
||||||
|
|
|
k |
°3 |
[V а , — МХ2]3 |
|
|
С помощью характеристики Sk получим некоторое |
предположе- |
ние о смещении кривой распределения случайной величины вправо или влево от исходного закона распределения. Если Sk > 0, то кривая распределения будет сдвинута в левую сторону (положи тельная асимметрия), если S k < 0, то кривая распределения будет смещена вправо (отрицательная асимметрия);
г) коэффициентом крутизны закона распределения случайных величин или эксцесс:
Ех = - |
3 |
^4 _ |
а„ — 4а3Мх + 7М2х — ЗМ* |
|
■3. |
||
|
|
|
[Ѵа2- М 2х]* |
С помощью этого показателя можно получить представление о крутизне или плосковершинности закона распределения случайных величин. Если Ех > 0, то кривая распределения будет иметь острую вершину, а если Ех < 0, то она будет иметь плоскую вершину. На конец, если Ех = 0, то мы получим подтверждение того, что закон распределения исследуемых случайных величин подчиняется нор мальному закону распределения.
Счедовательно, имея в распоряжении все эти показатели, легко делать выводы о законе распределения исследуемых случайных ве личин. Далее остается принять критерий близости законов распре деления рассматриваемых случайных величин на сопоставление
52
с законом случайных процессов, описываемых соотношением
Ф — fi (х>а) (/ = 1, е),
аналитический выход которых задается последовательно. В мате матике разработано множество критериев для оценки близости зако нов распределений исследуемых случайных величин тому или иному закону распределения. Для такой оценки мы принимаем критерий «хи-квадрат» X2.
Процесс оценки близости законов распределений с помощью X2 выполняется в следующей последовательности.
1.Все численные характеристики исследуемой случайной ве личины вносятся в предъявляемый закон, т. е. Ф = / (ах, Мх, х, і).
2.Вычисляется показатель согласованности по формулам
X2 = я У |
(Р/ — Р{) |
(2.18) |
2_j |
|
|
і=і |
|
|
или более проще |
|
|
( т і — n p t ) |
|
|
|
^P: |
|
где Pi = ~ — статистическая |
вероятность; |
mt — благоприятст |
вующий случай, т. е. количество повторений данной величины при я опытах; я — общее число наблюдений.
Как известно показатель согласованности X2 зависит от степеней свободы рассматриваемой случайной величины г. Это число г равно числу разрядов минус число связей, наложенных на частоту рр
1)
|
t=i |
|
|
|
2) |
М [X] = 2 хіРі = |
Мх |
|
|
|
і=і |
|
|
|
3) |
М[х — Мх]2 = |
2 |
(Xt - Мху Pi = Dx, |
|
|
|
i = |
1 |
|
где Xi — дискретные значения |
случайных |
величин; Мх — матема |
||
тическое ожидание случайных |
величин; |
Dx — дисперсия случай |
ных величин, характеризующая разброс дискретных значений слу чайной величины вокруг своего математического ожидания.
Исходя из численных значений X2, с помощью формул (2.18) вычислен показатель согласованности и построены специальные
таблицы. При |
помощи этих таблиц для |
каждого значения |
и чис |
||
ла степени свободы можно выбрать вероятность того, что |
закон |
||||
распределения |
исследуемой случайной |
величины |
близок |
к |
за |
кону распределения ф/ = //(х ) Таким |
образом, с |
помощью |
X2 |
можно оценить степень согласованности закона распределения ис следуемой случайной величины с законом, предъявляемым в виде
53