Добавил:

Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.

Вуз:

Национальный университет водного хозяйства и природопользования

Предмет:

[НЕСОРТИРОВАННОЕ]

Файл:

Форма.docx

Скачиваний:

Добавлен:

12.02.2016

Размер:

275.07 Кб

Скачать

☆

<<< < Предыдущая 1 2 3 4 5 6 78 / 88

Висновок

Таким чином поставлена задача про створення програми, що дозволяє змоделювати багатоканальну систему масового обслуговування з обмеженою довжиною черги, вирішена. Створена програма дозволяє зробити розрахунок таких параметрів системи як:

P_i - імовірності станів;
P_від - імовірність відмовлення;
a - абсолютну пропускну здатність системи;
q - відносну пропускну здатність системи;
t_оч- середній час чекання в черзі;
t_сист - середній час перебування в системі;
N_реаліз - число реалізацій для досягнення імовірності відмовлення відмінної від теоретичної на величину  ;

при заданих початкових умовах :

n - число каналів обслуговування;
m - число місць у черзі;
 - інтенсивність надходження заявок;
 - інтенсивність обслуговування;

 - точність для розрахунку числа реалізацій;

 - рівень довіри.

А також візуально простежити роботу системи і порівняти теоретичні та практичні результати розрахунків. Програма може підлягати подальшій доробці , наприклад можна задати ще декілька різних законів розподілу для інтенсивностей  і  чи зробити так, щоб система сама вибирала оптимальний час роботи, тобто час коли система ввійде в сталий режим.

Список використаної літератури

1. Советов, Б. Я. Моделирование систем [Текст]: навчальний посібник / Б. Я. Советов,С. А. Яковлев – М.: Высшая школа, 2001. – 343c. – ISBN 966- 574-022-9

2. Советов, Б. Я. Моделирование систем [Текст] / Б. Я. Советов – Практикум. – М.: Высшая школа, 2003. – 295c. – ISBN 966-346-160-4

3. Наконечний, С. І. Стохастичні процеси та моделі в економіці, соціології, екології [Текст] / С. І. Наконечний, В. І. Жлуктенко, С. С. Савіна: Навч. посібник. – К.: КНЕУ, 2002. – 226 с. – ISBN: 966-574-422-4

4. Ситник, В. Ф. Імітаційне моделювання [Текст] / В. Ф. Ситник, Н. С. Орленко: Навч. посібник. – К.: КНЕУ, 1998. –232 с. – ІSBN 966-574- 404-6

5. Афанасьев, М. Ю. Исследование операций в экономике: модели, задачи, решения [Текст] / М. Ю. Афанасьев, Б. П. Суворов: Учеб. пособие. – М.: ИНФРА-М, 2003. – 444 с. – ISBN 966-683-024-8

Додаток

Лістинг програми в С++:

//**********************СМО з обмеженою довжиною черги*********************

#include <conio.h>

#include <iostream.h>

#include <math.h>

#include <stdlib.h>

#include <dos.h>

#include <stdio.h>

#include <string.h>

#define size 10

#define DOWN 80

#define UP 72

#define ENTER 13

float lam=7;

float mu=2;

int n=5;

int m=7;

float Q=0.97;

float eps=0.03;

struct Request

{

int Num;

unsigned long Time;

unsigned long Time_format;

};

struct Channel

{

int Busy;

int NumReq;

unsigned long TimeBegin;

unsigned long TimeBegin_format;

float TimeToProcess;

};

Channel ChanArray[size];

Request ReqArray[size];

float ProcTime[size];

int ReqInArray;

long unsigned LastRequestTime;

float LastRivnVar;

int CurReqNum=1;

void Phon(void);

//-------------------------------Степiнь числа------------------------------

unsigned long Stepen(int a,int b)

{

unsigned long s=a;

for(int i=1;i<b;i++)

s*=a;

return s;

}

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

//------------------------------Факторiал числа-----------------------------

int Faktor(int c)

{

int i=1;

int f=1;

f*=i++;

while (i<=c);

return f;

}

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

//------------------------------Системний час-------------------------------

unsigned long GetSTime(int format)

{

struct time t;

gettime(&t);

if (!format)

return ((unsigned long)t.ti_hour*1000000+(unsigned long)t.ti_min*10000+(unsigned long)t.ti_sec*100+(unsigned long)t.ti_hund);

else

return ((unsigned long)t.ti_hour*100*60*60+(unsigned long)t.ti_min*100*60+(unsigned long)t.ti_sec*100+(unsigned long)t.ti_hund);

}

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

//-----Генератор псевдо-виподкових чисел з рiвномiрним розподiлом-----------

float RivnRoz()

{

int i;

unsigned long n;

float nn;

i=random(100);

const int M=1236;

const int L=3;

const int n0=225;

const int m=300;

n=Stepen(L,i)*n0+((M*(Stepen(L,i)-1)/(L-1)))%m;

nn=(n-(n/1000)*1000)/1000.0;

return (nn*2*lam);

}

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

//--------------------------Формування заявок-------------------------------

int GetRequest(Request &rq)

{

int f=0;

unsigned long SystemTime;

SystemTime=GetSTime(1);

if ((SystemTime-LastRequestTime)>=LastRivnVar)

{

delay(1);

LastRivnVar=RivnRoz();

rq.Num=CurReqNum;

rq.Time=SystemTime;

rq.Time_format=GetSTime(0);

LastRequestTime=SystemTime;

if (CurReqNum>=100) CurReqNum=1;

else CurReqNum++;

f=1;

}

// else cout<<"\n????";

// getch();

return f;

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

//--------------------------Звiльнення каналiв------------------------------

}

void CheckChannels(int &SumBusyChan)

{

unsigned long SystemTime;

unsigned long t;

int th,tx,ts,tm;

SystemTime=GetSTime(1);

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

if (ChanArray[i].Busy)

{

SumBusyChan++;

if (SystemTime-ChanArray[i].TimeBegin>=ChanArray[i].TimeToProcess)

{

ChanArray[i].Busy=0;

t=GetSTime(0);

tm=t%100;

ts=(t%10000)/100;

tx=(t%1000000)/10000;

th=(t%100000000)/1000000;

gotoxy(16,22);

cprintf("Заявка №%d покинула систему о %d:%d:%d.%d ",ChanArray[i].NumReq,th,tx,ts,tm);

}

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

//--------------------------Обслуговування заявок---------------------------

int ProcessRequest(Request &rq)

{

int f=0;

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

if (!ChanArray[i].Busy)

{

ChanArray[i].TimeBegin=GetSTime(1);

ChanArray[i].TimeBegin_format=GetSTime(0);

ChanArray[i].TimeToProcess=ProcTime[i];

ChanArray[i].Busy=1;

ChanArray[i].NumReq=rq.Num;

f=1;

break;

}

return f;

}

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

//--------------------------Формування черги--------------------------------

int InsertIntoOchered(Request &rq,long &vid)

{

int f=0;

if (ReqInArray<m)

{

ReqArray[ReqInArray++]=rq;

f=1;

}

else

{

gotoxy(16,22);

cout<<"Заявцi було вiдмовила в обслуговуваннi ";

vid++;

}

return f;

}

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

//------------------------Видiлення заявок з черги--------------------------

int GetFromOchered(Request &rq)

{

int f=0;

if (ReqInArray>0)

{

rq=ReqArray[0];

for(int i=0;i<ReqInArray-1;i++)

{

ReqArray[i].Time=ReqArray[i+1].Time;

ReqArray[i].Time_format=ReqArray[i+1].Time_format;

ReqArray[i].Num=ReqArray[i+1].Num;

}

ReqInArray--;

f=1;

}

return f;

}

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

//----------------------------Теоретичнi обрахунки--------------------------

void TeoretObrahynok(void)

{

float Tobsl,Toch,Tsystem,psi,P0,P1,x;

float Pvid,q,a,z;

float r=0;

Tobsl=1/mu;

psi=lam/mu;

P0=0;

x=psi/n;

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

P0+=pow(psi,i)/Faktor(i);

for(i=1;i<=m;i++)

P0+=pow(psi,n+i)/(Faktor(n)*pow(n,i));

P0=1/P0;

Pvid=(pow(psi,(n+m))/(pow(n,m)*Faktor(n)))*P0;

q=1-Pvid;

a=lam*q;

z=a/mu;

for(i=0;i<m;i++)

r+=(i+1)*pow(x,i);

r*=((pow(psi,(n+1))*P0)/(Faktor(n)))/n;

Toch=r/lam;

Tsystem=Toch+q*Tobsl;

//вивод результатiв

Phon();

textcolor(BLACK);

textbackground(BLACK);

gotoxy(14,22);

cprintf(" ");

for(i=4;i<22;i++)

{

gotoxy(67,i);

putch(' ');

}

textcolor(WHITE);

textbackground(BLUE);

gotoxy(13,3);

cprintf("г================Теоретичнi розрахунки===============¬");

gotoxy(13,4);

cprintf("¦ lam=%.2f mu=%.2f n=%d m=%d Q=%.2f eps=%.4f ¦",lam,mu,n,m,Q,eps);

gotoxy(13,5);

cprintf("¦=============================================T======¦");

gotoxy(13,6);

cprintf("¦ Iмовiрнiсть простою каналу обслуговування Po¦%.4f¦",P0);

gotoxy(13,7);