数学建模社区-数学中国

标题: 【求问大牛】这道题用lingo怎么编出来啊 [打印本页]

作者: 夏默默 时间: 2013-8-12 11:37
标题: 【求问大牛】这道题用lingo怎么编出来啊

作者: wujianjack2 时间: 2013-8-12 14:40
  楼主你好！问题不错！
  在我看来，这是一个排队论问题，可惜我对此类问题不太熟悉，但还是尝试给你解决了一下，得到了结果，程序代码如下：
MODEL:
SETS:
LINE:NN,P,N_FAC,RP_N,S_RP_N,C;
ENDSETS
DATA:
N=10;
LAMBDA=3;MU=2;
LINE=1..N;
ENDDATA
R=LAMBDA/MU;
@FOR(LINE(I):NN(I)=I);
@FOR(LINE(I):N_FAC(I)=@SMAX(@PROD(LINE(J)|J#LE#I:NN(J))));
@FOR(LINE(I):RP_N(I)=@POW(R,NN(I))/N_FAC(I));
@FOR(LINE(I):S_RP_N(I)=@SUM(LINE(J)|J#LE#NN(I):RP_N(J))+1);
@FOR(LINE(I):P(I)=RP_N(I)/S_RP_N(I));
@FOR(LINE(I):C(I)=R*(1-P(I)));
END

在LINGO 11上运行得到如下结果（仅列出比较重要的）：
Feasible solution found.
Total solver iterations:                            0
Variable          Value

P( 1)    0.6000000
P( 2)    0.3103448
P( 3)    0.1343284
P( 4)    0.4795737E-01
P( 5)    0.1418316E-01
P( 6)    0.3533261E-02
P( 7)    0.7565545E-03
P( 8)    0.1418338E-03
P( 9)    0.2363842E-04
P( 10)    0.3545750E-05

C( 1)    0.6000000
C( 2)       1.034483
C( 3)       1.298507
C( 4)       1.428064
C( 5)       1.478725
C( 6)       1.494700
C( 7)       1.498865
C( 8)       1.499787
C( 9)       1.499965
C( 10)    1.499995

由结果可知，当N=4时，P=0.04795737,对应的C=1.428064，与原题的结果是一致的。

但是，本程序实际上还有改进之处，但是我目前只能做到这样子，抱歉！以上也仅是个人的一点见解，望笑纳！

欢迎光临数学建模社区-数学中国 (http://www.madio.net/)