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

Метод Монте-Карло (вероятность попадания при стрельбе) - C++

Восстановить пароль Регистрация
 
Shelby_ph
0 / 0 / 0
Регистрация: 31.03.2014
Сообщений: 26
30.05.2014, 17:52     Метод Монте-Карло (вероятность попадания при стрельбе) #1
В общем задача такова: Добавить Метод Монтекарло, что бы показывало всевозможные попытки расчета вероятности.
Сама задача: Производится стрельба двумя снарядами по k бакам с го-
рючим (k>2), расположенным рядом друг с другом в одну ли-
нию. Каждый снаряд независимо от других попадает в первый
бак с вероятностью p1 во второй – с вероятностью р2 и т. д. Для
воспламенения баков требуется два попадания в один и тот же
бак или два попадания в соседние баки. Оценить вероятность
воспламенения баков.
Код
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
#include <iostream>
#include <vector>
#include <algorithm>
#include <numeric>
 
int main()
{
using namespace std;
int num_of_tanks;
cout « "Enter number of tanks:\n";
cin » num_of_tanks;
vector<double> tanks(num_of_tanks);
cout « "enter posibilities for each tank:\n";
for (int i = 0; i < tanks.size(); i++)
{
cout « "posibility for " « i + 1 « ") tank: ";
cin » tanks[i];
}
if (accumulate(tanks.begin(), tanks.end(), 0) > 1)
cout « "To correct ansver sum of pisibilities must be no more, than 1\n";
 
//расчет вероятности двойного попадания
double double_hit = 0;
for (int i = 0; i < tanks.size(); i++)
double_hit += (tanks[i] * tanks[i]);
 
//вероятность попадания в соседние баки
double neibour_hir = 0;
for (int i = 0; i < tanks.size() - 1; i++)
neibour_hir += tanks[i] * tanks[i + 1];
for (int i = tanks.size()-1; i > 0; i--)
neibour_hir += tanks[i] * tanks[i - 1];
 
cout « "inflammation posibility = "
« neibour_hir + double_hit « "%" « endl;
 
cin.get();
cin.get();
return 0;
}
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
30.05.2014, 17:52     Метод Монте-Карло (вероятность попадания при стрельбе)
Посмотрите здесь:

метод Монте-Карло C++
C++ метод Монте-Карло
C++ Метод Монте-Карло
C++ Метод монте Карло
C++ вроде метод монте карло
После регистрации реклама в сообщениях будет скрыта и будут доступны все возможности форума.
Kuzia domovenok
 Аватар для Kuzia domovenok
1883 / 1738 / 116
Регистрация: 25.03.2012
Сообщений: 5,907
Записей в блоге: 1
30.05.2014, 18:00     Метод Монте-Карло (вероятность попадания при стрельбе) #2
Во-первых, ты где писал этот код? У тебя IDE (студия)-то установлена?
Если установлена, откуда в коде эта мерзость?
Цитата Сообщение от Shelby_ph Посмотреть сообщение
cout « "Enter number of tanks:\n";
Если бы ты мог писать хотя бы Hello world, ты бы знал как записывается оператор , применяемый к cout! Это << а не «
Ты точно хоть одну программу, хоть хеллоуворлд сам написал?
Далее, зачем ты считаешь эту муть типа
C++
1
double_hit += (tanks[i] * tanks[i]);
тебя же чётко просили! методом Монте-Карло! По-сути, надо многократно провести численный эксперимент со стрельбой по бакам, а не вычислять эмпирически по законам теорвера!
Arkaniy
30.05.2014, 18:03
  #3

Не по теме:

Kuzia domovenok, могу предположить, что код был скинут ТСу в Вк. Там в сообщениях заменяется << на «. Сам с этим мучаюсь, когда одногруппники кидают код с просьбами о помощи...

Shelby_ph
0 / 0 / 0
Регистрация: 31.03.2014
Сообщений: 26
30.05.2014, 18:03  [ТС]     Метод Монте-Карло (вероятность попадания при стрельбе) #4
Прошу прощения. Из контакта копиривал)
Arkaniy
30.05.2014, 18:04
  #5

Не по теме:

Как в воду глядел

Shelby_ph
0 / 0 / 0
Регистрация: 31.03.2014
Сообщений: 26
30.05.2014, 18:04  [ТС]     Метод Монте-Карло (вероятность попадания при стрельбе) #6
Покажите пожалуйста как надо)
Kuzia domovenok
 Аватар для Kuzia domovenok
1883 / 1738 / 116
Регистрация: 25.03.2012
Сообщений: 5,907
Записей в блоге: 1
30.05.2014, 18:43     Метод Монте-Карло (вероятность попадания при стрельбе) #7
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
#include <vector>
#include <ctime>
#include <cstdlib>  
 
int main(){   
    int attempts, num_of_tanks;     
    std::cout<<"Input number of iterations (please more): ";
    std::cin>>attempts; 
    std::cout<<"Input number of tanks: ";
    std::cin>>num_of_tanks;
    std::vector<double> probs(num_of_tanks);
    std::vector<int>    Fbounds(num_of_tanks);
    std::cout << "enter hit probability for each tank:\n";
    for (int i = 0; i < probs.size(); ++i)
    {
        std::cout << "hit probability for (" << i + 1 << ") tank: ";
        std::cin >> probs[i];
        if (i)
            Fbounds[i]=Fbounds[i-1]+ probs[i]*RAND_MAX;
        else
            Fbounds[i]=probs[i]*RAND_MAX;
    }
    Fbounds[num_of_tanks-1]=RAND_MAX;
    srand(time(NULL));
    int triggered_counter=0;
    for (int i=0; i<attempts; ++i)
    {
        int random_key, random_hit1, random_hit2;
 
        random_key=rand();
        for(random_hit1=0; random_hit1<num_of_tanks && random_key>Fbounds[random_hit1]; ++random_hit1);  
        random_key=rand();
        for(random_hit2=0; random_hit2<num_of_tanks && random_key>Fbounds[random_hit2]; ++random_hit2);  
        if (abs(random_hit1-random_hit2)<=1)//если баки соседние
            ++triggered_counter;
    }
    double explosion= triggered_counter*1.0/attempts;
    std::cout<<"explosion probability "<< explosion*100<<"%"<<std::endl;
 
    std::cin.get();
    return 0;
}
MoreAnswers
Эксперт
37091 / 29110 / 5898
Регистрация: 17.06.2006
Сообщений: 43,301
05.06.2014, 23:04     Метод Монте-Карло (вероятность попадания при стрельбе)
Еще ссылки по теме:

Метод Монте-Карло. Объем сферы C++
C++ Метод Монте-Карло
C++ Метод Монте Карло (неправильные значения)

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

Или воспользуйтесь поиском по форуму:
Shelby_ph
0 / 0 / 0
Регистрация: 31.03.2014
Сообщений: 26
05.06.2014, 23:04  [ТС]     Метод Монте-Карло (вероятность попадания при стрельбе) #8
Пожалуйста! Поясните вот эту часть кода. 23-40. Помогите начинающему программисту)
Yandex
Объявления
05.06.2014, 23:04     Метод Монте-Карло (вероятность попадания при стрельбе)
Ответ Создать тему
Опции темы

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