Форум программистов, компьютерный форум, киберфорум
C# для начинающих
Войти
Регистрация
Восстановить пароль
Блоги Сообщество Поиск Заказать работу  
 
Рейтинг 4.86/29: Рейтинг темы: голосов - 29, средняя оценка - 4.86
1 / 1 / 0
Регистрация: 02.05.2015
Сообщений: 105
.NET Core

Написать программу, рассчитывающую площадь фигуры методом Монте-Карло

17.06.2020, 18:33. Показов 5823. Ответов 3

Студворк — интернет-сервис помощи студентам
Фигура ограничена линиями y=x*x + 1, y=(x-1)*(x-1), y=2/x; x = [0,2].
Найти площадь фигуры и относительную погрешность (eps = |S-S*|/S*; S* - реальная площадь фигуры)
Реальная площадь (по онлайн калькулятору через интеграл) = ~2
Написал программу, но получаю площадь сильно отличную от "реальной". При увеличении n точность меняется, но значение всё ещё далеко.
По примеру должно выполняться условие S-eps = S*, у меня оно не выполняется

Код
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
static void Main(string[] args)
{
            var rand = new Random();
            int n = 1000;
            double
                S = 0, S2 = 2,           // Площадь и истинная площадь 
                a = 0, b = 2, c = 2,    // границы и пик 
                x = 0, y = 0,              // точка
                eps,                           // относительная погрешность
                y1 = x * x + 1,           // первый график
                y2 = (x - 1) * (x - 1),  // второй график
                y3 = 2 / x;                 // третий график
            for (int i = 1; i < 5; i++)
            {
                x = 0; y = 0;
                for (int j = 1; j < n; j++)
                {
                    x = b * rand.NextDouble();
                    y = c * rand.NextDouble();
 
                    //Console.WriteLine("x: "+x + "; y: "+y);
 
                    if (x < 1 && y > 1) // Верхняя левая четверть.
                    {
                        if (y < y1) S++; 
                    }
                    else if (x < 1 && y < 1) // Нижняя левая четверть
                    {
                        if (y < y2) S++; 
                    }
                    else if (x > 1 && y > 1) // Верхняя правая часть 
                    {
                        if (y < y3) S++; 
                    }
                    else if (x > 1 && y < 1)     // Нижняя правая часть
                    {
                        if (y < y2) S++; 
                    }
                }
            }
            S = S/n;
            eps = Math.Abs(S - S2) / S2;
            Console.WriteLine("Число наблюдений: " + n);
            Console.WriteLine("Площадь методом Монте-Карло: " + S);
            Console.WriteLine("Точная площадь: " + S2);
            Console.WriteLine("eps: " + eps);
            Console.WriteLine("S-eps:" + (S - eps));
            return 0;
}


Добавлено через 1 час 15 минут
Ой, чет совсем поплыл с матаном, из-за этого стал воспринимать переменные как функции
Исправлено:
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
static void Main(string[] args)
{
            var rand = new Random();
            int n = 1000;
            double
                S = 0, S2 = 2,            // Площадь и истинная площадь 
                a = 0, b = 2, c = 2,    // границы и пик 
                x = 0, y = 0,              // точка
                eps;                         // относительная погрешность
            for (int i = 1; i < 5; i++)
            {
                x = 0; y = 0;
                for (int j = 1; j < n; j++)
                {
                    x = b * rand.NextDouble();
                    y = c * rand.NextDouble();
 
                    //Console.WriteLine("x: "+x + "; y: "+y);
 
                    if (x < 1 && y > 1) // Верхняя левая четверть.
                    {
                        if (y <  x * x + 1) S++; 
                    }
                    else if (x < 1 && y < 1) // Нижняя левая четверть
                    {
                        if (y < (x - 1) * (x - 1)) S++; 
                    }
                    else if (x > 1 && y > 1) // Верхняя правая часть 
                    {
                        if (y < 2 / x;) S++; 
                    }
                    else if (x > 1 && y < 1)     // Нижняя правая часть
                    {
                        if (y < (x - 1) * (x - 1)) S++; 
                    }
                }
            }
            S = S/n;
            eps = Math.Abs(S - S2) / S2;
            Console.WriteLine("Число наблюдений: " + n);
            Console.WriteLine("Площадь методом Монте-Карло: " + S);
            Console.WriteLine("Точная площадь: " + S2);
            Console.WriteLine("eps: " + eps);
            Console.WriteLine("S-eps:" + (S - eps));
            return 0;
}
Добавлено через 5 минут
В этом и была проблема *facepalm*

Добавлено через 8 минут
А, нет, не решилась проблема
0
Лучшие ответы (1)
IT_Exp
Эксперт
34794 / 4073 / 2104
Регистрация: 17.06.2006
Сообщений: 32,602
Блог
17.06.2020, 18:33
Ответы с готовыми решениями:

Вычисление площади фигуры на фотоснимке методом Монте-Карло
Дан фотоснимок электронного микроскопа с изображением некоторой фигуры. Фотоснимок предварительно обрабатывается и переводится в монохром...

Площадь фигуры методом Монте-Карло в Matlab
Мне нужен код для матлаба, для вычисления площади параллелограмма методом Монте-Карло с указанием двух сторон и углом между ними. Если...

Определить методом Монте Карло площадь фигуры
Здравствуйте пожалуйста помогите по решение этой задачке. Определить Определить методом Монте Карло площадь фигуры.

3
172 / 91 / 36
Регистрация: 22.05.2020
Сообщений: 308
17.06.2020, 18:50
magdake, мне кажется, у вас ошибка при сравнении координат в строках 26 и 34. Написано y < ..., а должно быть y > ....
0
1 / 1 / 0
Регистрация: 02.05.2015
Сообщений: 105
17.06.2020, 19:59  [ТС]
valker, Да, благодарю, заработало
Но не могу понять, почему, когда убираю внешний цикл for (int i = 1; i < 5; i++), программа начинает считать сильно заниженное значение, ведь, по-сути, цикл увеличивает количество точек в 4 раза (остался от примера с матлаба, там это используется для вычисления сразу для разного кол-ва шагов, как я понял)
0
172 / 91 / 36
Регистрация: 22.05.2020
Сообщений: 308
17.06.2020, 23:08
Лучший ответ Сообщение было отмечено magdake как решение

Решение

magdake, дело в том, что S / n это отношение площади фигуры к площади окружающего прямоугольника. Окружающий прямоугольник у вас 2*2 = 4. То есть можно было убрать внешний цикл и просто написать S = 2 * 2 * S / n;
Кстати, j в цикле надо инициализировать нулём, а не единицей.
1
Надоела реклама? Зарегистрируйтесь и она исчезнет полностью.
BasicMan
Эксперт
29316 / 5623 / 2384
Регистрация: 17.02.2009
Сообщений: 30,364
Блог
17.06.2020, 23:08
Помогаю со студенческими работами здесь

Определить площадь фигуры (методом Монте-Карло)
Здраво народ! Стоит задача: с помощью метода Монте-Карло найти площадь фигуры. Для этого я генерирую случайные числа в диапазоне от 0 до 8....

Определить площадь фигуры методом Монте-Карло
Здравствуйте, эксперты. Есть задача определения площади фигуры(по имеющимся графическим данным)- график не создаётся! Нужно...

Методом Монте-Карло определить площадь фигуры в Python
Вычислить площадь заштрихованной фигуры методом Монте-Карло, не используя случайные числа.(на языке Python)

Методом Монте-Карло определить площадь заштрихованой фигуры
Нарисовать линиями в Chart for (int i=0; i&lt;=100; i++) { x=i; ...

Определить площадь представленной фигуры методом Монте-Карло
Пожалуйста, помогите составить программу определения площади заданной фигуры методом Монте-Карло.


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

Или воспользуйтесь поиском по форуму:
4
Ответ Создать тему
Новые блоги и статьи
Отправка уведомления на почту при изменении наименования справочника
Maks 24.03.2026
Программная отправка письма электронной почты на примере изменения наименования справочника "Склады" в конфигурации БП3. Перед реализацией необходимо выполнить настройку системной учетной записи. . .
модель ЗдравоСохранения 5. Меньше увольнений- больше дохода!
anaschu 24.03.2026
Теперь система здравосохранения уменьшает количество увольнений. 9TO2GP2bpX4 a42b81fb172ffc12ca589c7898261ccb/ https:/ / rutube. ru/ video/ a42b81fb172ffc12ca589c7898261ccb/ Слева синяя линия -. . .
Midnight Chicago Blues
kumehtar 24.03.2026
Такой Midnight Chicago Blues, знаешь?. . Когда вечерние улицы становятся ночными, а ты не можешь уснуть. Ты идёшь в любимый старый бар, и бармен наливает тебе виски. Ты смотришь на пролетающие. . .
Контроль уникальности заводского номера - вариант №2
Maks 24.03.2026
В отличие от предыдущего варианта добавлено прерывание циклов, также добавлены новые переменные для сохранения контекста ошибки перед прерыванием цикла: Процедура ПередЗаписью(Отказ, РежимЗаписи,. . .
SDL3 для Desktop (MinGW): Вывод текста со шрифтом TTF с помощью библиотеки SDL3_ttf на Си и C++
8Observer8 24.03.2026
Содержание блога Финальные проекты на Си и на C++: finish-text-sdl3-c. zip finish-text-sdl3-cpp. zip
Жизнь в неопределённости
kumehtar 23.03.2026
Жизнь — это постоянное существование в неопределённости. Например, даже если у тебя есть список дел, невозможно дойти до точки, где всё окончательно завершено и больше ничего не осталось. В принципе,. . .
Модель здравоСохранения: работники работают быстрее после её введения.
anaschu 23.03.2026
geJalZw1fLo Корпорация до введения программа здравоохранения имела много невыполненных работниками заданий, после введения программы количество заданий выросло. Но на выплатах по больничным это. . .
Контроль уникальности заводского номера - вариант №1
Maks 23.03.2026
Алгоритм контроля уникальности заводского (или серийного) номера на примере документа выдачи шин для спецтехники с табличной частью в конфигурации КА2. Данные берутся из регистра сведений, по. . .
КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin
Copyright ©2000 - 2026, CyberForum.ru