Форум программистов, компьютерный форум CyberForum.ru
Наши страницы

С++ для начинающих

Войти
Регистрация
Восстановить пароль
 
litviak
0 / 0 / 0
Регистрация: 01.11.2012
Сообщений: 6
#1

Метод-Монте карло. Не работает генератор случайных чисел - C++

01.11.2012, 01:25. Просмотров 800. Ответов 4
Метки нет (Все метки)

Аглоритм простой:
Координаты - случайные величины.
При каждой итерации проверяется пододает ли случайная точка из куба 1*1*1 в заданную область или нет.
Объем области считается как процент попавших случайных величин(чем больше итераций тем точнее соответственно).
Собственно код и вывод:
C++
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
#include <math.h>
#include <iostream>
#include <vector>
#include <stdlib.h>
#include <time.h>
#include <string>
 
using namespace std;
 
class monteCalroMethod
{
 
    vector <int> Tests;
    int numberOfTest;
    
    bool arc1(float x,float y,float z)
    {
        if((x*x+y*y+z*z)<=1)
            return 1;
        else
            return 0;
    };
    bool arc2(float x,float y,float z)
    {
        if(z>=0.6)
            return 1;
        return 0;
    };
    bool onceTest(float x,float y,float z)
    {
        if(arc1(x,y,z)&&(arc2(x,y,z))) return 1;
        else return 0;
    };
    int posSize()
    {
        int num;
        for(int i=0;i<numberOfTest;i++)
            if(Tests[i])num++;
        return num;
            
    }
    
public:
    monteCalroMethod(int n=1)
    {
        numberOfTest=pow(10,n);
        Tests.resize(numberOfTest,0);
    };
    void calculatePlanArea()
    {
        for(int i=0;i<numberOfTest;i++)
        {
            float x,y,z;
            srand(time(NULL));
            x=rand()/RAND_MAX;
            y=rand()/RAND_MAX;
            z=rand()/RAND_MAX;
            if(onceTest(x,y,z))
                Tests[i] = 1;
            cout<<z<<endl;
        }
        float arc = posSize()/numberOfTest*4;
        cout<<"Area = "<<arc<<endl;      
    };
    void Change()
    {
        cout<<"not today";
    };
    
};
 
int main()
{
    int num=3;
    char answ;
    cout<<"Calculator Volume"<<endl;
    monteCalroMethod mont;
    cout<<"Now functions are: x^2+y^2+z^2=1"<<endl;
    cout<<"                 z>=0"<<endl;
    cout<<"Your want change?              (y,n)"<<endl;
    cin>>answ;
    if(answ=='y'){
        mont.Change();
    };
    mont.calculatePlanArea();
    return 0;
}
Вывод:
Код
Calculator Volume
Now functions are: x^2+y^2+z^2=1
                 z>=0
Your want change?              (y,n)
n
0
0
0
0
0
0
0
0
0
0
Area = -7.61152e+08
Kirills-MacBook-Air:Documents duckie$
Компилирую через g++
0
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
01.11.2012, 01:25
Здравствуйте! Я подобрал для вас темы с ответами на вопрос Метод-Монте карло. Не работает генератор случайных чисел (C++):

метод Монте-Карло - C++
Трехмерное тело образовано объединением нескольких сфер произвольного размера и взаимного расположения. Найти объем этого тела, используя...

Метод Монте-Карло - C++
Помогите написать программу для вычисления определенного интеграла методом Монте-Карло.

метод Монте-Карло - C++
всем привет, у меня вопрос по методу Монте - Карло, у меня есть код, #include &lt;stdlib.h&gt; #include &lt;iostream&gt; #include &lt;time.h&gt; ...

Метод Монте-Карло - C++
как мне перевести этот код на с++..... помогите пожалуста(( program MonteKarlo; uses crt; Label l1,l2; var ...

Метод Монте-Карло - C++
Помогите написать программу на С++. Из множества 1, 2, ..., n, случайным образом извлекается с возвращением по одному числу S раз. Пусть...

Метод монте Карло - C++
Здравствуйте! Посмотрите, пожалуйста, почему программа неправильно считает определенный интеграл методом Монте Карло. Считаю интеграл от -3...

Надоела реклама? Зарегистрируйтесь и она исчезнет полностью.
Kuzia domovenok
1891 / 1746 / 118
Регистрация: 25.03.2012
Сообщений: 5,925
Записей в блоге: 1
01.11.2012, 03:12 #2
C++
1
2
3
4
for(int i=0;i<numberOfTest;i++)
       {
            float x,y,z;
            srand(time(NULL));
это изначально так было? зачем вызывать srand() перед каждой генерацией сл. чисел? Возможно, ошибка здесь. вызови srand один раз в начале программы.
0
litviak
0 / 0 / 0
Регистрация: 01.11.2012
Сообщений: 6
01.11.2012, 11:40  [ТС] #3
Я так сделал что бы делать каждое испытание незаисимым друг от друга, какая разница сколько переинициализировать генератор, я же просто меняю его начальную точку?
Проверил, разницы никакой...
0
kravam
быдлокодер
1695 / 882 / 45
Регистрация: 04.06.2008
Сообщений: 5,459
01.11.2012, 12:12 #4
Цитата Сообщение от litviak Посмотреть сообщение
Проверил, разницы никакой...
потому что дело не в этом. ПОКА не в этом.
Сначала напиши
C++
1
            z=rand()/(float)RAND_MAX;
потом включи srand ОДИН раз, вот так:
C++
1
2
3
int main()
{
    srand(time(NULL));
И всё будет круто.
0
litviak
0 / 0 / 0
Регистрация: 01.11.2012
Сообщений: 6
01.11.2012, 12:56  [ТС] #5
Я понял, спасибо, круто, даже и не подумал бы что оно вот так.
0
MoreAnswers
Эксперт
37091 / 29110 / 5898
Регистрация: 17.06.2006
Сообщений: 43,301
01.11.2012, 12:56
Привет! Вот еще темы с ответами:

вроде метод монте карло - C++
распишите пожалуйста что делает это программа?? int i,a,b,n,k,c,d,e,f; double s1,s,x,y,z; int _tmain(int argc, _TCHAR* argv) { ...

Не работает генератор случайных чисел - C++
#include &lt;iostream&gt; #include &lt;stdlib.h&gt; #include &lt;time.h&gt; using namespace std; int main(){ int a; srand(time(0));...

Метод Монте-Карло. Объем сферы - C++
Задан радиус сферы и количество итераций. Найти объём сферы методом Монте-Карло и определить точность

Метод Монте Карло (неправильные значения) - C++
Пишу программу для подсчета определенного интеграла 3мя способами: прямоугольников, трапеций и Монте Карло. Первые 2 реализованы правильно,...


Искать еще темы с ответами

Или воспользуйтесь поиском по форуму:
Yandex
Объявления
01.11.2012, 12:56
Ответ Создать тему
Опции темы

КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2017, vBulletin Solutions, Inc.
Рейтинг@Mail.ru