TCP throughput and timeout steady state and time-varying dynamics

Size: px
Start display at page:

Download "TCP throughput and timeout steady state and time-varying dynamics"

Transcription

1 TCP throughput an timeout steay state an time-varying ynamics Stephan Bohacek an Khushboo Shah Dept. of Electrical an Computer Engineering Dept. of Electrical Engineering University of Delaware University of Southern California Abstract While many moels of the TCP s ynamics have been evelope, few focus on the effects of timeout an high loss probability. Active queue management (AQM) is an important application of these ynamic moels. However, recent work has shown that AQM provies little performance benefit over rop-tail queueing for HTTP traffic, except possibly at high utilizations. It is at these utilizations that the ynamic moels of TCP are the least accurate. This paper presents a ynamic moel of TCP that accurately moels timeout. This moel is also applicable to the static case. This paper also presents a moel of the variance an the istribution of the congestion winow. It is shown that, while the ynamics of the mean value of the congestion winow are rather mil, the ynamics of timeout isplay large oscillations that take several secons to ecay. These oscillations cause the average bit-rate to also willy oscillate. Finally, this paper inclues results from several million simulations proviing a etaile view of the ynamics of timeout. I. Introuction s of TCP are extensively use throughout networking research. For example, TCP moels are use in fiels such as AQM [], [], [3], [4], esign of TCP frienly transport protocols [5], planning an provisioning networks [6], preicting file transfer time [7], [8], etc. In this way, TCP moels are founation for research. The accuracy of these moels is critical to the accuracy an relevance of the work that rests upon them. s for TCP have evolve over time. Initially, a simple relationship between the sening rate an the loss probability was evelope. The well-known "/sqrt(p)" formula, or more specifically MSS p 3// T T p, was foun an verifie [9], []. While this formula is useful in many settings, it was foun to not be accurate in situations where the packet loss probability is high. The reason for this inaccuracy is that this simple moel oes not account for TCP s timeout mechanism. Some stuies have shown that timeout, as oppose to triple uplicate acknowlegement, is a significant way that TCP etects packet loss []. In orer to inclue the impact of timeout, the simple This work was supporte in part by NSF grant Simulations for the work escribe in this paper were supporte by the University of Southern California Center for High Performance Computing an Communications. "square root of p" formula was extene to inclue the effect of timeout []. This formula was shown to be more accurate for large loss probabilities where timeout is significant. Sikar et al. [] evelope a moel for TCP throughput that inclues time-out. While this moel is ifferentfrom[],theyreportthatitgivessimilarnumerical results. However, we have performe extensive simulations, over 5 million in total, an have foun that when applie to the ns- implementation of TCP-SACK an when loss is probabilistic, the moel given in [] results in errors as large as a factor of two. Consiering the importance of the relationship between sening rate an loss probability, these ifferences are significant an require closer examination. This paper provies an accurate metho to estimate TCP s sening rate even for large loss probabilities. Another important extension to the TCP moels was to make them ynamic. That is, the initial moels assume that the loss probability was fixe. Inthecasethat the loss probability varie, ynamic moels were require. The first moel in this irection was given in [4], an several closely relate moels have followe [], [3], [7]. The principle application of these moels is for esign an analysis of AQM. These moels are extensively use to prove stability an to inicate the performance of AQM. Since AQM was first introuce, there have been a large number of researchers active in this area. ecently, several groups have investigate what, if any, benefit AQM has on the performance of the network (see [3] an references therein). A recent paper along these lines showe that the AQM, be it rop tail or any more elaborate version, has little impact on the performance of the network for utilization below 9%. However, for utilization above 9%, some AQM techniques were shown to improve the performance [4]. This work is significant as it inicates the environment where AQM might have a positive impact, specifically, in networks with high utilization. On the other han, this conclusion is troubling since the moels of TCP use to esign AQM neglect timeout. That is, in the environment where AQM might have the largest impact on performance, the analysis an esign of AQM is least accurate.

2 Because of this inaccuracy, the possible impact of AQM on network performance at high utilizations has never been carefully explore. One reason that AQM esign has not consiere high utilization networks is that ynamics of TCP in such an environment are not known. By ynamics, we mean the variation of the throughput in the face of a time-varying loss probability. This paper provies a simple moel of the ynamics of TCP s throughput that inclues timeout. Perhaps one of the most significant contributions of this paper is that the ynamics of timeout is quite complicate an can result in willy varying bit-rates that can take several secons to ecay. For example, the lower right of Figure 4 shows such oscillations. This ynamics is far more complicate than the ynamics of the mean congestion winow use by most AQM esign. As a result, this work calls into question the proof of stability of AQM at high loss probabilities. As a by-prouct of this work, the variance an the istribution of the congestion winow is also foun. Specifically, it is shown that the winow size is accurately moele with a negative binomial istribution. An finally the huge number of simulations use to valiate this moel provie a etaile view of the behavior of the ynamics of TCP s timeout mechanism. The paper procees as follows: The next section iscusses the simulations use for verification. In Section III the istribution of the congestion winow is foun. Next, in Section IV, a moel of the probability of a flow being in timeout is evelope. The results from Sections III an IV yiel a static moel of TCP throughput. In Section V, this moel is then extene to the ynamic case. The paper closes with some concluing remarks in Section VI. II. The simulations The goal of this effort is to closely examine the behavior of TCP. There is little oubt that in toay s networks the sening rate of a TCP connection is complicate by interacting flows, router inuce packet reorering, faulty loa balancing inuce packet reorering, server stalls, implementation iiosyncrasies, etc. While these effects are critical an a better unerstaning of them is necessary for a complete unerstaning of TCP, this paper focuses strictly on the behavior of TCP. The rationale for this is that in orer to unerstan the behavior of TCP in the wil, it is necessary to thoroughly unerstan TCP in a controlle environment. Perhaps the most significant ifference between these simulations an the environment foun in the wil is that in these simulations packet losses are entirely ranom an never ue to queue overflow. In the case of AQM, such simulations are appropriate as the goal of AQM is to rop packets in a controlle fashion, not when the queue fills. The simulations presente in this paper are for a single flow over a single bottleneck topology where rops at the bottleneck are ranom. The ns- error moel was use to prouce the rops. The results presente here are for a fixe roun-trip time of 3ms. Other roun-trip times have been explore an yiel the same results as the ones presente here. These simulations use ns- implementation of TCP-SACK in version.b8a. Delaye acknowlegement an a maximum receiver winow were not use. Our investigation into TCP inclue over 5 million simulations. This large number of simulations is require to estimate the ynamics of timeout. These simulations were performe on the University of Southern California s Linux cluster [5] an the large simulation results were transferre via the Internet to the University of Delaware for post-processing an storage. There is little oubt that without such high-performance computing, this work woul not be possible. III. The istribution of the congestion winow This section presents a simple moel for the istribution of the congestion winow. While the istribution may be of interest in its own right, we are particularly intereste in this istribution since it is use to etermine the probability of a flow entering timeout. We justify the moel briefly with some analysis an then present simulations results. Misra [4] introuce the iea of moeling TCP as a stochastic ifferential equation. Specifically, he suggeste that the winow size varies accoring to W t = T T W tn t () where N is a Cox process that counts the number of packet losses. This moel is for the evolution of the congestion winow when the flow is not in timeout. Thus, all probabilities shoul be conitione on the flow not in timeout. We enote this conition as TO. Thismoelwasfurther investigate in [7] where the partial ifferential equation of the winow size was foun. Specifically, if p w, t TO is the probability ensity of the congestion winow taking the value w at time t, thenp satisfies p w, t TO = Ã p w, t TO () t T T w +w δ (t T T ) 4p w, t TO p w, t TO, where δ is the loss probability. If δ is constant, then in p( w,t TO) steay state, i.e., t =, the istribution of the winow solves p w TO = w δ 4p w TO p w TO. (3) w From (3) it is straightforwar to show that E w m TO = c m δ m/,

3 for some constants c m. If m =, then this is the wellknown "square root of p" formula where c has been foun to be between. an.3. Consiering m =, we fin that the variance is of the form Var w TO := E w γ TO E w TO = δ, (4) where γ := c c. This relationship is borne out in simulations that inicates that γ.3. A close form solution has not been foun for (), whereas a complicate close form solution has been foun for (3). However, simple yet accurate approximate solution to (3) is the negative binomial istribution, i.e., p (w) = Γ (N + w ) Γ (N)(w )! ( q)n q w. Since the mean of this negative binomial ranom variable is Nq/( q)+an the variance is Nq/( q), q = E w TO Var w TO an N = q q E w TO (5) an from mean an variance of the w given above we get q = δ γ µ c + δ an N = ( q) q µ c. (6) δ In orer to etermine the istribution of the winow size, the parameters q an N must be etermine. This can be one by plugging the mean an variance into (5). However, if only the loss probability is known, then (6) must be use, in which case the parameters c an γ must beetermine. Wehavefounhatiftheobjectiveisto use the istribution to etermine the probability of timeout, the c = p 3/ an γ =.3 is sufficient. On the other han, if one is intereste in the istribution of the congestion winow (the objective of this section), then the selection of these parameters, especially c requires some more analysis. However, to keep this paper focuse, this issue is not iscusse. Once the parameters are etermine, the istribution of the congestion winow can be estimate. Figure shows the istribution of the congestion winow given the flow is not in the timeout state, i.e., p w TO. These figures shows the three types of curves, the observe istribution (histogram), the istribution with c = p 3/, anheistribution where the parameters are etermine by plugging in the observe mean an variance into (5). This last istribution coincies with the istribution given by (6) but with c =.7 for δ = 4 an c =.4 for δ =.5. For δ =., thevalueofc playslessofarole,c =.4, c = p 3/, anc =.7, or by using the observe mean an variance, all yiel nearly the same istribution. In all cases, γ =.3. probability. observe. observe.9 / /.9 / / (3/).8 /δ (3/).8 /δ.7.7/ δ / or.7 observe m.6 observe mean δ= -4.3 δ= w w probability Figure. Distribution of the congestion winow for loss probability of δ = 4 (left), δ =. (mile), an δ =.5 (right). IV. A moel of timeout In this section, we present a moel of the probability of a flow being in the timeout state. We say that a flow is in the timeout state if the next packet will be sent only after the retransmission timer expires. To unerstan the evelopment of this moel, one shoul consier a large collection of TCP experiments all running in parallel on ifferent networks with ientical network characteristics (e.g., link spees, propagation elays, etc.). Our goal is to etermine the fraction of these flows that are in timeout. A flow can enter the timeout state in the following three ways: IIfaflow experiences so many losses that triple uplicate acknowlegements are not receive, i.e., if at least max (w, ) losses occur in one winow. II In the case of the ns- implementation of TCP-SACK, if more than w/ packets are roppe. III If a retransmitte packet is roppe. Consiering I an II, we see that a flow will timeout if a single rop is followe by at least max (min (bw/c,w 3), ) rops out of the next w packets. Given that rops occur at a rate w δ,ropsthat lea to timeout occur at a rate ρ (w, δ) = w w δ X k=max(min(bw/c,w 3),) µ w δ k ( δ) w k. k Let λ (δ) betheratethataflow moves congestion avoiance to timeout ue to I an II. This rate is foun taking the expecte value of ρ (w, δ), i.e., λ (δ) =E ρ (w, δ) TO = X w ρ (w, δ) p δ (w), where p δ (w) is given in Section III, however, we have shown the epenence on δ. This epenence is This way of entering time-out is iscusse in []. FC-357 eliminates this way to enter time-out. It is straightforwar to ajust the evelopment below to reflect FC-357.

4 mae explicit by (6). Note that the approximation E ρ (w, δ) TO ρ E w TO,δ will result in large errors. Consier the situation where δ =.9 an, using c = p 3/, the mean value of the congestion winow is 4. Inthiscase,E ρ (w, δ) TO =.3 while ρ E w TO,δ =.38, aifference aroun a factor of three. Next we etermine λ, the rate that a flow moves from congestion avoiance to timeout because a retransmitte packet is roppe. The rate that a packet is first roppe is δ w. The rop will lea to a retransmission only if triple uplicate acknowlegements are receive, i.e., if less than max (min (bw/c,w 3), ) packets are roppe out of the next w packets. In this case, the retransmission is roppe an timeout is entere with probability δ. Thus ρ (w, δ) :=δ w δ ( w X k=max(min(bw/c,w 3),) µ w δ k ( δ) w k. k The average rate that a flow enters timeout ue to roppe retransmissions is λ (δ) :=E (ρ (w, δ)) = δ µ (δ) Finally, the rate that a TCP flow moves into timeout is λ + λ. Up to this point we consiere the rate that a single flow moves into timeout. Now we will etermine the fraction of flows entering timeout. Let I (t), enotetherate that flows enter timeout at time t. That is, in small time interval t, the probability that some flow enters timeout is I (t) t. Let T O bethetimeinwhichflows remain inthetimeoutstateuponthefirst timeout. If a flow exits timeout, but immeiately experiences another rop, then this flow again enters timeout, but for T O secons. We enote the rate that flows enter timeout for this secon time with I (t). The fraction of flows in timeout are t t T O I (τ) τ + t t T O I (τ) τ. The rate that flows enter timeout is the prouct of the fraction of flows not in timeout an the rate that such flows enter timeout, i.e., I (t) =(λ (δ (t)) + λ (δ (t))) (7) µ µz t Z t I (τ) τ + I (t) τ. t T O t T O Whiletheratethatflows exit timeout only to reenter timeout is I (t) =δi (t T O). (8) We on t consier the case where a flow enters time-out a thir time for a time interval 4 T O. While it is straightforwar to exten the calculations, it oes not appear to improve the accuracy of the moel. In steay state, I (t) an I (t) are constant. So (7) an (8) reuce to I =( I T O I T O)(λ (δ)+λ (δ)) I = δi. This can be solve λ (δ)+λ (δ) I = +T O (λ (δ)+λ (δ)) ( + δ) δ (λ (δ)+λ (δ)) I =. +T O (λ (δ)+λ (δ)) ( + δ) Thus, the fraction of flows in timeout is Z t Z t P (TO)= I (τ) τ + I (τ) τ (9) t T O t T O = ( + δ) T O (λ (δ)+λ (δ)) +(+δ) T O (λ (δ)+λ (δ)). If the loss probability varies, then (7), (8), an the relationship P (in timeout at time t) = t t T O I (τ) τ + t t T O I (τ) τ must be use to etermine the probability that a flow is in timeout. To verify (9), a large number of ns- simulations were performe so that the probability of timeout coul be etermine even for small loss probabilities. Figure shows the observe relationship of P (in timeout) an the loss probability as well as the estimate of this probability given by (9). These calculations use c = p 3/. Note that the fit is quite goo, the observe probability an the moel are essentially the same. Pahye et al. [] provies an estimate of the throughput that accounts for timeout. This formula also provies an estimate of the fraction of time a flow spens in timeout, i.e., P (TO). Specifically, [] fins P (timeout) where Q (W (δ),δ) G (δ) T o δ 8 (W (δ)+)+q(w (δ),δ) G (δ) T o, δ () ³+( 3 ³ w 3 δ) ( δ) Q (w, δ) =min, ( ( δ) w ) / ³ 3 ( δ) r 8( δ) W (δ) =+ + 3δ G (δ) =+δ +δ +4δ 3 +8δ 4 +6δ 5 +3δ 6. Figure inclues this estimate. We see that this moel, while proviing a qualitative fit in that the shape of the relationship is correct, it over estimates the probability of being in timeout. For example, for δ =%, the observe probability of being in timeout is 5%, while () gives an

5 probability of timeout observe / /.4 (3/) /δ.3 Pahye et al δ(loss probability) probability of timeout δ(loss probability) Figure. Probability of Timeout. The soli blue line represents the observe probability of a flow being in timeout. The ashe green line is the moel s estimate of this probability. Note that these two lines nearly coincie. The thir curve is the relationship erive from []. elative Error in Throughput (fraction).5 Pahye et al..4 FC δ(loss probability) Figure 3. elative Error of Preicte Throughput. error fraction of.5 implies a 5% error. Note that an estimate of 43%. However, this comparison must be qualifie. This paper focuses on TCP-SACK while the focus of [] is on TCP-ENO. Furthermore, this paper examines how TCP behaves uner ranom losses such as arises in AQM. In [], the queueing iscipline is rop tail. Furthermore, the loss probability in [] was not exactly the probability of a packet being roppe, but more along the lines of the probability of a rop event occurring where a rop event leas to the rest of the packets in the "roun" being roppe. Thus, the δ use in () is less than the packet loss probability. However, this moel is useful as a benchmark as it was not known how accurate this moel woul be for TCP-SACK. From the probability of timeout, the average throughput can be foun via T = c ( P (TO)) MSS. δ Figure 3 shows the relative error of this moel when compare to the observe throughput. Figure 3 also shows the relative error of the moel in [] as well as the moel suggeste in FC-3448 for a TCP frienly sening rate. The figure shows only the relative error, an examination of the error shows that the moel of [] an FC-3448 unerestimate the sening rate (this is ue to the over estimation of the timeout probability). As mentione, these other moels are for benchmarking purposes only. Figure shows that for even moerately large loss probabilities, the probability of being in timeout is quite high. E(W(t)) Prob(timeout at time) variance(w(t)) bit-rate Mean Only - 34 Figure 4. Simulation an moeling of a jump in loss probability from % to %. The jump occurre at t =. The soli blue curves show the observe an the ashe re curves shows the moel. In the lower right, the otte green inicates the bit-rate if only the mean size of the congestion winow is use (i.e., the ashe re curve in the upper left). Note the units of the bit-rate is in bps. Specifically, for δ>%, timeout plays an important role in the sening rate of TCP. For δ>5%, aflow will spen 6% of the time in timeout. Thus, timeout ominates the sening rate of TCP. There are many implications of this. ecently there have been some propose moifications to TCP such as limite transmit an ECN. These moifications significantly affect the behavior of timeout, specifically, they enter timeout much less frequently. Thus, for δ>%, one can expect that TCP implementations with these moifications will be much more aggressive than implementations without these moifications. There has been extensive work in ynamic moeling of TCP, but most has neglecte timeout. For example, many approaches to AQM utilize a moel of TCP that neglects timeout. We see that such an approach is only reasonable if the loss probability is less than %. However, in the area of AQM, it is hope that these moels are applicable to higher loss probabilities. Next we evelop a moel for the ynamic behavior of TCP that inclues timeout. V. Dynamics of TCP Here we exten the ynamics that is commonly use to moel the ynamics of TCP. The commonly use moel is for only the mean sening rate of TCP. Here we present a moel for the variance. Once the variance is known, the istribution is foun as in Section III. With the istribution, we can fin the probability of being in timeout as was one in the previous section. Of course, the paper to this point has only consiere a constant loss probability. Here we allow the loss probability to vary with time. Thus, this moel is appropriate for AQM esign.

6 8 6 5 variance(w(t)) E(W(t)) bit-rate Prob(timeout at time) Mean Only variance(w(t)) Mean Only bit-rate Prob(timeout at time) E(W(t)) Figure 5. These figures are similar to the ones shown in Figure 4. However, here two sets of ynamics are shown, the loss probability jumping from.% to % an falling from % to.%. In both cases, the changes occur at t = Figure 6. Valiation of the. These figures are near ientical to Figure 5, however, here the loss probability jumps between % an 5%. Let w (t) = E w (t) TO an w (t) = E w (t) TO. An approximation of the ynamics of the mean value of the congestion winow is often given by w (t) = δ (t ) w (t ) w (t). c () This is often approximate as w (t) = δ (t ) w (t). c () Interestingly, this formula is often erive from (). However, if proper stochastic calculus is applie, the correct ynamics for the mean are w (t) = δ (t ) E (w (t ) w (t)), which can be approximate as w (t) = δ (t ) w (t). (3) Note that in (3) the secon moment is utilize, while () have c multiplie by the first moment square. For tunately, we have foun that () an () give goo approximation to the ynamics of the mean. It seems that the mean varies slowly enough that this ifference between w (t) an c w (t) is not significant. From () it is possible to etermine the ynamics of the secon moment of the congestion winow, 3 w (t) = w (t) δ (t ) E w3 (t). 4 Using the same iea as above, we approximate 3/ ³ an arrive at E w3 (t) 83 β 3/ w (t) (c.3) w (t) = w (t) (4) ³ 3/ 38 β. δ (t ) w (t) 3/ 4 3 (c.3) As in the case of (), this approximation is exact in steay state. Equations () an (4) form a system of orinary ifferential equations (ODEs). With the mean an secon moment, the istribution of the congestion winow can be foun as iscusse in Section III. With this approximation of the istribution, the probability of being in timeout can be foun as in Section IV, i.e., P (T O at time t) = t t I (τ ) τ + t T O I (τ ) τ, where I (τ ) is given t T O by (7) an I (τ ) is given by (8). Figures 4-6 provie some valiation of this ynamic moel of³ TCP. Through extensive simulation, E w (t) T O, E w (t) T O, an P (T O at time t) were foun uner ifferent time variations of the loss probability. Specifically, in each simulation, the loss probability took two values δ an δ. The simulation began an ran for secons with loss probability set to δ. In these simulations, secons was enough to ensure that the system was in steay state (i.e., E w (t) T O was constant). Then, secons after the beginning of the simulation, the loss probability was switche to δ. In the figures, this moment is labele as t =. This simulation was carrie out, times. With these simulations it is possi ble³ to etermine the value of the statistics E w (t) T O, E w (t) T O, an P (T O at time t) at each time point.

7 For example, to etermine P (TO at time 5 ms), we foun the fraction of the, flows that were in the timeout state at time t = 5 ms. This large number of simulations allows for high confience in the observations an allows the etails of the ynamics to be observe. Figures 4-6 inicate that the ynamic moel is quite accurate. They also show, as expecte, that using the ynamics of the mean congestion winow given by () is not aequate for large loss probability, but is aequate for small loss probability. We also see that the ynamics for the mean an the variance of the congestion winow perform rather poorly for very large loss probability. However, in this case, the behavior of TCP is completely ominate by the timeout, which is well moele. As a result, the bit-rate is accurately moele. Perhaps the most striking aspect of Figures 5-6 is the ynamics of the bit-rate at large loss probabilities. Most significantly, we see that the bit-rate experiences wil oscillations that can take several secons to ecay. In all cases, the moel for the mean congestion winow () oes not prouce any oscillations. To the best of the authors knowlege, this behavior has not been previously observe, much less moele. Oscillating "step" responses are often a sign of instability, thus it is likely that these ynamics will impact the stability results of previous AQM research. Figure 6 shows some unexpecte behavior of the mean an the variance of the congestion winow. Specifically, we see that the mean has two istinct phases of growth with a transition at t = ms while the variance shows non-monotonic growth. The ODEs () an (4) o not moel this behavior. This behavior is ue to the feeback of the fraction of flows in timeout to the mean an the variance of the congestion winow. Notice that the nonmonotonic growth of w occurs just as flow exit timeout. Since the probability of a flow being in the timeout state is so high, a slight error in the mean congestion winow has a minor impact on the average bit-rate. VI. Conclusions We have evelope a new moel for the sening rate of TCP that inclues the effect of timeout. Previous moels have inclue the effect of timeout, but we have foun these moels to greatly over estimate the probability of being in timeout. The moel presente is easily extene to a ynamic moel of timeout. This ynamic moel was shown to agree with simulations quite well. Furthermore, it was shown that timeout can cause the bit-rate to oscillate substantially for several secons after a change in the loss probability. It was also shown that for large loss probabilities, the behavior of TCP is ominate by the behavior of timeout. This moel is significantly ifferent from those use in previous esign of AQM. Future work will examine the possibility of esigning an AQM that utilizes these moels for TCP. Perhaps, more substantial performance gains will be achieve for large loss probability. eferences [] S.Athuraliya,V.H.Li,S.H.Low,anQ.Yin, EM:Active queue management, IEEE Network, vol. 5, pp ,. [] C. V. Hollot, V. Misra, D. F. Towsley, an W. Gong, On esigning improve controllers for AQM routers supporting TCP flows, in INFOCOM,. [3] S. Kunniyur an. Srikant, Analysis an esign of an aaptive virtual queue (AVQ) algorithm for active queue management, in Proceeings of ACM SIGCOMM,. [4] V. Misra, W. Gong, an D. Towsley, Stochastic ifferential equation moeling an analysis of tcp-winowsize behavior, in Proceeings of PEFOMANCE99, 999. [5] M.Hanley,S.Floy,J.Pahye,anJ.Wimer, FC TCP frienly rate control (TFC): Protocol specification, 3. [6] J.Beckers,I.Henrawan,.Kooij,an.V.erMei, Generalize processor sharing performance moels for internet access lines, in 9th IFIP CONFEENCE on PEFOMANCE MODELLING AND EVALUATION OF ATM an IP NET- WOKS, (Buapest),. [7] S. Bohacek, A stochastic moel of tcp an fair vieo transmission, in Infocom, 3. [8] N. Carwell, S. Savage, an T. Anerson, ing TCP latency, in INFOCOM (3), pp ,. [9] T.Ott,J.H.B.Kemperman,anM.Mathis, Winowsize behavior in TCP/IP with constant loss probability, in Proc. Of the DIMACS Workshop on Performance of ealtime Applications on the Internet, 996. [] M. Mathis, J. Semke, J. Mahavi, an T. Ott, The macroscopic behavior of the TCP congestion avoiance algorithm, ACM Comput. Comm. eview, vol. 7, no. 3, 997. [] J. Pahye, V. Firoiu, D. Towsley, an J. Krusoe, ing TCP throughput: A simple moel an its empirical valiation, Proceeings of the ACM SIGCOMM 98 conference on Applications, technologies, architectures, an protocols for computer communication, pp , 998. [] K. S. V. Biplab Sikar, Shivkumar Kalyanaraman, Analytic moels for the latency an steay-state throughput of TCP tahoe, reno, an SACK, IEEE/ACM Transactions on Networking, vol.,no.6,pp ,3. [3] S. Bohacek, K. Shah, G. Arce, an M. Davis, A review of central issues - signal processing challenges in active queue management, IEEE Signal Processing Magazine, vol., pp , 4. [4] L. Le, J. Aikat, K. Jeffay, an F. D. Smith, The effects of active queue management on web performance, in SIGCOMM, 3. [5]

TCP throughput and timeout steady state and time-varying dynamics

TCP throughput and timeout steady state and time-varying dynamics TCP throughput an timeout steay state an time-varying ynamics Stephan Bohacek an Khushboo Shah Dept. of Electrical an Computer Engineering Dept. of Electrical Engineering University of Delaare University

More information

Modeling RED with Idealized TCP Sources

Modeling RED with Idealized TCP Sources Moeling RED with Iealize TCP Sources P. Kuusela, P. Lassila, J. Virtamo an P. Key Networking Laboratory Helsinki University of Technology (HUT) P.O.Box 3000, FIN-02015 HUT, Finlan Email: {Pirkko.Kuusela,

More information

Survey Sampling. 1 Design-based Inference. Kosuke Imai Department of Politics, Princeton University. February 19, 2013

Survey Sampling. 1 Design-based Inference. Kosuke Imai Department of Politics, Princeton University. February 19, 2013 Survey Sampling Kosuke Imai Department of Politics, Princeton University February 19, 2013 Survey sampling is one of the most commonly use ata collection methos for social scientists. We begin by escribing

More information

On Characterizing the Delay-Performance of Wireless Scheduling Algorithms

On Characterizing the Delay-Performance of Wireless Scheduling Algorithms On Characterizing the Delay-Performance of Wireless Scheuling Algorithms Xiaojun Lin Center for Wireless Systems an Applications School of Electrical an Computer Engineering, Purue University West Lafayette,

More information

Math 342 Partial Differential Equations «Viktor Grigoryan

Math 342 Partial Differential Equations «Viktor Grigoryan Math 342 Partial Differential Equations «Viktor Grigoryan 6 Wave equation: solution In this lecture we will solve the wave equation on the entire real line x R. This correspons to a string of infinite

More information

Homework 7 Due 18 November at 6:00 pm

Homework 7 Due 18 November at 6:00 pm Homework 7 Due 18 November at 6:00 pm 1. Maxwell s Equations Quasi-statics o a An air core, N turn, cylinrical solenoi of length an raius a, carries a current I Io cos t. a. Using Ampere s Law, etermine

More information

Improved Rate-Based Pull and Push Strategies in Large Distributed Networks

Improved Rate-Based Pull and Push Strategies in Large Distributed Networks Improve Rate-Base Pull an Push Strategies in Large Distribute Networks Wouter Minnebo an Benny Van Hout Department of Mathematics an Computer Science University of Antwerp - imins Mielheimlaan, B-00 Antwerp,

More information

Lecture 2 Lagrangian formulation of classical mechanics Mechanics

Lecture 2 Lagrangian formulation of classical mechanics Mechanics Lecture Lagrangian formulation of classical mechanics 70.00 Mechanics Principle of stationary action MATH-GA To specify a motion uniquely in classical mechanics, it suffices to give, at some time t 0,

More information

Time-of-Arrival Estimation in Non-Line-Of-Sight Environments

Time-of-Arrival Estimation in Non-Line-Of-Sight Environments 2 Conference on Information Sciences an Systems, The Johns Hopkins University, March 2, 2 Time-of-Arrival Estimation in Non-Line-Of-Sight Environments Sinan Gezici, Hisashi Kobayashi an H. Vincent Poor

More information

Linear First-Order Equations

Linear First-Order Equations 5 Linear First-Orer Equations Linear first-orer ifferential equations make up another important class of ifferential equations that commonly arise in applications an are relatively easy to solve (in theory)

More information

Numerical Key Results for Studying the Performance of TCP Traffic under Class-Based Weighted Fair Queuing System

Numerical Key Results for Studying the Performance of TCP Traffic under Class-Based Weighted Fair Queuing System Journal of ommunication an omputer 13 (2016) 195-201 oi:10.17265/1548-7709/2016.04.007 D DAVID PUBLISHING Numerical Key Results for Stuying the Performance of TP Traffic uner lass-base Weighte Fair Queuing

More information

'HVLJQ &RQVLGHUDWLRQ LQ 0DWHULDO 6HOHFWLRQ 'HVLJQ 6HQVLWLYLW\,1752'8&7,21

'HVLJQ &RQVLGHUDWLRQ LQ 0DWHULDO 6HOHFWLRQ 'HVLJQ 6HQVLWLYLW\,1752'8&7,21 Large amping in a structural material may be either esirable or unesirable, epening on the engineering application at han. For example, amping is a esirable property to the esigner concerne with limiting

More information

Thermal conductivity of graded composites: Numerical simulations and an effective medium approximation

Thermal conductivity of graded composites: Numerical simulations and an effective medium approximation JOURNAL OF MATERIALS SCIENCE 34 (999)5497 5503 Thermal conuctivity of grae composites: Numerical simulations an an effective meium approximation P. M. HUI Department of Physics, The Chinese University

More information

ensembles When working with density operators, we can use this connection to define a generalized Bloch vector: v x Tr x, v y Tr y

ensembles When working with density operators, we can use this connection to define a generalized Bloch vector: v x Tr x, v y Tr y Ph195a lecture notes, 1/3/01 Density operators for spin- 1 ensembles So far in our iscussion of spin- 1 systems, we have restricte our attention to the case of pure states an Hamiltonian evolution. Toay

More information

Sliding mode approach to congestion control in connection-oriented communication networks

Sliding mode approach to congestion control in connection-oriented communication networks JOURNAL OF APPLIED COMPUTER SCIENCE Vol. xx. No xx (200x), pp. xx-xx Sliing moe approach to congestion control in connection-oriente communication networks Anrzej Bartoszewicz, Justyna Żuk Technical University

More information

Lectures - Week 10 Introduction to Ordinary Differential Equations (ODES) First Order Linear ODEs

Lectures - Week 10 Introduction to Ordinary Differential Equations (ODES) First Order Linear ODEs Lectures - Week 10 Introuction to Orinary Differential Equations (ODES) First Orer Linear ODEs When stuying ODEs we are consiering functions of one inepenent variable, e.g., f(x), where x is the inepenent

More information

Schrödinger s equation.

Schrödinger s equation. Physics 342 Lecture 5 Schröinger s Equation Lecture 5 Physics 342 Quantum Mechanics I Wenesay, February 3r, 2010 Toay we iscuss Schröinger s equation an show that it supports the basic interpretation of

More information

THE VAN KAMPEN EXPANSION FOR LINKED DUFFING LINEAR OSCILLATORS EXCITED BY COLORED NOISE

THE VAN KAMPEN EXPANSION FOR LINKED DUFFING LINEAR OSCILLATORS EXCITED BY COLORED NOISE Journal of Soun an Vibration (1996) 191(3), 397 414 THE VAN KAMPEN EXPANSION FOR LINKED DUFFING LINEAR OSCILLATORS EXCITED BY COLORED NOISE E. M. WEINSTEIN Galaxy Scientific Corporation, 2500 English Creek

More information

Generalization of the persistent random walk to dimensions greater than 1

Generalization of the persistent random walk to dimensions greater than 1 PHYSICAL REVIEW E VOLUME 58, NUMBER 6 DECEMBER 1998 Generalization of the persistent ranom walk to imensions greater than 1 Marián Boguñá, Josep M. Porrà, an Jaume Masoliver Departament e Física Fonamental,

More information

TCP has been widely successful on the wired Internet. Stability and Delay Consideration for Flow Control over Wireless Networks

TCP has been widely successful on the wired Internet. Stability and Delay Consideration for Flow Control over Wireless Networks Stability an Delay Consieration for Flow Control over Wireless Networks Minghua Chen, Alessanro Abate, Avieh Zakhor, Shankar Sastry Department of Electrical Engineering an Computer Sciences University

More information

IPA Derivatives for Make-to-Stock Production-Inventory Systems With Backorders Under the (R,r) Policy

IPA Derivatives for Make-to-Stock Production-Inventory Systems With Backorders Under the (R,r) Policy IPA Derivatives for Make-to-Stock Prouction-Inventory Systems With Backorers Uner the (Rr) Policy Yihong Fan a Benamin Melame b Yao Zhao c Yorai Wari Abstract This paper aresses Infinitesimal Perturbation

More information

Switching Time Optimization in Discretized Hybrid Dynamical Systems

Switching Time Optimization in Discretized Hybrid Dynamical Systems Switching Time Optimization in Discretize Hybri Dynamical Systems Kathrin Flaßkamp, To Murphey, an Sina Ober-Blöbaum Abstract Switching time optimization (STO) arises in systems that have a finite set

More information

State observers and recursive filters in classical feedback control theory

State observers and recursive filters in classical feedback control theory State observers an recursive filters in classical feeback control theory State-feeback control example: secon-orer system Consier the riven secon-orer system q q q u x q x q x x x x Here u coul represent

More information

Fluid model for a data network with α-fair bandwidth sharing and general document size distributions: two examples of stability

Fluid model for a data network with α-fair bandwidth sharing and general document size distributions: two examples of stability IMS Lecture Notes Monograph Series??? Vol. 0 (0000) 1 c Institute of Mathematical Statistics, 0000 Flui moel for a ata network with α-fair banwith sharing an general ocument size istributions: two examples

More information

Transmission Line Matrix (TLM) network analogues of reversible trapping processes Part B: scaling and consistency

Transmission Line Matrix (TLM) network analogues of reversible trapping processes Part B: scaling and consistency Transmission Line Matrix (TLM network analogues of reversible trapping processes Part B: scaling an consistency Donar e Cogan * ANC Eucation, 308-310.A. De Mel Mawatha, Colombo 3, Sri Lanka * onarecogan@gmail.com

More information

Implicit Differentiation

Implicit Differentiation Implicit Differentiation Thus far, the functions we have been concerne with have been efine explicitly. A function is efine explicitly if the output is given irectly in terms of the input. For instance,

More information

Perturbation Analysis and Optimization of Stochastic Flow Networks

Perturbation Analysis and Optimization of Stochastic Flow Networks IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL. XX, NO. Y, MMM 2004 1 Perturbation Analysis an Optimization of Stochastic Flow Networks Gang Sun, Christos G. Cassanras, Yorai Wari, Christos G. Panayiotou,

More information

APPROXIMATE SOLUTION FOR TRANSIENT HEAT TRANSFER IN STATIC TURBULENT HE II. B. Baudouy. CEA/Saclay, DSM/DAPNIA/STCM Gif-sur-Yvette Cedex, France

APPROXIMATE SOLUTION FOR TRANSIENT HEAT TRANSFER IN STATIC TURBULENT HE II. B. Baudouy. CEA/Saclay, DSM/DAPNIA/STCM Gif-sur-Yvette Cedex, France APPROXIMAE SOLUION FOR RANSIEN HEA RANSFER IN SAIC URBULEN HE II B. Bauouy CEA/Saclay, DSM/DAPNIA/SCM 91191 Gif-sur-Yvette Ceex, France ABSRAC Analytical solution in one imension of the heat iffusion equation

More information

inflow outflow Part I. Regular tasks for MAE598/494 Task 1

inflow outflow Part I. Regular tasks for MAE598/494 Task 1 MAE 494/598, Fall 2016 Project #1 (Regular tasks = 20 points) Har copy of report is ue at the start of class on the ue ate. The rules on collaboration will be release separately. Please always follow the

More information

SYNCHRONOUS SEQUENTIAL CIRCUITS

SYNCHRONOUS SEQUENTIAL CIRCUITS CHAPTER SYNCHRONOUS SEUENTIAL CIRCUITS Registers an counters, two very common synchronous sequential circuits, are introuce in this chapter. Register is a igital circuit for storing information. Contents

More information

The Principle of Least Action

The Principle of Least Action Chapter 7. The Principle of Least Action 7.1 Force Methos vs. Energy Methos We have so far stuie two istinct ways of analyzing physics problems: force methos, basically consisting of the application of

More information

Calculus in the AP Physics C Course The Derivative

Calculus in the AP Physics C Course The Derivative Limits an Derivatives Calculus in the AP Physics C Course The Derivative In physics, the ieas of the rate change of a quantity (along with the slope of a tangent line) an the area uner a curve are essential.

More information

The influence of the equivalent hydraulic diameter on the pressure drop prediction of annular test section

The influence of the equivalent hydraulic diameter on the pressure drop prediction of annular test section IOP Conference Series: Materials Science an Engineering PAPER OPEN ACCESS The influence of the equivalent hyraulic iameter on the pressure rop preiction of annular test section To cite this article: A

More information

A Better Model for Job Redundancy: Decoupling Server Slowdown and Job Size

A Better Model for Job Redundancy: Decoupling Server Slowdown and Job Size A Better Moel for Job Reunancy: Decoupling Server Slowown an Job Size Kristen Garner 1 Mor Harchol-Balter 1 Alan Scheller-Wolf Benny Van Hout 3 October 1 CMU-CS-1-131 School of Computer Science Carnegie

More information

A Fair Comparison of Pull and Push Strategies in Large Distributed Networks

A Fair Comparison of Pull and Push Strategies in Large Distributed Networks 1 A Fair Comparison of Pull an Push Strategies in Large Distribute Networks Wouter Minnebo an Benny Van Hout, Member, IEEE Abstract In this paper we compare the performance of the pull an push strategy

More information

Table of Common Derivatives By David Abraham

Table of Common Derivatives By David Abraham Prouct an Quotient Rules: Table of Common Derivatives By Davi Abraham [ f ( g( ] = [ f ( ] g( + f ( [ g( ] f ( = g( [ f ( ] g( g( f ( [ g( ] Trigonometric Functions: sin( = cos( cos( = sin( tan( = sec

More information

An M/G/1 Retrial Queue with Priority, Balking and Feedback Customers

An M/G/1 Retrial Queue with Priority, Balking and Feedback Customers Journal of Convergence Information Technology Volume 5 Number April 1 An M/G/1 Retrial Queue with Priority Balking an Feeback Customers Peishu Chen * 1 Yiuan Zhu 1 * 1 Faculty of Science Jiangsu University

More information

Closed and Open Loop Optimal Control of Buffer and Energy of a Wireless Device

Closed and Open Loop Optimal Control of Buffer and Energy of a Wireless Device Close an Open Loop Optimal Control of Buffer an Energy of a Wireless Device V. S. Borkar School of Technology an Computer Science TIFR, umbai, Inia. borkar@tifr.res.in A. A. Kherani B. J. Prabhu INRIA

More information

d dx But have you ever seen a derivation of these results? We ll prove the first result below. cos h 1

d dx But have you ever seen a derivation of these results? We ll prove the first result below. cos h 1 Lecture 5 Some ifferentiation rules Trigonometric functions (Relevant section from Stewart, Seventh Eition: Section 3.3) You all know that sin = cos cos = sin. () But have you ever seen a erivation of

More information

05 The Continuum Limit and the Wave Equation

05 The Continuum Limit and the Wave Equation Utah State University DigitalCommons@USU Founations of Wave Phenomena Physics, Department of 1-1-2004 05 The Continuum Limit an the Wave Equation Charles G. Torre Department of Physics, Utah State University,

More information

Math Notes on differentials, the Chain Rule, gradients, directional derivative, and normal vectors

Math Notes on differentials, the Chain Rule, gradients, directional derivative, and normal vectors Math 18.02 Notes on ifferentials, the Chain Rule, graients, irectional erivative, an normal vectors Tangent plane an linear approximation We efine the partial erivatives of f( xy, ) as follows: f f( x+

More information

ELECTRON DIFFRACTION

ELECTRON DIFFRACTION ELECTRON DIFFRACTION Electrons : wave or quanta? Measurement of wavelength an momentum of electrons. Introuction Electrons isplay both wave an particle properties. What is the relationship between the

More information

Role of parameters in the stochastic dynamics of a stick-slip oscillator

Role of parameters in the stochastic dynamics of a stick-slip oscillator Proceeing Series of the Brazilian Society of Applie an Computational Mathematics, v. 6, n. 1, 218. Trabalho apresentao no XXXVII CNMAC, S.J. os Campos - SP, 217. Proceeing Series of the Brazilian Society

More information

Construction of the Electronic Radial Wave Functions and Probability Distributions of Hydrogen-like Systems

Construction of the Electronic Radial Wave Functions and Probability Distributions of Hydrogen-like Systems Construction of the Electronic Raial Wave Functions an Probability Distributions of Hyrogen-like Systems Thomas S. Kuntzleman, Department of Chemistry Spring Arbor University, Spring Arbor MI 498 tkuntzle@arbor.eu

More information

Linear analysis of a natural circulation driven supercritical water loop

Linear analysis of a natural circulation driven supercritical water loop TU Delft Bachelor Thesis Linear analysis of a natural circulation riven supercritical water loop D J van er Ham 4285816 supervise by Dr. Ir. M. Rohe July 3, 216 Nomenclature Symbol Units Description A

More information

Chapter 6: Energy-Momentum Tensors

Chapter 6: Energy-Momentum Tensors 49 Chapter 6: Energy-Momentum Tensors This chapter outlines the general theory of energy an momentum conservation in terms of energy-momentum tensors, then applies these ieas to the case of Bohm's moel.

More information

IERCU. Institute of Economic Research, Chuo University 50th Anniversary Special Issues. Discussion Paper No.210

IERCU. Institute of Economic Research, Chuo University 50th Anniversary Special Issues. Discussion Paper No.210 IERCU Institute of Economic Research, Chuo University 50th Anniversary Special Issues Discussion Paper No.210 Discrete an Continuous Dynamics in Nonlinear Monopolies Akio Matsumoto Chuo University Ferenc

More information

Calculus Class Notes for the Combined Calculus and Physics Course Semester I

Calculus Class Notes for the Combined Calculus and Physics Course Semester I Calculus Class Notes for the Combine Calculus an Physics Course Semester I Kelly Black December 14, 2001 Support provie by the National Science Founation - NSF-DUE-9752485 1 Section 0 2 Contents 1 Average

More information

Analytic Scaling Formulas for Crossed Laser Acceleration in Vacuum

Analytic Scaling Formulas for Crossed Laser Acceleration in Vacuum October 6, 4 ARDB Note Analytic Scaling Formulas for Crosse Laser Acceleration in Vacuum Robert J. Noble Stanfor Linear Accelerator Center, Stanfor University 575 San Hill Roa, Menlo Park, California 945

More information

Computing Exact Confidence Coefficients of Simultaneous Confidence Intervals for Multinomial Proportions and their Functions

Computing Exact Confidence Coefficients of Simultaneous Confidence Intervals for Multinomial Proportions and their Functions Working Paper 2013:5 Department of Statistics Computing Exact Confience Coefficients of Simultaneous Confience Intervals for Multinomial Proportions an their Functions Shaobo Jin Working Paper 2013:5

More information

Similarity Measures for Categorical Data A Comparative Study. Technical Report

Similarity Measures for Categorical Data A Comparative Study. Technical Report Similarity Measures for Categorical Data A Comparative Stuy Technical Report Department of Computer Science an Engineering University of Minnesota 4-92 EECS Builing 200 Union Street SE Minneapolis, MN

More information

Introduction to the Vlasov-Poisson system

Introduction to the Vlasov-Poisson system Introuction to the Vlasov-Poisson system Simone Calogero 1 The Vlasov equation Consier a particle with mass m > 0. Let x(t) R 3 enote the position of the particle at time t R an v(t) = ẋ(t) = x(t)/t its

More information

2 Mathematical moelling of the Internet ing that, with an appropriate formulation of an overall optimization problem, the network's implicit objective

2 Mathematical moelling of the Internet ing that, with an appropriate formulation of an overall optimization problem, the network's implicit objective 1 Mathematical moelling of the Internet Frank Kelly Statistical Laboratory, University of Cambrige, 16 Mill Lane, Cambrige CB2 1SB, UK Abstract Moern communication networks are able to respon to ranomly

More information

Math 1B, lecture 8: Integration by parts

Math 1B, lecture 8: Integration by parts Math B, lecture 8: Integration by parts Nathan Pflueger 23 September 2 Introuction Integration by parts, similarly to integration by substitution, reverses a well-known technique of ifferentiation an explores

More information

ON THE OPTIMALITY SYSTEM FOR A 1 D EULER FLOW PROBLEM

ON THE OPTIMALITY SYSTEM FOR A 1 D EULER FLOW PROBLEM ON THE OPTIMALITY SYSTEM FOR A D EULER FLOW PROBLEM Eugene M. Cliff Matthias Heinkenschloss y Ajit R. Shenoy z Interisciplinary Center for Applie Mathematics Virginia Tech Blacksburg, Virginia 46 Abstract

More information

A Better Model for Job Redundancy: Decoupling Server Slowdown and Job Size

A Better Model for Job Redundancy: Decoupling Server Slowdown and Job Size A Better Moel for Job Reunancy: Decoupling Server Slowown an Job Size Kristen Garner, Mor Harchol-Balter, Alan Scheller-Wolf, an Benny Van Hout 1 Abstract Recent computer systems research has propose using

More information

arxiv:physics/ v2 [physics.ed-ph] 23 Sep 2003

arxiv:physics/ v2 [physics.ed-ph] 23 Sep 2003 Mass reistribution in variable mass systems Célia A. e Sousa an Vítor H. Rorigues Departamento e Física a Universiae e Coimbra, P-3004-516 Coimbra, Portugal arxiv:physics/0211075v2 [physics.e-ph] 23 Sep

More information

Equilibrium in Queues Under Unknown Service Times and Service Value

Equilibrium in Queues Under Unknown Service Times and Service Value University of Pennsylvania ScholarlyCommons Finance Papers Wharton Faculty Research 1-2014 Equilibrium in Queues Uner Unknown Service Times an Service Value Laurens Debo Senthil K. Veeraraghavan University

More information

The total derivative. Chapter Lagrangian and Eulerian approaches

The total derivative. Chapter Lagrangian and Eulerian approaches Chapter 5 The total erivative 51 Lagrangian an Eulerian approaches The representation of a flui through scalar or vector fiels means that each physical quantity uner consieration is escribe as a function

More information

Systems & Control Letters

Systems & Control Letters Systems & ontrol Letters ( ) ontents lists available at ScienceDirect Systems & ontrol Letters journal homepage: www.elsevier.com/locate/sysconle A converse to the eterministic separation principle Jochen

More information

PD Controller for Car-Following Models Based on Real Data

PD Controller for Car-Following Models Based on Real Data PD Controller for Car-Following Moels Base on Real Data Xiaopeng Fang, Hung A. Pham an Minoru Kobayashi Department of Mechanical Engineering Iowa State University, Ames, IA 5 Hona R&D The car following

More information

A Modification of the Jarque-Bera Test. for Normality

A Modification of the Jarque-Bera Test. for Normality Int. J. Contemp. Math. Sciences, Vol. 8, 01, no. 17, 84-85 HIKARI Lt, www.m-hikari.com http://x.oi.org/10.1988/ijcms.01.9106 A Moification of the Jarque-Bera Test for Normality Moawa El-Fallah Ab El-Salam

More information

Least-Squares Regression on Sparse Spaces

Least-Squares Regression on Sparse Spaces Least-Squares Regression on Sparse Spaces Yuri Grinberg, Mahi Milani Far, Joelle Pineau School of Computer Science McGill University Montreal, Canaa {ygrinb,mmilan1,jpineau}@cs.mcgill.ca 1 Introuction

More information

The Exact Form and General Integrating Factors

The Exact Form and General Integrating Factors 7 The Exact Form an General Integrating Factors In the previous chapters, we ve seen how separable an linear ifferential equations can be solve using methos for converting them to forms that can be easily

More information

Qubit channels that achieve capacity with two states

Qubit channels that achieve capacity with two states Qubit channels that achieve capacity with two states Dominic W. Berry Department of Physics, The University of Queenslan, Brisbane, Queenslan 4072, Australia Receive 22 December 2004; publishe 22 March

More information

Quantum optics of a Bose-Einstein condensate coupled to a quantized light field

Quantum optics of a Bose-Einstein condensate coupled to a quantized light field PHYSICAL REVIEW A VOLUME 60, NUMBER 2 AUGUST 1999 Quantum optics of a Bose-Einstein conensate couple to a quantize light fiel M. G. Moore, O. Zobay, an P. Meystre Optical Sciences Center an Department

More information

Chapter 4. Electrostatics of Macroscopic Media

Chapter 4. Electrostatics of Macroscopic Media Chapter 4. Electrostatics of Macroscopic Meia 4.1 Multipole Expansion Approximate potentials at large istances 3 x' x' (x') x x' x x Fig 4.1 We consier the potential in the far-fiel region (see Fig. 4.1

More information

θ x = f ( x,t) could be written as

θ x = f ( x,t) could be written as 9. Higher orer PDEs as systems of first-orer PDEs. Hyperbolic systems. For PDEs, as for ODEs, we may reuce the orer by efining new epenent variables. For example, in the case of the wave equation, (1)

More information

Optimal Cooperative Spectrum Sensing in Cognitive Sensor Networks

Optimal Cooperative Spectrum Sensing in Cognitive Sensor Networks Optimal Cooperative Spectrum Sensing in Cognitive Sensor Networks Hai Ngoc Pham, an Zhang, Paal E. Engelsta,,3, Tor Skeie,, Frank Eliassen, Department of Informatics, University of Oslo, Norway Simula

More information

How the potentials in different gauges yield the same retarded electric and magnetic fields

How the potentials in different gauges yield the same retarded electric and magnetic fields How the potentials in ifferent gauges yiel the same retare electric an magnetic fiels José A. Heras a Departamento e Física, E. S. F. M., Instituto Politécnico Nacional, México D. F. México an Department

More information

Lower Bounds for the Smoothed Number of Pareto optimal Solutions

Lower Bounds for the Smoothed Number of Pareto optimal Solutions Lower Bouns for the Smoothe Number of Pareto optimal Solutions Tobias Brunsch an Heiko Röglin Department of Computer Science, University of Bonn, Germany brunsch@cs.uni-bonn.e, heiko@roeglin.org Abstract.

More information

Throughput Optimal Control of Cooperative Relay Networks

Throughput Optimal Control of Cooperative Relay Networks hroughput Optimal Control of Cooperative Relay Networks Emun Yeh Dept. of Electrical Engineering Yale University New Haven, C 06520, USA E-mail: emun.yeh@yale.eu Ranall Berry Dept. of Electrical an Computer

More information

TRACKING CONTROL OF MULTIPLE MOBILE ROBOTS: A CASE STUDY OF INTER-ROBOT COLLISION-FREE PROBLEM

TRACKING CONTROL OF MULTIPLE MOBILE ROBOTS: A CASE STUDY OF INTER-ROBOT COLLISION-FREE PROBLEM 265 Asian Journal of Control, Vol. 4, No. 3, pp. 265-273, September 22 TRACKING CONTROL OF MULTIPLE MOBILE ROBOTS: A CASE STUDY OF INTER-ROBOT COLLISION-FREE PROBLEM Jurachart Jongusuk an Tsutomu Mita

More information

Delay Limited Capacity of Ad hoc Networks: Asymptotically Optimal Transmission and Relaying Strategy

Delay Limited Capacity of Ad hoc Networks: Asymptotically Optimal Transmission and Relaying Strategy Delay Limite Capacity of A hoc Networks: Asymptotically Optimal Transmission an Relaying Strategy Eugene Perevalov Lehigh University Bethlehem, PA 85 Email: eup2@lehigh.eu Rick Blum Lehigh University Bethlehem,

More information

A simple model for the small-strain behaviour of soils

A simple model for the small-strain behaviour of soils A simple moel for the small-strain behaviour of soils José Jorge Naer Department of Structural an Geotechnical ngineering, Polytechnic School, University of São Paulo 05508-900, São Paulo, Brazil, e-mail:

More information

Non-deterministic Social Laws

Non-deterministic Social Laws Non-eterministic Social Laws Michael H. Coen MIT Artificial Intelligence Lab 55 Technology Square Cambrige, MA 09 mhcoen@ai.mit.eu Abstract The paper generalizes the notion of a social law, the founation

More information

Lecture 6: Control of Three-Phase Inverters

Lecture 6: Control of Three-Phase Inverters Yoash Levron The Anrew an Erna Viterbi Faculty of Electrical Engineering, Technion Israel Institute of Technology, Haifa 323, Israel yoashl@ee.technion.ac.il Juri Belikov Department of Computer Systems,

More information

The Second Order Contribution to Wave Crest Amplitude Random Simulations and NewWave

The Second Order Contribution to Wave Crest Amplitude Random Simulations and NewWave The Secon Orer Contribution to Wave Crest Amplitue Ranom Simulations an NewWave Thomas A.A. Acock Department of Engineering Science University of Oxfor Parks Roa Oxfor Unite Kingom Scott Draper School

More information

II. First variation of functionals

II. First variation of functionals II. First variation of functionals The erivative of a function being zero is a necessary conition for the etremum of that function in orinary calculus. Let us now tackle the question of the equivalent

More information

On the Aloha throughput-fairness tradeoff

On the Aloha throughput-fairness tradeoff On the Aloha throughput-fairness traeoff 1 Nan Xie, Member, IEEE, an Steven Weber, Senior Member, IEEE Abstract arxiv:1605.01557v1 [cs.it] 5 May 2016 A well-known inner boun of the stability region of

More information

The Press-Schechter mass function

The Press-Schechter mass function The Press-Schechter mass function To state the obvious: It is important to relate our theories to what we can observe. We have looke at linear perturbation theory, an we have consiere a simple moel for

More information

Web-Based Technical Appendix: Multi-Product Firms and Trade Liberalization

Web-Based Technical Appendix: Multi-Product Firms and Trade Liberalization Web-Base Technical Appeni: Multi-Prouct Firms an Trae Liberalization Anrew B. Bernar Tuck School of Business at Dartmouth & NBER Stephen J. Reing LSE, Yale School of Management & CEPR Peter K. Schott Yale

More information

MATHEMATICS BONUS FILES for faculty and students

MATHEMATICS BONUS FILES for faculty and students MATHMATI BONU FIL for faculty an stuents http://www.onu.eu/~mcaragiu1/bonus_files.html RIVD: May 15, 9 PUBLIHD: May 5, 9 toffel 1 Maxwell s quations through the Major Vector Theorems Joshua toffel Department

More information

Section 7.1: Integration by Parts

Section 7.1: Integration by Parts Section 7.1: Integration by Parts 1. Introuction to Integration Techniques Unlike ifferentiation where there are a large number of rules which allow you (in principle) to ifferentiate any function, the

More information

12.11 Laplace s Equation in Cylindrical and

12.11 Laplace s Equation in Cylindrical and SEC. 2. Laplace s Equation in Cylinrical an Spherical Coorinates. Potential 593 2. Laplace s Equation in Cylinrical an Spherical Coorinates. Potential One of the most important PDEs in physics an engineering

More information

Reduction of chemical systems by delayed quasi-steady state assumptions

Reduction of chemical systems by delayed quasi-steady state assumptions Reuction of chemical systems by elaye quasi-steay state assumptions Tomáš Vejchoský a,b,, Raek Erban a, Philip K. Maini a a Wolfson Centre for Mathematical Biology, Mathematical Institute, University of

More information

Cascaded redundancy reduction

Cascaded redundancy reduction Network: Comput. Neural Syst. 9 (1998) 73 84. Printe in the UK PII: S0954-898X(98)88342-5 Cascae reunancy reuction Virginia R e Sa an Geoffrey E Hinton Department of Computer Science, University of Toronto,

More information

Optimization of Geometries by Energy Minimization

Optimization of Geometries by Energy Minimization Optimization of Geometries by Energy Minimization by Tracy P. Hamilton Department of Chemistry University of Alabama at Birmingham Birmingham, AL 3594-140 hamilton@uab.eu Copyright Tracy P. Hamilton, 1997.

More information

Estimation of the Maximum Domination Value in Multi-Dimensional Data Sets

Estimation of the Maximum Domination Value in Multi-Dimensional Data Sets Proceeings of the 4th East-European Conference on Avances in Databases an Information Systems ADBIS) 200 Estimation of the Maximum Domination Value in Multi-Dimensional Data Sets Eleftherios Tiakas, Apostolos.

More information

Design A Robust Power System Stabilizer on SMIB Using Lyapunov Theory

Design A Robust Power System Stabilizer on SMIB Using Lyapunov Theory Design A Robust Power System Stabilizer on SMIB Using Lyapunov Theory Yin Li, Stuent Member, IEEE, Lingling Fan, Senior Member, IEEE Abstract This paper proposes a robust power system stabilizer (PSS)

More information

6 General properties of an autonomous system of two first order ODE

6 General properties of an autonomous system of two first order ODE 6 General properties of an autonomous system of two first orer ODE Here we embark on stuying the autonomous system of two first orer ifferential equations of the form ẋ 1 = f 1 (, x 2 ), ẋ 2 = f 2 (, x

More information

Thermal runaway during blocking

Thermal runaway during blocking Thermal runaway uring blocking CES_stable CES ICES_stable ICES k 6.5 ma 13 6. 12 5.5 11 5. 1 4.5 9 4. 8 3.5 7 3. 6 2.5 5 2. 4 1.5 3 1. 2.5 1. 6 12 18 24 3 36 s Thermal runaway uring blocking Application

More information

. Using a multinomial model gives us the following equation for P d. , with respect to same length term sequences.

. Using a multinomial model gives us the following equation for P d. , with respect to same length term sequences. S 63 Lecture 8 2/2/26 Lecturer Lillian Lee Scribes Peter Babinski, Davi Lin Basic Language Moeling Approach I. Special ase of LM-base Approach a. Recap of Formulas an Terms b. Fixing θ? c. About that Multinomial

More information

Vectors in two dimensions

Vectors in two dimensions Vectors in two imensions Until now, we have been working in one imension only The main reason for this is to become familiar with the main physical ieas like Newton s secon law, without the aitional complication

More information

Calculus of Variations

Calculus of Variations Calculus of Variations Lagrangian formalism is the main tool of theoretical classical mechanics. Calculus of Variations is a part of Mathematics which Lagrangian formalism is base on. In this section,

More information

Generalizing Kronecker Graphs in order to Model Searchable Networks

Generalizing Kronecker Graphs in order to Model Searchable Networks Generalizing Kronecker Graphs in orer to Moel Searchable Networks Elizabeth Boine, Babak Hassibi, Aam Wierman California Institute of Technology Pasaena, CA 925 Email: {eaboine, hassibi, aamw}@caltecheu

More information

3.7 Implicit Differentiation -- A Brief Introduction -- Student Notes

3.7 Implicit Differentiation -- A Brief Introduction -- Student Notes Fin these erivatives of these functions: y.7 Implicit Differentiation -- A Brief Introuction -- Stuent Notes tan y sin tan = sin y e = e = Write the inverses of these functions: y tan y sin How woul we

More information

Research Article When Inflation Causes No Increase in Claim Amounts

Research Article When Inflation Causes No Increase in Claim Amounts Probability an Statistics Volume 2009, Article ID 943926, 10 pages oi:10.1155/2009/943926 Research Article When Inflation Causes No Increase in Claim Amounts Vytaras Brazauskas, 1 Bruce L. Jones, 2 an

More information

Harmonic Modelling of Thyristor Bridges using a Simplified Time Domain Method

Harmonic Modelling of Thyristor Bridges using a Simplified Time Domain Method 1 Harmonic Moelling of Thyristor Briges using a Simplifie Time Domain Metho P. W. Lehn, Senior Member IEEE, an G. Ebner Abstract The paper presents time omain methos for harmonic analysis of a 6-pulse

More information

Angles-Only Orbit Determination Copyright 2006 Michel Santos Page 1

Angles-Only Orbit Determination Copyright 2006 Michel Santos Page 1 Angles-Only Orbit Determination Copyright 6 Michel Santos Page 1 Abstract This ocument presents a re-erivation of the Gauss an Laplace Angles-Only Methos for Initial Orbit Determination. It keeps close

More information