PZMS
.pdf
41
фозы, или можно определить варьированием, выбрав значение, при котором данные линеаризуются на большем интервале.
По коэффициенту наклона определяется значение параметра γ=0.4. Анаморфоза воспроизводит исходное значение параметра γ.
2.9. Гамма-распределение
Плотность вероятности гамма-распределения:
x
xk 1 e C f (x) Г(k) Ck .
Рассмотрим применение анаморфозы на модельных данных, полученных функцией rgamma программы Mathcad (рис.2.37), генерирующей случайные числа с плотностью гаммараспределения, параметры k = 2, C = 1.
Анаморфоза для гамма-распределения:
f (x) (k 1) 1 . f (x) x C
Построение данных, соответствующих нормальному распределе-
нию, дает спрямление в координатах |
f (x) |
~ |
1 |
(рис.2.38). |
|
|
|||
|
f (x) x |
|
||
Рис.2.37 Модельные данные, |
Рис.2.38 Анаморфоза для данных, |
полученные функцией rgamma. |
полученные функцией rgamma. |
42
По коэффициенту наклона прямой, равному (k-1), определяется значение параметра k ≈ 2. Анаморфоза воспроизводит исходное значение параметров k.
2.10. Логистическое распределение
Плотность вероятности логистического распределения:
f (x) k F(x) (1 F(x)) или |
f (x) |
k e |
k (x x0) |
|
(1 e k (x x0))2 |
||||
|
|
|||
i
где F(x) – функция распределения. F(xi) f (xi) x
j 1
Рассмотрим применение анаморфозы на модельных данных, полученных функцией rlogis программы Mathcad (рис.2.39), генерирующей случайные числа с логистической плотностью распределения, параметр k = 0,5.
Анаморфоза для логистического распределения:
f (x) k (1 F(x)).
F(x)
Построение данных, соответствующих логистическому распреде-
лению, дает спрямление в координатах f (x) ~ F(x) (рис.2.40).
F(x)
Рис.2.39 Модельные данные, |
Рис.2.40 Анаморфоза для данных, |
полученные функцией rlogis. |
полученных функцией rlogis. |
43
По коэффициенту наклона прямой значение параметра k 0,51. При построении этой анаморфозы некоторое количество начальных точек не используется из-за того, что сумма накопленных значений F(x) в начале имеет большую ошибку.
2.11. Распределение Гомперца
Плотность вероятности распределения Гомперца:
f (x) k F(x) ln(F(x)) или |
|
A |
e k x |
, |
|
||||
f (x) A e k xe k |
||||
|
i |
|
||
где F(x) – функция распределения. F(xi) f (xi) |
x |
|||
|
j 1 |
|
||
Рассмотрим применение анаморфоз на модельных данных для распределения Гомперца (табл. 6), A = 4 , k = 0,6 (рис.2.41).
Таблица 2.7.
Модельные данные для распределения Гомперца, A=4, k=0.6
xi |
f(xi) |
xi |
f(xi) |
xi |
f(xi) |
xi |
f(xi) |
0 |
0,00509 |
2,5 |
0,20165 |
5 |
0,14290 |
7,5 |
0,04126 |
0,5 |
0,02123 |
3 |
0,21965 |
5,5 |
0,11537 |
8 |
0,03116 |
1 |
0,05656 |
3,5 |
0,21652 |
6 |
0,09109 |
8,5 |
0,02342 |
1,5 |
0,10816 |
4 |
0,19820 |
6,5 |
0,07075 |
9 |
0,01753 |
2 |
0,16176 |
4,5 |
0,17175 |
7 |
0,05428 |
9,5 |
0,01309 |
Анаморфоза для распределения Гомперца: |
|
|
|
|
|
||||
|
f (x |
x) f (x x) |
|
2 e k |
x ek |
x |
A |
k x. |
|
ln(ln( |
|
|
)) ln |
|
|
||||
|
f 2(x) |
|
|||||||
|
|
|
|
|
|
k |
|
||
Построение данных, соответствующих распределению Гомперца,
дает спрямление в координатах ln(ln( |
f (x |
x) f (x x) |
)) ~ x. |
|
|
||
|
|
f 2(x) |
|
По коэффициенту наклона прямой значение параметра k 0.6, а
0,6 e 0,5033 |
x 4 (рис.2.42). |
A 2 e k x ek |
Следующая анаморфоза для распределения Гомперца:
|
44 |
|
|
f (x) |
|
ln |
|
ln(A) k x |
|
||
|
|
|
|
F(x) |
|
Рис.2.41 Функция плотности рас- |
Рис.2.42 Анаморфоза для дан- |
пределения для данных табл.2.7 |
ных табл.2.7 |
Построение данных, соответствующих распределению Гомперца,
f (x) |
|
~ x |
(рис.2.43). |
дает спрямления в координатах ln |
|
||
|
F(x) |
|
|
Рис.2.43 Анаморфоза для данных табл.2.7
По коэффициенту наклона прямой значение параметра k 0.6, а A e1.2559 3.5. При построении этой анаморфозы, некоторое количество начальных точек не используется из-за того, что сумма накопленных значений F(x) в начале имеет некоторую ошибку.
45
2.12. Распределение Вейбулла
Плотность вероятности распределения Вейбулла:
f (x) |
k |
x k 1 |
||
|
|
|
|
|
λ |
|
|||
|
|
λ |
||
x k e λ .
Рассмотрим применение анаморфоз на модельных данных для распределения Вейбулла (табл.2.7), k = 3 , λ = 5 (рис.2.44).
Таблица 2.7.
Модельные данные для распределения Вейбулла, k=3, λ = 5
xi |
f(xi) |
xi |
f(xi) |
xi |
f(xi) |
xi |
f(xi) |
0,5 |
0,00599 |
3 |
0,17404 |
5,5 |
0,19182 |
8 |
0,02556 |
1 |
0,02381 |
3,5 |
0,20863 |
6 |
0,15348 |
8,5 |
0,01275 |
1,5 |
0,05256 |
4 |
0,23013 |
6,5 |
0,11269 |
9 |
0,00570 |
2 |
0,09005 |
4,5 |
0,23444 |
7 |
0,07563 |
9,5 |
0,00227 |
2,5 |
0,13237 |
5 |
0,22073 |
7,5 |
0,04619 |
10 |
0,00081 |
Анаморфоза для распределения Вейбулла:
|
f (x) |
k |
x |
|||||
ln |
|
ln |
|
|
(k 1)ln |
|
. |
|
|
|
|
||||||
|
|
|
λ |
λ |
||||
|
1 F(x) |
|||||||
Построение данных, соответствующих распределению Вейбулла,
|
f (x) |
|
дает спрямления в координатах ln |
|
~ ln x (рис.2.45). |
|
|
|
|
1 F(x) |
|
Рис.2.44 Функция плотности рас- |
Рис.2.45 Анаморфоза для данных |
пределения для данных табл.2.8. |
табл. 2.8. |
46
По коэффициенту наклона прямой (рис.2.45) определяется значение параметра k 2.09 1 3.09. При построении этой анаморфозы, некоторое количество начальных точек не используется изза того, что сумма накопленных значений F(x) вначале имеет большую ошибку.
Таким образом, использование анаморфоз позволяет определить принадлежность (либо отсутствие принадлежности) эмпирических данных к распределению заданного вида. В связи с тем, что число наиболее часто используемых распределений ограничено, перебор анаморфоз для этих распределений позволяет определить вид распределения наилучшим образом соответствующего анализируемым эмпирическим данным.
3. ПРЕДЕЛЫ ПРИМЕНИМОСТИ ОСНОВНЫХ ФУНКЦИЙ РАСПРЕДЕЛЕНИЯ МАТЕМАТИЧЕСКОЙ СТАТИСТИКИ
Использования статистических выборок в целях определения действительного поведения сложной системы и управления ее развитием ставит перед необходимостью иметь действенные методы обработки данных. По результатам измерений мы имеем информацию либо о динамике фиксированной системы, когда аргументом является время (правая часть рис. 3.1), либо для фиксированного момента времени характеристики состояния совокупности объектов определенного класса (левая часть рис.3.1).
Рис.3.1 Схема соответствий структуры данных статического разреза и динамики
47
Связь между этими методами представления информации дается эргодической теоремой (иногда ее называют эргодической гипотезой), в соответствии с которой средние по времени и по фазовому пространству равны между собой [9]. Отсюда появляется часто используемая в практике научных и практических исследований возможность по динамике конкретного объекта судить о характеристиках совокупности таких объектов, либо по данным о характеристиках статистической совокупности в данный момент времени характеризовать динамику индивидуального объекта.
Например, посадив малька рыбки в аквариум, можно наблюдать во времени рост ее длины и тем определять кривую роста. С другой стороны, можно спустить воду из пруда, собрать всех рыбок данного вида и построить функцию плотности распределения их по размерам. Как в одном, так и в другом случае мы получим информацию о том, какие рыбки бывают на разных стадиях развития.
На кривых роста мы наблюдали наличие критических точек, в которых принципиально меняется характер развития. В связи с тем, что в соответствии с эргодической теоремой динамика и данные статического разреза взаимосвязаны, естественно ожидать, что и в законах распределения математической статистики, представляющих информацию о данных статического разреза , будут реализованы критические состояния, являющиеся проекциями соответствующих критических возрастов из динамики развития.
Этим ставится проблема определения границ диапазонов однородности при статистической обработке экспериментальных данных. Знание точных значений границ сохранения свойств статистических совокупностей диктуется необходимостью решения широкого круга прикладных задач.
Связь между входом системы x(t) и выходом y(t) определяется разнесением по времени на интервал τ. В результате сигнал на выходе системы в данный момент времени характеризуется входным сигналом в момент времени (t-τ) [10-13].
Тогда характеристики линейной системы представляются
|
|
|
|
|
|
|
|
|
48 |
|
|
|
|
уравнением: |
|
d(m) y(t) |
|
|
|
|
(m 1) y(t) |
|
|
|
|||
am |
|
am 1 |
d |
... a0 y |
|
||||||||
|
dtm |
|
dt(m 1) |
(3.1) |
|||||||||
|
|
|
|
|
|
|
|
|
|||||
|
|
d(p)x(t τ) |
|
|
|
|
d(p 1)x(t τ) |
|
|||||
b |
p |
b |
p 1 |
... b |
x(t τ) |
||||||||
|
dt(p 1) |
||||||||||||
|
|
dtp |
|
|
0 |
|
|||||||
Аналитическое выражение взаимосвязей между динамическим описанием и данных статического разреза в соответствии с эргодической гипотезой можно проследить, рассматривая характеристическое уравнение для представленного выше уравнения.
Будем искать решение в экспоненциальной форме x ezt, y H(z) ezt .
При этом H(z) y, т.е. определяет соотношение сигналов на x
входе и выходе системы.
Тогда характеристическое уравнение для (3.1) записывается в виде:
H(z)(amzm am 1amzm 1 ... a0)
(3.2)
(bpzp bp 1zp 1 ... b0) e zτ 0
Левая часть этого уравнения называется характеристическим квазиполиномом.
Это уравнение не содержит времени в отличие от (3.1) и определяет только соотношение между параметрами рассматриваемой системы, которые могут быть определены из анализа данных для группы объектов в фиксированный момент времени. В связи с этим квазиполином (3.2) можно рассматривать как некоторую функцию плотности распределения математической статистики. Такие функции, как известно, используются при обработке экспериментальных данных, относящихся к определенному моменту времени. Покажем, что основные практически используемые функции плотности вероятности являются частными случаями квазиполинома (3.2).
Из характеристического квазиполинома соотношение сигналов на входе и выходе системы определяется соотношением
49
|
(b |
p |
zp |
b |
p 1 |
zp 1 |
... b ) |
|
|||||||
H(z) |
|
|
|
|
|
|
|
|
0 |
|
e zτ |
(3.3) |
|||
(a |
|
zm a |
|
|
a |
|
|
|
|
||||||
|
m |
m |
1 |
m |
zm 1 ... a ) |
|
|||||||||
|
|
|
|
|
|
|
|
|
0 |
|
|
||||
которое, при z=iv, называется частотной характеристикой. Условие устойчивости при этом определяется как равенство
нулю знаменателя частотной характеристики
amzm am 1amzm 1 ... a0 0
Однако для общего случая, когда корни характеристического уравнения являются комплексными, они могут быть и псевдоположительными, т.е. иметь положительную действительную часть при ненулевом коэффициенте при мнимой части, что также соответствует потере устойчивости.
Назовем соотношение (3.2) комплексной характеристикой для случая z=u+iv.
При нахождении решений уравнения (3.2) часто пользуются оценкой весомостей членов, входящих в его состав. Выбрав члены с наибольшей весомостью, определяют характер решения. При этом оказывается, что при различных диапазонах значений переменных и параметров, превалирующими являются различные члены. Так, при z>>0 доминируют члены с максимальными значениями показателей степени m и p: Тогда в случае одного доминирующего члена в ряду функций текущего значения времени и одного в ряду функций отклоняющегося значения аргумента
H(z) |
bpzp |
e zτ |
|
bp |
zp m e zτ |
(3.4) |
|
amzm |
am |
||||||
|
|
|
|
|
Заметим, что основные законы распределения, используемые в математической статистике, являются частными случаями характеристического квазиполинома уравнения (3.1), соответствующими комплексной характеристике (3.3).
Все эти распределения представляют собой показательные функции, либо их комбинацию со степенными функциями. В результате, в качестве основного требования к моделям обработки данных о надежности сложных систем требуется ввести такое описание, которое будет включать класс этих функций в качестве
50
частных случаев. В табл.3.1 эти функции сопоставлены с параметрами уравнений (3.4).
Таблица 3.1
Параметры уравнения (3.4) и соответствующие им функции плотности вероятностей
Распределение |
|
|
|
bp |
|
|
|
p m |
|
|
τ |
|
|
|
|
|
|
|
H(z) |
|
|
|
|||||||||||||
|
am |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
Нормальное |
|
|
|
1 |
|
|
|
|
|
|
0 |
|
|
z |
|
|
|
|
1 |
|
|
|
|
|
|
|
z2 |
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
e |
2σ2 |
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
||||||||||||
σ |
|
|
2π |
|
2σ |
|
|
|
|
||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
σ |
|
|
2π |
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
Экспоненци- |
|
|
τ |
|
|
0 |
|
|
τ |
|
|
|
|
|
|
|
τ e τz |
|
|
|
|||||||||||||||
альное |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
Пуассона |
|
|
|
|
1 |
|
|
|
|
|
|
k |
|
1 |
|
|
|
|
|
1 |
z |
k |
e |
z |
|
||||||||||
|
|
k! |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
k! |
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
1 |
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
k 1 |
|
z |
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
C |
|
|
|
|||||||||||||
Гамма |
|
|
|
|
|
|
|
|
|
|
|
k 1 |
|
|
|
|
|
|
|
z |
|
|
|
|
e |
|
|
|
|
|
|
||||
(k 1)!C |
k |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||
|
C |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||
|
|
|
|
|
(k 1)!C |
k |
|
||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
Исследованию корней характеристических уравнений посвящено много работ, в которых получены теоремы сравнения для практически важных случаев при переменных значениях коэффициентов и запаздываний. Эти результаты позволяют выделить в пространстве параметров области, где сохраняется характер решений уравнения (3.1).
Умножим обе части уравнения (3.3) на zezt . Тогда
|
|
|
z H(z) ezt |
|
bp |
zp m ez(t τ) |
|||
|
|
|
|
|
|||||
|
|
|
|
|
|
am |
|
||
dy(t) |
bp |
|
|
или |
|
||||
zp m ez(t τ) |
z ezτ H(z) ez(t τ) |
z ezτ y(t τ)(3.5) |
|||||||
|
|
|
|||||||
|
am |
||||||||
dt |
|
|
|
|
|
|
|||
Таким образом, сигнал на выходе системы определяется линейным дифференциальным уравнением с запаздывающим аргу-
