待ち行列シュミレーションのプログラム

締切済

質問者：ahoodesu
質問日時：2008/06/18 13:29
回答数：2件

今現在待ち行列シュミレーションを勉強しているのですが、M/M/2モデルのプログラムをそれぞれ窓口数１(M/M/1)、窓口数３(M/M/3)のプログラムにしたいのですが上手くできません。
窓口数２のプログラムを記載します(窓口数の増減に関係ないと思われる後半のプログラムは記載してません)ので、どなたかわかる方がいましたらご教授のほうお願いします。

M/M/2　待ち行列のモデル
-- ポアソン到着 => 到着時間間隔は指数分布
-- 指数サービス => サービス時間は指数分布
*/

#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include <time.h>

#define N_SERVICE2/*　窓口数 */
#define N_USER 10000/*　総到着人数 */
#define ARRIVAL_INTERVAL 30.0/*　平均到着間隔*/
#define SERVICE_TIME 70.0/*　平均サービス時間 */

double uni_random(void); /* 区間[0,1]の一様乱数 */
double exp_random(double); /* 平均値を引数にとる指数乱数*/

int main(void){
int i; /* 客番号 */
int j; /*窓口番号*/
double arrival[N_USER]={0.0}　/*N人の到着時刻*/
int counter[N_USER]={0}; /*N人の利用窓口*/
int current_user[N_SERVICE]={0};/*最近の窓口利用者 */
double srv_start[N_USER]={0.0}/*N人のサービス開始時刻 */
double wait[N_USER]={0.0};/*N人の待ち時間 */
double srv_time[N_USER]={0.0}/* N人のサービス時間 */
double srv_finish[N_USER]={0.0}; /* N人のサービス終了時刻*/
double avr_wait; /* 平均待ち時間*/
double avr_srvtime; /* 平均サービス時間　 */
double avr_srvrate[N_SERVICE]={0.0}; /* 窓口の稼働率 */

/* 乱数の「種」を初期化*/
srand((unsigned)time(NULL));

/* 最初の客が到着した時刻から開始. 窓口０に入ってすぐサービス開始 */
arrival[0]=0.0;
counter[0]=0;
current_user[0]=0;
srv_start[0]=0.0;
wait[0]=0.0;
srv_time[0]=exp_random(SERVICE_TIME);
srv_finish[0]=srv_start[0] + srv_time[0];

/* 2人目の客が到着. 窓口０が空いていれば窓口０に、塞がっていれば窓口１に入ってすぐサービス開始 */
arrival[1]=arrival[0]+exp_random(ARRIVAL_INTERVAL);
if(srv_finish[0]>arrival[1]){
counter[1]=1;
current_user[1]=1;
}
else{
counter[1]=0;
current_user[0]=1;
}
srv_start[1]=arrival[1];
wait[1]=0.0;
srv_time[1]=exp_random(SERVICE_TIME);
srv_finish[1]=srv_start[1]+srv_time[1];

for(i=2; i<N_USER; i++){
/* 客の到着時刻　＝　前の客の到着時刻　＋　到着間隔　*/
arrival[i]=arrival[i-1]+exp_random(ARRIVAL_INTERVAL);

/* 利用窓口　＝　2つに窓口で早く空くほう　*/
if(srv_finish[current_user[1]]<srv_finish[current_user[0]]){
counter[i]=1;
}
else{
counter[i]=0;
}