Дискуссионный математический форумМатематический форум

Математический форум Math Help Planet

Обсуждение и решение задач по математике, физике, химии, экономике

Теоретический раздел
Часовой пояс: UTC + 4 часа [ Летнее время ]
MathHelpPlanet.com RSS-лента Математического форума

Часовой пояс: UTC + 4 часа [ Летнее время ]




Начать новую тему Ответить на тему  [ Сообщений: 2 ] 
Автор Сообщение
 Заголовок сообщения: Пределы интегрирования в методе Монте Карло
СообщениеДобавлено: 29 фев 2016, 09:37 
Не в сети
Начинающий
Зарегистрирован:
25 фев 2016, 22:22
Сообщений: 4
Cпасибо сказано: 1
Спасибо получено:
0 раз в 0 сообщении
Очков репутации: 1

Добавить очки репутацииУменьшить очки репутации
Доброго дня.
Помогите разобраться с пределами интегрирования в кратном интеграле при решении методом Монте Карло. В интеграле [math]\boldsymbol{d} \boldsymbol{r}[/math] зависит от [math]\boldsymbol{d} \xi[/math] (рис.1.).
Пробовал генерировать числа в интервале в зависимости [math]\boldsymbol{d} \boldsymbol{r}[/math] (см. код программы), но это не помогло. Результат получается, как без зависимости [math]\boldsymbol{d} \boldsymbol{r}[/math] зависит от [math]\boldsymbol{d} \xi[/math] (рис.2,3).

Подскажите в чем ошибка?
Спасибо. Удачной рабочей недели.

Изображение
▼ "Код программы на C#"
z0 = 10;
ks0 = 0.142;
r0 = 1;
h = 0.01;

n = 10000;

dzn = 0; //пределы интегрирования dz
dzv = z0;

dksn = (3 * Math.PI / 2) - ks0; //пределы интегрирования dksi
dksv = (3 * Math.PI / 2) + ks0;

drn = r0 - h; //пределы интегрирования dr
drv = r0;

for (i = 1; i <= n; i++)
{
dks = dksn + (dksv - dksn) * Convert.ToDouble(rnd.Next(RAND_MAX)) / RAND_MAX; //генерация dksi

drn = (r0 - h) / Math.Cos((3 * Math.PI / 2) - dks); // генерация dr в интервале,
dr = drn + (drv - drn) * Convert.ToDouble(rnd.Next(RAND_MAX)) / RAND_MAX; // зависимом от dksi

Sum += Math.Exp((dr - ((r0 - h) / (Math.Cos((3 * Math.PI) / 2 - dks)))));
}

drn = r0 - h;
P0 = (dzv - dzn) * (dksv - dksn) * (drv - drn) * Sum / n; // итоговое нахождение результата

Вернуться к началу
 Профиль  
Cпасибо сказано 
 Заголовок сообщения: Re: Пределы интегрирования в методе Монте Карло
СообщениеДобавлено: 29 фев 2016, 13:40 
Не в сети
Начинающий
Зарегистрирован:
25 фев 2016, 22:22
Сообщений: 4
Cпасибо сказано: 1
Спасибо получено:
0 раз в 0 сообщении
Очков репутации: 1

Добавить очки репутацииУменьшить очки репутации
Не надо было генерировать в зависимости от [math]\boldsymbol{d} \boldsymbol{\xi}[/math].
Надо было отсеивать те значения, которые не попали в интервал.
▼ Код
for (i = 1; i <= n; i++)
{
dks = dksn + (dksv - dksn) * Convert.ToDouble(rnd.Next(RAND_MAX)) / RAND_MAX;
dr = drn + (drv - drn) * Convert.ToDouble(rnd.Next(RAND_MAX)) / RAND_MAX;
if (dr >= (r0 - h) / Math.Cos((3 * Math.PI / 2) - dks) && dr <= r0)
{
Sum += Math.Exp((dr - ((r0 - h) / (Math.Cos((3 * Math.PI) / 2 - dks)))));
}
}
P0 = (dzv - dzn) * (dksv - dksn) * (drv - drn) * Sum / n;

Изображение

Вернуться к началу
 Профиль  
Cпасибо сказано 
Показать сообщения за:  Поле сортировки  
Начать новую тему Ответить на тему  [ Сообщений: 2 ] 

 Похожие темы   Автор   Ответы   Просмотры   Последнее сообщение 
Метод Монте-Карло

в форуме Microsoft Excel

mariya1509

1

1005

24 апр 2013, 16:41

Метод Монте-Карло

в форуме Численные методы

galachel

4

213

24 фев 2016, 21:25

Метод Монте-Карло

в форуме Математическая статистика и Эконометрика

Yrii

2

177

13 сен 2015, 14:58

Задача по методу Монте-Карло

в форуме Исследование операций и Задачи оптимизации

swivelin

2

378

06 ноя 2014, 09:54

Метод Монте-Карло, регрессионная модель

в форуме Математическая статистика и Эконометрика

EcoFace

16

229

27 окт 2017, 02:54

Монте-Карло: решение любых систем

в форуме Размышления по поводу и без

Avgust

9

435

26 янв 2016, 06:29

Метод Монте-Карло для двойного интеграла

в форуме Численные методы

Jexio

2

101

14 фев 2018, 17:38

Нахождение площадей полигонов Вороного методом Монте Карло

в форуме Аналитическая геометрия и Векторная алгебра

kuzo

0

162

24 сен 2015, 02:37

Метод Монте Карло, разыграть значения нормально распр Сл.Вел

в форуме Математическая статистика и Эконометрика

olga21

6

932

22 ноя 2012, 14:43

Пределы интегрирования

в форуме Интегральное исчисление

Qller

4

87

05 дек 2017, 13:31


Часовой пояс: UTC + 4 часа [ Летнее время ]



Кто сейчас на конференции

Сейчас этот форум просматривают: нет зарегистрированных пользователей и гости: 1


Вы не можете начинать темы
Вы не можете отвечать на сообщения
Вы не можете редактировать свои сообщения
Вы не можете удалять свои сообщения
Вы не можете добавлять вложения

Найти:
Перейти:  

Яндекс.Метрика

Copyright © 2010-2016 MathHelpPlanet.com. All rights reserved