close

Вход

Забыли?

вход по аккаунту

?

Аппроксимация законов распределения ортогональными полиномами.

код для вставкиСкачать
УДК 681.518.3, 514:681.323/043.3/
С.А. Прохоров, И.А. Лезин, И.В. Солдатова
АППРОКСИМАЦИЯ ЗАКОНОВ РАСПРЕДЕЛЕНИЯ ОРТОГОНАЛЬНЫМИ
ПОЛИНОМАМИ
Приводится методика и алгоритмы аппроксимативного анализа законов распределения случайных
процессов ортогональными полиномами Лежандра, Чебышева, Лагерра Эрмита, положенными в
основу разработанной автоматизированной системы
При обработке статистического материала часто приходится решать вопрос о том, как подобрать для данного статистического ряда теоретическую кривую распределения, выражающую лишь существенные черты статистического материала. Такая задача, относящаяся к классу аппроксимативных задач, в данном случае называется задачей выравнивания (сглаживания)
статистических рядов [1].
Авторами разработана автоматизированная систеНачало
ма аппроксимативного анализа плотностей распределения вероятностей, являющаяся дальнейшим развитием разработанных ранее систем, используемая как при
подготовке специалистов в области обработки информации и управления, так и при проведении научных
Ввод исходных данных
исследований [2, 3].
Система включает в себя следующие подсистемы:
задания
входных воздействий;
Построение гистограммы,
генерирования временных рядов с заданным закооценка статистических
ном распределения;
характеристик
первичной статистической обработки временных
рядов;
идентификации законов распределения по фазоИдентификация законов
вым портретам;
распределения
аппроксимации плотностей распределения ортогональными полиномами Лежандра, Чебышева, Лагерра, Эрмита [4];
сплайн-интерполяции;
Выбор метода аппроксиработы с базой данных.
мации, критерия приближения
В основу работы системы положены как традиционные алгоритмы генерирования временных рядов с
заданными законами распределения, первичной обработки информации, связанной с оценкой числовых хаСоставление системы
рактеристик, построения гистограммы, так и оригиуравнений (уравнения)
нальные алгоритмы аппроксимации плотностей распределения вероятностей классическими ортогональными полиномами.
Выбор численного метода
Последовательность этапов решения задачи сглаи решение системы уравживания представлена на рис. 1 [2].
нений
После ввода исходных данных важным и не формализованным шагом решения этой задачи является
построение гистограммы наблюдаемого статистичеN
ского ряда {xi }i =1 . Для этого необходимо выполнить
Критерии согласия
нижеследующие этапы.
1. Определить диапазон изменения статистического ряда x min - x max .
2. При неизвестном законе распределения опредеКонец
лить ширину дифференциального коридора:
x - x min
Р и с. 1. Блок-схема сглаживания
D x = max
,
(1)
M
статистических рядов
где М - количество дифференциальных коридоров [1].
128
Из выражения (1) следует, что значение величины M не определено и, следовательно, исследователь должен задать её самостоятельно, ориентируясь на рекомендации специалистов.
Как следует из литературных источников, значение М рекомендуется выбирать в диапазоне от
10 до 20 [1].
При известном или предполагаемом законе распределения при линейной интерполяции
функции распределения величина Dx определяется выражением
8d
Dx=
,
(2)
F ¢¢(x ) max
где d - погрешность аппроксимации функции распределения; F ¢¢(x ) max - максимум модуля вто-
рой производной функции распределения [2];
é x - x min ù
M = ent ê max
(3)
ú.
Dx
ë
û
3. Определить частоту попадания анализируемой случайной величины в j-тый дифференциальный коридор:
^
1 N
p j = ådi j ,
(4)
N i =1
где
ì
é xi - x min ù
ï 1, если ent ê
ú + 1 = j Ù xi = x max ;
ë Dx û
ï
ï1
(5)
d i j = í , если xi = j D x ;
ï2
ï 0 , иначе ,
ï
î
- индикатор состояния.
Следует отметить, что di, j+1=1/2, если xi =jDx Ù x¹ xmax, т.е. в этом случае в j и j+1 коридоры
добавляется по 1/2.
4. Если частота попадания в какой-либо k-тый дифференциальный коридор мала
(pj<0,01¸0,02), то для уменьшения влияния случайности его объединяют с k+1 коридором. Эта
операция может быть применена неоднократно.
Исходя из условия нормировки
M ^
h j × Dx = 1 ,
å
j =1
(6)
вычисляем высоту столбцов:
^
^
hj =
pj
Dx
.
(7)
Далее на интервале [x min , x max ] строится линейная сплайн-модель. Для этого введем массив
{
}
^
пар x j , y j , j = 1, M , где xj – середина j-того интервала гистограммы, а y j = h j – соответственно высота j-того столбца. Значения сплайн-модели рассчитываются следующим образом:
ì y1 + ( y 2 - y1 ) / ( x2 - x1 ) × (x - x1 ), x < x1 ,
ï
f s , M (x ) = í y j + y j +1 - y j / x j +1 - x j × x - x j , x j < x < x j +1 ,
(8)
ï
î y M -1 + ( y M - y M -1 ) / ( x M - x M -1 ) × ( x - x M -1 ), x > x M .
Здесь f s, M (x ) – линейный сплайн, построенный на M дифференциальных коридорах, который
(
)(
)(
)
представляет собой модель аппроксимируемой плотности вероятности.
Очевидно, что внешний вид сплайн-модели сильно зависит от количества дифференциальных коридоров, причем утверждать однозначно, при каком именно значении M сплайн-модель
наилучшим образом соответствует истинной плотности вероятности, невозможно. Поэтому зачастую целесообразнее применять механизм “усреднения” модели плотности вероятности по
нескольким сплайн-моделям.
129
С использованием (8) для расчета сплайн-модели плотности вероятности усредненная модель определяется по формуле
M +l
f s , M (x ) =
å f (x )
s, m
m = M -l
2l + 1
.
(9)
0,5
1
0,4
0,8
0,3
0,6
f(x)
f(x)
В данном случае в точке x суммируются значения 2l+1 сплайн-моделей, построенных на
одних и тех же исходных выборках для различного числа коридоров M (от M-l до M+l), затем
сумма делится на количество просуммированных значений 2l+1, тем самым определяется
среднее значение для указанного ряда сплайнов.
Следующим этапом решения задачи сглаживания является идентификация законов распределения [2]. В том случае, когда закон распределения в виде функций определенного вида с
неизвестными параметрами подобрать не удается, единственным способом является применение ортогональных разложений в том или ином базисе. Отметим, что плотность распределения
вероятности может находиться в следующих диапазонах (- ¥ , ¥ ) , (- ¥ , a ] , [b , ¥ ) (см. рис. 2):
0,2
0,1
0,4
0,2
0
0
-4
-3
-2
-1
0
1
2
3
4
0
1
2
x
3
4
5
x
а
б
Р и с. 2. Типовые законы распределения:
а- нормальный закон; б- экспоненциальный закон
Таким образом, при построении модели плотности распределения вероятности необходимы ортогональные функции или полиномы, ортогональные как на конечном, так и на бесконечном (полубесконечном) интервалах: Лежандра, Чебышева, Лагерра, Эрмита [4].
Рассмотрим возможность аппроксимации плотности вероятности произвольного вида ортогональными полиномами. Дана функция вида f ( x ) , определенная на интервале [L , R ] . Она
может быть разложена в абсолютно сходящийся ряд вида [4]:
^
f (x ) =
¥
åb
k
× y k (x ) ,
(10)
k =0
где b k – коэффициенты Фурье, y k ( x ) – семейство базисных функций, ортогональных на определенном интервале с весом m (x ) . Это семейство характеризуется интегралом:
ìï y 2 , k = n,
y k (x )y n (x )m (x )dx = í n
ïî0, k ¹ n,
L
R
ò
(11)
где [L , R ] – интервал, на котором выполняется условие ортогональности указанных функций.
Следует подчеркнуть, что на практике приходится ограничиваться конечным числом членов ряда (10). Это приводит к появлению методической погрешности, значение которой зависит в немалой степени от способа оценки параметров модели.
Поэтому для модели аппроксимирующей функции
^
f (x ) =
m
åb
k
× y k (x ) ,
(12)
k =0
имеющей ограниченное число параметров, коэффициенты разложения, обеспечивающие минимум среднеквадратической погрешности аппроксимации:
130
R
2
n
é
ù
D = ê f (x ) b k ×y k (x )ú m (x )dx = min ,
k =0
û
Lë
определяются по формуле [3]
å
ò
(13)
R
b k = ò f (x )y k (x )m (x )dx .
(14)
L
При таком способе определения коэффициентов разложения погрешность аппроксимации, с
учетом свойств ортогональных полиномов, равна
R
D=
ò
f 2 (x )m ( x )dx -
m
åb
2
k
.
(15)
k =0
L
Так как выбранная система ортогональных полиномов определена на конкретном интервале, а интервал существования аппроксимируемой функции [x min , xmax ] произволен, то они могут
не совпадать. Поэтому приходится вводить коэффициенты линейного переноса, и формула (14)
примет вид
x max
1
f (x )y k ((x - b ) / a )m ((x - b )/ a )dx .
bk =
ax
ò
(16)
min
Примеры аппроксимации различными системами полиномов и значения коэффициентов
линейного переноса a и b представлены в табл. 1.
При численном нахождении значения интеграла (16) в силу того, что на границах интервала аппроксимации трудно добиться хорошего приближения, имеет смысл предварительно преобразовать аппроксимируемую функцию так, чтобы ее значение на одном из концов интервала
(или сразу на обоих) стало равно нулю. В таком случае удается точнее посчитать коэффициенты и избежать на концах аппроксимируемой функции явления Гиббса [5].
Преобразование осуществляется по формуле
(17)
f (0 ) ( x ) = f ( x ) - d × x - c .
Значения коэффициентов d и c рассчитываются в зависимости от того, сколько концов функции
требуется свести к нулю. Если к нулю сводится левый конец, то коэффициенты рассчитываются по следующим формулам:
ìd = 0,
(18)
í
îc = f (xmin ).
Если сводится только правый конец, то
ìd = 0,
í
îc = f (xmax ).
Для обоих концов сразу коэффициенты принимают значения:
ìd = ( f ( xmax ) - f (xmin )) (x max - xmin ),
í
îc = f (xmin ) - d × xmin .
Далее вычисляются коэффициенты разложения для функции f (0 ) (x ) :
(0 )
bk
1
=
a
(19)
(20)
x max
( )
ò f (x )y ((x - b )/ a )m ((x - b)/ a )dx .
0
k
(21)
x min
Полученные коэффициенты необходимо пересчитать в b k , которые используются в выражении (12). Массивы коэффициентов отличаются друг от друга только первыми двумя элементами. Примеры расчета коэффициентов для различных ортогональных систем представлены в
табл. 2.
Следующим приемом аппроксимации плотности распределения вероятностей, позволяющим уменьшить погрешность, является двухсторонняя аппроксимация. Суть этого метода заключается в том, что какая-то определенная точка x (например, точка экстремума) разбивает
интервал аппроксимации функции [xmin , x max ] на два независимых интервала [x min , x ] и [x , xmax ] .
На каждом из них функцию аппроксимируют отдельно, что в ряде случаев позволяет значи-
131
тельно улучшить результат, уменьшив погрешность и сократив количество членов в разложе^
^
нии функции f ( x ) . Тогда функция f ( x ) записывается с использованием единичной функции:
ì0, x < 0,
1(x ) = í
î1, x ³ 0.
(22)
Таблица1
Примеры расчета коэффициентов
Полиномы
Лежандра
[- 1;1]
Основные соотношения
m(x) = 1 ;
Pk (x ) =
s
å as × x k -2 s ,
2
×
2
a s = (- 1) ×
s
p
x max - x min
,
2
+ x min
x
b = max
2
a=
(2k - 2s)!
2k × s!×(k - s)!×(k - 2s)!
1
;
1 - x2
1
T0 (x ) =
Tk (x ) =
min
s =0
m (x ) =
[- 1;1]
x
1 max
b k = ò f (x )Pk ((x - b) / a )dx ,
ax
k div 2
2×k +1
as = (- 1) ×
Чебышева
Коэффициенты
p
x
1 max
bk = ò f ( x)Tk (( x -b) / a) 1-((x -b) / a)2 dx
ax
,
min
,
k div 2
å as × x k -2 s
s =0
× 2 k - 2 s -1 ×
- x min
x
a = max
,
2
x
+ x min
b = max
2
;
k (k - s - 1)!
s!×(k - 2s )!
m (x ) = e - x
L0 ( x ) = 1,
Лагерра
Lk ( x ) =
[ 0; ¥ )
k
Эрмита
(- ¥ ;+¥ )
a = 1,
k!
b = x min
(s!)2 × (k - s )!
å
s =0
as × x k - 2s ,
k!
2k × p
×
-(x -b) / a dx
,
s =0
k div 2
a s = (- 1)s
ò f (x)Lk ((x - b) / a)e
xmin
å as × x s ;
a s = (- 1)s ×
H k (x ) =
bk =
xmax
bk =
xmax
ò
f (x)Hk ((x -b) / a)e-(( x-b) / a) dx,
2
xmin
2 k - 2s
s!×(k - 2 s )!
a = 1,
b = (x max + x min ) 2
Выражение (12) с учетом (22) представим в виде
mЛ
^
f (x ) =
åb
k ,Л
×y k , Л ((x - bЛ )/ a Л ) ×1(x - x ) +
k =0
+
mП
åb
k ,П
(23)
×y k , П ((x - bП ) / a П ) × 1(x - x ).
i =0
Так как границы интегрирования изменились, то теперь коэффициенты рассчитываются по
модифицированным формулам:
b k ,Л
132
1
=
aЛ
x
ò f (x )y ((x - b )/ a )m ((x - b )/ a )dx ;
k ,Л
x min
Л
Л
Л
Л
Л
(24)
b k ,П =
1
aП
x max
ò f (x )y ((x - b )/ a )m ((x - b )/ a )dx
k ,П
П
П
П
П
П
(25)
x
для левой и правой ветви соответственно.
Таблица 2
Пересчет коэффициентов при сведении к нулю
Полиномы
Пересчет коэффициентов
Лежандра
b 1 = b 1(0 ) + d × a ×
b 0 = b 0(0 ) + (c + d × b) × 2 ,
2
,
3
b k = b k(0 ) , k = 2,...,m
b 0 = b 0(0 ) + (c + d × b ) × p ,
b 1 = b 1(0 ) + d a ×
Чебышева
p
2
,
b k = b k(0 ) , k = 2 ,...,m
b 0 = b 0(0 ) + c + d × (b + a ),
b 1 = b 1(0 ) - d × c ,
Лагерра
b k = b k(0 ) , k = 2 ,...,m
b 0 = b 0(0 ) + (c + d × b) × 4 p ,
b 1 = b 1(0 ) + d × a ×
Эрмита
p
2
,
b k = b k(0 ) , k = 2 ,..., m
Коэффициенты линейного переноса также изменяются соответственно изменению интервалов интегрирования.
Зачастую при двусторонней аппроксимации используют операцию совмещения максимумов. Смысл данной операции состоит в том, что после вычисления коэффициентов для выражения (23) их нормируют, то есть совмещают в точке x значения аппроксимирующих выражений для обеих ветвей со значением аппроксимируемой функции. Этот подход дает возможность использовать комбинированные модели, т.е. применять для аппроксимации левой и правой ветвей различные ортогональные полиномы, соответствующие наилучшим образом области существования плотности распределения вероятностей.
Для этого требуется выполнение условий
mЛ
ì^
f
x
(
)
=
b k , Л ×y k , Л ((x - b Л ) a Л );
ï
ï
k =0
(26)
í
mП
ï^
b k , П ×y k , П ((x - b П ) a П ).
ï f (x ) =
k =0
î
Так как при конечных значениях mЛ и mП условие (26) не выполняется, для его обеспече-
å
å
^
ния необходимо искать аналитическое выражение для f ( x ) в виде
^
f (x) =
mЛ
å
k =0
bk , Л ×y k , Л ((x - b Л ) a Л ) × 1(x - x ) +
mП
åb
k,П
×y k , П ((x - b П ) a П ) × 1(x - x ).
(27)
k =0
Для того, чтобы квадратическая погрешность аппроксимации была минимальной при дополнительных условиях (26), необходимо минимизировать ΔЛ и ΔП по bn,Л и bn,П соответственно:
133
2
x
mЛ
mЛ
ì
é
ù
ïD Л =
bk , Л ×y k , Л ((x - b Л ) / a Л ) ;
ê f (x ) - bk , Л ×y k , Л ((x - b Л ) / a Л )ú dx + l Л
ïï
ê
k =0
k =0
ûú
x min ë
í
2
x max
mЛ
mП
ï
é
ù
bk , П ×y k , П ((x - b П ) / a Л )ú dx + l П
bk , П ×y k , П ((x - b П ) / a П ) .
ïD П =
ê f (x ) úû
k =0
k =0
ë
x ê
îï
ò
ò
å
å
å
(28)
å
Найдем частные производные
¶D Л ¶D П
,
и приравняем их к нулю:
¶bn, Л ¶bn, П
x
mЛ
ì ¶D
é
ù
ï Л = -2 ê f (x ) - bk ,Л ×y k , Л ((x - b Л ) / a Л )úy k , Л ((x - b Л ) / a Л )dx +
ï ¶bn , Л
ê
k =0
ûú
x min ë
ï
+ l Л y k , Л ((x - b Л ) / a Л ) = 0;
ï
(29)
í
x max
mП
é
ù
ï ¶D П
= -2 ê f (x ) - bk , П ×y k , П (( x - b П ) / a Л )úy k , П ((x - b П ) / a Л )dx +
ï
ê
k =0
ï ¶bn , П
ûú
x ë
ï
+ l Пy k , П (( x - b П ) / a П ).
î
Выполнив промежуточные преобразования, получим
lЛ
ì
ïïbn , Л = b n , Л - 2 y n , Л (( x - b Л ) a Л );
(30)
í
lП
ïb = b
y
(
(
)
)
.
x
b
a
n, П
n, П
n, П
П
П
2
îï
С учетом (26) определим
ì
ö
æ mЛ
2ç b k ,Л ×y k , Л ((x - b Л ) a Л ) - f (x )÷
ï
÷
ç
ï
ø;
è k =0
l
=
Л
ï
mЛ
ï
y k2, Л ((x - b Л ) a Л )
ïï
k =0
(31)
í
mП
ö
æ
ï
2ç b k ,П ×y k , П ((x - b П ) a П ) - f (x )÷
ï
÷
ç
ø.
ïl = è k =0
П
mП
ï
ï
y k2, П ((x - b П ) a П )
ïî
k =0
Подставив это выражение в выражение для оценки коэффициентов разложения (30), окончательно получим
mЛ
ì
f (x ) b k , Л ×y k , Л ((x - b Л ) a Л )
ï
k =0
ïb = b
+
y n , Л ((x - b Л ) a Л );
n, Л
mЛ
ï n, Л
2
ï
y k , Л ((x - b Л ) a Л )
ï
k =0
(32)
í
mП
ï
f (x ) b k , П ×y k , П ((x - b П ) a П )
ï
k =0
ïbn , П = b n , П +
y n , П ((x - b П ) a П ).
mП
ï
2
y k , П ((x - b П ) a П )
ï
k =0
î
В табл, 3 приведены значения коэффициентов bn,Л и bn,П для различных систем ортогональных полиномов.
Так как полиномы являются степенными функциями, то от представления (12) можно перейти к степенному представлению аппроксимирующей функции, пересчитав коэффициенты.
В таком случае общий вид аппроксимирующей функции выглядит так:
ò
å
ò
å
å
å
å
å
å
å
å
å
^
f (x) =
m
åb
k =0
134
× ((x - b )/ a ) , x Î [xmin , xmax ] .
k
k, f
(33)
В табл. 4 приведены формулы расчета коэффициентов b k , f в формуле (32) для различных
систем полиномов. Применяемые в табл. 4 обозначения a s ,s + 2 m – это коэффициенты as полинома степени s+2m, используемые в табл. 1.
Таблица3
Коэффициенты разложения при двусторонней аппроксимации
Полиномы
Коэффициенты
f (x ) bn , Л = b n , Л +
mЛ
åb
× Pk (1)
k ,Л
k =0
mЛ
å (1)
Pn (1),
Pk2
k =0
Лежандра
f (x ) bn , П = b n , П +
mП
åb
k ,П
× Pk (- 1)
k =0
mП
å (- 1)
Pn (- 1)
Pk2
k =0
bn , Л = b n , Л
æ 1
2
b +
f (x ) - ç
ç p 0, Л
p
è
+
1 + 2 × mЛ
mЛ
åb
k ,Л
k =1
ö
÷
÷
ø T (1),
n
p
Чебышева
bn , П = b n , П
m
ö
æ 1
2 П
k
b 0, П +
(
- 1) × b k , П ÷
f (x ) - ç
÷
ç p
p k =1
ø T (- 1)
è
+
n
1 + 2 × mП
å
p
f (x ) bn , Л = b n , Л +
mЛ
åb
k ,Л
k =0
,
mЛ + 1
Лагерра
f (x ) bn , П = b n , П +
mП
åb
k ,П
k =0
mП + 1
Таблица4
Коэффициенты степенного ряда
Полиномы
Лежандра,
Чебышева,
Эрмита
bk , f =
Коэффициенты
(m -k )div 2
åb
k +2 s
× a s,k +2 s , k = 0,..., m
s =0
n
b k , f = å b s × a k ,s , k = 0,..., m
Лагерра
s= k
Заметим, что по полученному аналитическому выражению плотности распределения вероятностей легко находится аналитическое выражение и функции распределения:
ì0, x £ xmin ,
ï x
^
ï ^
F ( x ) = í f ((x - b ) / a )dx, xmin < x < xmax ,
(34)
ï xmin
ï
î1, x ³ xmax .
ò
135
Коэффициенты восстановленной функции распределения b i,F определяются следующим
образом:
ì b i , F = b i -1, f × a / i, i = 1, n + 1;
ï
n +1
(35)
í
i
ï b 0 ,F = - b i , F × ((x min - b ) / a ) .
i =1
î
Соответственно выражение для функции распределения с учетом (35) запишется в виде
å
^
F (x ) =
n +1
å b ((x - b) a ) .
i
(36)
i ,F
i= 0
На рис. 3 приведен пример реализации случайного процесса – оборотов газотурбинного
двигателя при пуске. В качестве примера для выборки, приведенной на рис. 3, построена аппроксимационная аналитическая модель плотности вероятностей, график которой представлен
на рис. 4.
5975
Обороты
5970
5965
5960
5955
5950
0
50
100
150
200
250
300
350
400
450
500
550
600
Отсчеты
Р и с. 3. Обороты газотурбинного двигателя при пуске
0,4
f(x)
0,3
0,2
0,1
0
5950
5953
5955
5958
5960
5963
5965
5967
5970
5972
5975
x
Р и с. 4. Аппроксимированная плотность вероятностей
Таким образом, в настоящей работе приведена и реализована методика построения аналитических выражений для законов распределения случайных процессов, заданных графически.
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1.
2.
3.
4.
5.
136
Вентцель Е.С. Теория вероятностей. М.: Наука, 1969. 576 с.
Прохоров С.А. Аппроксимативный анализ случайных процессов / 2-е изд., перераб. и доп. Самара: Самар. гос.
аэрокосм. ун-т. 2001. 380 с.
Прохоров С.А. Моделирование и анализ случайных процессов. Лабораторный практикум / 2-е изд., перераб. и
доп. Самара: СНЦ РАН, 2001. 380 с.
Д. Джексон. Ряды Фурье и ортогональные полиномы. М.: ГИИЛ, 1948. 260 с.
Андре Анго. Математика для электро- и радиоинженеров. М.: Наука, 1965. 780 с.
Поступила 19.05.2004 г.
Документ
Категория
Без категории
Просмотров
9
Размер файла
249 Кб
Теги
законов, полиномами, распределение, ортогональных, аппроксимация
1/--страниц
Пожаловаться на содержимое документа