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

Моделирование - C++

Восстановить пароль Регистрация
 
Yska
Сообщений: n/a
15.01.2014, 01:07     Моделирование #1
Преподаватель дал написать задачку из статейки, чтоб сверить результаты полученные программой и полученые во время исследование системы обслуживания.

Постановка задачи: В настоящей работе рассматривается одноканальная система обслуживания M/D/1 с временем обслуживания tau > 0, временем пребывания на орбите T = 2*tau, входящим потоком Пуассона заявок с параметром lambda > 0. По правилам обслуживания в системе может быть не более двух заявок. Заявки принимаются в систему только в следующих случаях: 1) система свободна, и поступившая заявка начинает немедленно обслуживаться линией в течении времени tau > 0; 2) в системе находится одна заявка, обслуживаемая линией, при этом поступающая заявка отправляется на орбиту, и после пребывания на орбите длительностью T = 2*tau ей гарантированон немедленно начало обслуживания; 3) в системе находится одна заявка на орбите, и если в момент поступления очередной заявки оставшееся время время пребывания на орбите не меньше tau, то она немедленно принимается к обслуживанию, в противном случае она покидает систему (теряется).
Если в системе находятся две заявки (одна на орбите, а другая на линии обслуживания), то очередная заявка входящего потока также покидает систему (теряется.)

Я написал, но результаты не совпадают.

Вот хочу спросить совета, по поводу правильно ли я реализовал алгоритм? Потому что есть вероятность что формула полученная в статье не верна, а алгоритм как раз таки работает на отлично.

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
#include "stdafx.h"
#include <iostream>
#include <math.h>
#include <time.h>
using namespace std;
 
int main()
{
    int number_of_experiments, countY, countN, summ = 0, summN = 0; // количество опытов, количество обслуженных заявок, количество потеряных заявок
    double  time_to_work, arrival_intensity, tau;// время на работу, интенсивность поступления, тау(время обслуживания)
    double release_time_line, time_of_admission, orbit_time, orbit_for_while; // время освобождения линии, время поступления, орбитное время, орбитное для цикла while
    cout << "Plese enter the number of experiments: " << endl;
    cin >> number_of_experiments;
    cout << "Please enter time to work: " << endl;
    cin >> time_to_work;
    cout << "Please enter arrival intensity :" << endl;
    cin >> arrival_intensity;
    cout << "Please enter tau(service time): " << endl;
    cin >> tau;
    srand(time(NULL));
    for (int j = 1; j <= number_of_experiments; j++){
        release_time_line = 0; // время освобождения линии
        time_of_admission = 0; // время поступления
        orbit_time = 0; // орбитное время
        orbit_for_while = 0;
        countY = 0; countN = 0; // количество обслуженных заявок, количество потеряных заявок
        while(time_of_admission < time_to_work) {
            if(orbit_time == 0) {
                if(release_time_line <= time_of_admission) {
                    release_time_line = time_of_admission + tau;
                    countY = countY + 1;
                    time_of_admission = time_of_admission - (1/arrival_intensity)*log(0.0001 * ((rand() % 10000)+1));
                } else {
                    orbit_time = time_of_admission + 2*tau;
                    time_of_admission = time_of_admission - (1/arrival_intensity)*log(0.0001 * ((rand() % 10000)+1));
                }
            } else {
                if(release_time_line <= time_of_admission) {
                    if((orbit_time - time_of_admission) >= tau) {
                        release_time_line = time_of_admission + tau;
                        countY = countY + 1;
                        time_of_admission = time_of_admission - (1/arrival_intensity)*log(0.0001 * ((rand() % 10000)+1));
                    } else {
                        release_time_line = orbit_time + tau;
                        countY = countY + 1;
                        orbit_for_while = orbit_time;
                        orbit_time = 0;
                        while(orbit_for_while > time_of_admission) {
                            countN = countN + 1;
                            time_of_admission = time_of_admission - (1/arrival_intensity)*log(0.0001 * ((rand() % 10000)+1));
                        }
                    }
                } else {
                    countN = countN + 1;
                    time_of_admission = time_of_admission - (1/arrival_intensity)*log(0.0001 * ((rand() % 10000)+1));
                }
            }
        }
        summN = summN + countN;
        summ = summ + countY + countN;
    }
    cout << "result program: " << double(summN)/summ << endl;
    cout << "result formula: " << (2*tau*arrival_intensity*exp(arrival_intensity*tau)+exp(arrival_intensity*tau)+0.5*exp(arrival_intensity*tau)-0.5*exp(-(arrival_intensity*tau))-1-2*arrival_intensity*tau)/(3*exp(arrival_intensity*tau)+2*arrival_intensity*tau*exp(arrival_intensity*tau)-2-2*arrival_intensity*tau) << endl;
    system("pause");
    return 0;
}
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
15.01.2014, 01:07     Моделирование
Посмотрите здесь:

Моделирование C++
C++ Нейрокомпьютерное моделирование
C++ Моделирование случайных величин
C++ Моделирование случайных величин
C++ Моделирование систем
C++ мат.моделирование
C++ Моделирование СМО
C++ Моделирование COM

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

Или воспользуйтесь поиском по форуму:
После регистрации реклама в сообщениях будет скрыта и будут доступны все возможности форума.
Ответ Создать тему
Опции темы

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