51124

Моделирование непрерывно-стохастической системы массового обслуживания

Лабораторная работа

Информатика, кибернетика и программирование

На вход n-канальной СМО с отказами поступает поток заявок с интенсивностью l = 6 заявок в час. Среднее время обслуживания одной заявки 0.8 часа. Каждая обслуженная заявка приносит доход 4у.е. Содержание одного канала обходится 2 у.е./час. Определить экономически целесообразное количество каналов.

Русский

2014-02-06

106.86 KB

3 чел.

Министерство образования Республики Беларусь

Учреждение образования

«Белорусский государственный университет

информатики и радиоэлектроники»

Кафедра программного обеспечения

информационных технологий

ОТЧЕТ

по лабораторной работе №4

“ Моделирование непрерывно-стохастической системы массового обслуживания ”

Выполнил:         Проверила:

студент гр. 750503       Костюк Елена

Макаревич Владимир

Минск 2009


Задание

7)  На  вход n-канальной  СМО  с  отказами  поступает  поток  заявок  с

интенсивностью l = 6  заявок  в  час.  Среднее  время  обслуживания  одной

заявки 0.8  часа.  Каждая  обслуженная  заявка  приносит  доход 4у.е.

Содержание одного канала обходится 2 у.е./час. Определить экономически

целесообразное количество каналов.

Теоретическое моделирование

#include <conio.h>

#include <stdio.h>

double fact(long n)

{

 if ( n==0 ) return 1.0;

 else return n * fact(n-1);

}

double extent(double x, int n)

{

 if ( n==0 ) return 1.0;

 return x * extent(x, n-1);

}

//--------------------------------

double S(int n, double w)

{

 if ( n==0 ) return 1.0;

 return S(n-1,w) + extent(w, n) / fact(n);

}

double P0(int n, double w)

{

 double s = S(n, w);

 return 1 / s;

}

int main()

{

 double w, h = 6.0, u = 1 / 0.8, c = 4.0, d = 2.0, Pk, p, profit, max_p = 0;

 int max_ch;

w = h/u;

printf("Chanels | Profit\n");

printf("========|========\n");

 for (int k=0; k<16; k++)

{

 p = P0(k, w);

 Pk = p * extent(w, k) / fact(k);

 profit = h * (1-Pk) * c - k * d;

 if ( max_p<profit ) max_p = profit, max_ch = k;

 printf(" %6d | %2.3f -> %f\n", k, profit, Pk);

}

printf(" ----------------\n");

printf(" %6d | %2.3f  <- Optimal\n", max_ch, max_p);

 return _getch();

}

Практическое моделирование

#include <conio.h>

#include <stdio.h>

#include <stdlib.h>

#include <math.h>

double RandExp(double h)

{

 double R = (double)rand()/RAND_MAX;

 return -log(1-R) / h;

}

class Channel

{

private: bool busy;

public: Channel()

 { busy = false; }

public: bool Busy()

 { return busy; }

public: bool Empty()

 { return !busy; }

public: void Occupy()

 { busy = true; }

public: void Liberate()

 { busy = false; }

};

class System

{

private: Channel **channel;

private: int n, state;

private: System() {}

public: System(int N)

 {

  state = 0;

  n = N;

  channel = new Channel*[n];

  for(int i=0; i<n; i++)

   channel[i] = new Channel();

 }

public: int State()

 { return state; }

public: bool InRequest()

 {

  for(int i=0; i<n; i++)

  {

   if ( channel[i]->Empty() )

   {

    channel[i]->Occupy();

    state++;

    return true;

   }

  }

  return false;

 }

public: bool OutRequest()

 {

  for(int i=0; i<n; i++)

  {

   if ( channel[i]->Busy() )

   {

    state--;

    channel[i]->Liberate();

    return true;

   }

  }

  return false;

 }

};

int main()

{

 double w, h = 6.0, u = 1 / 0.8, c = 4.0, d = 2.0, profit;

 double t1 = RandExp(h), t2 = RandExp(u), eps = 0.1;

 long k, N = 6, TotalTime = 1000000;

System *sys = new System(N);

k = N+1;

 double *P = new double[k];

 for(int i=0; i<k; i++) P[i] = 0;

 for(int i=0; i<TotalTime; i++)

{

 if( t1<=0 )

 {

  t1 = RandExp(h);

  sys->InRequest();

 }

 if( t2<=0 )

 {

  t2 = RandExp(u);

  sys->OutRequest();

 }

 t1-=eps, t2-=eps;

 P[sys->State()]++;

}

 for(int i=0; i<k; i++) P[i] = P[i] / TotalTime;

 for(int i=0; i<k; i++)

{

 printf("P[%d] = %2.8f\n", i, P[i]);

}

profit = h * (1-P[N]) * c - N * d;

printf("\n %6d | %2.8f\n", N, profit);

 return _getch();

}


 

А также другие работы, которые могут Вас заинтересовать

7717. Модели происхождения Вселенной 77.5 KB
  Модели происхождения Вселенной. 1. Ученик А.А. Фридмана Георг Гамов (эмигрировал в 1934г. из СССР в США) в 1948 г. разработал модель горячей Вселенной под названием Космология Большого взрыва» Радиус Вселенной в первоначальном состоянии был р...
7718. Схема анализа плана воспитательной работы классного руководителя 26 KB
  Схема анализа плана воспитательной работы классного руководителя наличие общей характеристики класса, детального анализа уровня его воспитанности, выделение положительных сторон и недостатков в поведении и деятельности учащихся ха...
7720. Теоретичні основи безпеки життєдіяльності 108.25 KB
  Теоретичні основи безпеки життєдіяльності Програмна анотація Предмет безпеки життєдіяльності. Теоретичні основи безпеки життєдіяльності. Науки, на яких базується безпека життєдіяльності. Теорія небезпеки та ризику. ОСНОВНІ ПОНЯТТЯ Безпек...
7722. Середовище життєдіяльності людини 122.98 KB
  Середовищежиттєдіяльності людини Програмна анотація Навколишнє середовище та середовище життєдіяльності людини. Класифікація негативних факторів. Негативні фактори активної групи та способи захисту від них. ОСНОВНІ ПОНЯТТЯ Навколишнє сер...
7723. Безпека життєдіяльності у законах та підзаконних актах 105.71 KB
  Безпека життєдіяльності у законах та підзаконних актах Програмна анотація Основні законодавчі акти з безпеки життєдіяльності. Екологічне законодавство. Основні положення Закону України Про охорону праці. ОСНОВНІ ПОНЯТТЯ Законодав...
7724. Організація і функціонування служби охорони праці на виробництві 112.5 KB
  Організація і функціонування служби охорони праці на виробництві Програмна анотація Організація служби охорони праці на підприємстві. Завдання інженерно-технічних працівників щодо забезпечення нешкідливих та безпечних умов праці. Захист працівників ...
7725. Створення комфортних умов праці на виробництві 163.99 KB
  Створення комфортних умов праці на виробництві Програмна анотація Нормалізація параметрів виробничого середовища. Методи і засоби контролю повітряного середовища. Освітлення. ОСНОВНІ ПОНЯТТЯ Параметри виробничого середовища Метеоро...