Research Article  Open Access
Adaptive Maximums of Random Variables for Network Simulations
Abstract
In order to enhance the precision of network simulations, the paper proposes an approach to adaptively decide the maximum of random variables that create the discrete probabilities to generate nodal traffic on simulated networks. In this paper, a statistical model is first suggested to manifest the bound of statistical errors. Then, according to the minimum probability that generates nodal traffic, a formula is proposed to decide the maximum. In the formula, a precision parameter is used to present the degree of simulative accuracy. Meanwhile, the maximum adaptively varies with the traffic distribution among nodes because the decision depends on the minimum probability generating nodal traffic. In order to verify the effect of the adaptive maximum on simulative precision, an optical network is introduced. After simulating the optical network, the theoretic average waiting time of nodes on the optical network is exploited to validate the exactness of the simulation. The proposed formula deciding the adaptive maximum can be generally exploited in the simulations of various networks. Based on the precision parameter K, a recursive procedure will be developed to automatically produce the adaptive maximum for network simulations in the future.
1. Intoduction
Simulations are an important technique for the design of systems, the estimation of performance, and the maintenance of systems [1, 2]. It is widely used in various fields. For the simulation of complex systems, how to save computing time is an important topic. Moreover, it is also worthwhile to discuss how to reach acceptable simulative precision. In general, to promote simulative precision will lower simulative efficiency. Therefore, in order to enhance the simulative precision of complex systems, it is very important to take appropriate tradeoff between precision and efficiency.
For example, DQDB networks are systems with asynchronous transfer mode (ATM) and timedivision multiple access (TDMA) [3]. Its medium access control (MAC) protocol is so complex that the performance analysis of the network is very difficult [4, 5]. To make an exact analysis on the performance of the network is almost impossible [5]. Most papers estimate the performance of the network by simulations [6โ20]. Thus, it only depends on the precision of simulations to exactly comprehend the behavior of the complex system. However, if simulative precision is overpromoted, simulative efficiency will be suppressed. Therefore, how to reach the acceptable simulative precision and meanwhile take good tradeoff between precision and efficiency is an important problem for simulating complex systems. It is worthwhile to explore, but there are few papers discussing the problem.
In simulations, random variables are used to create various distribution functions of probabilities. These probabilities are applied to direct the input amplitude of signals or noise as simulating communication systems [1, 2]. For network simulations, probabilities are exploited to control the generation of nodal traffic. The probabilities assigned to represent nodal traffic are continuous. However, the probabilities created from random variables are discrete. Due to the inherent difference between the continuous and discrete probabilities, simulative errors will take place absolutely.
In order to take precise simulations, this paper suggests a statistical model. Based on the model, it is obvious that simulative precision only depends on the maximum of random variables controlling the generation of network traffic. Based on the perception, a simple formula is proposed to decide the feasible lower bound of maximums of random variables. In the formula, the feasible lower bound is dependent on a precision parameter, denoted by , and the minimum probability generating nodal traffic. The larger the precision parameter, the more the simulative precision. Then, a prime number that is slightly greater than the lower bound can be chosen as the maximum of random variables. The chosen maximum can adapt to the traffic distribution of simulated networks. So, the adaptive maximum cannot only result in acceptable simulative precision but also take desired efficiency. In practice, due to the approach deciding the adaptive maximum, simulative systems can make optimal tradeoff between reaching high precision and saving computing time.
So as to understand the effect of the adaptive maximum on simulative precision, an optical TDMA network is introduced. The MAC protocol of the optical TDMA network implements traffic control. The average waiting time of a node on the network is in inverse proportion to the traffic of the node [21]. Due to the quantitative analysis of the optical TDMA network, the rootmeansquare (rms) difference between the simulative and theoretic average waiting times of nodes is calculated to validate the exactness of simulations.
In Section 2, the suggested statistical model exhibits the relationship between the statistical error and the maximum of random variables. The discussion for deciding the adaptive maximum of random variables will also be shown in this section. The MAC protocol performing traffic control and some working conditions assumed for simulating optical TDMA networks are presented in Section 3. Section 4 illustrates the effect of the adaptive maximum on the performance estimated by simulations. The validation of simulations will also be shown in this section. Section 5 includes conclusions.
2. The Decision of Adaptive Maximums
Before simulations, a set of continuous probabilities is assigned to predefine the distribution of nodal traffic. In simulations, a set of discrete probabilities, which corresponds to the set of continuous probabilities, controls the generation of nodal traffic. The difference between continuous and discrete probabilities is used to manifest the influence of the maximum of integral random variables on simulative precision. Then the minimum probability in the set of continuous probabilities is exploited to decide the adaptive maximum of integral random variables.
In simulated networks, every node has one queue. Queues consist of cell (packet) buffers. Queues provide firstinfirstout (FIFO) service. The first cell buffer in a queue is attached to the transmission system of the simulated network. When an available slot on the transmission system is passed through, the contents of the first cell buffer will be written into the available slot.
The number of cells temporarily storing in queues is dependent on the traffic generated by nodes. The heavier the nodal traffic, the longer the queuing delay. According to the complex MAC protocol of most networks, an available slot on transmission systems appears for some node randomly. Hence, the prediction of the queuing delay of a specified cell is difficult. In order to estimate performance, most simulations assume that networks are with heavy load [4, 6, 11, 21]. This assumption will lengthen queues so that the theoretic analysis of the queuing delay of a specified cell will become more difficult. Therefore, how to enhance the accuracy of queuing delays is a key topic for network simulations. In order for taking precise simulations, a statistical model must be first introduced. Based on the model, the data resulted from simulations can be applied to calculate the nodal mean of queuing delays. The queuing delay is defined as the period for which a cell stays in queues.
For a node within networks, every cell generated by disassembling procedures is first stored in queues. Before the cell is transferred to transmission systems, it must be sequentially shifted into the first cell buffer of nodal queues. When the cell is within the first cell buffer, the MAC protocol will be exploited to decide the moment after that the node can write the cell out. In above operations, there are two moments relative to the queuing delay. The first moment is the instant that a cell enters queues. The second one is the flash that the cell is sent onto transmission systems. Let and in sequence denote the first and second moments of the th cell generated by the th node, where i and j are the ordinal number of cells and nodes, respectively. The queuing delay of the th cell generated by the th node is designated by . Then can be represented as where Max() is the ordinal number of the last cell generated by the th node, is the maximum of the integral random variable used to create the set of discrete probabilities, and is the number of nodes within networks.
Let denote the probability that the th node generates traffic. These probabilities, s, are chosen in accordance with interested simulative scenarios. The chosen probabilities are continuous, but their corresponding probabilities generated by an integral random variable are discrete. When a discrete probability corresponds to a continuous probability, both probabilities must be equal theoretically. But, in practice, the discrete probability is inherently different from its corresponding continuous probability. If the difference between them is large, statistical delays could not converge on the precise level that is acceptable.
Let represent an integral random variable that uniformly distributes over [0, ]. Its probability density function can be presented by
The random variable is applied to create s for simulations. Let denote the discrete probability corresponding to . The difference between and can be represented as
From (3), if theoretically approaches infinite, the difference will become zero. Therefore, (3) can be rearranged as
If a simulation is with a preciseprobability set , the statistical delays of network simulations will also be precise. In a word, when theoretically approaches infinite, the statistical queuing delay will converge precisely.
After discussing the enhancement of statistical precision, statistics of the queuing delay are presented. According to (1), the is a set of positive random numbers. Let () denote the mean of . Then () can be calculated by The mean () represents the average queuing delay of the th node.
If is theoretically infinite, the computing time taken in simulations is also infinite. In practice, can be chosen in accordance with s. If s is small, must be large enough to ensure statistical precision. In other words, must adaptively vary with s to promote simulative precision.
Let represent the minimum probability among s. must be greater than the inverse of for minimum precise requirement in simulations, that is,
In order to guarantee that the precise degree is acceptable, must be chosen so that where is a positive real number.
Then, , the adaptive maximum, can be represented by the equation
In (8), the value is a precision parameter. To enlarge will result in higher precision. On the other hand, if simulative procedures are with overlarge , they will consume more computing time but promote little precision. Thus, the decision of the value is dependent on the acceptable degree of precision.
Because the recursive formula of the powerresidue method is computationally very efficient [1, 22], it is widely adopted for generating random variables with uniform distribution. As the powerresidue method is exploited to generate the uniformdistribution random sequence, the maximum of the integral random variable must be a prime number. Therefore, the lower bound of , indicated by , must be first taken with the following equation:
Then, a prime number that is slightly greater than can be assigned as the adaptive maximum . The adaptive maximum can result in the optimal tradeoff between simulative precision and efficiency.
3. Optical TDMA Networks
In order for comprehending the effect of adaptive maximums of random variables on the precision of network simulations, an optical TDMA networks is introduced. Before depicting the structure of the network, the deduction for the average waiting time of nodes on the network (the waiting mean) is presented. The waiting time of a cell is the queuing delay that the cell waits in the first cell buffer of queues for an available slot on transmission systems. Based on the structure of the network, several working conditions are assumed for simulations. Due to these working conditions, an MAC protocol implementing traffic control is described.
For TDMA networks, a node must send requests to preserve empty slots when it is going to transmit messages. More requests preserve more slots. As the number of preserved slots of a node becomes large, the waiting mean of the node will be reduced. Therefore, if a node has more traffic, its waiting mean will be decreased. The relationship between the waiting mean and the traffic of the th node [21] can be presented as where () and () are the waiting mean and traffic of the th node, respectively, and is the slot rate of the optical TDMA networks.
Because optical TDMA networks are highspeed networks, the slot rate can approach infinite. Therefore, () can be rearranged as
In (11), () is a function of (). It exhibits that the waiting mean of a node on an optical TDMA network is in inverse proportion to the traffic of the node. The theoretic waiting mean will be exploited to validate simulations in Section 4. The structure of the optical TDMA network is shown in Figure 1.
In Figure 1, the medium between the slot generator and the slot terminator is an optical fiber. The slot flow on optical fibers is sent by slot generators and sinks into slot terminators. The number of nodes within the network is . Nodes are numbered from 0 to () in sequent order. The ordinal number of every node also relates to the nodal position in the topology. The period that the slot generator just completely sends a slot onto optical fibers is called a slot time. A slot length is the distance that a slot spreads on the optical fiber. Other working conditions concerning the space between adjacent nodes, the length of messages, and the traffic distribution among nodes are described as follows.
For all interested simulative scenarios, the space between adjacent nodes is one slot length. Messages are similar in length. Every message can be contained in the payload of a slot.
The traffic distribution among nodes affects the operation of traffic control in the MAC protocol. For the benefit of easily performing traffic control, a basic traffic denoted by is introduced. The amount of is dependent on the defined distribution of traffic. In a scenario, the traffic of some node can be several times the amount of . In order to obviously present the influence of the adaptive maximum on simulative precision, it is assumed that traffic is uniform distribution among nodes in every simulative scenario. Hence, the traffic of every node in a simulative scenario is equal to one . Because the optical fiber is a oneway bus, and all messages are not transmitted out of the network, the ()th node does not generate any traffic. Let denote the network traffic. Then can be presented by Therefore, the in simulations can be shown as Based on the introduction of , the approach of traffic control can be described as follows.
In this paper, slot frames are used to implement traffic control. The slot flow on optical fibers is partitioned into frames. There are 1/ slots in a frame. When a frame is passed to the th node, the node can only write one message into an empty slot within the frame when its queue is not empty. After the operation, the node must immediately stop writing messages out regardless of whether its queue is empty or not. After the moment, the node waits for the arrival of the next frame to restart the controlling process.
4. Simulations
The simulation is applied to present the influence of the adaptive maximum of a random variable on simulative precision. The waiting mean of optical TDMA nodes manifests the precise degree caused by different adaptive maximums. The simulative efficiency is dependent on the size of the adaptive maximum. The larger the maximum of the random variable, the lower the efficiency of the simulative system. On the other hand, the theoretic waiting mean calculated by (11) will be used to validate simulative data. The rms difference between the simulative and theoretic data, denoted by , is defined as where () and () are the simulative and theoretic waiting means of the th node, respectively.
Two parameters must be chosen before simulations. The parameters are the number of nodes and the network traffic . The chosen s and s are based on interested scenarios. After the choice of s and s, (13) can be exploited to calculate the basic traffic . Because messages are one slot in length, the basic traffic can be regarded as the minimum probability . In order for clearly distinguishing between precise degrees corresponding to different adaptive maximums, small is necessary. So, the chosen in all interested scenarios is equal to 0.25.
In simulations, two scenarios are interested. In order to certainly manifest the effect of adaptive maximums, the s in two scenarios are 40 and 50, respectively. On account of sufficiently exhibiting the influence of precision parameters on simulative precise degrees, three s are assigned in every scenario. The three s in every scenario are 300, 3000, and 30โ000, respectively. Then, (9) is used to calculate the corresponding of every parameter. In accordance with the powerresidue method, the adaptive maximum , which is a prime number and slightly greater than its corresponding , can be finally found. According to the description above, those relative parameters derived from the , s, and s are listed in Table 1.

For traffic control, the number of slots in a frame, which is equal to 1/, must vary with scenarios. Due to (13), the number of slots in a frame, denoted by , can be presented as Hence, the of the scenario with 40 nodes is 156 and that of the scenario with 50 nodes is 196. Because uniform traffic distribution is assumed, the traffic of every node is equal to the in each scenario. Consequently, the theoretic waiting mean calculated by (11) is the same as the in each scenario.
In the following figures that show simulative results, the horizontal axis is the ordinal number of nodes. Because the ordinal number of nodes is discrete, all curves in figures consist of piecewise lines. The waiting mean of nodes on the vertical axis is expressed in slot times.
After simulations, (5) is used to calculate the simulative waiting mean of nodes. Figures 2 and 3 show the variation of waiting means corresponding to two scenarios, respectively. In every figure, the horizontal solid line represents the theoretic waiting mean. Other three curves are relative to simulative data corresponding to three precision parameters. In two figures, it is obviously exhibited that curves will become smoother when the precision parameter is enlarged.
Based on the theoretic waiting mean, s calculated by (14) is used to validate the simulations. Table 2 presents these s. Every represents the difference between the horizontal solid curve and a simulative curve in each figure. Observing Table 2, larger precision parameter and smaller will result in smaller . Therefore, a simulative curve will approach the horizontal solid line if precision parameter is consecutively enlarged. However, the decrease of s does not linearly correspond to the increase of parameters. Consequently, simulative procedures with overlarge will consume more computing time but promote little precision. To adjust to a proper value can take an adaptive maximum to result in an acceptable precise level with relative high efficiency. In a word, network simulations with appropriate adaptive maximums can take optimal tradeoff between simulative precision and efficiency.

5. Conclusions
The paper discusses tradeoff between simulative precision and efficiency. Based on the statistical model of queuing delays, the difference between continuous and discrete probabilities is used to manifest the effect of the maximum of random variables on the statistical error. Then a simple method with precision parameter is proposed to decide the maximum of random variables. The maximum adapts to the minimum probability among s. The adaptive maximum must be enlarged as the minimum probability of s becomes smaller.
In order to manifest the effect of the adaptive maximum on the simulative precision, an optical TDMA network whose MAC protocol performs traffic control is simulated. The average waiting time of some optical TDMA node is in inverse proportion to the traffic of the node. The theoretic average waiting time is exploited to calculate s to validate the exactness of simulations. Simulative results exhibit that the adaptive maximum can take optimal tradeoff between simulative precision and efficiency.
In network simulations, the adaptive maximum not only results in the acceptable degree of precision but also suitably saves computing time. Based on the precision parameter , a recursive procedure will be developed to automatically generate the adaptive maximum in future.
References
 M. C. Jeruchim, P. B. Balaban, and K. S. Shanmugan, Simulation of Communication Systems: Modeling, Methodology, and Techniques, Kluwer Academic Publishers, New York, NY, USA, 2000.
 C. B. Rorabaugh, Simulating Wireless Communication Systems, Pearson Education, Upper Saddle River, NJ, USA, 2004.
 IEEE Std 802.6, โDistributed Queue Dual Bus (DQDB) Subnetwork of a Metropolitan Area Network (MAN),โ December 1990. View at: Google Scholar
 C. Y. R. Chen, G. A. Makhoul, and D. S. Meliksetian, โA queueing analysis of the performance of DQDB,โ IEEE/ACM Transactions on Networking, vol. 3, no. 6, pp. 872โ881, 1995. View at: Publisher Site  Google Scholar
 B. Mukherjee and S. Banerjee, โAlternative strategies for improving the fairness in an analytical model of DQDB networks,โ in Proceedings of the 10th Annual Joint Conference of the IEEE Computer and Communications Societies (INFOCOM '91), vol. 2, pp. 879โ888, Bal Harbour, Fla, USA, April 1991. View at: Publisher Site  Google Scholar
 R. N. Chiou and J. K. Guo, โThe effect of bursty lengths on DQDB networks,โ in Proceedings of the 1st International Conference on Complex, Intelligent and Software Intensive Systems (CISIS '07), pp. 85โ91, Vienna, Austria, April 2007. View at: Publisher Site  Google Scholar
 R. N. Chiou and J. K. Guo, โThe characteristics of queuing delays on DQDB networks,โ in Proceedings of the International Symposium of Photonic Devices and System Application, pp. 25โ26, Taipei, Taiwan, April 2006. View at: Google Scholar
 S. M. Banik, S. Radhakrishnan, T. Zheng, and C. N. Sekharan, โDistributed floor control protocols for computer collaborative applications on overlay networks,โ in Proceedings of the International Conference on Collaborative Computing: Networking, Applications and Worksharing, pp. 1โ10, San Jose, Calif, USA, December 2005. View at: Publisher Site  Google Scholar
 N. A. ELFishawy and L. M. Zawra, โPerformance evaluation of DQDB+/$\u0e42\x88\x92$ metropolitan area networks supporting isochronous and nonisochronous traffic,โ in Proceedings of the 20th National Radio Science Conference (NRSC '03), pp. 1โ8, Cairo, Egypt, March 2003. View at: Publisher Site  Google Scholar
 N. Funabiki, A. Sugano, and T. Higashino, โA proposal of a minimalstate processing search algorithm for isochronous channel reuse problems in DQDB networks,โ in Proceedings of the 17th IEEE International Conference on Advanced Information Networking and Applications (AINA '03), pp. 229โ232, Xian, China, March 2003. View at: Publisher Site  Google Scholar
 T.J. Kim and D.H. Cho, โBandwidth tuning for fairness of DQDB in clientserver traffic environments,โ IEEE Communications Letters, vol. 4, no. 6, pp. 208โ210, 2000. View at: Publisher Site  Google Scholar
 R. G. de Silva and W. J. Dewar, โCongestion control in DQDB/ATM interworking using PWM,โ in Proceedings of the IEEE Global Telecommunications Conference (GLOBECOM '98), vol. 3, pp. 1527โ1532, Sydney, Australia, November 1998. View at: Publisher Site  Google Scholar
 G. B. Brewster and M. K. Vernon, โThe fairness of DQDB networks with slot reuse,โ in Proceedings of the 14h Annual Joint Conference of the IEEE Computer and Communication Societies (INFOCOM '95), vol. 3, pp. 1154โ1163, Boston, Mass, USA, April 1995. View at: Publisher Site  Google Scholar
 W. M. Moh, Y.J. Chien, I. Zhang, and T. Moh, โDelay performance evaluation of high speed protocols for multimedia communications,โ in Proceedings of the 4th International Conference on Computer Communications and Networks (ICCCN '95), pp. 352โ355, Las Vegas, Nev, USA, September 1995. View at: Publisher Site  Google Scholar
 L. N. Kumar and A. D. Bovopoulos, โAn access protection solution for heavy load unfairness in DQDB,โ in Proceedings of the 11th Annual Conference of the IEEE Computer and Communications Societies (INFOCOM '92), vol. 1, pp. 190โ199, Florence, Italy, May 1992. View at: Publisher Site  Google Scholar
 S. Y. Cheung, โControlled request DQDB: achieving fairness and maximum throughput in the DQDB network,โ in Proceedings of the 11th Annual Conference of the IEEE Computer and Communications Societies (INFOCOM '92), vol. 1, pp. 180โ189, Florence, Italy, May 1992. View at: Publisher Site  Google Scholar
 M. Sumita, P. C. Fetterolf, and A. D. Little, โEffect of bandwidth balancing mechanism on fairness and performance of DQDB MANs,โ in Proceedings of the 11th Annual Conference of the IEEE Computer and Communications Societies (INFOCOM '92), vol. 3, pp. 1096โ2005, Florence, Italy, May 1992. View at: Publisher Site  Google Scholar
 K. H. Liang and Y. C. Chen, โAn approach for fairness improvement in DQDB networks,โ in Proceedings of the 3rd Workshop on Future Trends of Distributed Computing Systems, pp. 165โ171, Taipei, Taiwan, April 1992. View at: Publisher Site  Google Scholar
 E. L. Hahne, A. K. Choudhury, and N. F. Maxemchuk, โDQDB networks with and without bandwidth balancing,โ IEEE Transactions on Communications, vol. 40, no. 7, pp. 1192โ1204, 1992. View at: Publisher Site  Google Scholar
 J. Filipiak, โAccess protection for fairness in a distributed queue dual bus metropolitan area network,โ in Proceedings of the IEEE International Conference on Communications (ICC '89), vol. 2, pp. 635โ639, Boston, Mass, USA, June 1989. View at: Publisher Site  Google Scholar
 R. N. Chiou and C. N. Shyi, โThe performance of optical networks with TDMA,โ in Proceedings of the 2nd International MultiConference on Engineering and Technological Innovation (IMETI '09), Orlando, Fla, USA, July 2009. View at: Google Scholar
 R. Y. Rubenstein, Simulation and the MonteCarlo Method, John Wiley & Sons, New York, NY, USA, 1981.
Copyright
Copyright © 2009 Rong Nan Chiou and ChiaNian Shyi. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.