Асимптотические разложения в центральной предельной теореме в многомерных пространствах тема автореферата и диссертации по математике, 01.01.05 ВАК РФ

Осмоловский, Игорь Юрьевич АВТОР
кандидата физико-математических наук УЧЕНАЯ СТЕПЕНЬ
Москва МЕСТО ЗАЩИТЫ
2009 ГОД ЗАЩИТЫ
   
01.01.05 КОД ВАК РФ
Диссертация по математике на тему «Асимптотические разложения в центральной предельной теореме в многомерных пространствах»
 
Автореферат диссертации на тему "Асимптотические разложения в центральной предельной теореме в многомерных пространствах"

Московский Государственный Университет имени М.В.Ломоносова Механико-математический факультет

На правах рукописи УДК 519.214.4

Осмоловский Игорь Юрьевич

АСИМПТОТИЧЕСКИЕ РАЗЛОЖЕНИЯ В ЦЕНТРАЛЬНОЙ ПРЕДЕЛЬНОЙ ТЕОРЕМЕ В МНОГОМЕРНЫХ ПРОСТРАНСТВАХ

Специальность: 01.01.05 - Теория вероятностей и математическая

статистика

АВТОРЕФЕРАТ

диссертации на соискание ученой степени кандидата физико-математических наук

11&7СГ

Москва - 2009 г.

003470708

Работа выполнена на кафедре теории вероятностей механико-математического факультета в Московском Государственном Университете имени М. В. Ломоносова.

Научный руководитель: доктор физико-математических наук

Сенатов Владимир Васильевич.

Официальные оппоненты: доктор физико-математических наук

профессор

Круглов Виктор Макарович,

доктор физико-математических наук профессор

Хохлов Юрий Степанович.

Ведущая организация: Санкт-Петербургское отделение мате-

матического института.

Защита диссертации состоится 29 мая 2009 года в 16 часов 40 минут на заседании диссертационного совета Д.501.001.85 при Московском Государственном Университете имени М. В. Ломоносова по адресу: 119991, ГСП-1, Москва, Ленинские Горы, Главное Здание МГУ, механико-математический факультет, аудитория 16-24.

С диссертацией можно ознакомиться в библиотеке механико-математического факультета МГУ (Главное здание, 14 этаж).

Автореферат разослан 28 апреля 2009 года.

Ученый секретарь диссертационного совета Д.501.001.85 при МГУ, доктор физико-математических наук, профессор И.Н.Сергеев

1 Общая характеристика работы

Актуальность темы

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

Пусть ... - н.о.р.с.в. с ЕХх = 0 и = 1. Обозначим че-

рез Р общую функцию распределения (ф.р.) этих случайных величин и Рп(х) = < х) - функцию распределения нормированной

у I , у"

суммы ^ " первых п из этих случайных величин. ЦПТ утверждает, что

ад -» Ф(х) при п оо

X

равномерно по —оо < х < оо, где Ф(ж) = ^ / е~и !2йи - функ-

—00

ция распределения стандартного нормального закона. При выполнении некоторых дополнительных условий у ф.р. Рп(х) существует плотность рп(х) и рп{х) —при п -> оо равномерно по —оо < х < оо, где <р(х) = ^ ~ плотность стандартного нормального закона.

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

Р*п{х)= I ... I Р(Х-У1-----Т/я-1)^(У1)..-^(Уп-1), —ОО < X < ОО,

—оо —со

4—V—'

и-1

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

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

г=0

при х > 0. Прямые расчеты по этой формуле при больших п невозможны хотя бы из-за того, что 70! > 10100. Как уже отмечалось, ЦПТ позволяет заменять многократные свертки нормальными законами, работа с которыми не вызывает трудностей. Однако, при такой замене мы всякий раз (за исключением тривиального случая, когда F - нормальная функция распределения) совершаем некоторую ошибку, и возникает естественный вопрос о величине этой ошибки, или, как иногда говорят, о точности аппроксимации в ЦПТ.

Одним из самых известных результатов в этом направлении является теорема БерриМЭссена2, которая гарантирует, что

p{Fn, Ф) = sup (1)

—ооСКоо у"

где с - некоторая константа. Эта оценка является неулучшаемой с точностью до значения константы с, для которой известна как верхняя оценка3 с < 0,7056, так и нижняя оценка с ^ ^/^г = 0,409..Л К сожалению, точность оценки Берри-Эссена невелика. Если мы захотим гарантировать с помощью (1) справедливость неравенства p(Fn, Ф) ^ 10_3, то в силу того, что Е|Хх|3 ^ 1 (это следует из неравенства Ляпунова), величина п должна быть более (103с)2 > 160000. В случае, когда нормированная сумма состоит из нескольких десятков слагаемых, оценка теоремы Берри-Эссена, по существу, бессодержательна.

'Berry A. The accuracy of the Gaussian approximation to the sum of independent variates. Trans. Amer. Math. Soc. 1941. V.49. № 1. P. 122-136.

2Esseen C.-G. On the Liapunoff limit of error in the theory of probability. Ark. Mat. Astr. F>s. 1942. V.28A. № 9. P. 1-19.

3Шевцова И.Г., Уточнение верхней оценки абсолютной постоянной в неравенстве Берри-Эссена, Теория вероятностей и ее применения, 2006, т. 51, в. 3, с.с. 622 - 626.

4Esseen, Moment inequalities with application to the Central Limit Theorem, Scand. Actuarietidskrist. N 34.1956. P.160-170.

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

(2)

а пример оценки, содержащей псевдомоменты, дает неравенство

p(Fn, Ф) ^ при п > 4, (3)

уп

где С - некоторые постоянные, а 1/3 - метрика на множестве функций распределения, которая называется вариацией с весом (здесь вес равен |а;|3) и определяется равенством5

00 оо

V3{V,W) = J |x|3|d(V - И0(®)| = sup J x3w{x)d{V-W){x)

—00 (-00

где верхняя грань берется по множеству таких измеримых функций w, что |го(а;)[ ^ 1, —00 < х < <х>, a.V и W - произвольные функции распределения.

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

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

По-видимому, малая точность аппроксимации, которую гарантируют приведенные оценки, связана с тем, что они применимы для очень широкого класса распределений F: они справедливы для любой ф.р. F с конечным третьим моментом. Существенное продвижение в оценках точности аппроксимации в ЦПТ можно получить за счет сильного сужения множества функций распределения исходных случайных величин. Так известно, что если для некоторого натурального m ^ 2 абсолютный момент ¡Зт+2 функции распределения F конечен, моменты F совпадают с

5Ульянов В.В. К уточнению оценок скорости сходимости в центральной предельной теореме. Теория вероятностей и ее применения. 1978. Т23. Вьш. 3. С. 684-687.

моментами Ф вплоть до порядка т + 1, F обладает некоторой гладкостью, то

+ (4)

где с зависит лишь от тп. Известны явные оценки величин о (^д)- Это неравенство приводилось в спецкурсе "Дополнительные главы теории вероятностей", прочитанном В.В. Сенатовым на механико-математическом факультете МГУ в 2002 г.

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

Одним из условий на гладкость F, которое гарантирует выполнение (4), является условие Крамера6

lim sup \f(t)\ < 1, |f|—»00

где / - характеристическая функция распределения F. Выполнение этого условия гарантирует существование у F непрерывной компоненты в ее лебеговом разложении. При выполнении условия Крамера и упомянутых выше условий на моменты справедливо неравенство (4), в котором для 0 {'^mß) можно получить явную оценку, в которую входит величина

.a(T) = sup{|/(i)|:i>r}<l

для некоторого Т > 0. При этом о (-¿д) из неравенства (4) при росте п убывает экспоненциально быстро.

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

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

6Крамер Г. Случайные величины и распределения вероятностей. М.: Изд-во иностр. литер., 1947.

7Ибрагимов И.А. О точности аппроксимации распределения сумм независимых величин нормальным распределением, Теория вероятностей и ее применения. 1966. Т.Н. Вып.4. С.632-655.

совпадением моментов F и Ф вплоть до порядка m + 1, очень сильно сужает область применения неравенства (4).

Хорошо известен еще один подход к аппроксимации распределений Fn, он связан с так называемыми асимптотическими разложениями. В этом подходе нормальный закон Ф рассматривается только как первое приближение распределения Fn и аппроксимация F„ ищется в виде суммы Ф и некоторых слагаемых, стремящихся к нулю при росте п. Асимптотические разложения в ЦПТ появились в работах Грама8 1883 года, Шарлье9 1913-1914 годов и в работе Эджворта10 1905 года. В 1920-х годах асимптотические разложения интенсивно изучались Г. Крамером11, а затем и другими исследователями, которыми были получены важные результаты, однако подавляющее большинство этих результатов давало оценки точности для асимптотических разложений в терминах О (£) и о где а > 0 - некоторое число, зависящее от количества моментов, которые существуют у распределения F. По-видимому, первые результаты с явными оценками точности для асимптотических разложений появились в конце 20-го века в работах Шимицу12, Добрич, Гош13. В 1990-х годах появились результаты В. Сенатова, которые были получены с использованием сопровождающих зарядов.

Основные результаты по исследованию явных оценок асимптотических разложений в ЦПТ были получены В. В. Сенатовым. Им были получены явные равномерные оценки остаточных частей разложений.

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

Цель работы

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

8Gram J.P., J. reine und angew. Math, 1883, Bd 94, S. 41-43.

'Charlier C.V.L., Arkiv. Mat. Astr. Fys., 1913/1914, Bd 9, N 25, S. 1-17.

10Edgeworth, The law of error, Camb. Phil. Soc. Proc. 20,1905, 36-141.

"Крамер Г. Случайные величины и распределения вероятностей. М.: Изд-во иностр. литер., 1947.

12Shimizu R_, On the remainder term for the central limit theorem, Ann. Inst. Stat. Math. 1974. V.26. P.195-201.

13Dobric V., Ghosh B.K. Some analogs of the Berry-Esseen bound for first order Chebychev-Edgeworth expansions, Stat. Decis. 1996. V.14.

Научная новизна

Основные результаты диссертации являются новыми и состоят в следующем:

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

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

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

• Для всех асимптотических разложений получены явные оценки остаточных членов.

Методы исследования

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

Теоретическая и практическая значимость

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

Апробация работы

Основные результаты диссертации докладывались и обсуждались на Большом семинаре кафедры теории вероятностей механико-математического факультета МГУ (руководитель - член-корр. РАН А.Н.

Ширяев) в 2008 году и на семинаре кафедры статистики ВМК МГУ (руководитель - академик Ю. В. Прохоров) в 2008 году. Был сделан доклад на международной конференции "Ломоносов-2009". Результаты диссертации были представлены на семинаре "Прикладные аспекты теории вероятностей и математической статистики "кафедры теории вероятностей и математической статистики РУДН (руководитель - д.ф.-м.н. Ю.С.Хохлов) в 2009 г.

Публикации

По теме диссертации опубликованы 3 работы. Список приведен в конце автореферата [1] - [3].

Структура диссертации

Диссертация состоит из введения, пяти глав, приложения и списка литературы, насчитывающего 33 наименования и организованного в алфавитном порядке. Результаты, полученные автором диссертации, оформлены в виде Теорем и Лемм. Нумерация лемм и теорем состоит из двух чисел. Первое число относится к номеру главы, второе к номеру утверждения (леммы или теоремы). Нумерация формул - сквозная. Общий объем работы составляет 78 страниц.

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

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

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

Обозначим Ей - евклидово пространство размерности й, (•, •) - скаляр-

14Сенатов В.В. Об одном многомерном аналоге разложения Чебышева, Теория вероятностей и ее применения, 2007, т. 52, в. 3, с.с. 603 - 610.

ное произведение в Ed, | • | - норму в Ed, ат - т-й момент одномерного стандартного нормального распределения, <р(х) = ^ ~ плот~

ность нормального закона с нулевым средним и единичной матрицей ко-вариаций в Ed. Многомерные аналоги многочленов Чебышева-Эрмита -это множество многочленов l,j — 0,1,2,..., которое являются последовательностью серий, состоящих из счетного числа элементов. Первый элемент каждой серии является функцией на Ed, второй элемент каждой серии - линейный функционал, третий элемент - билинейный функционал и т.д.

Эти многочлены можно свести в таблицу

о-(О) о-(О) гт(0) rr(0) гг(0) гг(0)

я0 ! -"1 > 2 > п3 > п4 > пЪ >••■

и(2) гг(2) тт( 2) гг(2)

О ' 1 ' 2 ' -"3 »• • •

Я( 4) я(4)

Многочлены Чебышева-Эрмита из первой серии (ей соответствует ] = 0), точнее, действие многочленов (ж) на векторы }ц,..., Ы 6 можно получить по формуле

... л) = (-1)' Ы)/Ф)

(производная понимается в смысле Фреше), которая аналогична равенству Н\{х) = (—1)'(<р(х))^ /¡р(х), справедливому в одномерном случае. Величина Н1°\х)(к1,..., Лг) является симметрической функцией переменных Лх,..., /г;, поэтому мы иногда будем записывать ее в виде Щ (х)({/1х,..., Ы}), где {/11,..., 1ц} означает множество, состоящее из векторов Ъ,х,...,

Первый элемент первой серии - функция на Етождественно равная 1, второй элемент первой серии - линейный функционал (х, •), третий элемент этой серии - билинейный функционал (ж, -)(а;, ■) — (•, •), х Е Еа, и т.д.

Справедлива формула обращения

I ^ . . . ^ ^ = И?\х)(Нъ . ■ • , Н,)ф),

Е*

где г - мнимая единица.

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

равенствами

Н^\х){кь ..., ЫЫх) = ¿)'(*, Ь)... (<, /ц)сй.

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

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

Н[(х) = Ш^х), 1 = 1,2,... Аналог этого свойства для многочленов Я® (х) имеет вид

(я/0) (*)(&!,..., Ы))х(Ы+х) = кг)н1%(х)({Ьъ ...,Ы}\ {кт}).

г=1

Также выполнено рекуррентное соотношение Я® (х)^,..., км) = (х, Ы+1)н1°\х)(к1г..., /ц)—

I

-^(Ь+и Лг)Я,®({!ч,...,Л1}\{Лг}). (5)

Г=1

Эти равенства получены В.В. Сенатовым в упомянутой работе14.

Для многомерных аналогов многочленов Чебышева-Эрмита произвольной строки в первой главе получены соотношения

I

(я/зд(х)(/11,..., Ы))х (Л/+1) = £>г+1, йг)ЯЙ)(®)({ЛьЧ \ {Лг})+

Г—1

и рекуррентная формула

я£?(®)(АьЫ+1) = (х, /гг+1)ЯР(х)(Л1,..., Ы)-

I

- /гг)я/!{)(х)({^,..., ад \ {Лг}) - 2щМ-^хЖЬ,..., ы+1}).

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

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

Пусть Х\,Х2,... - независимые одинаково распределенные величины в Е* с нулевым средним и единичным ковариационным оператором и Р - распределение Х\. Пусть для характеристической функции $(1) распределения Р выполнено условие

< оо (а)

Е*

для некоторого г/>0, и для некоторой пары (ц,Т), где функция е-1<1 /2 ^ ^ 1 и число Т > 0, выполняется неравенство |/(£)| ^ /¿(г) при всех |г| < Т.

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

Определим моментные характеристики распределения Р, которые нам понадобятся. Для в £ N и вектора е е Еа, (е, е) = 1, положим

аДе) = с*8(е,Р):= ^е^ГР^и), &(е) = &(е, Р) ^ |(е, «)|'Р(йи),

Е* Е"

= вир Ме)|, Д, = вир Д,(е), 5 = 0,1,...

|е|=1 |е| = 1

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

"Сенатов В.В. Несколько асимптотических разложений в ЦПТ в многомерном случае. Теория вероятностей и ее применения, 2008, т. 53, в. 2, с.с. 293 - 306.

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

Мы будем использовать многомерные моменты Чебышева-Эрмита в8, з = 0,1,..., которые определяются следующим образом:

вв _ в6(е) (~1)Ч-2,(е) Щ У

si

! ,ч!

у. (—l)Jg3_2j(e) у. (-1)' [(е иу-ЪрЫи) J—и J-и fid

а также величины

|6>s| = sup \9s{e)\, s = 0,1,...,

|e|=l

i^2 K-2jl

II л и E -Ч7Й7-^лТ > Для четных s > 4

JIM = J i=0 — 2j)!

В оценках остаточных частей асимптотических разложений будут использоваться величины ||^т)||, которые определяются формулами для ||03||, в которых нужно положить |а;| — 0 при I ^ S + 1.

Отметим, что в${е) - момент Чебышева-Эрмита s-ro порядка проекции Р на направление вектора е.

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

Д,(е) = 1(е,иУ(Р-Ф)(<1и), 6S = sup |Аа(е)|, 0,(е) = J |(е,1х)|'Ф(Л).

Ed Ed

Отметим, что для любого е 6 Ed, |е| = 1, величина /З3(е) совпадает с s-м абсолютным моментом одномерного стандартного нормального закона. Введем следующие величины

Е* |t|<TVn

1

Величина Bs асимптотически (при d —> оо) равна ^ ¿/2 ' Для распределений с конечным четвертым моментом пару (/л, Т) можно подобрать так, чтобы Вя>п —У Bs при п —» оо для любого фиксированного s > 0.

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

¿•-Е^+та (6)

в=0 4 '

где к - произвольное неотрицательное целое число и 7 - непрерывная комплекснозначная функция такая, что | ^ 1. Итак, пусть

рп(х) - плотность нормированной суммы

дп(х) - плотность нормированной п-кратной свертки заряда.

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

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

• оценка близости плотностей нормированных сумм \рп(х) — <7„(ж)|;

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

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

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

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

|041 N(¿-1) , N<¿(¿+2)

9 6 Зп 72 п

(при (З4 < 9 величина р < 1 для всех п, начиная с некоторого), при всех х

€ Еа для плотности рп(х) справедливо равенство

Рп(х) = ф) + |М у Н^{х){и,и,и)Р{йи)+

Е6

да

Е*

+ 111 Н^^щи^УЛУ^ы^Р^Р^Р^ + К,

Е* Е* Е*

где 11 = О (¿щ) , п —)■ оо, - остаточная часть, для которой выписывается явная оценка.

Явный вид многочленов I = 3,4,6,7,9, Яд4' и Н^ указан в главе 1. В частности,

Н^(х)(и,и,и) = (х,и)3 -3(х,и){и,и),

поэтому слагаемое в асимптотическом разложении рп(х), связанное с можно записать в виде

01Е((х,Х1)3-3(Х,Х1)(Х1,Х1)) =

= 01 (|х|3Е(Хье1)3 - ЗМЕрГьезХХьХО),

где ех = - с.в. распределения Р. Величина ех)3 являет-

ся третьим моментом проекции с.в. Х\ на направление вектора ех, а Е(Х, ех) (X, X) можно рассматривать как некий смешанный момент третьего порядка с.в. X.

Все интегралы в приведенном асимптотическом разложении являются величинами, зависящими от |ж| и от моментов исходных случайных величин Х1,Х2,... (среди этих моментов - моменты проекций с.в. Х\ на направление вектора ех, четвертый момент нормы Е|Х1|4, смешанные моменты, моменты ЕрГ^Хг)3 и т.д.)

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

Теорема 2.6 Пусть распределение Р таково, что < оо и ¡3^ < 9. Тогда при всех х £ Еа для всех п ^ тах(3, и) таких, что

Р = Зп + 72п <Х'

для плотности рп{х) справедливо равенство

рп{х) = „»(¡с) + IН^(х)(и,и,и)Р(йи)+

Е*

+^^Н[?\х){ч,ч,и,и)Р{<1ч) - ЗЯ^! + 11 н1°\х)(и,и,и,у,ь,ь)Р(ди)Р(^)+

ЕЛ Е*

+ 1I ^\х){щщч,щи)Р{ди) - ^ I Н?\х)(щи,и)Р(с1и)+ ЕЛ Е*

Е* Е<1

/Н?\х)(и,и,и)Р(с1и)+ Е*

Е* ЕЛ Е*

где В. = О(^), п —> оо, - остаточная часть, для которой выписывается явная оценка.

Главная часть разложения из теоремы 2.6 отличается от главной части разложения теоремы 2.4 тем, что в нее добавились величины, связанные с Я® и за счет этого улучшается скорость стремления к нулю остаточной части.

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

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

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

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

Д.(е) = J(e,u)s(P — Ф)(с?и);

Е"

• слагаемые из £, убывающие при росте п как (9(n;m)1)/2), включаем в остаточную часть разложения, оставшиеся слагаемые - в главную.

Представление разности плотностей с помощью As(et) дает

Лемма 3.1 Во введенных предположениях и обозначениях при п ^ тах(^, 2т) существует плотность рп(х) и для любого х € Ed справедливо равенство

Рп{х) - <р{х) =

(2тг)' ]е Тл Ч^з 51

где et = ий- остаточная часть, для которой выписывается явная оценка.

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

t л/п

dt + R,

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

Рп(я) = ¥>(¡0+

т • 1

П! 1

<р( х)

+£ £

(п — £)!&з! - • • &т+1! (З!)''3... ((тп +1)!)*™+1 п(3&+-+(т+1)*:т+,)/2

ууИ У (ГУ [ ... [ н(2^(и{3] и(3) и(3) ,/.(3) и{3] и(3) и(т+1) 3=° Е* Е*

• • ■. с^т^) • • • • • • ол.<1и{г1)) ■ ■ ■ я(<1и{^)+д

где <3 = Р — Ф и внутреннее суммирование ведется по всем наборам ..., таким, что

к3 Н-----1- = А;, З^з Н-----1- (т + 1)&т+1 ^т + 2к,

"т к3 Н-----\-(т- 1)кт+1

2

«7 =

|Л| сГ1Вз(т+1^т+1) вг1*

П 2

Рт+2 + /?т+2 Г1 ^3(1-1)+т+2,п-1

"Т / . ■ / у °П 3(1-1)+т+2 ^Г "Г"

(т + 2)! ^

п

№Т

С1

С1

41

е г

О

1=1

1=1

П 2

+ Сп-^у1ЬЗт+1-ЛТ^/п),

где для действительных а, Ь, Т\ обозначено Ьа,ь{т1) = ГЛ, п / е" а г <1г а для натурального I и положительного действительного Т

т+З г п»+1 г

з=3

в=3

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

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

3 Благодарности

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

4 Список публикаций автора по теме диссертации

[1] Осмоловский И.Ю. "О некоторых свойствах многомерных аналогов многочленов Чебышева-Эрмита", Теория вероятностей и ее применения, 2008, т. 53, в. 2, сс. 373 - 378.

[2] Осмоловский И.Ю. "Об оценке точности аппроксимации для асимптотических разложений в многомерном случае", Теория вероятностей и ее применения, 2009, т. 54, в. 1, сс. 152 - 158.

[3] Осмоловский И.Ю. "Об оценке точности аппроксимации для асимптотических разложений в многомерном случае", Деп. В ВИНИТИ 24.03.2009 № 152-B2Û09.

Подписано в печать 2. ?■ 09 Формат 60x90 1/16. Усл. печ. л. /.25 Тираж 100 экз. Заказ 23

Отпечатано с оригинал-макета на типографском оборудовании механико-математического факультета МГУ имени М. В. Ломоносова

 
Содержание диссертации автор исследовательской работы: кандидата физико-математических наук, Осмоловский, Игорь Юрьевич

Введение

1 Многомерные аналоги многочленов Чебышева-Эрмита и их свойства

1.1 Многомерные аналоги многочленов Чебышева-Эрмита. Свойства и примеры аналогов многочленов Чебышева-Эрмита.

1.2 Вспомогательные утверждения.

1.3 Доказательство рекуррентной формулы для многочленов Н^\з = 1,2,.

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

2.1 Основные обозначения и предположения.

2.2 Постановка задачи.

2.3 Две леммы.

2.4 Разложения для плотностей при конечности моментов порядков 5 и

3 Асимптотические разложения для плотностей в общем случае

3.1 Асимптотические разложения для плотностей в общем случае при конечности момента порядка га ^ 3.

3.2 Формулировка и доказательство леммы

4 Асимптотические разложения для решетчатых распределений

4.1 Асимптотические разложения в локальной ЦПТ.

5 Асимптотические разложения для вероятностей

5.1 Получение разложений для вероятностей с помощью разложений для плотностей.

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

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

Пусть - н.о.р.с.в. с 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 > 160000. В случае, когда нормированная сумма состоит из нескольких десятков слагаемых, оценка теоремы Берри-Эссена, по существу, бессодержательна.

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

2) а пример оценки, содержащей псевдомоменты, дает неравенство p(Fn, Ф) < при 4, (3) у/П где С - некоторые постоянные, a v^ - метрика на множестве функций распределения, которая называется вариацией с весом (здесь вес равен |я|3) и определяется равенством оо оо u3(V, W) = J \x\3\d(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{|/(í)|:í>T}<l для некоторого Т > 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) = <p(x) + J2 -^-ЩхЫх) 1=3 L и числа 9i(qn) достаточно просто вычисляются через числа 61

ЧЯп) у^ п\ ( в3 \кз (Jh ^ l\ ^ k0\k3\.ki\ \3\п3/2 J "\l\n1/2, для 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)\ ~ j\2?(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 j\2i(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^3/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 Н^0)(х)(и, и, и, v, v, и)р(сь)р(е^)+ I н(°\х)(и,щщи,и)Р(аи)-^^ I нР(х)(и,и,и)Р(с1и)+

120п3/2 р{х)

УIН$°\х)(и,и,и,у,у,у,у)Р{ёи)Р№)-1^ I Н<£\х){и,и,и)Р{<1и)+

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 пришлось ограничиться лишь одним результатом, тем более, что обобщения этого результата на общий случай достаточно очевидны.

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