Научная статья на тему 'Об эффективности оценок при параметрической идентификации нелинейных статических объектов'

Об эффективности оценок при параметрической идентификации нелинейных статических объектов Текст научной статьи по специальности «Математика»

CC BY
60
13
i Надоели баннеры? Вы всегда можете отключить рекламу.

Аннотация научной статьи по математике, автор научной работы — О А. Кацюба, М Б. Линеенко

Рассматривается задача параметрического оценивания стохастических статических моделей, нелинейных относительно параметров. В качестве метода оценивания применяется метод квазиправдоподобных оценок; дано доказательство состоятельности и асимптотической нормальности оценок максимального квазиправдоподобия. Показано, что применение линейно-комбинированной оценки позволяет получить эффективность оценивания как угодно близкую к эффективности оценок максимального правдоподобия в условиях отсутствия априорной информации о законах распределения помех наблюдений.

i Надоели баннеры? Вы всегда можете отключить рекламу.
iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

ABOUT AN EFFICIENCY OF ESTIMATIONS WITH REFERENCE TO PARAMETRICAL IDENTIFICATION OF NONLINEAR STATIC OBJECTS

It is considered the problem of a parametrical estimation of stochastic static models, nonlinear concerning parameters. The method of quasi-likelihood estimates is applied as a method of estimation; it is adduced proofs of a consistency and of an asymptotic normality of maximum quasi-likelihood estimates in conditions of independent heterogeneous observations. It is shown, that application of linearly-combined estimations allows receiving effectiveness of estimation somehow close to effectiveness of maximum likelihood estimations in conditions of lack of the information of distribution functions of parasites of observations.

Текст научной работы на тему «Об эффективности оценок при параметрической идентификации нелинейных статических объектов»

УДК 519.254

ОБ ЭФФЕКТИВНОСТИ ОЦЕНОК ПРИ ПАРАМЕТРИЧЕСКОЙ ИДЕНТИФИКАЦИИ НЕЛИНЕЙНЫХ СТАТИЧЕСКИХ ОБЪЕКТОВ

©2004 О. А. Кацюба, М. Б. Линеенко

Самарская государственная академия путей сообщения

Рассматривается задача параметрического оценивания стохастических статических моделей, нелинейных относительно параметров. В качестве метода оценивания применяется метод квазиправдоподобных оценок; дано доказательство состоятельности и асимптотической нормальности оценок максимального квазиправдоподобия. Показано, что применение линейно-комбинированной оценки позволяет получить эффективность оценивания как угодно близкую к эффективности оценок максимального правдоподобия в условиях отсутствия априорной информации о законах распределения помех наблюдений.

Пусть имеет место стандартная задача несмещенного оценивания регрессионных

функций; наблюдаемые величины у,

(i = 1, n) представляется в виде

yt = u (xt, a0) + s(i), где u( xt, a0) - функция,

вид которой известен, dim a0 = p, s(i) - независимые случайные величины, удовлетворяющие условиям M (s(i)) = 0,

M(s 2(i)) < a 2, (m - оператор математического ожидания).

Требуется определить состоятельные оценки a(n), не зная априорно закона распределения шума

f,(s(i)) = f, (yi, a0) = f (У, " u(X, a0)).

В условиях априорной неопределенности существуют методы состоятельного оценивания (метод эмпирического риска, М-оценка и т.д.), однако очевидно, что эффективность этих оценок в общем случае далека от эффективности оценки максимального правдоподобия. В связи с этим возникает следующая задача: создание методов и алгоритмов определения таких состоятельных несмещенных оценок на основе множества оценок (линейно-комбинированных оценок), эффективность которых в широком смысле была бы в принципе достаточно близка к эффективности оценок максимального правдоподобия. Наиболее просто вопрос получения множества состоятельных оценок, из которых мо-

жет быть сконструирована достаточно эффективная оценка, решается, если воспользоваться методом эмпирического риска, когда эмпирический функционал конструируется на основе некоторых законов распределения, и тогда указанный метод можно трактовать как непосредственное расширение метода максимального правдоподобия на случай априорной неопределенности при модификации некоторых его условий (такие оценки названы квазиправдоподобными), что позволяет получить более конструктивные условия состоятельности и асимптотической нормальности оценок, причем, этот метод распространен на случай неоднородных наблюдений, а также на задачу оценивания параметров линейных разностных уравнений при наличии помех в выходных переменных.

Отдельные вопросы состоятельности квазиправдоподобных оценок и теории линейно-комбинированных оценок рассмотрены в [1] и [2]. Предлагаемая статья дает возможность рассмотреть вопросы эффективности оценивания параметров нелинейных моделей в целом.

Пусть эксперименты {х(п), Л(и п), Ра(п)} порождаются наблюдениями У = (у^,...,уп)т со значениями в {х, Ли} и распределениями

{ра}, пусть семейство р(п) доминируется некоторой мерой и и существует функция

cFg(n)(yi,...,yn; a) = 5F(n)(YT, a) сьо^^у) du

= f(n) (YT, a).

Определение 1. Функцией квазиправдоподобия параметров а, отвечающей

{х(п), А(ип),F{")} и у, называется неотрицательная вещественная функция

ф(п)(УТ,а) * /{п)(УГ,а), а е А,

где

фп (уТ, а) = С..р„ (У " ^ а\...Уп - и(хп, а)).

Оценкой максимального квазиправдоподобия (МКП) а0 для заданной функции квазиправдоподобия ф(П)(УТ, а) по наблюдениям у называется а(п, ф), определяемая из

ф(п) (уТ , ^ ф)) = ^Р ф(п) (уТ , а) .

аеА

В частности, если

ф( п)(УТ, а) = /(п)(УТ, а), то оценка будет правдоподобной.

Оценку, при которой функция ф(п) (У Т, а) принимает абсолютно наибольшее значение на А, назовем оценкой МКП в строгом смысле. Пусть функция фг. (у., а) имеет частные

производные по всем аргументам ( j = 1, p ), тогда функция

n

ln ф( n)(XT, a) = jln ф,(yt, a)

,(j)

,=1

j дln ф,(yi, a) = 0

(1)

на оценке МКП в строгом смысле. В дальнейшем будем предполагать, что для всех п, за исключением конечного числа значений, существует самое большее одна оценка МКП.

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

Известная ограниченность доказательства состоятельности, основанного на "локальных" суждениях, не связанных с понятием точной верхней грани функции МКП заключается в том, что, во-первых, имеет место требование единственности корня, во-вторых, требуется выполнение довольно жестких ограничений на поведение функции квазиправдоподобия (существование производных и т.д.), поэтому далее рассмотрено доказательство состоятельности оценок МКП, не требующее существования производных и основанное на понятии точной верхней грани.

Введем следующие обозначения: $г (а) - замкнутый шар в А с центром в точке а и радиусом г ;

ф* (У., а г) = эир ф. (у., а)

( а) '

ф* Ч У,,= suP ф, (y,, a)

и> 0;

(в случае независимых случайных величин ) также имеет частные производные. Если

для некоторого а е А эта функция достигает максимума, то в этой точке

ф*( У,, ar) =

(ф* (y,, a г х ф* (y,, a г) >1

.=1 даи) Система уравнений (1) называется системой уравнения МКП. Ясно, что всякая оценка МКП в строгом смысле удовлетворяет системе уравнений МКП, поэтому любой корень системы уравнений (1), лежащий в А, называется оценкой МКП в широком смысле. Из нижеизложенного следует, что оценка МКП, как в строгом смысле, так и в широком смысле не единственна, причем оценка МКП в широком смысле не обязательно тождествен-

[1, впротивном случае

Утверждение 1. Пусть выполняются следующие условия:

10. При каждом уТ функция ф(п)(УТ,a)

непрерывна относительно a в A (A - замкнутое множество).

20. Вектора xi - неслучайные. 30. Для каждого малого r > 0 выполняются M [log ф * (y,, a, r) J - конечно (2)

и

1 ¿M [log ф*( y,, a, r )J

n i=1

сходится к

lim — ¿M [log ф*( yt, a, r )J

n—^ro n ^^ , =1

равномерно отно-

сительно r .

40. Для каждого a0 е A существует шись (4), получаем

1 п

lim -YM[logф,(вг)].

и^-вд n ^Т " г=1

50. Для всех у и каждого a0 е A , a ^ a0

lim lim - Y M[log ф* (у,, a, г)] =

г^-0 и^-вд и ^^ " 1=1

lim- Y M [log ф, (у,, a)]

<

< lim- Y M [log ф, (в,)]; " ,=1

lim lim - ¿M[log ф** (у,, и)]

= —вд

= lim1 ¿M [log ф* (у,, a, r )] = " i=1

= lim1 ¿M[log ф, (у,, a)].

" 1=1

Далее выберем и = и 0 так, чтобы

Um1 ¿M |log ф* *( yt, Uo)]<

я^-да и ^^ " i=1

< lim1 ¿M[log ф, (уг, a)].

(5)

(6)

60. При каждом a е A и r , и для дис-

i =1

персии величины logф (V-,a,r) , logф,(s,), „ г0

от^./^ /з o-i^v Существование такого вытекает из 50.

log ф**(yt, и) выполняются условия теоремы Обозначим через A" е A' множество, состо-Маркова [3, с.209].

ящее из всех точек a из A, для которых вы-

Тогда для любого замкнутого множества полняется неравенство ||а|| < и0, тогда в со-

А'е А и ао ^ А' и любого С > 0, 5 > 0 име- ответствии с (5), (6) для каждой точки а е А" ет место

q

a0

SUP ф(и)(уГ , a)

aeA'

ф( и)(УГ, a0)

> - — 5

поставим в соответствие значение г , что

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

если п > п(С, 5), где qa - значение вероятности.

Доказательство утверждения 1. Из предположения 10 следует, что

М logф*(у.,a,г) = logф*(у,,а) ; (3)

lim - YM [log ф* (у,, a, г)] <

и^-вд и ^^

< lim-¿M[logф, (в,)].

(7)

Множество A'' - компактно, поэтому [6, с.!93] в найдется конечное число точек

ф,*( У г, a) =

Гф, (у,, a), если ф, (у,, a) > - • ,

f-1

что A" ^ U (k)(ak ). Очевидно,

k=

Ц впротивном случае.

Так как log ф* (у{, a, г) - невозрастающая функция г , тогда из (2) и (3) [4, с.202]

lim M[logф*(у,,a,г)] = M[logф*(у,,a)]. Аналогично доказательство и для случая Г- впротивном случае;

[ф, (у,, a) если ф, (у,, a) < 1,

тогда

lim M[log ф*(у,, a, г)] = M[log ф*(у,, a)]. (4)

Из п. 30 следует возможность почленного перехода к пределу [5] и, воспользовав-

что для каждого уТ имеем

ц и тТ '--

УТ, a )<2Пф* (у,, aj, гк ) +

aeA'

k=! i=-

ф*( у,, a) =

+ Пф*Чу, , U0).

,=-

Следовательно, достаточно доказать, что для любых Q, 5', 5''

(

q

a0

Пф*(у,, ak, гк)

ф (и)(УТ, a0)

Л

> - — 5'

, (8')

n > n(C,5',ak), k = 1, | ;

q

a0

Пф * *( у,-, u 0)

i=1

Ф(n)(YT,a0)

<С'

> 1 -5'

, (8'')

i=i

i=i

1 ^Г l1 n Г 1

- YMM* (у, ak, rk )J(у,, M;

,=1

,=1

Z1^* (y,, ak, rk) —Zlog^ (y,, ao -

-»—да •

Х^Гсу, «0) —Zlogi (y,, a0)

—^

откуда

л

Zlog)*(у,,ak,rk) —Ylogk (y,,a0) <lo§?

VV ,=1 ,=1

>

> 1 — 5';

(i

л

л

Y^gfe,u0) —(y,,a0) <logc'

,=1

>

п > п(С", 8м, и0), так как д не зависит от п .

Из условия теоремы и теоремы Слуцкого [7, с.282] следует, что

1 n 1 n

- Z1^* (у,, ak, rk)—Zlog^ (у,, ao -

ция а(п, ф) такая, что для всех Ут е %(п) выполняется равенство

^Р ф(п) (уТ , а) = ф(п) (уТ , ф)) .

аеЛ

Тогда {а(п, ф)} - состоятельная последовательность оценивающих функций для а е л .

Доказательство следствия 1. По определению а(п, ф)

ф^, a(n, ф)х

> 1

1 n 1 n

- (у,, u0)—Zlog^ (у,, оО -

n ,=1 n

- ¿М [logф; *( у,, ио)]--¿М [хоёф, (у,, ао)].

п ,=1 п ,=1

Однако по условиям (6), (7) эти величины отрицательны, поэтому

Ф(n)(YT,

для всех n и всех у1, ..., ун. Предположим,

что {a(n, ф)} несостоятельная последовательность, тогда можно выбрать последовательность целых чисел 0 < n1 < n2 <... и указать

два положительных числа С 0, 50, таких, что

qa, (¡|a(n,, ф—<С 0 )<1—50, , = 1,2,...,

т.е. имеют место такие множества х(n), во всех точках Y^ в которых с определенной вероятностью ||<a(n, ф) — a0|| >С 0, откуда следует YT е х(n):

sup _ ф(„,) (YT, a) > ф(И,) , a(n, ф))

a— a0 >s0,aeA

и

sup

ф (n )(YT, a)

> 1

||>?0,а^Л ф(п,)(Уп, , ао)

для неограниченно больших . Так как множество Л п {а: ||а - а01| > £0} - замкнуто, то по утверждению 1 множества %(п) должны иметь произвольную малую qao меру, если п1

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

> 1 -8'',

что доказывает (8') и (8'').

Аналогично доказывается неравенство

г достаточно велико, то это приводит к проти-

для любых k = 1, д . Так как д не зависит от воречию и доказательству утверждения. Тог-

п, то можно найти такое п*(8х, Сх), что од- да {а(п, ф)} - состоятельная п°след°ватель-

новременно выполняются все неравенства, из ность оценивающих функций для а е л . чего следует доказательство теоремы. Пусть имеется возможность получить

Следствие 1 Пусть выполняются все последовательность состоятельных оценок с

предположения теоремы и пусть для каждо- разной асимптотической эффективностью, а

го п > 1 определена на %(п) измеримая функ- для сравнения эффективности этих оценок с

эффективностью оценок максимального асимптотически несмещенной и асимптоти-правдоподобия необходимо знать дисперси- чески нормальной с ковариационной матри-

онные матрицы (хотя бы при п ^ да ).

Как известно, оценки максимального правдоподобия при определенных условиях регулярности [7] являются асимптотически эффективными оценками, поэтому, определив условия асимптотической нормальности оценки МКП, можно сравнить их эффективность с эффективностью оценки максимального правдоподобия.

Утверждение 2. Пусть выполняются следующие условия.

10. Все условия состоятельности (утверждение 1).

цей (при п ^ да )

Е М

,=1

^дМп^Чу.а)

да(т) да(1)

а = а,

> х

х Wn (а,)

М

' д 21п ф(г )(у1, а)

ЕЕ да(т)даи)

,=1

а = а,

20. При всех а е Л и любых п существуют и единственны при а = а0 первые и вто-

Доказательствоутверждения 2. Пусть для всех а е Л и каждого у{ существует положительная дважды дифференцируемая

рые производные 1п ф, (у,, а) и соответствен- , . Л , , Л

^ ^ функция б(а) и такая функция иф (у,), что

но матрицы н (УТ, а) =

1Е д 1п ф, (уг, а)

п Е да(1)

д

2 (

да{й) да(т)

б (а)

д 1п ф, (У,, а)

да

(у)

< иф(У,)

J (УТ, а) =

1Е д21п ф, (у,, а)

п Е да(т)да(1)

для всех d, т , у . Разложив функцию

1 б(а)Е д 1п ф,(у, ,а) Т ~

_ б (а)^ в ряд Тейлора относи-

п

1=1

да

М-

Е н (у,, а) Нт (у,, а)

г =1

а=а

= Ж.

тельно а = а0 , получаем

б (а)н (Ут, а) = б а )Н (Ут, а0) + + б (¿0) J (Ут, Ь0)~ +

у, т = 1, р , - положительно определенная матрица.

30. Для любых а е Л существует момент

+

1 ЕЕ дб(а0) д 1п ф,(уг, а0)

п ^ да(т) да(1)

а +

М-

д 1п ф,(у,, а)

да

а = а,

2+5

+1 Р(а, а0, Ут )а та , где Р(1 )(а, а0,Ут ) =

1

(9)

= "ЕР(1)(^,а,(у,)Р(1)(у,,а,«0^ < 1.

,=1

и

lim Е М-

,=1

д 1п ф, (у,, а) Ж ]-^ да

а=а

2+5

Подставим а(п, ф) в (9) и приравняем к 0:

б а )н (Ут, а0) + б а ) J (Ут, а0) X

= 0

х (а(п, ф) - а0) +

1 ЕЕ дб(а0) д 1п ф, (у,, а)

п^ да(т) да( 1

(условие Ляпунова), где 5 - любое положительное число, тогда оценка МКП является

1

х (а(п, ф) - а0) + — Р(а(п, ф), а0, У ) х

х

(а(п, ф) - а0)Т (а(п, ф) - а0), откуда

g а ) н (УТ, а0) = - ^ (а0) J (УТ, а0) +

+

1 Х дg (а0) д 1П ф, (У г, а)

п ¿1 да(т) да(1)

1

х (а(п, ф) - а0) + — Р(а(п, ф), а0, У ) х

х (а(п, ф) - а0 )Т (а(п, ф) - а0) = = (а0) J (УТ ,а0) +

+

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

1 Х дg (а0) д 1П фг (У г,а)

п ^ да(т) да(1)

+

+1 Р(а(п, ф), а0, УТ )(а(п, ф) - а0 )Т | х

х (а(п, ф) - а0) . Окончательно

(а(^ ф) - а 0) = (а0) J (уТ, а0) +

+

1 Х дg(а0) д 1П фг (Уг,а)

п Х да(т) да(1)

+

форме единичной матрицы, откуда V д 1П фг (У г ,а)

Х р имеет нормальный закон рас-

г=1 да

пределения с нулевым математическим ожиданием и дисперсионной матрицей Wn, далее,

1Х д 1пф,(у,,а) п г=1 да

нормальный закон распре-

Wn

деления с дисперсионной матрицей —^.

Таким образом, числитель выражения (10) распределен по р -нормальному закону с нулевым математическим ожиданием и дис-

Wn

персионной матрицей , откуда [7, с.281],

получаем, что распределение величины (10) - р -нормальное с нулевым математическим ожиданием и дисперсионной матрицей

1

х м

г=1

^д 21п )(Уг,а)

да(т) да(1)

г х

+ |р(а(п,ф),а0,уТ )(а(В,ф) -а0)Т | х

х g (а0) н (УТ ,а0).

По теореме Слуцкого 1 Р(а(п, ф), а0, УТ )(а(п, ф) - а0 )Т р х р).

Как уже было показано выше,

(10)

->0

->0

ный вектор

д

асимпто-

х Wn К)

м

С д21п ф(г)(Уг ,а)

да(т) да(1)

х

г=1

а = а,

1 у дg(а0) д 1П фг (Уг,а)

п у да(т) да(1)

Таким образом, знаменатель выражения (10) сходится к некоторой постоянной матрице -gа)J(УТа).

Исследуем числитель выражения (10): при выполнении условий Ляпунова случай" №п ]-у2 ХХ д 1П фг (Уг ,а)

тически нормален с нулевым математическим ожиданием и ковариационной матрицей в

Если же оценка правдоподобная, то эта матрица равна (при п ^го ) Wl-1(a0) (при условии регулярности в смысле вторых частных производных по а для каждого " г").

Пусть имеется возможность получить последовательность состоятельных оценок с разной асимптотической эффективностью, при этом возникает задача повышения эффективности оценок путем комбинации их.

Определение 2. Любая линейная комбинация конечного числа квазиправдоподобных

оценок а(п, ф(/)) (I = 1, k ) относительно разных функций квазиправдоподобия ф(/) с определенными весами называется линейно-комбинированной квазиправдоподобной а ъ.

Покажем, что определенных условиях,

х

1

накладываемых на веса линейно-комбинированной оценки, имеет место следующее неравенство: 1ейТак_х/а0 < 1ейТак/а0, причем,

если в последовательности {а(п, ф(1))} имеет

место а(п, ф(1)) = а(п, f) ( f (у,, а) - закон распределения помех наблюдений), то а к =а f (под значением 1еАТ а ка0 понимается отношение значений обобщенной дисперсии оценки максимального правдоподобия и исследуемой оценки).

Для последовательности квазиправдоподобных оценок имеет место утверждение 3.

Утверждение 3. Пусть имеет место последовательность состоятельных, асимптотически несмещенных оценок а(п, ф(1))

... а(п,ф(к)) , при всех а е А имеет место равенство

Г д 1п ф(1)(уг, а)

да(1)

-У (у,, а)Ф, = 0

1

0'

0'

т

01

1 0!... 0

0 1 !■■■ 0

"Г • 1 • 1

0 0!... 1

1' 0~ 0^ ... ¡у ?0Г" -1--- а1[п, ф (1)]

01

0" —'-1— 0' I ■•• -1---Т ! 1' а р [п, ф(1)]

будет состоятельной, асимптотически несмещенной с асимптотической эффективностью

1ейТ а ка0 < 1ейТ а к/а0 , причем, если имеет место а(п,ф(5)) = а(п,У) ( ^ = 1, к ), то а к = а(п, У), и 1ейТ а к/а0 = 1, где

а 1 [п, ф(1)] =

а( 1 )(п, ф(1)) а( 1 )(п, ф(2))

а( 1 )(п, ф(к))

0 =

соу(а(1) (п, ф(1)), а(1)(п, ф(5))) соу(й(2)(и, ф(1)), а(1){«,~ф(5))}

соу(а( р)(п, ф(1)), а(1)(п, ф(5))) 1

соу(а(1) (п, ф(1)), а( р)(п, ф)))

(1)

( р )

соу(а(2)(п, ф(1)), а( р)(п, ф(5)))

соу(0?( р )(п, ф(1)), а( р)(п, ф(5)))

и выполняется условие дифференцирования по параметру под знаком интеграла (в частном случае, когда ф(1) = имеет место условие регулярности в смысле вторых производных по ), то линейно-комбинированная оценка

- ковариационная матрица квазиправдоподобных оценок а(п, ф(1)) (положительно определенная, невырожденная, рк х рк ), порядок матрицы: 1: к х 1; 0: к х1; а 1 (п,ф(1)): к х1.

Следствие утверждения 3. Ковариационная матрица линейно-комбинированных оценок вычисляется по формуле

Дак] =

1'

0'

0'

0

Г

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

-1— 1 : 1 0!

01 0 . 1! 1 • 1

•! 1' 0 0!

Доказательство основных свойств полученных линейно-комбинированных оценок, исходя из свойств составляющих их квазиправдоподобных оценок, в том числе состоятельность, асимптотическая несмещенность, значение меры эффективности

1еАТ ак/а0 относительно наблюдений п приведено в [8].

Достаточно сложным вопросом является выбор функции ф для получения квазиправдоподобных оценок, при этом функции ф должны обладать свойствами, описанными в утверждении 2. Было бы желательно, чтобы последовательность |ф(1)} включала в

1

'

'

0

1

1

0

1

0

а к =

х

0

X

себя наиболее распространенные критерии, соответствующие законам распределения Лапласа, Гаусса, которые обладают определенными экстремальными свойствами при робастном оценивании параметров.

Одним из возможных примеров такой последовательности {ф(/) | является последовательность квазиправдоподобных оценок, соответствующая / -обобщенному нормальному закону

ФСу, X,J) =

1

J+1

2(2X2 )1J j

exp

^ У -u(x, ai

(2X2 ),

min

I

1

\yt - u (xt, a)|

(

ln ф(' )(YT, X, J) = ln

J + 1

~T

П^2 ).

Y

-I

1 syl - u(xt, a)|'

(2X2 ),

где г( +у1) - гамма-функция; 2Х2^ - величина, являющаяся функцией дисперсии; / -любое положительное число: при / = 1 - закон Лапласа, при / = 2 - закон Гаусса. Тогда в качестве формальных критериев при целых / > 0 для получения квазиправдоподобных оценок примем последовательно

и 1 ^

min V--— (у, - u (xf, a))

что соответствует методу наименьших модулей, квадратов и т.д., так как логарифмическая функция квазиправдоподобия имеет вид

На основе разработанных алгоритмов оценивания создано программное обеспечение, которое нашло применение при расчете прогноза концентрации вредных веществ в атмосфере.

СПИСОК ЛИТЕРАТУРЫ

1. Кацюба О. А. О методе квазиправдоподобных оценок в задачах идентификации нелинейных объектов // Сиб. отд-ние АН СССР. Автометрия. 1986. №6.

2. Кацюба О. А., Хакимов Б. Б. Алгоритм нелинейного параметрического оценивания в многомерных задачах статистической обработки // Сиб. отд-ние АН СССР. Автометрия. 1984. №2.

3. Гнеденко Б. В. Курс теории вероятностей. М., 2001.

4. Ширяев А. Н. Вероятность. М.: Наука. 1980.

5. Будак Б. М., Фомин С. В. Кратные интегралы и ряды. М.: Наука. 1965.

6. Александров П. С. Введение в теорию множеств и общую топологию. М.: Наука, 1977.

7. Крамер Г. Математические методы статистики. М.: Мир, 1975.

8. Кацюба О. А., ЛинеенкоМ. Б. Параметрическое оценивание нелинейных статических объектов в условиях априорной неопределенности // Труды III межд. конф. "Идентификация систем и задачи управления" SICPRO'04. М., 2004 г.

ABOUT AN EFFICIENCY OF ESTIMATIONS WITH REFERENCE TO PARAMETRICAL IDENTIFICATION OF NONLINEAR STATIC OBJECTS

©2004 O.A. Katsyuba, M.B. Lineenko

Samara State Academy of Ways of Communications, Samara

It is considered the problem of a parametrical estimation of stochastic static models, nonlinear concerning parameters. The method of quasi-likelihood estimates is applied as a method of estimation; it is adduced proofs of a consistency and of an asymptotic normality of maximum quasi-likelihood estimates in conditions of independent heterogeneous observations. It is shown, that application of linearly-combined estimations allows receiving effectiveness of estimation somehow close to effectiveness of maximum likelihood estimations in conditions of lack of the information of distribution functions of parasites of observations.

n

i Надоели баннеры? Вы всегда можете отключить рекламу.