Научная статья на тему 'О случайных пуассоновских заполнениях ячеек'

О случайных пуассоновских заполнениях ячеек Текст научной статьи по специальности «Математика»

CC BY
99
22
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
РАЗМЕЩЕНИЕ ЧАСТИЦ ПО ЯЧЕЙКАМ / ПРЕДЕЛЬНОЕ РАСПРЕДЕЛЕНИЕ / ЧИСЛО ЯЧЕЕК ЗАДАННОГО ОБЪЕМА / PARTICLE ALLOCATION TO CELLS / LIMIT DISTRIBUTION / NUMBER OF CELLS WITH A GIVEN SIZE

Аннотация научной статьи по математике, автор научной работы — Хворостянская Елена Владимировна

Рассматривается модель заполнения N различных ячеек частицами, в которой случайные величины η 1,..., η N, равные числу частиц в ячейках, независимы и имеют распределение Пуассона. Для подмножества реализаций такой схемы, удовлетворяющих условию η 1 +... + η N ≤ n, при N → ∞ то получены предельные распределения числа ячеек μ r, содержащих ровно r частиц.

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

ON RANDOM POISSON ALLOCATION TO CELLS

We study a model for particle allocation to N different cells in which the random variables η 1,..., η N equal to the number of particles in cells are independent and have a Poisson distribution. For a subset of realizations of such a scheme satisfying η 1 +... + η N ≤ n limit distributions of cell number μ r with exactly r particles are obtained as N → ∞.

Текст научной работы на тему «О случайных пуассоновских заполнениях ячеек»

Труды Карельского научного центра РАН № 1. 2013. С. 112-116

УДК 519.2

О СЛУЧАЙНЫХ ПУАССОНОВСКИХ ЗАПОЛНЕНИЯХ ЯЧЕЕК

Е. В. Хворостянская

Институт прикладных математических исследовании Карельского научного центра РАН Петрозаводский государственный университет

Рассматривается модель заполнения N различных ячеек частицами, в которой случайные величины ni,... ,Vn, равные числу частиц в ячейках, независимы и имеют распределение Пуассона. Для подмножества реализаций такой схемы, удовлетворяющих условию ni +... + Vn ^ n, при N ^ то получены предельные распределения числа ячеек лr, содержащих ровно r частиц.

Ключевые слова: размещение частиц по ячейкам, предельное распределение, число ячеек заданного объема.

E. V. Khvorostyanskaya. ON RANDOM POISSON ALLOCATION TO CELLS

We study a model for particle allocation to N different cells in which the random variables пъ ..., nN equal to the number of particles in cells are independent and have a Poisson distribution. For a subset of realizations of such a scheme satisfying ni + ... + nN ^ n limit distributions of cell number лr with exactly r particles are obtained as N ^ то.

Key words: particle allocation to cells, limit distribution, number of cells with a given size.

Пусть п различных частиц независимо одна от другой размещаются по N ячейкам, занумерованным числами от 1 до N, так, что каждая частица с вероятностью N-1 попадает в одну из ячеек. Такая схема размещения частиц по ячейкам носит название классической схемы размещения, и ее изучению посвящено множество работ (см., например, [2, 3]). Хорошо известно, что совместное распределение заполнений ячеек в классической схеме совпадает с условным распределением N независимых пуассоновских случайных величин при условии, что их сумма равна п. Используя это свойство, при N ^ те и всех возможных значениях п, г были, в частности,

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

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

Рк = Р{п = к} = “к! е-\ к = 0,1,2,...,

г = 1,...,^ Очевидно, что число всех

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

0

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

Р{Пі = кі,... ,Пм = км} = (

дкі+...+к^ 0—ЛМ

к1!... км!

і

\

Е Е

Лге

і=0 іь ..., іи ^ 0

Іі + ... + Іи = і

їі! ...їм!

/

Введем независимые одинаково распределенные по закону Пуассона с параметром Л случайные величины Сі,..., См• Легко видеть, что

Р{Сі = кі,..., См = км| Сі + ... + См ^ п}

Лкі+...+к^е—лм ^ ™ (лж)е—лм'

Ї!

і

кі!... км!

чг=0

Отсюда и из полиномиальной формулы следует, что выполнено соотношение

Р{П1 = к1,... ,пм = км}

= Р{6 = к1,..., См = км| £1 + ... + См ^ п}.

Эта модель является частным случаем схемы, рассмотренной в [4] и основанной на идеях обобщенной схемы размещения частиц по ячейкам [2]. Заметим, что в [4] нет конкретных примеров, удовлетворяющих условиям этой схемы.

Обозначим через случайную величину,

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

п - А^ (А - г)2р

у = -ТАГ’ /г(А) = -^Т-■

Справедливы следующие утверждения.

Тогда

Р К = к} = е-мрг (1 + о(1))

равномерно относительно целых неотрицательных к, для которых (к — Npr) /у^рг лежит в любом конечном фиксированном интервале.

Теорема 2. Пусть N ^ те, 0 <А1 ^ А ^

А2 < те, г ^ 0 фиксировано и выполнено одно из следующих условий:

1) у ^ +те;

2) |у| ^ С< те, г — А ^ 0;

3) у ^ —те, г—А ^ 0 так, что у(г — А) ^ 0,

у2/л/Ж ^ 0.

Тогда

Р {^- = к} =

1 + о(1)

л/2пЖрг (1 - р-)

3—«2/2

равномерно относительно и = (к — Npr) / РУ в любом конечном фиксированном интервале. Утверждение теоремы остается справедливым, если выполнены условия теоремы 1 и Npr ^ те, а также в случае, когда г = 0, у ^ те, N ^ те, А ^ 0 так, что А^ ^ те.

Теорема 3. Пусть г = 1, п ^ 1 фиксировано, А ^ 0, N ^ те так, что NА ^ С < те. Тогда для целых неотрицательных к

Р {^- = к} =

(ЖЛ)к е

к „-мл

к!

Е

чг=0

(ЖЛ)* е

мЛ

і

(1 + о(1))

Теорема 1. Пусть N ^ те, у2\/Рг7^ + у2/г(А) ^ 0 при у ^ —те и выполнено одно из следующих условий:

1) г ^ 0, А ^ те;

2) г ^ те, 0 < А1 ^ А ^ А2 < те;

3) г ^ 2, А ^ 0;

4) г = 1, у ^ те, А ^ 0 так, что А2N ^ те.

при 0 ^ к ^ п и Р {^г = к} = 0 при к > п.

Замечание 1. Нетрудно видеть, что для нашей модели теорема 1 имеет значительно более общий характер, чем теорема 9 работы [4], поскольку в ней условия 1), 2), 4) не рассматривались, а условие 3) требует дополнительных ограничений. Теорема 2 о сходимости к нормальному закону также охватывает существенно более широкую зону изменения параметров по сравнению с теоремами 4 и 5 [4].

Согласно теореме 3 [4], для к = 0,1,..., N справедливо равенство

Р {^г = к}

(1)

"к (1 - г- )м—к

Р

(г)

м—к

113

X

X

где

См = £1 + ... + £м, См^-к = £1 ^ + ... + к,

/-(г) ^(г)

независимые случайные величины £1 ,..., £м имеют распределения

Р {С(г) = к} = Р {£1 = к| £1 = г} ,

г = 1,..., N. Для доказательства теорем 1, 2 с помощью (1) нам потребуются вспомогательные утверждения, которые приводятся ниже в леммах 1, 2.

Лемма 1. Пусть N ^ те и АN ^ те. Тогда

2/2^.

■ См-АЖ ^ 1 = 1 + о^ Г е—^2

л/т

Р

где

о>\/і

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

тг

Е Сі

\/2П

(г) Л грг

1 - Рг :

е ^2/2^г,

(2)

ог

О Сі

(г)

Л (1 - Рг - /г (Л))

(1 - рг)2

Доказательство. При Л — 0, Л -то, г ^

2 и 0 <Лі ^ Л ^ Л2 < то слабая сходимость распределения ^С«г) - тгі^ / (ог/і)

к стандартному нормальному закону получена в [1], [3]. Учитывая явный вид характеристической функции случайной величины (С«г) — тгв) /оу/в, можно показать, что утверждение леммы остается в силе и в оставшихся случаях. □

Докажем теорему 1. Поскольку рг ^ 0, выполнено соотношение

N

к)рг

рк (1 - Рг)м—к = е—мрг (1 + 0(1))

к!

Доказательство. Обозначим через ^>(£), ^(£) характеристические функции случайных величин £1, (См — АN) /\/А^. Легко видеть, что

<р(£) = ел(е*‘-1),

•0(£) = ехр | — г£/АЖ| ^м ^//ж).

Используя формулу Тейлора, при любом фиксированном £ получаем, что

1п ^) = —г£/АЖ + АN (ег*/^лм — 1)

= —£2/2 + о(1).

Отсюда и из теоремы непрерывности следует утверждение леммы. □

Лемма 2. Пусть в ^ те, Ав ^ те и при А ^ 0, г = 1 выполнено условие А2в ^ те. Тогда

(3)

равномерно относительно к, для которых и = (к — Npr) /у^Т лежит в любом конечном фиксированном интервале.

Пусть выполнено одно из условий теоремы и, кроме того, АN ^ те. С помощью лемм 1, 2 получаем, что

(г)

(4)

Р -І ОТ—к ^ п - кг Р {См < п}

/ е ^2/2^г

-----------(

у '

/ е—г2/2^

-(1 + 0(1)),

где

хг =

п - кг - тг (Ж - к)

оу \/ N - к ’

тг, оу определены в (2). Для выбранных значений к несложно показать, что

хг =

п - ЛЖ + и(Л - г)/Жрг/ (1 - рг)

/лж (1 - і—л)і/2( 1 - и/—2м)

При условиях 1)-3) теоремы /г(Л) — 0 и

Л , /г(Л) + и^Рг/Ж Хг = У 1 + --------^^--------------- + 0 (/г(Л))

2(1 - рг)

р7\' \ . -у/(Л) 8ИП(Л - г)

+0І,/^ 1 +

1 - рг

(5)

(1 + 0(1)).

Если выполнено условие 4) теоремы, то /і (Л) — 1 и

- •' ■ - /-Лі Г (■ - - €?

-1/2

(6)

-и^ /і(Л)(1 + о(1)) = /==(1 + о(1))

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

Р -І См—к ^ п - кг Р {(м < п}

= 1 + о(1).

(7)

Если |у| ^ С < то, то из (4) получаем, что

Р -І См"—к ^ п - кг Р {См < п}

1+

V

/ е—22/2^г

(1 + о(1)),

и, поскольку из (5) следует соотношение хг — у ^ 0, легко видеть, что выполнено (7).

Нетрудно проверить, что при любом а > 0

І е—^2/2^г ^ а2 | е 2/ = а 2 е—“2/2.

.] .] г2 1 + а2

а а

Используя эти неравенства и (5), получаем оценки

—І

— е—^2/2^г

Г е—22/2^г

|у|

^ 1 + у2 е—(хГ—у2)/2, УХг

где

/ е г2/2^

|Хг 1_________ > уХг е—(хГ—у2)/2

о° 1 + ™2 ,

Г е—*2/2^г 1 + х

|у|

1+^ = (1 + 1/у2) (1+ 0(1)),

^(! + «С1)),

1 + 1 + 1/у2

х2 — у2 = у2 (/г(А) + и^РТТж) (1 + о(1)).

Отсюда и из (4) следует, что при у ^ —те и

у2/г(А) + у2л/РгТ^ ^ 0 справедливо (7).

Пусть теперь АN ^ а, где а - некоторая неотрицательная постоянная. Очевидно,

что распределение суммы См сходится к закону Пуассона с параметром а. Для характеристической функции ^ (£) случайной величины >(г)

См-к выполнено соотношение 1п (£) = ^ — к) 1п

еЛ(е^* —і) - е^гРг 1 - Рг

= NА (ей — 1)(1 + о(1)),

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

и по теореме непрерывности распределение (г)

С(^_к сходится к закону Пуассона с параметром а. Поскольку Npr ^ 0, получаем, что к = 0 и справедливо равенство (7).

Утверждение теоремы 1 следует из (1), (3),

(7).

Для доказательства теоремы 2 положим к = + иу^рг (1 — рг). Пусть 0 < А1 ^

А ^ А2 < те, г ^ 0 фиксировано. С помощью лемм 1, 2 находим, что выполнено равенство (4), где

= (у (1 - Рг) + и^/г(Л) (1 - Рг) 8ИП(Л - г))

х (1 - Рг - /г(Л)) —і/2 .

у^1 - Рг - и^Рг (1 - Рг) /Ж

Аналогично доказательству теоремы 1 можно показать, что справедливо соотношение (7).

При Л — 0, г = 0 находим, что Л3(Ж - к) = Л4Ж(1 + о(1)) — то, и из лемм 1, 2 следует (4), где х0 = \/2у/Л(1 + о(1)). Очевидно, что если у — то, то Х0 — то и выполнено (7).

При Жрг (1 - Рг) — то имеет место равен-

ство

Ж

к

Рк (1 - Рг)

м—к

1+ о(1)

л/2пЖрг (1 - Рг)

е—“2/2.

Отсюда и из (1), (7) получаем утверждение теоремы 2.

Доказательство теоремы 3 аналогично доказательству теоремы 8 [4].

Работа выполнена при поддержке Программы стратегического развития Петрозаводского государственного университета на 2012-2016 гг.

ЛИТЕРАТУРА

1. Колчин А. В. Предельные теоремы для обобщенной схемы размещения // Дискретная математика. 2003. Т. 15, № 4. С. 148-157.

2. Колчин В. Ф. Случайные графы. М.: Физмат-лит, 2000. 256 с.

115

3. Колчин В. Ф., Севастьянов Б. А., Чистяков В. П. Случайные размещения. М.: Наука, 1976. 159 с.

СВЕДЕНИЯ ОБ АВТОРЕ:

Хворостянская Елена Владимировна

старший научный сотрудник, к. ф.-м. н.

Институт прикладных математических исследований Карельского научного центра РАН ул. Пушкинская, 11, Петрозаводск, Республика Карелия, Россия, 185910 эл. почта: [email protected] тел.: (8142) 781218

4. Чупрунов А. Н., Фазекаш И. Аналог обобщенной схемы размещения. Предельные теоремы для числа ячеек заданного объема // Дискретная математика. 2012. Т. 24, № 1. С. 140-158.

Khvorostyanskaya, Elena

Institute of Applied Mathematical Research, Karelian Research Centre, Russian Academy of Sciences 11 Pushkinskaya St., 185910 Petrozavodsk, Karelia, Russia

e-mail: [email protected] tel.: (8142) 781218

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