Research Article

# Estimation Some Parameters of K-station Series Model

M.M. El Genidy

ABSTRACT

This study deals with the problem of statistical point estimation of K-station series model with unlimited inter queue capacity. K-station series model all over the world, both public and private, are facing continuously increasing traffic loads. The maximum likelihood function is used to estimate the arrival rates, service rates and the utilization factor in steady state when the service discipline is general.

 Services Related Articles in ASCI Search in Google Scholar View Citation Report Citation

 How to cite this article: M.M. El Genidy , 2012. Estimation Some Parameters of K-station Series Model. Asian Journal of Applied Sciences, 5: 307-313. DOI: 10.3923/ajaps.2012.307.313 URL: https://scialert.net/abstract/?doi=ajaps.2012.307.313

Received: December 07, 2011; Accepted: January 25, 2012; Published: March 26, 2012

INTRODUCTION

This study considers a system with K-station in series, supposes that arrivals at station 1 are generated from an infinite population according to a Poisson distribution with mean serviced units will move successively from one station to the next until they are discharged from station K. Service time at each station I has exponential distribution. The problem of point estimation is formulated as follows: A random sample is drawn from a population whose distribution has a known mathematical form but involves a certain number of unknown parameters. On the basis of this random sample, it is required to estimate the value of unknown parameters of the population.

Research has done on parameter estimation of the extended generalized Gaussian family distributions using maximum likelihood scheme (Lee, 2010). Time domain DS-UWB channel maximum likelihood algorithm was in working by Yin et al. (2011) but the work of Shi et al. (2011) was in solving project scheduling problems using estimation of distribution algorithm with local simplex search. Desta et al. (2011) carried out review on selected channel estimation algorithms for orthogonal frequency division multiplexing system. Wang et al. (2011) deal with a novel joint estimation algorithm for multi-parameter of underwater acoustic channels. Abd El-Salam (2011) deals with an efficient estimation procedure for determining ridge regression parameter. Ojo (2010) has work done in on the estimation and performance of one-dimensional autoregressive integrated moving average bilinear time series models. Oyeyemi et al. (2010) consider on the estimation of power and sample size in test of independence. Ojo (2008) carried out on the performance and estimation of spectral and bispectral analysis of time series data.

DESCRIPTION OF THE PROBLEM AND ITS SOLVING

Under these conditions, the out put of station I follows Poisson process with arrival rate λ, where station I follows the model (M/Mi/1):(GD/∞/∞).

This means for the ith station, the steady state probabilities are given by:

For I = 1, 2,…, k, where ni is the number in the system consisting of station I only. Steady state results will exist only if utilization factor ρi = λ/μi<1.

To find the likelihood function, it is necessary to find three basic components, Namely:

 • The initial number of units ni with probability Pr (ni) = for I = 0, 1, 2,…, k • The inter-arrival time of lengths t1 for R arrival units with probability • The service time of durations qj for M units with probability:

It follows that the likelihood function is given by:

Hence:

Hence:

Where:

After differentiating we get:

 (1)

And:

 (2)

where, and are the maximum likelihood estimators of the parameters λ and μi, respectively. Then:

Where:

Using the Eq. 1 and 2, then:

 (3)

And:

 (4)

where, is the maximum likelihood estimator of E (ni).

From the Eq. 3 and 4 we get:

 (5)

And:

 (6)

Substituting for its value in the Eq. 6, we get the equation:

Where:

Hence:

 (7)

Suppose that M = R = ni and using the Eq. 5 and 7, it follows that:

And:

 Table 1: Application with numerical results* *Suppose that the following data access times of units were given

In Table 1, we can obtain the following results:

 • Arrival rate of all customers to server 1 = 4.5 • Arrival rate of all customers to server 2 = 9.5 • Arrival rate of all customers to server 3 = 19.5 • Service rate for server 1 = 5 • Service rate for server 2 = 10 • Service rate for server 3 = 20.6666666666667 • Number of servers = 3 • Sum of arrival times of customers to server 1 = 27 • Sum of arrival times of customers to server 2 = 57 • Sum of arrival times of customers to server 3 = 117 • Sum of arrival times of customers to all servers = 201 • Sum of service times of server 1 = 30 • Sum of service times of server 2 = 60 • Sum of service times of server 3 = 124 • Sum of service times of all servers = 214 • Number of arrival customers to the system = 6 • Number of arrival customers to any server = 6 • Number of departure customers from the system = 6 • Utilization factor (traffic intensity) = 0.868570344017331

CONCLUSION

The mentioned above results show that the system is stable because the coefficient of efficiency has ρ (traffic intensity) is less than or equal to 1 and if it is greater than 1, then that leads to instability system and losses, so the system must be stopped or modify the service times in this server or change the arrival times of the units. Therefore, we can control the system to obtain less expensive, rapid processes and stability of the system.

In Fig. 1-3, the following formula was obtained:

 y = ax4+bx3+cx2+dx+e

where, x represents the arrival times, y represents the times of service received by the server, thus we can forecast the length of service time required from the server if the arrival time of the unit was given.

Figure 4 shows the distribution of the service times. Therefore, we can obtain the average service time of any unit or client within the system from the equation:

 y = aebx

 Fig. 1: The relation between arrival times of all customers and service times of server 1

 Fig. 2: The relation between arrival times of all customers and service times of server 2

 Fig. 3: The relation between arrival times of all customers and service times of server 3

 Fig. 4: The distribution of service times of servers

In Fig. 4, the values of a and b depend on the unit within the system.

The fourth unit has the values a = 3.188 and b = 0.602, if the average of the arrival time was given, then we can obtain the average service time in the system. As a result, we can predict the time of service required from the system and know the number of servers that allows us to manage the financial cost to add or delete a server to the system, like as the cost of traffic loads and waiting the cargo ships in the port to get the service from several servers.

REFERENCES

1:  Lee, J.Y., 2010. Parameter estimation of the extended generalized gaussian family distributions using maximum likelihood scheme. Inform. Technol. J., 9: 61-66.

2:  Yin, Z., Z. Shi, G. Qi and Y. Zhou, 2011. Time domain DS-UWB channel estimation using maximum likelihood algorithm. Inform. Technol. J., 10: 127-133.

3:  Shi, Y.J., G.J. Shen and W. Chen, 2011. Solving project scheduling problems using estimation of distribution algorithm with local simplex search. Inform. Technol. J., 10: 1374-1380.
CrossRef  |

4:  Desta, Y.T., J. Tao and W. Zhang, 2011. Review on selected channel estimation algorithms for orthogonal frequency division multiplexing system. Inform. Technol. J., 10: 914-926.

5:  Wang, B., G. Yang and G. Liu, 2011. A novel joint estimation algorithm for multi-parameter of underwater acoustic channels. Inform. Technol. J., 10: 440-445.

6:  Abd El-Salam, M.E.F., 2011. An efficient estimation procedure for determining ridge regression parameter. Asian J. Math. Stat., 4: 90-97.