マンガでよめる痔のこと・薬のこと

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

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));

/* 最初の客が到着した時刻から開始. 窓口0に入ってすぐサービス開始 */
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人目の客が到着. 窓口0が空いていれば窓口0に、塞がっていれば窓口1に入ってすぐサービス開始 */
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;
}

A 回答 (2件)

> シュミレーション



シミュレーション(simulation)が正しいです。
正しい用語を使いましょう。
    • good
    • 0

「上手くできない」というのは, 何がどう「上手くできない」んでしょうか?

    • good
    • 0

お探しのQ&Aが見つからない時は、教えて!gooで質問しましょう!

このQ&Aを見た人が検索しているワード

このQ&Aと関連する良く見られている質問

Q待ち行列のシミュレーションのプログラムが書けないです。

現在待ち行列のシミュレーションを勉強しています。
M/M/1のプログラムは書けて、シミュレーションもしっかりとできました。
そこで、次はサーバーは1つのままで、複数の到着流をもつ待ち行列をシミュレーションしてみようと思ったのですが、
思うようにプログラムが書けなくて困っています。
どなたか要点だけでも良いので教えていただき、プログラムの手助けをしていただけるとうれしいです。

Aベストアンサー

アイデアだけ。(自信なし)
やりかたはイロイロあると思うけど・・・。
---

「到着流」というのを「客発生器」と呼ぶとして…

結果としては、たとえ2個の「客発生器」を作ったとしても、
「サービス処理側」からしてみれば、
「客発生器」が複数あるなんてことはつゆ知らず、
単に、
「ああ、『到着間隔が1/2になった客発生器』を使ったのね・・・、ああ忙しい」
としか思わないハズ。
つまり、ソースコード的には、
客発生器に変更が生ずる以外、
他に変更はないはず、ということ。

つまり、いわば
『到着間隔が1/2(くらい)になった客発生器』を1つ
作ってしまえばいいのだ。
---

たとえば、いままで、「次の客到着時間」を、
唯一の「客発生器」を使って
time = nextTime();
などと取得いていたとする。
この「元の」発生器を使って、
新しい発生器を作ってみる。

#define kINVAL = 0;//無効
long time1, time2;//待避用
/**---------------
*新たな「客生成器」
-----------------*/
long newNextTime() {
 long res=0;
 if (time1 == kINVAL) {
  time1 = nextTime();//1個目の客生成器
 }
 
 if (time2 == kINVAL) {
  time2 = nextTime2();//2個目の客生成器
  //(※一個目の発生器と同じであってはならない)
 }
 
 //この時点で、「返すべき時間の候補」が2つ存在する
 
 //2つの候補のうち、小さいほう(時刻の早いほう)を返す
 if(time1<time2) {
  res = time1;
  time1 = kINVAL;//使用したので無効
  return time1;
 } else {
  res = time2;
  time2 = kINVAL;//使用したので無効
  return time2;
 }
}

---
で、今までnextTime()を使っていた場所で
代わりにこのnewNextTime()を使う。

結果的に到着間隔が1/2になると思うが…。
実験してみてね。
(ちょっと頭の中だけで考えてソースコード書いたので、
間違ってたらごめん)
---
↓ネット上にあるソースコード。
(※僕が書いたアイデアとは直接の関係はない)
http://www.sist.ac.jp/~suganuma/cpp/3-bu/18-sho/18-sho.htm#e-18-6

参考URL:http://www.sist.ac.jp/~suganuma/cpp/3-bu/18-sho/18-sho.htm#e-18-6

アイデアだけ。(自信なし)
やりかたはイロイロあると思うけど・・・。
---

「到着流」というのを「客発生器」と呼ぶとして…

結果としては、たとえ2個の「客発生器」を作ったとしても、
「サービス処理側」からしてみれば、
「客発生器」が複数あるなんてことはつゆ知らず、
単に、
「ああ、『到着間隔が1/2になった客発生器』を使ったのね・・・、ああ忙しい」
としか思わないハズ。
つまり、ソースコード的には、
客発生器に変更が生ずる以外、
他に変更はないはず、ということ。

つまり、い...続きを読む


人気Q&Aランキング