PZMS
.pdf
61
Рис.4.6 Ограничения для дан-  | 
	Рис.4.7 Анаморфоза для дан-  | 
ных рис.2.17  | 
	ных рис. 4.6  | 
Таким образом, функции плотности распределения могут рассматриваться как квазиполиномы моделей с памятью, а моменты появления в решениях характеристического уравнения псевдоположительных корней - как точки фазовых переходов.
В математической статистике смысл параметров распределений не выясняется. Это параметры аппроксимации. В действительности это влияние релаксационных механизмов, в результате чего параметры распределений могут определяться непосредственными измерениями релаксационных характеристик систем.
Таким образом, в качестве принципа полноты описания предъявляется требование одновременного анализа данных статического разреза и динамики. Это позволяет восстанавливать динамику по законам распределения и наоборот, также решать задачи идентификации релаксационныхсвойств по законамраспределения.
5. ВЛИЯНИЕ ПОЛНОТЫ ВЫБОРКИ НА РЕЗУЛЬТАТЫ МАЛЫЕ И БОЛЬШИЕ ВЫБОРКИ
Основные результаты теории вероятностей и математической статистики, используемые при анализе надежности, требуют, чтобы выборка была большая и однородная. Эти требования, как правило, не выполняются. Реально в задачах оценки надежности приходится иметь дело с малыми выборками, которые еще
62
исодержат данные о качественно разных состояниях.
Вматематической статистике существует метод Р.А.Фишера определения доверительного интервала средней для малых выборок [19].
Ввыборке размера n в предположении о нормальном рас-
пределении количество μ M 
Sμ распределено в t распределе-
нии Стьюдента с (n-1) степенями свободы при стандартном отклонении среднейSn Sμ 

n.
Здесь μ- выборочная средняя, M - вероятная средняя, sμ - выборочное стандартное отклонение.
Доверительный интервал в 95% равен
μ t0.05Sn M μ t0.05Sn.
Величина t0,05 определяется по таблице t Фишера [19] (приложение 1).
Таким образом, если производится большое число выборок размера n и вычисляются доверительные интервалы, то полученный в результате интервал должен включать истинную величину вероятной средней приблизительно в 95% случаев.
Получаемая при этом оценка задает диапазон однородности измерений, выход за который связан с изменением качественного состояния системы.
Использование этого метода осложняется необходимостью обращения к таблице t Фишера. Кроме того, предположение о нормальном распределении отношения разности выборочной и вероятной средних к стандартному отклонению средней для малых выборок проблематично, т.к. нормальный закон справедлив для больших однородных выборок.
В принципе целесообразно иметь метод анализа данных, позволяющий оперативно делать оценки границ, классифицирующих диапазонов для выборок, содержащих данные, принадлежащие к качественно разным группам.
Как показано выше, постоянство относительной ошибки результатов измерений определяется линейной зависимостью между логарифмом характеристики состояния системы и аргументом.
Точки, в которых данные отклоняются от линейной зависи-
63
мости, называются критическими. Пара соседних критических точек определяет диапазон однородности эмпирических данных для анализируемой выборки.
Обработка большого количества малых статистических выборок с построением их ранговых распределений в полулогарифмическом масштабе и сопоставление полученных критических точек с результатами обработки тех же данных по Фишеру показала их полную идентичность по возможностям классификации однородных групп в малых статистических выборках.
В качестве примера рассмотрим данные о числе отклонений от конструкторской документации при изготовлении партии сложных технических систем, табл.5.1.
  | 
	
  | 
	Таблица 5.1  | 
Ранговое распределение числа отказов, N  | 
||
N  | 
	lnN  | 
	(N-m)2  | 
135  | 
	4.905275  | 
	2079.36  | 
127  | 
	4.844187  | 
	1413.76  | 
107  | 
	4.672829  | 
	309.76  | 
102  | 
	4.624973  | 
	158.76  | 
101  | 
	4.615121  | 
	134.56  | 
91  | 
	4.51086  | 
	2.56  | 
89  | 
	4.488636  | 
	0.16  | 
86  | 
	4.454347  | 
	11.56  | 
81  | 
	4.394449  | 
	70.56  | 
80  | 
	4.382027  | 
	88.36  | 
75  | 
	4.317488  | 
	207.36  | 
73  | 
	4.290459  | 
	268.96  | 
72  | 
	4.276666  | 
	302.76  | 
66  | 
	4.189655  | 
	547.56  | 
56  | 
	4.025352  | 
	1115.56  | 
Прежде всего, данные сортируются для построения рангового распределения в полулогарифмическом масштабе. По оси абсцисс откладывается порядковый номер точки, а по оси ординат – логарифм числа отклонений, рис.5.4.
64
В теории размерностей известно, что в критических точках значение безразмерной переменной порядка единицы [20]. Отсюда в связи с тем, что показатель степени у экспоненты величина безразмерная,
kt* 1 и y* e. y0
В результате можно ожидать, что изменение характеристик цепных реакций квантуется изменением функции в e=2,71828… раз. Это простое следствие теории размерностей.
На рис.5.1 приведены результаты первых экспериментов по исследованию динамики радиоактивного распада, проведенных в
1905 г. [21].
  | 
	
  | 
	Радиоактивный распад  | 
	
  | 
	
  | 
||
  | 
	4  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
  | 
	3.5  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
  | 
	3  | 
	
  | 
	
  | 
	1  | 
	
  | 
	
  | 
  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
|
  | 
	2.5  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
E  | 
	2  | 
	
  | 
	
  | 
	1  | 
	
  | 
	
  | 
ln  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
|
  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
|
  | 
	1.5  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
  | 
	1  | 
	
  | 
	
  | 
	1  | 
	
  | 
	
  | 
  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
|
  | 
	0.5  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
  | 
	0  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
  | 
	0  | 
	10  | 
	20  | 
	30  | 
	40  | 
	50  | 
t, сутки
Рис.5.1 Динамика радиоактивного распада, по оси ординат активность в относительных единицах, масштаб натуральный ло-
гарифмический
Здесь очевидны разрывы, характеризующие переходы через критические состояния. Расположены они регулярно через равные интервалы, равные для натуральных логарифмов активности единице, что соответствует изменению самой активности между разрывами в е раз.
При построении шкал отношений проблема состоит в опре-
65
делении таких чисел, которые единым образом определяют структуру природных систем. Оказалось, что такие числа связаны, в частности, с числом Непера е=2,71828..., которое является основанием натуральных логарифмов. При решении задачи о реализации умножения через сложение основоположник логарифмов Джон Непер сопоставил члены арифметической и геометрической прогрессий [22]. Он разбил единичный интервал на n равных частей и каждому члену полученной арифметической прогрессии поставил в соответствие члены геометрической прогрессии. Последний член геометрической прогрессии, соответствующий 1 в арифметической прогрессии, при n, стремящемся к бесконечности, и определяет значение числа Непера. Подобного рода механизмы взаимодействий между членами арифметической и геометрической прогрессий широко представлены как в природе, так и в рукотворных системах. Так, клавиши пианино представляют арифметическую прогрессию, а частоты соответствующих им нот - геометрическую.
На рис.5.2-5.7 построены ранговые распределения для данных величины сердечного ритма с разным количеством элементов в линейном масштабе и логарифмических масштабах по оси ординат. Пунктирными линиями обозначены границы, полученные методом Фишера.
Рис.5.2 Ранговое распределение  | 
	Рис.5.3 Логарифмический мас-  | 
для данных сердечного ритма,  | 
	штаб для данных рис.5.2  | 
n = 10  | 
	
  | 
66
Рис.5.4 Ранговое распределение  | 
	Рис.5.5 Логарифмический мас-  | 
для данных сердечного ритма,  | 
	штаб для данных рис.5.4  | 
n = 20  | 
	
  | 
Рис.5.6 Ранговое распределение  | 
	Рис.5.7 Логарифмический мас-  | 
для данных сердечного ритма,  | 
	штаб для данных рис.5.6  | 
n = 30  | 
	
  | 
Из представленных выше рисунков (рис.5.2-рис.5.6) видно совпадение интервалов однородности намечаемых методом Фишера и ранговым распределением. Метод Фишера отмечает один диапазон, ранговое распределение намечает несколько диапазонов.
6. ПОСТРОЕНИЕ СТРУКТУР СИСТЕМ ПО ДАННЫМ СТАТИСТИЧЕСКИХ ВЫБОРОК
Возможности теории химического строения А.М.Бутлерова
67
явно к настоящему времени не используются в других областях. Рассмотрим реализацию ее результатов на примере природной популяции локальных стад одного из видов лососевых Камчатки - нерки. Лососи Камчатки образуют практически полностью изолированные локальные стада. В ходе исследований генетической структуры популяций локальных стад нерки оз. Азабачьего ежегодно в 1971-1973 гг. на ряде нерестилищ определялось количество гомозиготных (ВВ и B’B’) и гетерозиготных (BB’) особей всего на 56 нерестилищах [23]. Рассмотрим, как из полученных в ходе исследований данных определить формулу генетической структуры этих стад.
Обработка данных показала наличие двух принципиально разных соотношений между частотами генов, в то же время стабильность соотношений генов внутри каждого из этих соотношений оказалась достаточно высокой.
  | 
	Зависимость между процентами противоположных  | 
	
  | 
||
  | 
	гомозиготных генотипов локальных стад нерки  | 
	
  | 
||
  | 
	
  | 
	
  | 
	
  | 
	0  | 
-2  | 
	-1.5  | 
	-1  | 
	-0.5  | 
	-0.5 0  | 
  | 
	
  | 
	
  | 
	
  | 
	-1  | 
  | 
	
  | 
	
  | 
	
  | 
	-1.5  | 
ln BB  | 
	
  | 
	
  | 
	
  | 
	-2  | 
  | 
	
  | 
	
  | 
	-2.5  | 
|
  | 
	
  | 
	
  | 
	
  | 
	-3  | 
  | 
	
  | 
	
  | 
	
  | 
	-3.5  | 
  | 
	
  | 
	
  | 
	
  | 
	-4  | 
  | 
	
  | 
	
  | 
	
  | 
	-4.5  | 
  | 
	
  | 
	ln B'B'  | 
	
  | 
	
  | 
Рис.6.1. Зависимость между процентами гомозиготных гено-  | 
||||
  | 
	типов в локальных стадах нерки  | 
	
  | 
||
На рис. 6.1 приведены экспериментальные данные, характеризующие соотношение между частотами встречаемости геноти-
  | 
	
  | 
	
  | 
	68  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
пов B’B’ и ВВ. Регрессионная зависимость между этими часто-  | 
||||||||
тами выражается соотношением  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
|||
  | 
	
  | 
	
  | 
	BB a  | 
	
  | 
	
  | 
	,  | 
	
  | 
	(6.1)  | 
где а - постоянный коэффициент.  | 
	B B  | 
	
  | 
	
  | 
	
  | 
||||
  | 
	
  | 
	
  | 
	
  | 
	
  | 
||||
Введем новую переменную, характеризующую линию регрессии,  | 
||||||||
  | 
	
  | 
	
  | 
	k1 BB  | 
	
  | 
	,  | 
	
  | 
	(6.2)  | 
|
  | 
	
  | 
	
  | 
	
  | 
	B B  | 
	
  | 
	
  | 
	
  | 
	
  | 
и в функции от этой переменной построим зависимость час-  | 
||||||||
тоты встречаемости третьего генотипа, рис.6.2.  | 
	
  | 
	
  | 
||||||
  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
	-1.5  | 
-1.4  | 
	-1.2  | 
	-1  | 
	-0.8  | 
	-0.6  | 
	
  | 
	-0.4  | 
	-0.2  | 
	-1.6 0  | 
ln [B'B']  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
	-1.7  | 
  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
	-1.8  | 
|
  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
	-1.9  | 
|
+(1/2)  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
	-2  | 
  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
	-2.1  | 
|
[BB]  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
|
  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
	-2.2  | 
|
ln  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
|
  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
	-2.3  | 
|
  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
|
  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
	-2.4  | 
  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
	-2.5  | 
ln BB'
Рис.6.2. Зависимость процента гетерозиготного генотипа локальных стад нерки от линии регрессии (6.1) как аргумента
Как видно из результатов, представленных на рис.6.2, большое количество точек попадает в совокупность, которая характеризуется линией регрессии вида
BB C  | 
	3  | 
	
  | 
	
  | 
	3  | 
	
  | 
	,  | 
	(6.3)  | 
  | 
	2  | 
	
  | 
|||||
BB  | 
	
  | 
||||||
(  | 
	
  | 
	
  | 
	
  | 
	B B  | 
	)  | 
	
  | 
где С - константа, тогда как большая часть данных отклоняется от этой зависимости. Возьмем линию регрессии (6.3) за новую
69
переменную.
В функции от наиболее редко встречающегося генотипа B’B’ для точек, принадлежащих линии регрессии (6.3), новая пере-
2
менная BB (BB
B B ) 3 должна давать постоянный уровень. Данные, приведенные на рис.6.3, показывают, что разбросы наблюдаются для редко встречающегося генотипа B’B’ при частоте их встречаемости меньшей 12%. Если частоты встречаемости генотипа B’B’ превышают эту величину, соотношения между генотипами этих стад соответствуют постоянству зависимо-
сти (6.3).
  | 
	-2  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
  | 
	0  | 
	0.1  | 
	0.2  | 
	0.3  | 
	0.4  | 
	0.5  | 
(ln B'B')/2]  | 
	-2.1  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
-2.2  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
|
-2.3  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
|
BB +  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
|
-2.4  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
|
+ (2/3[)ln  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
|
-2.5  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
|
-2.6  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
|
ln BB'  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
|
-2.7  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
|
  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
|
  | 
	-2.8  | 
	
  | 
	
  | 
	
  | 
	
  | 
	
  | 
  | 
	
  | 
	
  | 
	
  | 
	B'B'  | 
	
  | 
	
  | 
Рис.6.3. Оценка диапазона стабильности линии регрессии (6.2)  | 
||||||
Таким образом, выделяется группа локальных стад, для которых соотношение (6.3) определяет взаимосвязи между частотами генотипов, при этом в данную группу попадают стада с частотой встречаемости генотипа B’B’ 12%. Это условие отделяет речные нерестилища, на которых реализованы существенно более жесткие условия отбора по сравнению с озерными нерестилищами, где реализуется условие B’B’ 12%.
  | 
	
  | 
	70  | 
	
  | 
	
  | 
На основе соотношения (6.3) можно, используя неравенство  | 
||||
Коши, в соответствии с которым среднее геометрическое не пре-  | 
||||
восходит среднего арифметического, оценить элементарную  | 
||||
формулу такой структуры:  | 
	
  | 
	
  | 
	
  | 
|
  | 
	
  | 
	
  | 
	
  | 
	-4  | 
-2  | 
	-1.5  | 
	-1  | 
	-0.5  | 
	0  | 
  | 
	
  | 
	
  | 
	
  | 
	-5  | 
  | 
	
  | 
	
  | 
	
  | 
	-6  | 
BB'  | 
	
  | 
	
  | 
	
  | 
	-7  | 
+ 6 ln  | 
	
  | 
	
  | 
	
  | 
|
  | 
	
  | 
	
  | 
	-8  | 
|
B'B'  | 
	
  | 
	
  | 
	
  | 
	-9  | 
ln  | 
	
  | 
	
  | 
	
  | 
|
  | 
	
  | 
	
  | 
	
  | 
|
  | 
	
  | 
	
  | 
	
  | 
	-10  | 
  | 
	
  | 
	
  | 
	
  | 
	-11  | 
  | 
	
  | 
	
  | 
	
  | 
	-12  | 
  | 
	
  | 
	ln BB  | 
	
  | 
	
  | 
Рис.6.4. Стада с частотой встречаемости генотипа B’B’ 12%  | 
||||
[(B’B’)(BB’)3 (ВВ)2] 1/6 (1/6)[(ВВ)2(B’B’)(BB’)3], где в правой части запись соответствует элементарной формуле, принятой в химии, т.е. (ВВ)2(B’B’)(BB’)3 означает, что на единицу (B’B’) приходится 2 единицы (ВВ)) и 3 единицы (BB’).
Для построения структурной формулы в соответствии с теорией химического строения А.М.Бутлерова требуется знать валентности - число связей между элементами структур. Полученная элементарная формула соответствует формуле угольной кислоты Н2СО3, в которой углерод имеет 4 валентности, кислород - 2 и водород - 1. В результате структурная формула генетической структуры локальных стад нерки на речных нерестилищах имеет вид рис.6.5.
Рассмотрением данных о генетической структуре озерных локальных стад можно аналогичными преобразованиями полу-
