Главная » Просмотр файлов » Соболь И.М. Численные методы Монте-Карло (1973)

Соболь И.М. Численные методы Монте-Карло (1973) (1186217), страница 44

Файл №1186217 Соболь И.М. Численные методы Монте-Карло (1973) (Соболь И.М. Численные методы Монте-Карло (1973)) 44 страницаСоболь И.М. Численные методы Монте-Карло (1973) (1186217) страница 442020-08-26СтудИзба
Просмтор этого файла доступен только зарегистрированным пользователям. Но у нас супер быстрая регистрация: достаточно только электронной почты!

Текст из файла (страница 44)

..., хь ..., припадле'каших интервалу (О, !), называется вполне равнол!арно распределенной, если при каждом натуральном и последовательность точек (хг,..., х„), (х,..., х„г!), (хз,..., х„.!.7), ... (21) равномерно распределена в К". Из этого определения и п. 2.! следует, что если Ф(уь ..., у.) ннтегрируема по Риману в К", то, каково бы ни было натуральное а, А — ! !Ип — ~ Ф(хг-с1, °,хг+ ) = ) Ф(Р)г(Р. (22) Может показаться, что расчет по формуле (22) выгоднее, чем расчет по формуле (20), так как при этом затрачивается примерно в л раз меньше чисел хь Несомненно однако, чтс вычислители, инстинктивно стремясь к «независимости» псевдослучайных точек, отдадут предпочтение точкам (!9) и формуле (20).

Вполне вероятно, что они правы: может оказаться, что для точек (21) отклонение Р (при любом п?) растет быстрее, чем для точек (!9). Вопрос этот пе исследован. ') таким свойс!вом обладают почти все а)1, $ 31 УНИВЕРСАЛЬНЫЕ ПСЕВДОСЛУЧАЙНЫЕ ЧИСЛА 973 3.3. Асимптотически вполне равномерно распределенные последовательности чисел. Рассмотрим семейство последовательностей х>, ..., х„..., зависящее от нату. ральпого параметра д(д= 1, 2,...): х>(у), х,(д)...,, х,(д)..., (23)' Семейство (23) естественно назвать асимптотически вполне равномерно распределенным, если для любой функции Ф(уь ..., у„), интегрируемой по Рнману В К", при любом натуральном и справедливо равенство 1 д н †! 1пп Ип! — ~'. Ф (х>и, (д),..., х;„;, (д)) = ) Ф (Р) дР, Иех,а >О ю' 'и >св (24) илн равенство А — ! 1!!и 1пп — „, ~, Ф(х>,(у),...,.к>+в(у)) =-) Ф(Р) дР. (25) и н- ~ !о кп Термин этот введен И.

Франклином в !963 г. 1126). Ев!у же принадлежит второе утверждение следующей теоремы, первое утверждение которой доказано в 133). Т е о р е м а. Существует бесконечно многа чисел а, 0(а~1, таких, что семейство последовательностей х,(д) =Д(ад'), 1=1, 2, ... ву асимптотически вполне равномерно распределено как в смысле (24), так и в смысле (25)*). Интересно отметить, что при каждом фиксированном д последовательность х,=Д(ай') может быть равномерно распределенной в К!, но точки (хв>-ь хм) не могут быть равномерно распрел ! х ! деленными в Кв. В самом деле, Рис 69. точки (хь хт), (хм хи) ... можно записать в виде (х>, Д(дх!)), (хм Д(дхв)), Отсюда видно, что все они расположены на линии у=Д(ух), которая состоит из д отрезков прямых у=ух — й при * Таким свойством обладают сочти все и ив интервала (О, 1), 2?4 ' !ГЛ 7 НЕСЛУЧАИНЫЕ ТОЧКИ В АЛГОРИТМАХ й 4.

Проверка псевдослучайных чисел с детерминистической точки зрения 4.1. Интерпретация простейших критериев. В гл 1, 4 3 чы рассмотрели некоторые тесты, используемые лля проверки псевдослучайных чисел, и отметили, что все такие тссты только необходимы, ио не достаточны: они могут опровсргнуть гипотезу о том, чзо у„..., ул! — независимые значения случайной величины у, но ие могут доказать ее. Сузим теперь постановку вопроса: вместо «можно ли числа у„ ..., узг считать независимыми значениями у?» мададнм вопрос; можно ли числа у!, ° ° °, уз! использовать вмесго независимых значений у в формуле ! и ! Ф(х) !(х — ~ы Ф(у,) 'о 7=! (26) пра любых Ф(х) из достаточно широкого класса? При такой постановке вопроса рассмотренные в гл, 1 критерии оказываются и необходимыми н достаточными — каждый на своеи классе функций Ф(х).

Впервые такая интерпретация была изложена автором в 1968 г. 4.1.1. Критерий ы'. Из формулы (48) стр. 129 вытекает, что ! А' зпр ~ ) Ф (х) с(х — — 7, Ф(у.)! = й шее' !ц ! а 2 (2?) Следовательно, малость величины ыш нсобходпзза и достаточна лля 2 того, чтобы ошибка приближения (26) была мала для всех ф)акций класса (Р2 (1.). 4Л.2. Критерий Хе. Разобьем интервал (О, 1) на г интервалов А',+... + л, = (О, !), клипы которых равны рп ., ., р, так, что Р!+. ° .

+ Рг = ! и все р. )О. Обозначим черсз т. колнче. у ? а=6„„, Аг — ! (рис. 69). Тем более последовательность х, пе может быть вполне равномерно распределенной. Сформулированная теорема в какой-то мере объясняет, почему метод сравнений (п. 2.2.2 гл. !) позволяет строить «хорошие» псевдослучайные числа и почему приходится выбирать большие йч из (7) гл. ! вытекаст, что у!=Д(усу*). Однако никаких выводов о том, как выбирать «хорошие» Аг, М, то отсюда сделать нельзя. Как доказал С. М. Е!рмаков !'ЗЗ),скорость сходимости в формуле (24) (но не в (25)!) связана с у!ОФ(!")?Л?.

й 41 провеикл псевдослтчлиных чисел 275 ство значений у, прп 1(!<Ь', принадлежащих ХХ Критерий хз основан на асимптотвческом распределении величины !'=! l (28) которое не зависит от Хь ...,Х, а только от числа степеней свог' боды г — !. Введем множество функций А,(Е), постоянных на каждом из Х1, т, е, фуакций вида Ф(х) = =с. при х с Х, 1 < ! < г, l' п таких, что Г рс ! ! !'=! (чр) Ф(х)шА!(ь), то г г 1(егрздпо вычпсшпь, что если ! ~Фг(х — '~', Ф(у!) = Ят', Е й(гг 1! Используя известное неравенство(2!! л.)з < 'и!2о, получим оценку ч !! !' у.л! причем оценка зта гочная, так как прп последнее нсравсисзво обращается в равенство. )!таге, и зо р ~ ( Ф (х) г(х — — х~л Ф (у.) = й 1 св жгс1 'о !У,—.! Л~ (30) Отсюда видно, что малость величины 2мд необходима н достаточна для того, чтооы ошибка приближения (26) была мала для всех функций нласса А!(Ь), Мы.

знаем, олнако, что к формуле (26) сводятся лишь такие ал!орит!гы Моитс-Карло, у которых к.р.=!. Стало быть, расспотренныс критерии (как и критерий Колмогорова, см. упражнение 8 НЕСЛУЧА!ЧНЫЕ ТОЧКИ В АЛГОРИТМАХ 27Б !Гл т гл 1 и упражнение 6 гл 71 гарантируют применимость проверенных псездослучаииых чисел у!,..., Ум только для расчетов по алто. ритмам с к. Р.=!. 4.2. Миогомериые критерии.

Если мы хотим использовать зтч же числа для расчета Во алгоритмам с к. р.=л, то необходимо проверить вх с помощью кзкого-либо л.мерного теста. Например, разобьем куб К" иа г областей Х!+... + Х„=-К", объемы которых равны р!, ... р, так что рг+...+Р =1 и все р )О. Из чисел ! уь ° ° ° ум образуем А!„=-П(АГ7п) точек с координатами "' Тл)' (Уп+!' "" 12л)' '"' (Уп!мл — !Н.!' ' ' Улмп)' и пусть т — количество таких точек, принадлежащих Х . ! Тот же крвтерий 7!з с (г — 1)-й степенью спободй позволяет проверить распределение зтих точек вКп с помощью величины Обозначим через Ал!ь) множество функций Ф(х, ..., хп), постоянных на кап!дом Х7, т.

е. функций вида Ф(х,,...,х)мс прп (то...,х)цХ, 1<)(г, и таких, что выполнено условие (29). Повторяя вычисления п. 4.!.2, получиы, что Лп — ! зпР,~ Ф(Р) г1Р— А! л, Ф (Уж+2, ..., У, ) = й')гг:м Фщил!С! Ап (31) Такпч обРазом, малость Ух действительно гаРантиРУет пРиыенн. .2 л!ость чисел уь..., Ум лля расчета по некоторым алгоритмам с к. р.=л. Конечно, класс алгоритмов, которым отвечают функшш Ф (х„..., х ) из А„(ь) весьма узок. Однако если г велвко и все Х.

малы, то многие функции допускают хорошее приближение кусочно постояннымн функциями из А (Е). И тогда псевдослучайл пые числа, пригодные для интегрирования функций из А„1ь) окажутся пригодными для интсгрпровгиия гораздо более разнообраз. ных функппй. 4.3. О системах тестов. С точки зрения изложенного в пп. 4.1 и 4.2, система тестов, рекомендованная в гл. 1, представляется логичной, но довольно слабой.

В самом деле, она включает простейшие проверки одномерного и двумерного распределения у и заодно распределения некоторых попутно вычисляемых второстепенных величин. Лишь проверка серий моуксы рассматриваться как тест, упплжнкния к Гл. т Упражнения к главе 7 1. Рассмотрим пример п 1.2, где значения $ вычисляются ме.

тодом Неймана. Ограничим число проб величиной ! и условимся полагать Э=О, если су!)р!а~), ..., суг> р(а!). Для такого алго итма к. р.=2!. оказать, что для того, чтобы погрешность такого приближения не превосходила вероятной ошибки метода гм — 0,675Р')»/!4)//»' должно выполняться неравенство ! > 1п (г,//) 1п ' (1 — з). 2. Выберем в кубе /т" равпомсрну!о кубическую сетку, состоя- и ую нз !У вЂ” 31® точен с кооРлинатами — 1/2 ) л М [ !» — 1/2 где 1ь..., 1 =1, 2,..., 34.

/!оказать, что отклонение этой сетки равно 0„=05л/~ '" (при л=1 величпнабн 05 минималь. па, но прп больших л отношение Рм/й/ убывает крайне медленно) 3. Вычислить точку!3!з в 4-мерном кубе, Ответ: !/!з — — !!!Пб, 13/16, !3/16, 15/!6). 4. Используя равенства 1'1!' — — 2 ' при з= 1, 2, ..., доказать, что пеРвые кооРдинаты 4! ! точек !е! совпадают с числами Р»Я. рассчитанный на алгоритмы с к.р, еь, правда, алгоритмы весьма спепиального вила Если бы речь шла только об использовании алгоритмов с любыми к.

р.<л, то можяо было бы предложить очень хороший тест— расчет отклонения/)м или, другими словами, л-ыерный критерий Колмогорова !см. упра>кнсние 8 гл. 1 и упражнение 6 гл. 7). Правда, вычисление /тм при больших /У очень трудоемко. Обеспечить же знпирической проверкой «универсальность» псевдослучайных чисел, по-видимому, крайне трудно, Например, без теоретических результатов, подобных результатам 4 3, пришлось бы проверять л-мерное распределение не один раз, а л раз, так как (Уы ...,7„), (Ул)!, ..., У „)... хорошо рас. пределены, не следует, что точки(у/,...,7;+„) (у;+„+1, ..., у;+ „),... при 1<!<и также хороши.

Это видно пз следуюшего простого примера: в последовательности 11110000!11!0000 ... пары 11, 11, 00, 00. 11, 11, 00, 00 распределены неравномерно (так как комбинапнй О! и !О совсем нет); отбросив первую единицу последовательности, получим пары !1, 1О, 00, 01, 11, 10, 00, 01, ..., где каждая пз возможных комбинаппй встречается одинаково часто. Результаты этого параграфа содержатся в статье [34). 2ТЕ НИОЛУИАЙИЫН ТОЧКИ Н АЛГОРИТМАХ [ГЛ. т б.

Какова бы пи была последовательность чисел кь..., хг ... нз интервала (О, 1), равенство 1 ) 1' (х) г(х — — У„) (х.) = о ~ —,~ о — — (А~ длп 1 х и ) хз одновременно невозможно. (И. Н Ченпов). 6. Обозначим (Р', ((.) множество непрерывных функций 1'(х1, определенных при 0<к~! и пмг;гщпх нусочно нспрсрывные производные )'(х) такие, что 1 ) ()'(л)) ало Доказать, что лля любых чисел уь..., Тх 11 зпр ~~Ф(х) г(х — ~~ Ф(у.)~ = (.

ФаЗР (А) О й! " ФА)' где Км — — У)У(з — величина, Фигурирующая в критерии Колмогоро- ва (упражнение 8 гл 1) (И. б!. Соболь [84)). ГЛАВА 8 НЕКОТОРЫЕ ДРУ1'ИЕ ЗАДАЧИ 9 1. Интерполирование функций от большого числа переменных Рнс Уо В книге уже встречались задачи, в которых методы Монте-Карло оказываются эффективнее классических методов при большом числе переменных (см.упражнение 9 гл.

3 и п. 5.6 гл. 5). Здесь П,1/ изложена одна задача такого типа, рассмотренная впервые Д ж. Х з м м е р с л и 1134) . 1.1. Постановка задачи. Рассмотрим функцию )(хь - гх,„гА х„), значения которой -------+--- известны лишь в вершинах единичного и-мерного куба К". Требуется ироннтерполи- 4Ф ЙИ ровать значение втой функции в точке (хь ..., х„), расположенной внутри К" (рис. 70). Интерполяция линейная по каждому из переменных. В случае а= ! интерполяцпонная форзгула всем хоро. шо знакома: !(х1) = (1 — х1)!(О)+хи(1). !1етрудно проверить, что при а=2 получается аналогичная формула: )(хь лг) = (1 — х,) (1 — х,)!(О, О)+х,(! — хз))(1,0)+ +(1 — хДхз1(0, 1)+х1хз~(1, 1).

Характеристики

Тип файла
DJVU-файл
Размер
6,74 Mb
Тип материала
Высшее учебное заведение

Список файлов книги

Свежие статьи
Популярно сейчас
Зачем заказывать выполнение своего задания, если оно уже было выполнено много много раз? Его можно просто купить или даже скачать бесплатно на СтудИзбе. Найдите нужный учебный материал у нас!
Ответы на популярные вопросы
Да! Наши авторы собирают и выкладывают те работы, которые сдаются в Вашем учебном заведении ежегодно и уже проверены преподавателями.
Да! У нас любой человек может выложить любую учебную работу и зарабатывать на её продажах! Но каждый учебный материал публикуется только после тщательной проверки администрацией.
Вернём деньги! А если быть более точными, то автору даётся немного времени на исправление, а если не исправит или выйдет время, то вернём деньги в полном объёме!
Да! На равне с готовыми студенческими работами у нас продаются услуги. Цены на услуги видны сразу, то есть Вам нужно только указать параметры и сразу можно оплачивать.
Отзывы студентов
Ставлю 10/10
Все нравится, очень удобный сайт, помогает в учебе. Кроме этого, можно заработать самому, выставляя готовые учебные материалы на продажу здесь. Рейтинги и отзывы на преподавателей очень помогают сориентироваться в начале нового семестра. Спасибо за такую функцию. Ставлю максимальную оценку.
Лучшая платформа для успешной сдачи сессии
Познакомился со СтудИзбой благодаря своему другу, очень нравится интерфейс, количество доступных файлов, цена, в общем, все прекрасно. Даже сам продаю какие-то свои работы.
Студизба ван лав ❤
Очень офигенный сайт для студентов. Много полезных учебных материалов. Пользуюсь студизбой с октября 2021 года. Серьёзных нареканий нет. Хотелось бы, что бы ввели подписочную модель и сделали материалы дешевле 300 рублей в рамках подписки бесплатными.
Отличный сайт
Лично меня всё устраивает - и покупка, и продажа; и цены, и возможность предпросмотра куска файла, и обилие бесплатных файлов (в подборках по авторам, читай, ВУЗам и факультетам). Есть определённые баги, но всё решаемо, да и администраторы реагируют в течение суток.
Маленький отзыв о большом помощнике!
Студизба спасает в те моменты, когда сроки горят, а работ накопилось достаточно. Довольно удобный сайт с простой навигацией и огромным количеством материалов.
Студ. Изба как крупнейший сборник работ для студентов
Тут дофига бывает всего полезного. Печально, что бывают предметы по которым даже одного бесплатного решения нет, но это скорее вопрос к студентам. В остальном всё здорово.
Спасательный островок
Если уже не успеваешь разобраться или застрял на каком-то задание поможет тебе быстро и недорого решить твою проблему.
Всё и так отлично
Всё очень удобно. Особенно круто, что есть система бонусов и можно выводить остатки денег. Очень много качественных бесплатных файлов.
Отзыв о системе "Студизба"
Отличная платформа для распространения работ, востребованных студентами. Хорошо налаженная и качественная работа сайта, огромная база заданий и аудитория.
Отличный помощник
Отличный сайт с кучей полезных файлов, позволяющий найти много методичек / учебников / отзывов о вузах и преподователях.
Отлично помогает студентам в любой момент для решения трудных и незамедлительных задач
Хотелось бы больше конкретной информации о преподавателях. А так в принципе хороший сайт, всегда им пользуюсь и ни разу не было желания прекратить. Хороший сайт для помощи студентам, удобный и приятный интерфейс. Из недостатков можно выделить только отсутствия небольшого количества файлов.
Спасибо за шикарный сайт
Великолепный сайт на котором студент за не большие деньги может найти помощь с дз, проектами курсовыми, лабораторными, а также узнать отзывы на преподавателей и бесплатно скачать пособия.
Популярные преподаватели
Добавляйте материалы
и зарабатывайте!
Продажи идут автоматически
6384
Авторов
на СтудИзбе
308
Средний доход
с одного платного файла
Обучение Подробнее