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

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

Пусть — н.о.р.с.в. с EXi = 0 и DXi = 1. Обозначим через F общую функцию распределения (ф.р.) этих случайных величин и Fn (x) = < х) функцию распределения нормированной суммы Xl+^Xn первых п из этих случайных величин. ЦПТ утверждает, что.

Fn (x) —> Ф (ж) при п —> оо х равномерно по —оо < х < оо, где Ф (а-) = f е~и /2du — функция распределеоо ния стандартного нормального закона. При выполнении некоторых дополнительных условий у ф.р. Fn{х) существует плотность рп (х) и рп (х) —> (р (х) при п —> оо равномерно по —оо < х < оо, где р (х) = /2 — плотность стандартного нормального закона.

Важность ЦПТ объясняется тем, что она позволяет в практических расчетах заменять (при больших п) ф.р. Fn на ф.р. Ф, работа с которой не представляет трудностей. Функцию распределения Fn (x) можно записать в виде F*n (y/nx), где *" означает п—кратную свертку функции распределения F, точнее, оо оо.

F*n (x) = J. J F (x-yx——-Vn-x)dF{yi). dF (yn-i), -oo < a- < oo, oo —oo 4-v-" .

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

Например, в случае, когда F (x) является экспоненциальным распределением с параметром единица, то есть F (x) = 0 при х < 0 и F (x) — 1 — е~х при х ^ О га—1 т.

F™{x) = 1 — Р (X! + ¦ • • + xn > х) = 1 — е~Х 7[ г=0 при х > 0. Прямые расчеты по этой формуле при больших п невозможны хотя бы из-за того, что 70! > Ю100. Как уже отмечалось, ЦПТ позволяет заменять многократные свертки нормальными законами, работа с которыми не вызывает трудностей. Однако, при такой замене мы всякий раз (за исключением тривиального случая, когда F — нормальная функция распределения) совершаем некоторую ошибку, и возникает естественный вопрос о величине этой ошибки, или, как иногда говорят, о точности аппроксимации в ЦПТ.

Одним из самых известных результатов в этом направлении является теорема Берри-Эссена, которая гарантирует, что sup Fn (x) — Ф (х) ^ ¿-Щ^, (1) оо<�х<�оо Vn где с — некоторая константа. Эта оценка является неулучшаемой с точностью до значения константы с, для которой известна как верхняя оценка с < 0, 7056 ([32]), так и нижняя оценка с ^ = 0,409. (Эссен, 1956 год, [6]). К сожалению, точность оценки Берри-Эссена невелика. Если мы захотим гарантировать с помощью (1) справедливость неравенства p (Fn, Ф) ^ Ю-3, то в силу того, что E|Xt |3 ^ 1 (это следует из неравенства Ляпунова), величина п должна быть более (103с)2 > 160 000. В случае, когда нормированная сумма состоит из нескольких десятков слагаемых, оценка теоремы Берри-Эссена, по существу, бессодержательна.

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

2) а пример оценки, содержащей псевдомоменты, дает неравенство p (Fn, Ф) < при 4, (3) у/П где С — некоторые постоянные, a v^ - метрика на множестве функций распределения, которая называется вариацией с весом (здесь вес равен |я|3) и определяется равенством оо оо u3(V, W) = J x3d (V — W){x) = sup J x3w (x)d (V — W){x) oo l~oo где верхняя грань берется по множеству таких измеримых функций w, что |ги (ж)| ^ 1, -оо < i < оо, а У и Ж — произвольные функции распределения.

Значения постоянных в этих оценках оказывается существенно больше значения постоянной с в неравенстве (1).

Таким образом, оценка (2) не имеет преимуществ перед (1) при не очень больших |ж|, а (3) не имеет преимущество перед (1), если расстояние ^(.F, Ф) не очень мало.

По-видимому, малая точность аппроксимации, которую гарантируют приведенные оценки, связана с тем, что они применимы для очень широкого класса распределений F: они справедливы для любой ф.р. F с конечным третьим моментом. Существенное продвижение в оценках точности аппроксимации в ЦПТ можно получить за счет сильного сужения множества функций распределения исходных случайных величин. Так известно, что если для некоторого натурального т ^ 2 абсолютный момент (Зт+2 функции распределения F конечен, моменты F совпадают с моментами Ф вплоть до порядка т+1, F обладает некоторой гладкостью, то где с зависит лишь от т. Известны явные оценки величин о (^71) ¦ Это неравенство приводилось в спецкурсе «Дополнительные главы теории вероятностей», прочитанном В. В. Сенатовым на механико-математическом факультете МГУ в 2002 г.

Для справедливости (4) необходимо наложить на F ограничения, связанные с ее гладкостью. Эти условия, по крайней мере, должны гарантировать нерешетчатость F, поскольку для любой решетчатой ф.р. F величины p (Fn, Ф) ^ 0.125^= для всех п, начиная с некоторого, где h — шаг распределения F, при любых (сколь угодно сильных) ограничениях на моменты.

Одним из условий на гладкость F, которое гарантирует выполнение (4), является условие Крамера limsup|/(t)| < 1, оо где / - характеристическая функция распределения F. Выполнение этого условия гарантирует существование у F непрерывной компоненты в ее лебеговом разложении. При выполнении условия Крамера и упомянутых выше условий на моменты справедливо неравенство (4), в котором для о (—^77) можно получить явную оценку, в которую входит величина a (T)=sup{|/(i)|:i>T} 0. При этом о из неравенства (4) при росте п убывает экспоненциально быстро.

Условия, при которых доказывается неравенство (4) близки, а при четных m совпадают с условиями одной теоремы И. А. Ибрагимова, устанавливающей связь между скоростью стремления к нулю величины p (Fn, Ф) и значениями моментов F.

По-видимому, точности, которую гарантирует неравенство (4), достаточно для большинства практических расчетов уже при не очень больших ш, скажем, для ш, больших 3−5, однако ограничение, связанное с совпадением моментов F и Ф вплоть до порядка ш + 1, очень сильно сужает область применения неравенства (4).

Хорошо известен еще один подход к аппроксимации распределений Fn, он связан с так называемыми асимптотическими разложениями. В этом подходе нормальный закон Ф рассматривается только как первое приближение распределения Fn и аппроксимация Fn ищется в виде суммы Ф и некоторых слагаемых, стремящихся к нулю при росте п. Асимптотические разложения в ЦПТ появились в работах Грама [7] 1883 года, Шарлье [2] 1913;1914 годов и в работе Эджворта [4] 1905 года. В 1920;х годах асимптотические разложения интенсивно изучались Г. Крамером, а затем и другими исследователями, которыми были получены важные результаты, однако подавляющее большинство этих результатов давало оценки точности для асимптотических разложений в терминах О и о где, а > 0 — некоторое число, зависящее от количества моментов, которые существуют у распределения F. По-видимому, первые результаты с явными оценками точности для асимптотических разложений появились в конце 20-го века в работах Шимицу [8], Добрич, Гош [3]. В 1990;х годах появились результаты В. Сенатова. которые были получены с использованием сопровождающих зарядов. Кратко опишем содержание этого подхода, детальное описание которого можно найти в [27].

Для данного распределения .Р, у которого конечен момент /Зт+2, т ^ 2, и которое является достаточно гладким, попробуем подобрать гладкое распределение С с конечным моментом порядка т + 2 и первые т + 1 моментов которого совпадают с соответствующими моментамиР. Рассуждения, аналогичные тем, что использовались при доказательстве (4) показывают, что нормированные свертки Рп и (?" при росте п сближаются друг с другом быстрее, чем они сближаются с нормальным законом. Оценка величины р (Рп, Сп) имеет тот же вид, что (4) с явными оценками с и о (^72). Если нам удастся подобрать распределение С, для которого просто вычисляются функции Оп или просто получаются асимптотические разложения (?п, то в качестве асимптотического разложения Еп можно взять асимптотическое разложение бгп и тем самым задача построения асимптотического разложения для Рп будет решена.

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

Распределение (7 с указанными выше свойствами естественно назвать сопровождающим. В. Сенатов искал сопровождающие распределения в классе абсолютно непрерывных распределений, плотности д (х) которых допускают представление.

ОО д = Е тт^М*). (5) г=о где Щ (х) = (—1)1(р1х)/(р (х) — многочлены Чебышева-Эрмита, а — числа, называемые моментами Чебышева-Эрмита распределения С?

ОО.

01= J Н1{х)йС{х). оо.

Так как функции Н^х) являются многочленами 1-то порядка, то для вычисления 61 необходимо и достаточно знать моменты а^ = ЕХ (, 3 =2>,., 1.

Обратим внимание на то, что равенство (5) представляет собой разложение отношения по многочленам Чебышева-Эрмита, которые образуют полную ортогональную систему в пространстве Ь2{ф) действительных функций, квадрат которых интегрируем с весом <�р (х) по всей действительной прямой. Условия сходимости ряда (5) являются очень тяжелыми, но они оказываются ненужными, если еще сильнее сузить класс сопровождающих распределений, потребовав дополнительно, чтобы ряд из (5) состоял из конечного числа слагаемых. Этот класс распределений замкнут по отношению к переходу к нормированным сверткам. Действительно, если для некоторого распределения С?)ф), (6).

1=о то используя выражение для преобразования Фурье функций Нг (х)<�р (х), мы видим, что характеристическая функция.

1=0 характеристическая функция gn (t) = дп есть.

— t2/2 if' к i2 /О I 1 е это произведение можно представить в виде кп.

2,.

1=0 а поскольку обратные преобразования Фурье функций e~t2/2(it)s суть Hs (x)tp (x), то.

9"(«) = Е^Я/(хМХ)> (7) г=о L где 9i{qn) — моменты Чебышева-Эрмита плотности qn, то есть правая часть (7) имеет, с точностью до значений числовых параметров, тот же вид, что и правая часть (6).

Пусть F — гладкая функция распределения (условия гладкости будут уточнены ниже) с нулевым средним и единичной дисперсией, у которой существует момент ¡-Зт+2) тп ^ 2. Если мы при построении распределения G в (6) положим к = т + 1 и в качестве I = 1,., т + 1, возьмем моменты Чебышева-Эрмита распределения F, то первые т + 1 моментов F и G будут совпадать и можно показать, что рп (х) — qn (x)| = О (^kfi)) и Для величины О (—^72) можно получить явную оценку. Таким образом, G будет являться сопровождающим для F. Легко проверить, что момент Чебышева-Эрмита нулевого порядка равен единице для любого распределения, и для распределений с нулевым средним и единичной дисперсией моменты Чебышева-Эрмита первого и второго порядка равны нулю. Поэтому в данном случае (6) имеет вид m+l q q{x) = ч>{х) + Y^ lTHi (xMx)>

1=з m+l)n. qn (x) =.

ЧЯп) у^ п (в3 кз (Jh ^ l ^ k0k3. ki 3п3/2 J «ln½, для I ^ m + lпри I > т + 1 величины, связанные с 6j, j > т + 1, в правой части этого равенства следует опустить. Здесь суммирование производится по всем наборам ко, k3,., ki неотрицательных целых чисел, таким, что к0 + к3 +. + kt=n, 3 кг +. + Iki = I.

В правой части формулы для плотностей qn{x) участвует сумма растущего числа слагаемых. Естественно попробовать выделить из этой суммы фиксированное число слагаемых. Сумму этих слагаемых вместе с функцией <�р (х) мы будем называть главной частью разложения, а сумму оставшихся слагаемых мы будем называть остаточной частью. Задача построения асимптотического разложения для qn будет решена, если мы сможем оценить остаточную часть разложения. Для этого нужно задаться точностью аппроксимации, которая нас устраивает, и ответить на вопрос, как при п —> оо стремятся к нулю числа Oi{qn). Ответ на этот вопрос был получен Г. Крамером, из результатов которого следует, что.

8i (qn) = О ^ ПРИ.

Отсюда следует, что если мы хотим получить для qn (x) аппроксимацию точности О где целое число т ^ 2, то главная часть разложения qn (x) должна иметь вид.

3m-3 ", N ф)+? *ШН1(х)ф). (А) i=3,i#3m-4.

Для всех других I (в том числе для I = Зга — 4) О (фц). и при определенных условиях на значения 63,., вт+1 сумма этих слагаемых также есть O (-^jj). Выше говорилось, что рп{х) — qn (x)| = 0(^75-). Таким образом, в качестве главной части разложения рп (х) можно взять главную часть (А) асимптотического разложения для qn (x).

Следует отметить, что если плотность qn[x) существует и для нее ряд (5) сходится, то частичные суммы ряда из (5) не обязаны быть знакопостоянными функциями, то есть мы рассматриваем не только сопровождающие распределения, но и сопровождающие заряды.

Проиллюстрируем сказанное следующим примером. Пусть р — плотность распределения ф.р. F с единичной дисперсией, которое является симметричным и его шестой момент конечен. В этом случае момент о>2 = lj поэтому 0 = 62 = = = 0. Так как Нл (х) = гс4—6×2+3, то момент 04 распределения F равен а4—баг+З = «4—3, «4 = ЕХ±-. Нетрудно проверить, что при 0 ^ ^ 4 функция q (x), построенная с помощью равенства (6) с к = 5, то есть Q q (x) = ip (x) + -^ff4(x)ip (x), является плотностью распределения, а в иных случаях q (x) знакопеременна, но всеоо гда f q (x)dx = 1. Характеристическая функция для q (x) есть оо.

9(t) = e~^ (l + а характеристическая функция для дп (х) есть.

Ш-Н ШУ=§ Ше, 2/2(а)4 •.

Формула обращения для преобразования Фурье дает равенство дп (х) = ф) + С^Н4{х)ф) +? С* (-^У Нф) ф).

3=2 ^ ' '.

Можно показать, что при конечности шестого момента распределения ^ и выполоо нении условия / /(Ь)и (Ы < оо, где и > 0, а / - х.ф. ф.р. ^ (это условие гладкости оо.

Р), выполнено равенство.

Рп{х) ~ дп (х) = О прита-юо, и нетрудно получить явную оценку величины О Легко проверить, что абсолютная величина последней суммы по j = 2,., п не превосходит.

3=2 4 7 где В^ - абсолютные моменты нормального закона, деленные на /27г.

Из одной леммы, доказанной в [27], следует, что если 64 < 6, то при всех п таких, что т, «1 р=1т + т<1< последняя сумма не превосходит в^ 2.

2(1 — р) V41 / п2'.

Таким образом, для плотности рп получим асимптотическое разложение.

1п{х) = ф) +^Щ (х)ср (х) + О ^^, где для О можно получить явную оценку. При выполнении условия 64 < 6 эта оценка зависит лишь от шестого момента распределения ^ и от величины а (Т), которая упоминалась на странице 5. Условие 64 < 6 можно ослабить. Так, если число слагаемых в нормированных суммах четно, то от нормированных сумм можно перейти к нормированным суммам, где Х =, =. Момент Чебышева-Эрмита случайной величины равен у, и изложенная техника работает при 84 < 12.

Заряды с плотностями (6) — они называются зарядами типа 1 — не являются единственно возможными. В работе А. Е. Кондратенко [20] использовались заряды типа.

2, которые строились с помощью семиинвариантов, их характеристические функции суть где ае/, I = 3,., к, — семиинварианты распределения Р.

К сожалению, область применения зарядов типа 2 оказывается достаточно узкой. Так, если величина к в (В) равна 4 (то есть рассматривается распределение с конечным моментом /3 $), то необходимо налагать ограничение ае.1 < 0 или, что-то же самое, о:4 < 3. Это ограничение выполняется лишь для достаточно узкого класса распределений.

Стоит отметить, что в общем случае при работе с зарядами типа 1 для оценки остаточной части разложения плотности дп (х) приходится налагать ограничения на все используемые моменты Чебышева-Эрмита распределения Р. Одно из таких ограничений есть где Т > 0 — некоторое число, см. [27]. Это ограничение можно ослабить, разбивая суммы + • ¦ ¦ + Хп исходных случайных величин на блоки, состоящие из нескольких слагаемых. Применяя изложенную технику к блокам слагаемых, а затем возвращаясь к исходным величинам, можно получать требуемые результаты, но такой путь приводит к довольно громоздким рассуждениям.

Некоторые преимущества перед зарядами типов 1 и 2 имеют так называемые заряды типа 3, характеристические функции которых суть.

Заряд с х.ф. дк, к = 3,4,5, предназначен для получения разложений для распределений с конечным моментом порядка к +1. При использовании заряда дз ограничения на величину в3 не налагаются. При использовании зарядов с х.ф. д^ и д5 налагается только ограничение в4 < 6.

Можем продолжить рассуждения и по аналогии построить заряд с х.ф. де (£). Однако для этого заряда большая часть преимуществ зарядов типа 3 утрачивается. Кроме того, по-видимому, результатов, полученных с помощью зарядов с х.ф. дк, к = 3,4,5, достаточно для проведения большинства практических расчетов. При необходимости можно использовать технику без использования вспомогательных зарядов, которая помогает получать разложения сколь угодно высокой точности, если у исходного распределения существует достаточное количество моментов, однако при этом рассуждения окажутся более громоздкими.

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

3,(4) =.

Краткое содержание диссертации.

• Целью работы является обобщение результатов, о которых говорилось выше, на многомерный случай. При этом используются многомерные аналоги многочленов Чебышева-Эрмита, введенные в [30].

В первой главе дается определение многомерных аналогов многочленов Чебышева-Эрмита как полилинейных функционалов, полученных с помощью дифференцирования по Фреше функций, связанных с плотностью нормального закона, для некоторых из них указывается явный вид, перечисляются известные свойства этих функционалов, формулируются и доказываются утверждения, обобщающие свойства одномерных многочленов Чебышева-Эрмита'на эти функционалы.

Будем рассматривать обобщения многочленов Чебышева-Эрмита на многомерный случай. Обозначим Еа — евклидово пространство размерности с1, (-, •) — скалярное произведение в Еа, | • | - норму в Е'1. аш — т-й момент одномерного стандартного нормального распределения, <�р (х) = /2 — плотность нормального закона с нулевым средним и единичной матрицей ковариаций в Еа. Многомерные аналоги многочленов Чебышева-Эрмита — это множество многочленов 1,3= 0,1, 2,., которое являются последовательностью серий, состоящих из счетного числа элементов. Первый элемент каждой серии является функцией на Ей, второй элемент каждой серии — линейный функционал, третий элемент — билинейный функционал и т. д.

Эти многочлены можно свести в таблицу и-(О) и-(О) гг (0) И (0) гт (0) тт (0) л0) п 1 2 ! 173) 4) л5 ! ' ' «.

Н{2) и-(2) я (2) я (2) п0 1 п1 } л2 > п3) • • • Я (4) я (4).

Многочлены Чебышева-Эрмита из первой серии (ей соответствует ^ = 0), точнее, действие многочленов х) на векторы. , 1ц € Еа, можно получить по формуле н^хХЬ,., Ы) = (-1)1 (ф)){1) (нъ., Ы)/Ф) производная понимается в смысле Фреше), которая аналогична равенству Щ (х) = (—1)1 ((р (х))^ /1р (х), справедливому в одномерном случае. Величина Н1°х)(1г 1,., к{) является симметрической функцией переменных ., поэтому мы иногда будем записывать ее в виде х,., Ы}), где {/11,., Л,/} означает множество, состоящее из векторов .,.

Первый элемент первой серии — функция на Еа, тождественно равная 1, второй элемент первой серии — линейный функционал (х, •), третий элемент этой сериибилинейный функционал (х, -)(х, ¦) — (•,¦), х Е Ей, и т. д. Справедлива формула обращения I е-^е-1{12/2(*,. (?, й,)й = Н°х){1гъ ., Ы) ф),.

Е* где г — мнимая единица.

Для любого многочлены I = 0,1,2,., определяются равенствами г.

Н^х){Н ., ЫЫх) = щ~а] е-^е-1'12/2(?, *)''(*, Лх). (*, еа.

В первой главе указан явный вид некоторых из этих многочленов и рассматриваются их простейшие свойства.

В одномерном случае многочлены Чебышева-Эрмита являются многочленами Апреля, то есть для них справедливо равенство.

Н[(х) = 1Н1−1{х), 1 = 1,2,. Аналог этого свойства для многочленов Н^ (х) имеет вид я,(0)(аО (Л1, • • •, Л0) я №+0 = Х>1+1, Лг) я?>1(®)({'и" • • •' М (М)г=1.

Также выполнено рекуррентное соотношение Я$(аО (/ц,., Л/+1) = (*, Ы+1)Н}0)(х)(Ь, ., Л,)I.

— ]Г>+1, ., М {Л,}). (8).

Г=1.

Эти равенства получены в [30].

Для многомерных аналогов многочленов Чебышева-Эрмита произвольной строки в первой главе получены соотношения I Ы))я {Ы+1) = $>1+1, К) н2Л{х){{ни., А,} (М)+ г=1.

Лх,. Л+1″ и рекуррентная формула.

Ж)(/гь., Ы+1) = (х, Л,+1)Я,(!Ц)(*)(Ль ., Л,)I.

— Лг) Д1с^?(я:)({/ц, Л,} (М) — ЦН^^Шг,., Л^}).

7-=1.

Во.

ВТОрОИ главе в многомерном случае строятся новые асимптотические разложения плотностей нормированных сумм независимых случайных величин, у которых конечны моменты 5 и 6 порядков (асимптотические разложения в случае конечности моментов меньших порядков получены в [28]). Эти асимптотические разложения строятся с использованием вспомогательных сопровождающих зарядов, которые представляют собой функции, получаемые интегрированием многомерных аналогов многочленов Чебышева-Эрмита по распределению исходных случайных величин. Приводятся явные оценки остаточных частей разложений.

Здесь и далее мы будем действовать в рамках следующих обозначений и предположений.

Пусть ХЪХ2,. — независимые одинаково распределенные величины в Ed с нулевым средним и единичным ковариационным оператором и Р — распределение Х. Пусть для характеристической функции f (t) распределения Р выполнено условие.

У f{t)udt < оо (а) ed для некоторого и>0, и для некоторой пары (ц, Т), где функция е-'''2/2 ^ fi (t) ^ 1 и число Т > О, выполняется неравенство |/(t)| ^ ?i{t) при всех |?| ^ Т.

Через Рп обозначим распределение нормированной суммы. При выполнении условия (а) для всех п^ и существует плотность рп (х), х Е Ed, распределения Рп, которую можно вычислить по формуле обращения.

Определим моментные характеристики распределения Р, которые нам понадобятся. Для s е N и вектора е € Ed, (е, е) = 1, положим as (e) = as (e, P):= J{e, u) sP (du), &-(е) = &-(е, Р) := J (e, u) sP (du), ed ed as| = sup |ore (e)|, Ps = sup /3s (e), s = 0,1,. e|=l |e|=l.

Легко видеть, что аа (е) и Д3(е) — момент и абсолютный момент s-ro порядка проекции Р на направление вектора е.

Во всех леммах и теоремах мы будем предполагать, что для некоторого тп € N величина /Зт+2 конечна.

Мы будем использовать многомерные моменты Чебышева-Эрмита 6S, s = 0,1,., которые определяются следующим образом: si — si — j!2J (s — 2j) ~ j2?(s — 2j) J ^ ' >

J-U J-U Ed, а также величины в3 = sup|0B (e)|, 3 = 0,1,., e|=l.

JM:=< iff K2jl goi!2i{s-2j)V Д^ четных О 4.

7 Г + g j2i (S-2j) ' ДЛЯ НеЧеТНЫХ S ^ 5.

В оценках остаточных частей асимптотических разложений будут использоваться величины Ц^^Ц, которые определяются формулами для ||0Л||, в которых нужно положить |а:г| = 0 при / ^ й + 1.

Отметим, что в3(е) — момент Чебышева-Эрмита 5-го порядка проекции Р на направление вектора е.

Пусть Ф (:г) — нормальный закон с нулевым средним и единичным ковариационным оператором в Ed. Для е е Ed, e = 1, и натуральных s обозначим.

Дв (е) = 1(е, иУ (Р-Ф)((1и), 6S = sup |ДДе)|, = J |(е, и)|вФ (<�Ь). ed ed.

Отметим, что для любого е? Ed, е = 1, величина /3s (e) совпадает с s-м абсолютным моментом одномерного стандартного нормального закона. Введем следующие величины в- = (?F / W**"*, в"п — ^ / WV Л.

Величина Д, асимптотически (при d —> оо) равна ,.nds/2. Для распределений с.

27Гс{/2 конечным четвертым моментом пару (/i, Т) можно подобрать так, чтобы Д5]П —"¦ Д, при п —" оо для любого фиксированного s > 0.

Также отметим, что из формулы Тейлора с остаточным членом в интегральной форме следует, что для любого действительного, а к isas ak+1 э).

6 ?-J si {к +1)!' { } s=0 где к — произвольное неотрицательное целое число и 7 — непрерывная комплексно-значная функция такая, что |7| ^ 1. Итак, пусть рп (х) — плотность нормированной суммы {'Х1+^Хп), дп (х) — плотность нормированной п-кратной свертки заряда.

Мы можем записать плотность рп (х) в виде.

Рп (х) = Яп (х) + рп[х) — Яп (х)]. В случае использования заряда доказательство состоит из двух шагов:

• оценка близости плотностей нормированных сумм |рп (х) —.

• разложение плотности дп (х) по многочленам с некоторой погрешностью.

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

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

Теорема 2.4 Пусть распределение Р таково, что р5 < оо и /З4 < 9. Тогда в рамках введенных обозначений и предположений для всех п ^ тах (3, и) таких, что.

Р 6 Зп 72п при /?4 < 9 величина р < 1 для всех п, начиная с некоторого), при всех х Е Еа для плотности рп (х) справедливо равенство рп (х) = ф) + J ((х, и)3 — 3(х, и)(и, и)) Р (с1и)+ ф) 24п У (х) 72 п.

У Н? х)(и, и, и, и) Р (йи) — 3н? х) J ! Н^х)(и, и, и, у, у, у) Р (йи)Р (^)+.

ЕЛ ЕЛ.

Ф) I (x)(u, u, щv, v, v, v) P (du)P (dv)-^^ I Н^(х)(и, и, и)Р (аи)+.

144³/2 ел е<>

3!)%3/2 ел.

Щ Н?>(я)(и, «, и, «, «, V,», + Л, Еа где В, — остаточная часть, для которой выписывается явная оценка.

Теорема 2.5 Пусть распределение Р таково, что (3 $ < оо и /?4 < 9. Тогда при всех х б ЕЛ для всех п ^ тах (3, и) таких, что.

Р=Т + Зп + 72п <�Х' для плотности рп (х) справедливо равенство рп (х) = ф) + [ н!? х)(и, и, и)Р (с1и)+ J.

ЕЛ ф) ' 24? г.

J Н1°х)(и, и, и, и) Р{йи) — ЗН^х).

КЕ<1.

7^//- /нР (х)(и, и, и)Р (<1и)+ еа е<* еа.

УI Н? х)(и, и, и, у, у, у, у)Р^и)Р№) — ^^ I Н^х){и)Щи)Р{йи)Л.

144гг3/2 ел ел.

Ф).

3!)4п3/2.

J J J Нд°х)(и, и, и, у, V, у, w, и], w) P (du)P (dv)P (dw) + Я, ел ел ел где В, — остаточная часть, для которой выписывается явная оценка.

Теорема 2.6 Пусть распределение Р таково, что (Зе < оо и ?3± < 9. Тогда при всех х? ЕЛ для всех п ^ шах (3, и) таких, что Зп + 72п для плотности рп (х) справедливо равенство рп (х) = ф) + '0=1 н!? х)(щи, и) Р ((1и)+.

Ф).

24 п.

72 га.

J Н^х){и, и, и, и) Р (с1и)-ЗН^х) ! J Н⁰)(х)(и, и, и, v, v, и) р (сь)р (е^)+ I н (°х)(и, щщи, и) Р (аи)-^^ I нР (х)(и, и, и)Р (с1и)+.

120п3/2 р{х).

УIН$°х)(и, и, и, у, у, у, у)Р{ёи)Р№)-1^ I Н.

144 гг3/2 ел е* ел.

——————————;

3!)%3/2.

Е<1 где В. — остаточная часть, для которой выписывается явная оценка.

В приведенных теоремах оценки остаточных частей 7?, в каждой следующей теореме лучше, чем в предыдущей (во всяком случае, при больших п).

Явный вид оценок остаточных частей не приведен в этом введении ввиду их громоздкости.

При доказательстве теорем используется следующая.

Лемма 2.1 Пусть т^ 1, тбМ, и прг1 любом е? ЕЛ, (е, е) = 1, для распределения Р выполнено рт+2(р) = вир [ |(е, и) т+2р (и)(1и < ос, е|=1 3 е* тогда справедливо равенство ?9фапг + 7в^)/2Г+2 + где 7 — непрерывная комплекснозначная функция, такая, что |-~у| ^ 1 (здесь 7 может обозначать различные функции).

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

• с некоторой погрешностью представляем рп{х) — ср (х) в виде суммы Е. зависящей от разности моментов Р и Ф.

Д,(е) = У (в>и)в (Р-Ф)(<�г") — ел.

• слагаемые из Е, убывающие при росте п как 0(п (т|1)/2), включаем в остаточную часть разложения, оставшиеся слагаемые — в главную.

Представление разности плотностей с помощью ДДе^ дает.

Лемма 3.1 Во введенных предположениях и обозначениях при п ^ тах (г^, 2то) существует плотность рп (х) и для любого х Е ЕЛ справедливо равенство кш! г3Д5(е^).

1 Г 2 2 /.

Й.М — Р (Х) — ^?3 / е-'Ме-^С^ ?

I х I ^-«п /— &— 1 — 3 у/п.

Й+ Д,.

М^Ту/п где = щ и Я — остаточная часть, для которой выписывается явная оценка.

При помощи этой леммы получен следующий результат.

Теорема 3.1 В рамках введенных предполоо/сений и обозначений при всех п ^ тах (г/, 2т) у распределения Рп нормированной суммы существует плотность рп{х) и для всех х Е Ей рп (х) = 1р (х)+.

Е? к=1 кз,., кт+1.

П ср (х) п — к) к3 ¦ ¦ ¦ кт+1 (3!)**. ((т + 1)!)*-+1 п (зк3±+(™+1)кт+1)/2 X X.

V (-1)7 (Г у Г Г (2з), (3) (3) (3) (3) (3) (3) (т+1).

27 п) ] ] 1 ^ 1 1 1 ''''' кз ' кз ' кз ''''' ' «'' кт+1 '.

Е<1 Еа хсэ^и?) ¦ ¦ ¦ • • • • • • + я, кз / / т—-гего+1 где (5 = Р — Ф и внутреннее суммирование ведется по всем наборам ., кт+ таким, что кг ———Ь кт+1 = к, 3к3———1- (га + 1)/гт+1 < т + 2к,.

J=ш ^зН———{т-1)кт+1.

Д| < ^ш+1 -^3(т+1), п-(яг+1) дт+1 | П (2+ 5.

Рт+2 + ^3(г-1)+т+2,п-г г! т + 2)! г=1 п 2 т а/2.

2ТГ) тт т ¦> т I 1 13^- (л? +? {(¿-Я + + г=1 1=1 т И1 где для действительных а, Ъ, Т оо к11'2 Г «2 а г.

Т1 а для натурального I и положительного действительного Т ш+1 г т+1 г.

5 = 3 «я=3.

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

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

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