Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

PZMS

.pdf
Скачиваний:
1
Добавлен:
15.03.2025
Размер:
1.13 Mб
Скачать

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 (BBB 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.

Рассмотрением данных о генетической структуре озерных локальных стад можно аналогичными преобразованиями полу-