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 Delaare University of Southern California Abstract While many moels of the TCP s ynamics have been evelope, fe focus on the effects of timeout an high loss probability. Active queue management (AQM is an important application of these ynamic moels. Hoever, recent ork has shon 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 ino. It is shon that, hile the ynamics of the mean value of the congestion ino 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 illy oscillate. Finally, this paper inclues results from several million simulations proviing a etaile vie of the ynamics of timeout. I. INTRODUCTION s of TCP are extensively use throughout netorking research. For example, TCP moels are use in fiels such as AQM [], [], [3], [4], esign of TCP frienly transport protocols [5], planning an provisioning netorks [6], preicting file transfer time [7], [8], etc. In this ay, TCP moels are founation for research. The accuracy of these moels is critical to the accuracy an relevance of the ork that rests upon them. s for TCP have evolve over time. Initially, a simple relationship beteen the sening rate an the loss probability as evelope. The ell-knon /sqrt(p formula, or more specifically MSS 3// ( RT T p, as foun an verifie [9], []. While this formula is useful in many settings, it as foun to not be accurate in situations here 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 shon that timeout, as oppose to triple uplicate acknolegement, is a significant ay that TCP etects packet loss []. In orer to inclue the impact of timeout, the simple square root of p formula as extene to inclue the effect of timeout []. This formula as shon to be more accurate for large loss probabilities here timeout is significant. Sikar et al. [] evelope a moel for TCP throughput that inclues time-out. While this moel is ifferent from [], they report that it gives similar numerical results. Hoever, e have performe extensive simulations, over 5 million in total, an have foun that hen applie to the ns- implementation of TCP-SACK an hen loss is probabilistic, the moel given in [] results in errors as large as a factor of to. Consiering the importance of the relationship beteen 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 as to make them ynamic. That is, the initial moels assume that the loss probability as fixe. In the case that the loss probability varie, ynamic moels ere require. The first moel in this irection as given in [4], an several closely relate moels have folloe [], [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 as first introuce, there have been a large number of researchers active in this area. Recently, several groups have investigate hat, if any, benefit AQM has on the performance of the netork (see [3] an references therein. A recent paper along these lines shoe that the AQM, be it rop tail or any more elaborate version, has little impact on the performance of the netork for utilization belo 9%. Hoever, for utilization above 9%, some AQM techniques ere shon to improve the performance [4]. This ork is significant as it inicates the environment here AQM might have a positive impact, specifically, in netorks ith 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 here AQM might have the largest impact on performance, the analysis an esign of AQM is least accurate. Because of this inaccuracy, the possible impact of AQM on netork performance at high utilizations has never been carefully explore. One reason that AQM esign has not consiere high utilization netorks is that ynamics of TCP in such an environment are not knon. By ynamics, e 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 illy varying bit-rates that can take several secons to ecay. For example, the loer right of Figure 4 shos such oscillations. This ynamics is far more complicate than the Globecom /4/$. 4 IEEE

2 ynamics of the mean congestion ino use by most AQM esign. As a result, this ork calls into question the proof of stability of AQM at high loss probabilities. As a by-prouct of this ork, the variance an the istribution of the congestion ino is also foun. Specifically, it is shon that the ino size is accurately moele ith a negative binomial istribution. An finally the huge number of simulations use to valiate this moel provie a etaile vie of the behavior of the ynamics of TCP s timeout mechanism. The paper procees as follos: The next section iscusses the simulations use for verification. In Section III the istribution of the congestion ino is foun. Next, in Section IV, a moel of the probability of a flo 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 ith 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 netorks the sening rate of a TCP connection is complicate by interacting flos, 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 il, it is necessary to thoroughly unerstan TCP in a controlle environment. Perhaps the most significant ifference beteen these simulations an the environment foun in the il is that in these simulations packet losses are entirely ranom an never ue to queue overflo. In the case of AQM, such simulations are appropriate as the goal of AQM is to rop packets in a controlle fashion, not hen the queue fills. The simulations presente in this paper are for a single flo over a single bottleneck topology here rops at the bottleneck are ranom. The ns- error moel as 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 acknolegement an a maximum receiver ino ere 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 ere performe on the University of Southern California s Linux cluster [5] an the large simulation results ere transferre via the Internet to the University of Delaare for postprocessing an storage. There is little oubt that ithout such high-performance computing, this ork oul not be possible. III. THE DISTRIBUTION OF THE CONGESTION WINDOW This section presents a simple moel for the istribution of the congestion ino. While the istribution may be of interest in its on right, e are particularly intereste in this istribution since it is use to etermine the probability of a flo entering timeout. We justify the moel briefly ith 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 ino size varies accoring to W t = RT T t W tn t ( here N is a Cox process that counts the number of packet losses. This moel is for the evolution of the congestion ino hen the flo is not in timeout. Thus, all probabilities shoul be conitione on the flo not in timeout. We enote this conition as TO. This moel as further investigate in [7] here the partial ifferential equation of the ino size as foun. Specifically, if p (, t TO is the probability ensity of the congestion ino taking the value at time t, then p satisfies p (, t TO t = ( p(, t TO RT T + δ (t RT T ( 4p (, t TO p (, t TO, here δ is the loss probability. If δ is constant, then in steay state, i.e., solves p ( TO p(,t TO t ( =, the istribution of the ino = δ ( 4p ( TO p ( TO. (3 From (3 it is straightforar to sho that E ( m TO = c m δ, m/ for some constants c m.ifm =, then this is the ell-knon square root of p formula here c has been foun to be beteen. an.3. Consiering m =, e fin that the variance is of the form Var ( TO := E ( TO E ( TO = γ δ, (4 here γ := c c. This relationship is borne out in simulations that inicates that γ.3. A close form solution has not been foun for (, hereas a complicate close form solution has been foun for (3. Hoever, simple yet accurate approximate solution to (3 is the negative binomial istribution, i.e., p ( = Γ(N + Γ(N(! ( qn q. Since the mean of this negative binomial ranom variable is Nq/( q+ an the variance is Nq/( q, q = E ( TO Var ( TO an N = q q ( E ( TO (5 an from mean an variance of the given above e get q = δ ( ( c ( q c + an N =. (6 γ δ q δ Globecom /4/$. 4 IEEE

3 In orer to etermine the istribution of the ino size, the parameters q an N must be etermine. This can be one by plugging the mean an variance into (5. Hoever, if only the loss probability is knon, then (6 must be use, in hich case the parameters c an γ must be etermine. We have foun that if the objective is to use the istribution to etermine the probability of timeout, the c = 3/ an γ =.3 is sufficient. On the other han, if one is intereste in the istribution of the congestion ino (the objective of this section, then the selection of these parameters, especially c, requires some more analysis. Hoever, to keep this paper focuse, this issue is not iscusse. Once the parameters are etermine, the istribution of the congestion ino can be estimate. Figure shos the istribution of the congestion ino given the flo is not in the timeout state, i.e., p ( TO. These figures shos the three types of curves, the observe istribution (histogram, the istribution ith c = 3/, an the istribution here the parameters are etermine by plugging in the observe mean an variance into (5. This last istribution coincies ith the istribution given by (6 but ith c =.7 for δ = 4 an c =.4 for δ =.5. Forδ =., thevalueofc plays less of a role, c =.4, c = 3/, an c =.7, or by using the observe mean an variance, all yiel nearly the same istribution. In all cases, γ =.3. IV. A MODEL OF TIMEOUT In this section, e present a moel of the probability of a flo being in the timeout state. We say that a flo is in the timeout state if the next packet ill 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 netorks ith ientical netork characteristics (e.g., link spees, propagation elays, etc.. Our goal is to etermine the fraction of these flos that are in timeout. A flo can enter the timeout state in the folloing three ays: I If a flo experiences so many losses that triple uplicate acknolegements are not receive, i.e., if at least max (, losses occur in one ino. II In the case of the ns- implementation of TCP-SACK, if more than / packets are roppe. III If a retransmitte packet is roppe. Consiering I an II, e see that a flo ill timeout if a singleropisfolloebyatleastmax (min ( /, 3, rops out of the next packets. Given that rops occur at a rate Rδ, rops that lea to timeout occur at a rate ρ (, δ = R δ k=max(min( /, 3, ( k δ k ( δ k. This ay of entering time-out is iscusse in []. RFC-357 eliminates this ay to enter time-out. It is straightforar to ajust the evelopment belo to reflect RFC-357. Let λ (δ be the rate that a flo moves from congestion avoiance to timeout ue to I an II. This rate is foun taking the expecte value of ρ (, δ, i.e., λ (δ =E ( ρ (, δ TO = ρ (, δ p δ (, here p δ ( is given in Section III, hoever, here e have shon the epenence on δ. This epenence is mae explicit by (6. Note that the approximation E ( ρ (, δ TO ρ ( E ( TO,δ ill result in large errors. Consier the situation here δ =.9 an, using c = 3/, the mean value of the congestion ino is 4. In this case, E ( ρ (, δ TO =.3 R hile ρ ( E ( TO,δ =.38 R, a ifference aroun a factor of three. Next e etermine λ (i.e., III, the rate that a flo moves from congestion avoiance to timeout because a retransmitte packet is roppe. The rate that a packet is first roppe is δ R. The rop ill lea to a retransmission only if triple uplicate acknolegements are receive, i.e., if less than max (min ( /, 3, packets are roppe out of the next packets. In this case, the retransmission is roppe an timeout is entere ith probability δ. Thus ρ (, δ :=δ R δ ( ( δ k ( δ k. k k=max(min( /, 3, The average rate that a flo enters timeout ue to roppe retransmissions is λ (δ :=E (ρ (, δ Finally, the rate that a TCP flo moves into timeout is λ +λ. Up to this point e consiere the rate that a single flo moves into timeout. No e ill etermine the fraction of flos entering timeout. Let N I (t, enote the rate that flos enter timeout at time t, N is the total number of flos. That is, in small time interval t, the probability that some flo enters timeout is N I (t t. LetRT O be the time in hich flos remain in the timeout state upon the first timeout. If a flo exits timeout, but immeiately experiences another rop, then this flo again enters timeout, but for RT O secons. We enote the rate that flos enter timeout for this secon time ith N I (t. The fraction of flos in timeout are t RT O I (τ τ + t RT O I (τ τ. The rate that flos enter timeout is the prouct of the number of flos not in timeout an the rate that such flos enter timeout, i.e., I (t =(λ (δ (t + λ (δ (t (7 ( ( I (τ τ + I (t τ. t RT O t RT O We on t consier the case here a flo enters time-out a thir time for a time interval 4 RT O. While it is straightforar to exten the calculations, it oes not appear to improve the accuracy of the moel. Globecom /4/$. 4 IEEE

4 probability. observe. observe.9.9 (3/.8 /δ (3/.8 /δ.7.7/ δ or.7 observe mean.6 observe mean δ= -4.3 δ= probability probability. observe.8 (3/.6 /δ.4.4/ δ or. observe mean..8.6 δ= Figure. Distribution of the congestion ino for loss probability of δ = 4 (left, δ =. (mile, an δ =.5 (right. While the rate that flos exit timeout only to reenter timeout is I (t =δi (t RT O. (8 In steay state, I (t an I (t are constant. So (7 an (8 reuce to I =( I RT O I RT O(λ (δ+λ (δ I = δi. This can be solve λ (δ+λ (δ I = +RT O (λ (δ+λ (δ ( + δ δ (λ (δ+λ (δ I = +RT O (λ (δ+λ (δ ( + δ Thus, the fraction of flos in timeout is P (TO= t RT O I (τ τ + t RT O I (τ τ (9 = ( + δ RT O (λ (δ+λ (δ +(+δ RT O (λ (δ+λ (δ. If the loss probability varies, then (7, (8, an the relationship P (in timeout at time t = t RT O I (τ τ + t RT O I (τ τ must be use to etermine the probability that a flo is in timeout. To verify (9, a large number of ns- simulations ere performe so that the probability of timeout coul be etermine even for small loss probabilities. Figure shos the observe relationship of P (in timeout an the loss probability as ell as the estimate of this probability given by (9. These calculations use c = 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 flo spens in timeout, i.e., P (TO. Specifically, [] fins P (timeout Q (W (δ,δ G (δ T o δ R 8 (W (δ++q(w (δ,δ G (δ T o δ, ( 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 flo being in timeout. The ashe green line is the moel s estimate of this probability. Note that these to lines nearly coincie. The thir curve is the relationship erive from []. here (+( δ 3( ( δ 3 Q (, δ = min, ( ( δ / ( ( δ 3 8( δ W (δ =+ + 3δ G (δ =+δ +δ +4δ 3 +8δ 4 +6δ 5 +3δ 6. Figure inclues this estimate. We see that this moel, hile 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%, hile ( gives an estimate of 43%. Hoever, this comparison must be qualifie. This paper focuses on TCP- SACK hile the focus of [] is on TCP-RENO. Furthermore, this paper examines ho TCP behaves uner ranom losses such as arises in AQM. In [], the queueing iscipline is rop tail. Furthermore, the loss probability in [] as not exactly the probability of a packet being roppe, but more along the lines of the probability of a rop event occurring here 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. Hoever, this moel is useful as a benchmark as it as not knon ho accurate this moel oul be for TCP-SACK. From the probability of timeout, the average through-put Globecom /4/$. 4 IEEE

5 Relative Error in Throughput (fraction Pahye et al. RFC δ(loss probability Figure 3. Relative Error of Preicte Throughput. Note that an error fraction of.5 implies a 5% error. can be foun via T = c R ( P (TO MSS. δ Figure 3 shos the relative error of this moel hen compare to the observe throughput. Figure 3 also shos the relative error of the moel in [] as ell as the moel suggeste in RFC-3448 for a TCP frienly sening rate. The figure shos only the relative error, an examination of the error shos that the moel of [] an RFC-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 shos that for even moerately large loss probabilities, the probability of being in timeout is quite high. Specifically, for δ>%, timeout plays an important role in the sening rate of TCP. For δ>5%, a flo ill spen 6% of the time in timeout. Thus, timeout ominates the sening rate of TCP. There are many implications of this. Recently 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 ith these moifications ill be more aggressive than implementations ithout these moifications. There has been extensive ork 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 %. Hoever, in the area of AQM, it is hope that these moels are applicable to higher loss probabilities. Next e evelop a moel for the ynamic behavior of TCP that inclues timeout. V. DYNAMICS OF TCP Here e 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 e present a moel for the variance. Once the variance is knon, the istribution is foun as in Section III. With the istribution, e can fin the probability of being in timeout as as one in the previous section. Of course, the paper to this point has only consiere a constant loss probability. Here e allo the loss probability to vary ith time. Thus, this moel is appropriate for AQM esign. 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 sho the observe an the ashe re curves shos the moel. In the loer right, the otte green inicates the bit-rate if only the mean size of the congestion ino is use (i.e., the ashe re curve in the upper left. Note the units of the bit-rate is in bps. Let (t =E ( (t TO an (t =E ( (t TO. An approximation of the ynamics of the mean value of the congestion ino is often given by t (t = R c δ (t R (t R (t. ( R This is often approximate as t (t = R c R δ (t R (t. ( Interestingly, this formula is often erive from (. Hoever, if proper stochastic calculus is applie, the correct ynamics for the mean are t (t = R δ (t R E ( (t R (t, R hich can be approximate as t (t = R R δ (t R (t. (3 Note that in (3 the secon moment is utilize, hile ( have multiplie by the first moment square. Fortunately, c e have foun that ( an ( give goo approximation to the ynamics of the mean. It seems that the mean varies sloly enough that this ifference beteen (t an (t is not c significant. From ( it is possible to etermine the ynamics of the secon moment of the congestion ino, t (t = R (t 3 4 R δ (t R E ( 3 (t. Globecom /4/$. 4 IEEE

6 bit-rate Prob(timeout at time variance(w(t E(W(t 8 7 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 shon in Figure 4. Hoever, here to sets of ynamics are shon, 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, hoever, here the loss probability jumps beteen % an 5%. Using the same iea as above, e approximate E 3 (t 3/ β 8 an arrive at 3 (c.33/ (t (t = (t (4 t R 3/ 38 β δ (t R (t. 3/ 4 3 (c.3 R 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 ino can be foun as iscusse in Section III. With this approximation of the Globecom 4 istribution, the probability of being in timeout can be foun t as in Section IV, i.e., P (T O at time t = t RT O I (τ τ + t I (τ τ, here I (τ is given by (7 an I (τ is t RT O given by (8. Figures 4-6 provie some valiation of this ynamic moel of TCP. Through extensive simulation, E (t T O, E (t T O, an P (T O at time t ere foun uner ifferent time variations of the loss probability. Specifically, in each simulation, the loss probability took to values δ an δ. The simulation began an ran for secons ith loss probability set to δ. In these simulations, secons as enough to ensure that the system as in steay state (i.e., E (t T O as constant. Then, secons after the beginning of the simulation, the loss probability as sitche to δ. In the figures, this moment is labele as t =. This simulation as carrie out, times. With these simulations it is possible to etermine the value of the statistics E (t T O, E (t T O, an P (T O at time t at each time point. For example, to etermine P (T O at time 5 ms, e foun the fraction of the, flos that ere in the timeout state at time t = 5 ms. This large number of simulations allos for high confience in the observations an allos the etails of the ynamics to be observe. Figures 4-6 inicate that the ynamic moel is quite accurate. They also sho, as expecte, that using the ynamics of the mean congestion ino 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 ino perform rather poorly for very large loss probability. Hoever, in this case, the behavior of TCP is completely ominate by the timeout, hich is ell moele. As a result, the bit-rate is accurately moele. Perhaps the most striking aspect of Figures 4-6 is the ynamics of the bit-rate at large loss probabilities. Most significantly, e see that the bit-rate experiences il oscillations that can take several secons to ecay. In all cases, the moel for the mean congestion ino ( oes not prouce any oscillations. To the best of the authors knolege, 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 ill impact the stability results of previous AQM research. Figure 6 shos some unexpecte behavior of the mean an the variance of the congestion ino. Specifically, e see that the mean has to istinct phases of groth ith a transition at t = ms hile the variance shos nonmonotonic groth. The ODEs ( an (4 o not moel this behavior. This behavior is ue to the feeback of the fraction of flos in timeout to the mean an the variance of the congestion ino. Notice that the non-monotonic groth of occurs at RTO, just as flo exit timeout. Hoever, since the probability of a flo being in the timeout state is so high, a slight error in the mean congestion ino has a minor impact on the average bit-rate /4/$. 4 IEEE

7 VI. CONCLUSIONS We have evelope a ne moel for the sening rate of TCP that inclues the effect of timeout. Previous moels have inclue the effect of timeout, but e 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 as shon to agree ith simulations quite ell. Furthermore, it as shon that timeout can cause the bit-rate to oscillate substantially for several secons after a change in the loss probability. It as also shon 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 ork ill examine the possibility of esigning an AQM that utilizes these moels for TCP. Perhaps, more substantial performance gains ill be achieve for large loss probability. REFERENCES [] S. Athuraliya, V. H. Li, S. H. Lo, an Q. Yin, REM: Active queue management, IEEE Netork, vol. 5, pp ,. [] C. V. Hollot, V. Misra, D. F. Tosley, an W. Gong, On esigning improve controllers for AQM routers supporting TCP flos, in IN- FOCOM,. [3] S. Kunniyur an R. 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. Tosley, Stochastic ifferential equation moeling an analysis of tcp-inosize behavior, in Proceeings of PERFORMANCE99, 999. [5] M. Hanley, S. Floy, J. Pahye, an J. Wimer, RFC TCP frienly rate control (TFRC: Protocol specification, 3. [6] J. Beckers, I. Henraan, R. Kooij, an R. V. er Mei, Generalize processor sharing performance moels for internet access lines, in 9th IFIP CONFERENCE on PERFORMANCE MODELLING AND EVALUATION OF ATM an IP NETWORKS, (Buapest,. [7] S. Bohacek, A stochastic moel of tcp an fair vieo transmission, in Infocom, 3. [8] N. Carell, S. Savage, an T. Anerson, ing TCP latency, in INFOCOM (3, pp ,. [9] T. Ott, J. H. B. Kemperman, an M. Mathis, Wino size behavior in TCP/IP ith constant loss probability, in Proc. Of the DIMACS Workshop on Performance of Realtime 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. Revie, vol. 7, no. 3, 997. [] J. Pahye, V. Firoiu, D. Tosley, 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 Netorking, vol., no. 6, pp , 3. [3] S. Bohacek, K. Shah, G. Arce, an M. Davis, A revie 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 eb performance, in SIGCOMM, 3. [5] Globecom /4/$. 4 IEEE

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 Delaware 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

On Buffer Limited Congestion Window Dynamics and Packet Loss

On Buffer Limited Congestion Window Dynamics and Packet Loss On Buffer Limited Congestion Windo Dynamics and Packet Loss A. Fekete, G. Vattay Communication Netorks Laboratory, Eötvös University Pázmány P. sétány /A, Budapest, Hungary 7 Abstract The central result

More information

Pattern Propagation Speed in Synfire Chains with Excitatory- Inhibitory Couplings

Pattern Propagation Speed in Synfire Chains with Excitatory- Inhibitory Couplings Pattern Propagation Spee in Synfire Chains ith Excitatory- Inhibitory Couplings Baktash Babai School of Intelligent Systems, Institutes for Stuies in Theoretical Physics & Mathematics baktash@ipm.ir Abstract

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

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

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

'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

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

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

The Ritz Ballistic Theory & Adjusting the Speed of Light to c near the Earth and Other Celestial Bodies

The Ritz Ballistic Theory & Adjusting the Speed of Light to c near the Earth and Other Celestial Bodies College Park, MD 11 PROCDINGS of the NPA 1 The Ritz Ballistic Theory & Ajusting the Spee of Light to c near the arth an Other Celestial Boies Nina Sotina PhD in Physics, Mosco State University 448 Neptune

More information

Proof of Proposition 1

Proof of Proposition 1 A Proofs of Propositions,2,. Before e look at the MMD calculations in various cases, e prove the folloing useful characterization of MMD for translation invariant kernels like the Gaussian an Laplace kernels.

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

Can Punctured Rate-1/2 Turbo Codes Achieve a Lower Error Floor than their Rate-1/3 Parent Codes?

Can Punctured Rate-1/2 Turbo Codes Achieve a Lower Error Floor than their Rate-1/3 Parent Codes? Can Puncture Rate-1/2 Turbo Coes Achieve a Loer Error Floor than their Rate-1/3 Parent Coes? Ioannis Chatzigeorgiou, Miguel R. D. Rorigues, Ian J. Wassell Digital Technology Group, Computer Laboratory

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Dissipativity for dual linear differential inclusions through conjugate storage functions

Dissipativity for dual linear differential inclusions through conjugate storage functions 43r IEEE Conference on Decision an Control December 4-7, 24 Atlantis, Paraise Islan, Bahamas WeC24 Dissipativity for ual linear ifferential inclusions through conjugate storage functions Rafal Goebel,

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

θ 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

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

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

Conservation Laws. Chapter Conservation of Energy

Conservation Laws. Chapter Conservation of Energy 20 Chapter 3 Conservation Laws In orer to check the physical consistency of the above set of equations governing Maxwell-Lorentz electroynamics [(2.10) an (2.12) or (1.65) an (1.68)], we examine the action

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

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

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

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

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

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

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

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

Heat Transfer from Arbitrary Nonisothermal Surfaces in a Laminar Flow

Heat Transfer from Arbitrary Nonisothermal Surfaces in a Laminar Flow 3 Heat Transfer from Arbitrary Nonisothermal Surfaces in a Laminar Flo In essence, the conjugate heat transfer problem consiers the thermal interaction beteen a boy an a flui floing over or insie it. As

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

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

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

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

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

Resilient Modulus Prediction Model for Fine-Grained Soils in Ohio: Preliminary Study

Resilient Modulus Prediction Model for Fine-Grained Soils in Ohio: Preliminary Study Resilient Moulus Preiction Moel for Fine-Graine Soils in Ohio: Preliminary Stuy by Teruhisa Masaa: Associate Professor, Civil Engineering Department Ohio University, Athens, OH 4570 Tel: (740) 59-474 Fax:

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

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

Modelling TCP with a Discrete Time Markov Chain

Modelling TCP with a Discrete Time Markov Chain Modelling TCP with a Discrete Time Markov Chain José L Gil Motorola josegil@motorola.com ABSTRACT TCP is the most widely used transport protocol in the Internet. The end-to-end performance of most Internet

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

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

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

Using Partial Differential Equations to Model TCP Mice and Elephants in Large IP Networks

Using Partial Differential Equations to Model TCP Mice and Elephants in Large IP Networks Using Partial Differential Equations to Model TCP Mice and Elephants in Large IP Netorks M.Ajmone Marsan, M.Garetto, P.Giaccone, E.Leonardi, E.Schiattarella, A.Tarello Dipartimento di Elettronica, Politecnico

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

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

Adaptive Gain-Scheduled H Control of Linear Parameter-Varying Systems with Time-Delayed Elements

Adaptive Gain-Scheduled H Control of Linear Parameter-Varying Systems with Time-Delayed Elements Aaptive Gain-Scheule H Control of Linear Parameter-Varying Systems with ime-delaye Elements Yoshihiko Miyasato he Institute of Statistical Mathematics 4-6-7 Minami-Azabu, Minato-ku, okyo 6-8569, Japan

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

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

19 Eigenvalues, Eigenvectors, Ordinary Differential Equations, and Control

19 Eigenvalues, Eigenvectors, Ordinary Differential Equations, and Control 19 Eigenvalues, Eigenvectors, Orinary Differential Equations, an Control This section introuces eigenvalues an eigenvectors of a matrix, an iscusses the role of the eigenvalues in etermining the behavior

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

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

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

Yang Zhang, Xinmin Wang School of Automation, Northwestern Polytechnical University, Xi an , China

Yang Zhang, Xinmin Wang School of Automation, Northwestern Polytechnical University, Xi an , China oi:1.1311/1.39.6.48 Robust Aaptive Beamforming ith Null Broaening Yang Zhang, inmin ang School of Automation, Northestern Polytechnical University, i an 717, China Abstract Aaptive beamformers are sensitive

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

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

. 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

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

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

Spring 2016 Network Science

Spring 2016 Network Science Spring 206 Network Science Sample Problems for Quiz I Problem [The Application of An one-imensional Poisson Process] Suppose that the number of typographical errors in a new text is Poisson istribute with

More information

TMA 4195 Matematisk modellering Exam Tuesday December 16, :00 13:00 Problems and solution with additional comments

TMA 4195 Matematisk modellering Exam Tuesday December 16, :00 13:00 Problems and solution with additional comments Problem F U L W D g m 3 2 s 2 0 0 0 0 2 kg 0 0 0 0 0 0 Table : Dimension matrix TMA 495 Matematisk moellering Exam Tuesay December 6, 2008 09:00 3:00 Problems an solution with aitional comments The necessary

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

Two- and Three-Dimensional Validation of Icing Model

Two- and Three-Dimensional Validation of Icing Model APCOM & ISCM 11-14 th December, 13, Singapore To- an Three-Dimensional Valiation of Icing Moel *Ryosuke Hayashi¹ an Makoto Yamamoto 1 Grauate School of Mechanical Engineering, Tokyo University of Science

More information

4. Important theorems in quantum mechanics

4. Important theorems in quantum mechanics TFY4215 Kjemisk fysikk og kvantemekanikk - Tillegg 4 1 TILLEGG 4 4. Important theorems in quantum mechanics Before attacking three-imensional potentials in the next chapter, we shall in chapter 4 of this

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

Semiclassical analysis of long-wavelength multiphoton processes: The Rydberg atom

Semiclassical analysis of long-wavelength multiphoton processes: The Rydberg atom PHYSICAL REVIEW A 69, 063409 (2004) Semiclassical analysis of long-wavelength multiphoton processes: The Ryberg atom Luz V. Vela-Arevalo* an Ronal F. Fox Center for Nonlinear Sciences an School of Physics,

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

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

BEYOND THE CONSTRUCTION OF OPTIMAL SWITCHING SURFACES FOR AUTONOMOUS HYBRID SYSTEMS. Mauro Boccadoro Magnus Egerstedt Paolo Valigi Yorai Wardi

BEYOND THE CONSTRUCTION OF OPTIMAL SWITCHING SURFACES FOR AUTONOMOUS HYBRID SYSTEMS. Mauro Boccadoro Magnus Egerstedt Paolo Valigi Yorai Wardi BEYOND THE CONSTRUCTION OF OPTIMAL SWITCHING SURFACES FOR AUTONOMOUS HYBRID SYSTEMS Mauro Boccaoro Magnus Egerstet Paolo Valigi Yorai Wari {boccaoro,valigi}@iei.unipg.it Dipartimento i Ingegneria Elettronica

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