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

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

Восстановить пароль Регистрация
 
litviak
0 / 0 / 0
Регистрация: 01.11.2012
Сообщений: 6
01.11.2012, 01:25     Метод-Монте карло. Не работает генератор случайных чисел #1
Аглоритм простой:
Координаты - случайные величины.
При каждой итерации проверяется пододает ли случайная точка из куба 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++
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
01.11.2012, 01:25     Метод-Монте карло. Не работает генератор случайных чисел
Посмотрите здесь:

метод Монте-Карло C++
C++ метод Монте-Карло
C++ Метод Монте-Карло
C++ Неправильно работает генератор случайных чисел
C++ Метод монте Карло
После регистрации реклама в сообщениях будет скрыта и будут доступны все возможности форума.
Kuzia domovenok
 Аватар для Kuzia domovenok
1883 / 1738 / 116
Регистрация: 25.03.2012
Сообщений: 5,907
Записей в блоге: 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 один раз в начале программы.
litviak
0 / 0 / 0
Регистрация: 01.11.2012
Сообщений: 6
01.11.2012, 11:40  [ТС]     Метод-Монте карло. Не работает генератор случайных чисел #3
Я так сделал что бы делать каждое испытание незаисимым друг от друга, какая разница сколько переинициализировать генератор, я же просто меняю его начальную точку?
Проверил, разницы никакой...
kravam
быдлокодер
 Аватар для kravam
1512 / 872 / 44
Регистрация: 04.06.2008
Сообщений: 5,271
01.11.2012, 12:12     Метод-Монте карло. Не работает генератор случайных чисел #4
Цитата Сообщение от litviak Посмотреть сообщение
Проверил, разницы никакой...
потому что дело не в этом. ПОКА не в этом.
Сначала напиши
C++
1
            z=rand()/(float)RAND_MAX;
потом включи srand ОДИН раз, вот так:
C++
1
2
3
int main()
{
    srand(time(NULL));
И всё будет круто.
litviak
0 / 0 / 0
Регистрация: 01.11.2012
Сообщений: 6
01.11.2012, 12:56  [ТС]     Метод-Монте карло. Не работает генератор случайных чисел #5
Я понял, спасибо, круто, даже и не подумал бы что оно вот так.
Yandex
Объявления
01.11.2012, 12:56     Метод-Монте карло. Не работает генератор случайных чисел
Ответ Создать тему
Опции темы

Текущее время: 06:52. Часовой пояс GMT +3.
КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2016, vBulletin Solutions, Inc.
Рейтинг@Mail.ru