Introduction In this paper, we consider the probem of optimay scheduing a set of jobs on a singe machine which is susceptibe to occasiona breakdowns.

Size: px
Start display at page:

Download "Introduction In this paper, we consider the probem of optimay scheduing a set of jobs on a singe machine which is susceptibe to occasiona breakdowns."

Transcription

1 Stochastic Scheduing on a Repairabe Machine with Erang Uptime Distribution Wei Li, Chinese Academy of Sciences W. John Braun y, University of Winnipeg Yiqiang Q. Zhao y, University of Winnipeg Abstract A set of jobs is to be processed on a machine which is subject to breakdown and repair. When the processing of a job is interrupted by a machine breakdown, the processing ater resumes at the point at which the breakdown occurred. We assume that the machine uptime is Erang distributed and that processing and repair times foow genera distributions. Simpe permutation poicies on both machine parameters and the processing distributions are given which minimize the weighted number of tardy jobs, weighted ow times, and the weighted sum of the job deays. Key words: Stochastic scheduing, breakdowns, Erang uptimes, Nonpreemptive resume mode. Institute of Appied Mathematics, Chinese Academy of Sciences, Beijing, 8, China and aso with the Department of Mathematics and Statistics, University of Winnipeg, Winnipeg, Manitoba, Canada, R3B E9. y Department of Mathematics and Statistics, University of Winnipeg, Winnipeg, Manitoba, Canada, R3B E9.

2 Introduction In this paper, we consider the probem of optimay scheduing a set of jobs on a singe machine which is susceptibe to occasiona breakdowns. When a breakdown occurs, the job currenty being processed is hated unti the machine is repaired; then the job is resumed at the point at which it was interrupted. Scheduing probems of this sort have been studied thoroughy in the case of reiabe machines and in cases where job processing times and due dates are known in advance. However, such assumptions are not very reaistic; machines are often subject to unpredictabe breakdowns and possiby engthy repairs. Furthermore, processing times and deadines may be more adequatey modeed as random variabes. Singe-machine stochastic scheduing modes incorporating breakdowns and repairs were rst studied by Gazebrook [4]. Subsequenty, many authors have investigated stochastic probems on unreiabe or repairabe machines (see references []-[7], []-[5]). However, amost a of the resuts are based on one or more strong assumptions. Aahverdi and Mittentha [3, 4, 5], and Birge et a. [6] obtained some powerfu resuts in the case of deterministic processing times. In case of exponentia uptimes, it is aways possibe to transform a scheduing probem on an unreiabe machine to a probem on a reiabe machine by extending the processing time of a particuar job to the tota time that the job occupies the machine. For some exampes, see Gazebrook [4], Pinedo and Rammouz [3], Aahverdi and Mittentha [3], and Li and Cao [8]. Under the assumption of exponentia processing times, some simpe poicies have been obtained (see, for exampe, Pinedo and Rammouz [3], Chang et a. [7], Du and Pinedo [], and Li and Cao [7]), but these are usuay based on some stronger conditions, and the poicies obtained are not reated to the machine parameters. Few resuts have been obtained when specic distributiona assumptions have not been made about the machine uptime. In stochastic scheduing probems where specic distributiona assumptions have not been made about processing requirements, or machine uptimes and downtimes, optima processing strategies coud depend upon both the eapsed time since the ast repair period and the eapsed processing time of the current job. Therefore, it woud seem very dicut to nd optima index poicies which depend on uptime and downtime distributions as we as the processing time distribution. In view of this, we investigate a singe machine scheduing mode in which the machine is subject to breakdowns with Erang uptimes and genera repair times. Each job's processing time is aso a random variabe with a genera distribution. This mode woud seem to be as genera as can be

3 reaisticay considered, if one seeks a schedue which depends on the machine parameters as we as the processing distribution. Li and Gazebrook [9] have considered the competey genera case, but their resuts are not reated to the machine parameters. In order to proceed, we restrict our attention to the cass of simpe recourse strategies. Under such a strategy, the sequence of job competion times in a schedue is xed, but the competion times may be postponed as a resut of the machine breakdowns (see Aahverdi and Mittentha [4, 5] or Frenk []). Further assumptions and notation are as foows: A set of jobs J = f; ; ; Ng (N < ) is to be processed on a singe machine which is subject to breakdown and repair. Job j has weight h j, which is basicay a priority factor, denoting the importance of job j reative to the other jobs in the system. For exampe, this weight may represent the actua cost of keeping the job in the system. This cost coud be a hoding or inventory cost; it coud aso be the amount of vaue aready added to the job. A jobs have a common due date d which is a random variabe with an exponentia distribution and rate r; this index represents the committed shipping or competion date (the date the jobs are promised to the customer). The competion of a job after its due date is aowed, but a penaty is incurred. The jth job has a random processing requirement j with distribution function F j (), having hazard rate function j () and nite expectation = j. The machine is avaiabe for processing from time unti the rst breakdown occurs at time U. The machine then takes time D to be repaired, whie no processing takes pace. The repair having been competed, the machine is again avaiabe for processing from time U + D unti time U + D + U, and so on. The machine up times U ; U ; ; are independent random variabes with an Erang distribution function U() having m phases and where each phase is exponentiay distributed with parameter. That is, U(t) =? e?t m? k= (t) k : k! The repair times D ; D ; are aso independent and identicay distributed but with genera distribution function D(), having hazard rate () and nite expectation =. A of the above random variabes are assumed to be mutuay independent. Under the above assumptions, we sha consider the so-caed non-preemptive resume mode; that is, if a job is being processed when a breakdown occurs, the processing time is cumuative after the machine is repaired. Our goa is to nd schedues, depending on both machine parameters and the processing distributions, for the jobs which minimize one of severa possibe objectives. The objectives we consider are the weighted ow time, the weighted numbers of tardy jobs, and the weighted sum of job tardiness. Genera resuts are obtained which can be simpied for certain 3

4 specia cases aready in the iteratures [8], [9] and [3]. In addition, some new resuts are obtained for specia cases not yet considered in the iterature. A Reation Between Competion Time and System Parameters Competion times are important random variabes in scheduing theory. Since our goa is to optimize performance characteristics reated to job competion time, it is the purpose of this section to nd a usefu reationship between the competion time of job n and the parameters of the system. We rst set out some notation which wi be used throughout this paper. Let EY represent the expectation of a random variabe Y. C (n) represent the competion time of the nth job (n = ; ; ; N) under the scheduing poicy ; is a permutation of the vector (; ; ; N). p! i = exp? (i = ; ; ; m), the mth roots of unity. i m i (s) = + s?! i m p?sd repair time distribution.?i i (s) = s (i = ; ; ; m), where D is a random variabe having the (i = ; ; ; m? ) and m (s) = h?m s + (??sd ) i. e = (; ; ; ), e i = (; ; ; ; ; ; ) and e represent the transpose of e. (n;i;j) denote a scheduing poicy in which job i and job j have been assigned to the nth and (n + )st positions, respectivey, (n = ; ; ; N? ; i; j = ; ; ; N). (n;j;i) is the poicy which is identica to (n;i;j) except for a switch in the nth and (n + )st positions. For the purpose of considering the competion time of the nth job (n = ; ; ; N) in any given poicy, we dene S n (t) to be the state of the system at time t before the nth job has been competed. This stochastic process can be in any of the foowing states [ fg f(i; k) : i = ; ; ; ; m; k = ; ; ; ng where denotes the absorbing state corresponding to the competion of the nth job, and (i; k) denotes the state in which the machine is processing the kth job whie the uptime is in phase i (i = ; ; ; m); phase corresponds to machine downtime. That is, state (; k) corresponds to the machine being repaired after having broken down whie the kth job is being processed. 4

5 Ceary, the process S n (t) is not Markovian, but it can be extended to a Markov process with the use of suppementary variabes (e.g. Cox[9], Chaudhry and Tempeton[8]). To do this, we denote the eapsed processing time of the kth job at t by V k (t), and the eapsed repair time at t by V (t). Then f(s n (t); V k (t); V (t)); t ; k = ; ; ; ng is a Markov process. Set p i;k (t; x)dx = P fs n (t) = (i; k); x V k (t) < x + dxjs n () = (; ); V () = g p ;k (t; x; y)dy = P fs n (t) = (; k); V k (t) = x; y V (t) < y + dyjs n () = (; ); V () = g where i = ; ; ; m and k = ; ; ; n: Notice that the kth job in the schedue is the job whose processing time has hazard function (k). Lemma. The system is governed by the foowing (k)(x) p ;k (t; x) = p ;k (t; x; y)(y)dy; k = ; ; ; @x + + (k)(x) p i;k (t; x) = p i?;k (t; x) i = ; ; m; k = ; ; n () with boundary conditions and @y + (y) p ;k (t; x; y) = k = ; ; ; n (3) p i; (t; ) = (i = ; ; ; n) (4) p i;k (t; ) = Z p i;k? (t; x) (k?) (x)dx; i = ; ; ; m; k = ; ; n (5) p ;k (t; x; ) = p m;k (t; x); k = ; ; ; n (6) p i;k (; x) = i; k; (x); i = ; ; ; m; k = ; ; ; n (7) p ;k (; x; y) = k = ; ; ; n (8) where (x) denotes the Dirac distribution and i;k denotes the Kronecker deta. Some detais of the proof are suppied in the appendix to this paper. The previous emma aows us to nd a reation between the Lapace transform of the competion times and the Lapace transform of the processing times and downtimes, a resut which is crucia to most of the scheduing poicies provided ater. 5

6 Lemma. For n = ; ; ; N; and s,?sc (n) =? m m m i= The proof is suppied in the appendix. i (s)!?i k k (s) 4? E (s) n j= (j) A As a specia case, we can obtain the foowing resut, previousy obtained in [3]. 3 5 Coroary.3 When the uptime is exponentiay distributed with rate, then for n = ; ; ; N; and s,?sc (n) = E exp@?[s + (??sd )] n j= (j) A Proof: This foows directy from Lemma., since when m =, (s) = (s) = s + (??sd ). Remark. We note that when the uptime has a more genera form it is not aways possibe to obtain the Lapace transform of C (n) in terms of system parameters. Even if the Lapace transform of the competion times can be found, it is not aways convenient for obtaining optima scheduing poicies. In the case discussed, as we wi see in ater sections, the above resut gives an exceedingy usefu reation. 3 The Optima Schedue for the Weighted Fow Time In this section, we sha nd schedues which minimize the expectation of the tota weighted competion times, P k h (k) C (k). This quantity gives an indication of the tota hoding or inventory costs incurred by a schedue. The sum of the competion times is often referred to as the ow time whie the tota weighted competion time is referred to as the weighted ow time. Some simpe permutation poicies have been found, under the assumption that the machine's uptime is exponentiay distributed (see [3], [8] and [], for exampe). However, as far as we know, resuts depending on parameters of both processing and the machine for the case of nonexponentia uptime have not been reported. Here, we sha consider this often-studied scheduing probem in the sense of expectation. Frequenty, however, the poicies that minimize the objective in expectation minimize the objective stochasticay as we (see pp. 8, []). We have the foowing necessary and sucient condition: 6

7 Theorem 3. The scheduing poicy (n;i;j) is better than the poicy (n;j;i) if and ony if ( j h j + m + m? i! m = m? i h i ( + m + j m = [?? i ]E exp! [?? j ]E exp where = () = (?! ) ( = ; ; ; m? ). "? "? (k) + j!#)!#) (k) + i Proof: According to Lemma., we know that for any poicy the expected competion time of the nth job is: EC (n) = P where d in = n i? E exp? i N E =?h i ( + m h (n;i;j) (k) C (n;i;j) (k)!? + j m + n m i= + m? (i) m i=! i d in (k). Therefore, by direct cacuation, we obtain: m? = ( +h j + + m? m i m Thus, the Theorem is proved. = N h (n;j;i) (k) C (n;j;i) (k)!! [?? j ]E exp! [?? i ]E exp "? "? (k) + i!#)!#) (k) + j It shoud be pointed out that the above necessary and sucient condition is dependent on the sequence before the nth job in a schedue, so it is impossibe to get an index optima poicy reated to both the parameters of the processing times and the machine. However, since the above condition has a nice structure, some optima schedue poicies can sti be found as foows: Theorem 3. The scheduing poicy = (; ; ; N) is optima if for i < j and t where and the quantity H C (t; i; j) = h i i + j m h C (x) = me?x H C (t; i; j) H C (t; j; i) Z h C (x + t)p ( i x i + j )dx (x) km? (km? )! ; x is ony dependent on the parameters of the uptime distribution of the machine. 7

8 Proof: Let EW F T [] represent the expected sum of the weighted ow times under poicy. We sha prove that EW F T [ ] EW F T [] for every poicy. Without oss of generaity, suppose there exists a poicy and an n as we as i < j satisfying (n) = j and (n + ) = i. We next show that EW F T [ (n;i;j) ] EW F T [ (n;j;i) ]. In fact, by the identity: m? =! e? x =? + me?x (x) km? (km? )! and ( h i +!#) + m?! [?? j ]E exp "? m j m (k) + i = 8 < = h i : + Z! 9 = h C (x)p j m (k) + i x (k) + i + j dx ; = i j Z H C (t; i; j)dp (k) t! as we as the resut in Theorem 3., we see that EW F T [ (n;i;j) ] EW F T [ (n;j;i) ]. Thus, the theorem foows from a pairwise interchange argument (see [6]). The foowing consequence was previousy found in [3]. Coroary 3.3 The expected weighted ow time is minimized by the schedue where the jobs are arranged in decreasing order of h i i, provided one of the foowing conditions hods:. m = ;. the job processing times are a exponentiay distributed. Proof:. This foows directy from Theorem 3., since H C (t; i; j) = h i i + if m =.. This foows directy from Theorem 3.. However, we oer an aternative proof to show how to use the concusion in Theorem 3. to obtain the resut for this specia case. Note that, in this case, ( i h i + m + j m = i h i 8 < : + m Z " m = m? =!! e? x [?? j ]E exp # dp "? (k) + i + j x (k) + i!#)!9 = ; 8

9 using the corresponding equation with j and i interchanged together with m! e?x = me?x (x) km? (km? )! = We see that the condition in Theorem 3. hods if j h j i h i. A pairwise interchange argument competes the proof. If we are interested ony in optima poicies which are independent of the machine parameters we can get the foowing resut previousy obtained in [9]. Coroary 3.4 The poicy = (; ; ; N) minimizes the expected weighted ow time if hods for i < j and t. h j P ( j t i + j ) h i P ( i t i + j ) Proof: This is a straightforward appication of Theorem 3., since the above condition impies h i i h j j. 4 The Optima Schedue for the Weighted Number of Tardy Jobs In this section, we sha consider the probem of minimizing the expectation of the weighted number of tardy jobs, P k h (k)u (k). Here, U (k) = ( ; if C(k) > d ; if C (k) d where d, the common due date of the jobs, is an exponentia random variabe with rate r. This quantity is an often-used objective in practice since it is a measure that can be recorded very easiy. It is a more genera cost function than the one studied in Section 3. Now, the cost is discounted at a rate of r per unit time. That is, if the jth job is not competed by time t, an additiona expected cost h j re?rt dt is incurred over the period [t; t + dt]. If job j is competed at time t, the tota expected cost incurred over the period [; t] is h j (? e?rt ). In the present paper, we sha nd optima schedues which minimize the expected vaue of this objective. Some simpe permutation poicies have been found previousy in the case of exponentia uptime (see [3] and [8], for exampe), but in the nonexponentia case, no resuts have been reported. 9

10 Since EU (n;i;j) (k) =? E exp[?rc (n;i;j) (k) ] hods for k = ; ; ; N, a direct cacuation using Lemma. gives the foowing necessary and sucient condition: Theorem 4. The scheduing poicy (n;i;j) is better than the poicy (n;j;i) if and ony if h j m = h i m = m p= m p=!?p p (r)?? (r) i E exp (r)!?p p (r)?? (r) j E exp (r) " "? (r)? (r) (k) + j!#!# (k) + i (9) We note here that it is not possibe to obtain a genera index optima poicy since the above necessary and sucient condition depends on the sequence before the nth job. However, one optima scheduing poicy can sti be found : Theorem 4. The poicy = (; ; ; N) is optima, if for i < j and t, where H U (t; i; j) = h i R H U (t; i; j) H U (t; j; i) h U (x + t)p ( i x i + j )dx and the quantity h U (x) = me?(+r)x m p= k= is ony dependent on the parameters of the machine. x mp?rd km+p? p (r) (km + p? )! ; x Proof: The eft-hand-side of (9) can be written as h j m = = h j = Z m!?p p= Z 4 m m = p= H U (t; j; i)dp p (r) E (r) " p (r)!?p exp? (r)( j + e? (r)x (k) t 3 (k) )! 5 P (k) + j x!? exp? (r)( (k) + j + i ) (k) + i + j! dx The resut now foows from a pairwise interchange argument and Theorem 4.. The foowing specia cases can be found in [3].!#

11 Coroary 4.3 When m =, the schedue which minimizes the expected weighted number of tardy jobs is in decreasing order of h i? i?? i, where = r + (??rd ). Proof: This foows from Theorem 4. upon noting that H U (t; i; j) = e?t h i? i?? j or, more simpy, by using Theorem 4. directy and a pairwise interchange argument. Coroary 4.4 If the processing times are a exponentia, the schedue which minimizes the expected weighted number of tardy jobs is in decreasing order of h i i Proof: This is a straightforward appication of Theorem 4.. Aternativey, this resut can be obtained by using the same idea as in Coroary 3.3 and m =!?p e? (r)x ; x Theorem 4. yieds a simpe optima poicy reated ony to the parameters of the processing time distribution as in [9]. Coroary 4.5 The poicy = (; ; ; N) minimizes the expected weighted number of tardy jobs if h j P ( j t i + j ) h i P ( i t i + j ) hods for any i < j and t. 5 The Optima Schedue for the Sum of Weighted Tardiness In this section, we sha consider schedues which minimize the expectation of the sum of weighted tardiness, P k h (k) T (k), where T (k) = maxf; C (k)? dg. Many papers have deat with this scheduing probem in the deterministic case, but few have considered the stochastic case. Pinedo [] was the rst to consider this scheduing probem. He assumed that the processing time of job i is an exponentia random variabe with rate i, the weight associated with job i is h i and job i is due at time d i, which is a random variabe. By using a compatibiity condition, i.e., i h i k h k impies d i st d k, he proved that the poicy

12 that minimizes this objective is the one where the jobs are arranged in decreasing order of h i i. However, when the machine is unreiabe, the question is more compicated. The case in which uptimes are arbitrary has been considered by [3] in case of exponentia processing times. We now consider the case of arbitrary processing times with Erang uptimes. We have the foowing necessary and sucient condition: Theorem 5. The scheduing poicy (n;i;j) is better than the poicy (n;j;i) if and ony if ( h j + + m?! [?? i ]E exp "? m i m? mr m m = p= h i ( + m? mr m m = p=!?p p (r) (r) j + m!?p p (r) (r) =?? (r) i E exp m? =! " [?? j ]E exp?? (r) j E exp "? (r) (k) + j!#!#9 = (k) + j ;!# (k) + i "?? (r)!#9 = (k) + i ; Proof: The argument foows the same ines as the proofs of Theorems 3. and Theorem 4., after one notes that hods for k = ; ; ; N. ET (n;i;j) (k) = EC (n;i;j) (k)? r h??rc (n;i;j) (k) i Theorem 5. The poicy = (; ; ; N) is optima if, for i < j and t, where H T (t; i; j) = h i R h T (x) = + e?x H T (t; i; j) H T (t; j; i) h T (x + t)p ( i x i + j )dx and the quantity (x) km? (km? )!? r e?(+r)x m k= p= (x ) is ony dependent on the parameters of the machine. x mp?rd km+p? p (r) (km + p? )! Proof: This foows from Theorem 5. and + m m = + e?x! k e? kx? mr m m p= (x) km? (km? )!? r e?(+r)x p!?p k e? k(r)x m k= p= x mp?rd km+p? p (r) (km + p? )!

13 Theorem 5. aso yieds the foowing resut Coroary 5.3 When m =, the optima schedue arranges the jobs in the order (; ; ; N) if the compatibiity conditions h j j h i i and h j? j?? j h i?i?? i hod for i < j, where = r + (??rd ). Theorem 5. yieds three additiona resuts, of which the rst is new, and the others were ontained in [3] and [9]. Coroary 5.4 When m =, the schedue which minimizes the expected sum of the weighted tardiness of the jobs arranges the jobs in the order (; ; ; N) if hods for i < j and t. i h i + j h j + r? e?t j? i r? e?t i? j?? j?? i Proof: Here, we note that H T (t; i; j) = h i + r? e?t j? i?? j r j Coroary 5.5 If the processing times are a exponentia, the schedue which minimizes the expected sum of the weighted tardiness of the jobs arranges the jobs in decreasing order of h i i. Coroary 5.6 The poicy = (; ; ; N) minimizes the expected sum of the weighted tardiness if h j P ( j t i + j ) h i P ( i t i + j ) hods for any i < j and t. 3

14 6 Concusions and Further Research In this paper, we have studied a cass of singe unreiabe machine stochastic scheduing probems, with genera processing time, Erang uptime and genera downtime distributions. We have obtained some genera optima scheduing poicies which are not index poicies and which depend on both processing parameters and machine parameters, under the nonpreemptive resume mode assumption. Some of the specia cases provided are aso new to the iterature. The method of suppementary variabes has been used successfuy here to nd usefu representations for the competion times. It shoud be noted that, because objective functions are often functions of competion times, much research focuses on competion times, a probem for which the method of suppementary variabes seems we-suited. It shoud be emphasized that nding optima schedues which depend on machine parameters has proved dicut in cases where the uptime distribution is nonexponentia. Thus, the extension to Erang uptimes given here represents a signicant step forward. We aso wish to point out here that the optimaity conditions found in this paper are both necessary and sucient. This is in contrast to other papers in the ed which usuay provide ony sucient conditions. When conditions sucient for optima poicies of simpe structure fai, it is benecia to understand the cost impications of impementing simpe poicies nevertheess. Li and Gazebrook [9] estabished an upper bound on the oss incurred when a processing poicy is adopted under the simpifying assumption of an exponentia processing requirement. Severa suboptimaity bounds for the nonexponentia case sha be given in another paper. One question for future research that was pointed out by the referee concerns the case of deterministic uptime. By reaceing by m and etting m become innitey arge, we can obtain the specia cases of a our resuts for constant uptime. Obtaining simpe poicies wi require some eort and wi be the subject of a subsequent paper in which numerica impementation of our poicies wi aso be considered. We have restricted our attention to the nonpreemptive resume mode, but it woud be interesting to see if simiar resuts can be found for the nonpreemptive repeat mode. In that case, the makespan woud seem to be an important objective to study rst. So far, we have found that if the uptime has an Erang distribution, the Lapace transform of the competion time of the nth job in any 4

15 given poicy possesses a nice structure, but optimaity resuts are sti sought. Acknowedgements The authors woud ike to thank the referee for a carefu reading of the origina manuscript. His/Her comments and suggestions have ed to a much better presentation of the paper. The authors acknowedge that this work was supported by research grants from the Natura Sciences and Engineering Research Counci of Canada (NSERC). Dr. Wei Li aso acknowedges that this work was supported by research grants from the Nationa Natura Science Foundation of China (NNSFC) and from the University of Winnipeg, Manitoba, Canada. 5

16 Appendix Proof of Lemma.: The proof of these equations foows a standard probabiistic argument (for exampe, see [8] or [9]). Since these equations are the starting point for a other deveopments in the paper, we provide some detais of the proof for () and (7). The arguments for (), (3), (4), (5) and (6) are simiar to that for (), and equation (8) is immediate from the denition. To obtain equation (), consider the state of the system at t and at t + t. In order that S n (t + t) = (; k) and V k (t + t) = x + t (x > ), it is necessary that (a) at time t, S n (t) = (; k) and V k (t) = x. Within (t; t + t), the uptime of the machine is sti in phase and the kth job has not been competed; or (b) at time t, S n (t) = (; k), and V k (t) = x, V (t) = y for any y >. Within (t; t + t), the repair of the machine has been competed. Using a conditioning argument, we have Thus, Z p ;k (t + t; x + t) = p ;k (t; x)[? t? (k) (x)t] + p ;k (t; x; y)(y)tdy + o(t) p ;k (t + t; x + t)? p ;k (t; x + t) + p ;k (t; x + t)? p ;k (t; x) =?p ;k (t; x)[t + (k) (x)t] + Dividing by t and etting t! gives equation (). Z p ;k (t; x; y)(y)tdy + o(t) Most of the equations given by (7) foow immediatey from the mode assuptions. The ony exception needing more proof is p ; (; x) = (x). That is, we must show p ; (; x) = for a x >, whie Z p ; (; x)dx = for any " > : Z " p ; (; x)dx = The rst reation is obvious, because it is impossibe for the machine to have processed the rst job for a non-zero amount of time at time t=. Next, by the denition of p ; (; x), it foows that Z p ; (; x)dx = Z P fs n () = (; ); x V () < x + dxjs n () = (; ); V () = g = Thus, together with the rst reation, we can get the second reation. 6

17 R Remark: p ;(; x)dx = can be aso directy obtained from the genera reation in the equation (7) used in the foowing proof of the Lemma.. Proof of Lemma.: Soving the Lapace transformed version of (3) with the boundary condition (6) gives p ;k(s; x; y) = e?sy [? D(y)]p m;k(s; x); k = ; ; ; n () where p (s) denotes the Lapace transform of a function p(t). Using (), (), (7) and the corresponding denitions, we p k(s; x) = Ap k(s; x)? (k) (x)p k(s; x); + k (x)e (k = ; ; ; n) () where p k (t; x) = (p ;k (t; x); p ;k (t; x); ; p m;k (t; x)), and A = Noting that? (s);? (s); ;? m (s) are the m dierent eigenvaues of matrix A, we can write?i i? A = T f? i (s)g T? () C A mm where fx i g = diagfx i g, T = m?..... C. A and T? = m!?m!?m! m?m m?!! m?......! m! m m? C A Substituting () into () and soving, we have p (s; x) = [? F () (x)]?i T f? i (s)xgt? i? U(x)e (3) and p k(s; x) = [? F (k) (x)]?i T f? i (s)xgt? i? L k (4) where L k (k = ; ; n) is a coumn vector to be determined, and U(x) is one if x > and zero otherwise. Using boundary conditions (4), (5) and the above two resuts and dening L = e, vectors L k (k = ; ; n) can be recursivey determined as L k =?i T fe exp[? i (s) (k?) ]gt? i? L k? 7

18 = =?i T? T i? 8 < : Then substituting (5) into (4) gives us Since p k(s; x) = [? F (k) (x)] T? P (C (n) > t) = i? n 4 m e? i(s)x E exp 4?i (s) k? j= 9 3 = (j) 5 ; e (5)?i Z i= T 8 < : e? i(s)x E exp 4?i (s) k? j= 3 (j) 5 e (k = ; 3; ; n) (6) p i;k (t; x)dx + Z Z the resuts in equations (), (3) and (6) can be combined to give = = s Z e?st P (C (n) > t)dt = n 4 m Z i= e + n Z s (??sd )e m p k(s; x)dx h i se + (? E?sD ) T? i? = s e f i(s)g T :? i (s) = sm m m i= 8 < i (s)!?i k k (s) e?i p i;k(s; x)dx + T 8 < 4? E (s) 4? E (s) Z Z :? i (s) n j= (j) A (j) A 3 9 = ; p ;k (t; x; y)dxdy5 (7) p ;k(t; x; y)dxdy 3 5 4? E (s) 39 = 5 The ast equation above foows since it is a specia case of the formua n j= efy i gt fd i gt? fx i ge = m d k y i!?i k m i= when x = and x i = (i = ; 3; ; m). Now, by using the reationship: E[e?sC (n)] =? s Z 8 m 3 ; T? n j= i? (j) A e 39 = 5 ; 5 (8) e?st P (C (n) > m j= x j! j k A

19 the proof is readiy competed. References [] Adiri, I., Bruno, J., Frostig, E. and Rinnooy Kan, A.H.G. (989) Singe machine ow-time scheduing with a singe breakdown. Acta Informatica 6: 679{696. [] Aahverdi, Ai. (995) Two-stage production scheduing with separated setup times and stochastic breakdowns. The Journa of the Operationa Research Society 46: 896{94. [3] Aahverdi, Ai. and Mittentha, J. (994) Scheduing on M parae machines subject to random breakdowns to minimize expected mean ow time. Nava Research Logistics 4: 677{68. [4] Aahverdi, Ai. and Mittentha, J. (994) Two-machine ordered owshop scheduing under random breakdowns. Mathematica and Computer Modeing : 9{7. [5] Aahverdi, Ai. and Mittentha, J. (995) Scheduing on a two-machine owshop subject to random breakdowns with a makespan objective function. European Journa of Operationa Research 8: 376{387. [6] Birge, J., Frenk, J.B.G., Mittentha, J., and Rinnooy Kan, A.H.G. (99) Singe machine scheduing subject to stochastic breakdowns. Nava Research Logistics 37: 66{677. [7] Chang, C.S., Chao,.L., Pinedo, M. and Weber, R. (99) On the optimaity of LEPT and C Rues for machines in parae. J. App. Prob. 9: 667{68. [8] Chaudhry, M.L. and Tempetion, J.G.C. (983) A First Course in Buk Queues, Wiey, New York. [9] Cox, D.R. (955) The anaysis of non-markovian stochastic processes by the incusion of suppementary variabes. Proceedings of the Cambridge Phiosophica Society, 5: [] Du, Q. and Pinedo, M. (995) A note on minimizing the expected makespan in owshops subject to breakdowns. Nava Research Logistics, 5{6. [] Frenk, J.B.G. (987) Renewa theory and competey monotone functions. Report No. 8759/A, Econometric Institute, Erasmus University, Rotterdam. [] Frenk, J.B.G. (99) A genera framework for stochastic one-machine scheduing probems with zero reease times and no partia ordering. Probabiity in the Engineering and Informationa Sciences 5: 97{35. [3] Frostig, E. (99) A note on stochastic scheduing on a singe machine subject to breakdown the preemptive repeat mode. Probabiity in the Engineering and Informationa Sciences 5: 349{354. 9

20 [4] Gazebrook, K.D. (984) Scheduing stochastic jobs on a singe machine subject to breakdowns. Nava Research Logistics Quartery 3: 5{64. [5] Gazebrook, K.D. (987) Evauating the eects of machine breakdowns in stochastic scheduing probems. Nava Research Logistics Quartery 34: 39{335. [6] Gazebrook, K.D. (99) On non-preemptive poicies for stochastic singe machine scheduing with breakdown. Probabiity in the Engineering and Informationa Sciences 5: 77{87. [7] Li, W. and Cao, J.H. (994) Stochastic scheduing on an unreiabe machine with genera uptimes and genera setup times. Journa of System Engineering and System Sciences 3: 79{ 88. [8] Li, W. and Cao, J.H. (995) Stochastic scheduing on a singe machine subject to mutipe breakdowns according to dierent probabiities. Operations Research Letters 8: 8{9. [9] Li, W. and Gazebrook, K.D. (997) On stochastic machine scheduing with genera distributiona assumptions. To appear in European Journa of Operationa Research. [] Mittentha, I. (986) Scheduing on singe machine subject to breakdowns. Ph.D. dissertation. University of Michigan, Ann Arbor, MI. [] Pinedo, M. (983) Stochastic scheduing with reease dates and due dates. Operations Research 3: 559{57. [] Pinedo, M. (995) Scheduing: theory, agorithms, and systems. Prentice-Ha, Engewood Cis, NJ. [3] Pinedo, M. and Rammouz, E. (988) A note on stochastic scheduing on a singe machines subject to breakdown and repair. Probabiity in the Engineering and Informationa Sciences : 4{49. [4] Righter, R. (994) Scheduing. Stochastic Orders and Their Appications. ed. Shaked, M. and Shanthikumar, J.G., Academic Press, New York, 38{43. [5] Righter, R. and Shanthikumar,J.G. (989) Scheduing muticass singe server queueing systems to stochasticay maximise the number of successfu departures. Probabiity in the Engineering and Informationa Sciences 3: 33{333. [6] Ross, S.M. (983) Introduction to stochastic dynamic programming. Academic Press, New York.

Stochastic Complement Analysis of Multi-Server Threshold Queues. with Hysteresis. Abstract

Stochastic Complement Analysis of Multi-Server Threshold Queues. with Hysteresis. Abstract Stochastic Compement Anaysis of Muti-Server Threshod Queues with Hysteresis John C.S. Lui The Dept. of Computer Science & Engineering The Chinese University of Hong Kong Leana Goubchik Dept. of Computer

More information

A Brief Introduction to Markov Chains and Hidden Markov Models

A Brief Introduction to Markov Chains and Hidden Markov Models A Brief Introduction to Markov Chains and Hidden Markov Modes Aen B MacKenzie Notes for December 1, 3, &8, 2015 Discrete-Time Markov Chains You may reca that when we first introduced random processes,

More information

Minimizing Total Weighted Completion Time on Uniform Machines with Unbounded Batch

Minimizing Total Weighted Completion Time on Uniform Machines with Unbounded Batch The Eighth Internationa Symposium on Operations Research and Its Appications (ISORA 09) Zhangiaie, China, September 20 22, 2009 Copyright 2009 ORSC & APORC, pp. 402 408 Minimizing Tota Weighted Competion

More information

MARKOV CHAINS AND MARKOV DECISION THEORY. Contents

MARKOV CHAINS AND MARKOV DECISION THEORY. Contents MARKOV CHAINS AND MARKOV DECISION THEORY ARINDRIMA DATTA Abstract. In this paper, we begin with a forma introduction to probabiity and expain the concept of random variabes and stochastic processes. After

More information

<C 2 2. λ 2 l. λ 1 l 1 < C 1

<C 2 2. λ 2 l. λ 1 l 1 < C 1 Teecommunication Network Contro and Management (EE E694) Prof. A. A. Lazar Notes for the ecture of 7/Feb/95 by Huayan Wang (this document was ast LaT E X-ed on May 9,995) Queueing Primer for Muticass Optima

More information

Stochastic Automata Networks (SAN) - Modelling. and Evaluation. Paulo Fernandes 1. Brigitte Plateau 2. May 29, 1997

Stochastic Automata Networks (SAN) - Modelling. and Evaluation. Paulo Fernandes 1. Brigitte Plateau 2. May 29, 1997 Stochastic utomata etworks (S) - Modeing and Evauation Pauo Fernandes rigitte Pateau 2 May 29, 997 Institut ationa Poytechnique de Grenobe { IPG Ecoe ationae Superieure d'informatique et de Mathematiques

More information

Uniprocessor Feasibility of Sporadic Tasks with Constrained Deadlines is Strongly conp-complete

Uniprocessor Feasibility of Sporadic Tasks with Constrained Deadlines is Strongly conp-complete Uniprocessor Feasibiity of Sporadic Tasks with Constrained Deadines is Strongy conp-compete Pontus Ekberg and Wang Yi Uppsaa University, Sweden Emai: {pontus.ekberg yi}@it.uu.se Abstract Deciding the feasibiity

More information

Schedulability Analysis of Deferrable Scheduling Algorithms for Maintaining Real-Time Data Freshness

Schedulability Analysis of Deferrable Scheduling Algorithms for Maintaining Real-Time Data Freshness 1 Scheduabiity Anaysis of Deferrabe Scheduing Agorithms for Maintaining Rea-Time Data Freshness Song Han, Deji Chen, Ming Xiong, Kam-yiu Lam, Aoysius K. Mok, Krithi Ramamritham UT Austin, Emerson Process

More information

CS229 Lecture notes. Andrew Ng

CS229 Lecture notes. Andrew Ng CS229 Lecture notes Andrew Ng Part IX The EM agorithm In the previous set of notes, we taked about the EM agorithm as appied to fitting a mixture of Gaussians. In this set of notes, we give a broader view

More information

A NOTE ON QUASI-STATIONARY DISTRIBUTIONS OF BIRTH-DEATH PROCESSES AND THE SIS LOGISTIC EPIDEMIC

A NOTE ON QUASI-STATIONARY DISTRIBUTIONS OF BIRTH-DEATH PROCESSES AND THE SIS LOGISTIC EPIDEMIC (January 8, 2003) A NOTE ON QUASI-STATIONARY DISTRIBUTIONS OF BIRTH-DEATH PROCESSES AND THE SIS LOGISTIC EPIDEMIC DAMIAN CLANCY, University of Liverpoo PHILIP K. POLLETT, University of Queensand Abstract

More information

XSAT of linear CNF formulas

XSAT of linear CNF formulas XSAT of inear CN formuas Bernd R. Schuh Dr. Bernd Schuh, D-50968 Kön, Germany; bernd.schuh@netcoogne.de eywords: compexity, XSAT, exact inear formua, -reguarity, -uniformity, NPcompeteness Abstract. Open

More information

Schedulability Analysis of Deferrable Scheduling Algorithms for Maintaining Real-Time Data Freshness

Schedulability Analysis of Deferrable Scheduling Algorithms for Maintaining Real-Time Data Freshness 1 Scheduabiity Anaysis of Deferrabe Scheduing Agorithms for Maintaining Rea- Data Freshness Song Han, Deji Chen, Ming Xiong, Kam-yiu Lam, Aoysius K. Mok, Krithi Ramamritham UT Austin, Emerson Process Management,

More information

Optimal Control of Assembly Systems with Multiple Stages and Multiple Demand Classes 1

Optimal Control of Assembly Systems with Multiple Stages and Multiple Demand Classes 1 Optima Contro of Assemby Systems with Mutipe Stages and Mutipe Demand Casses Saif Benjaafar Mohsen EHafsi 2 Chung-Yee Lee 3 Weihua Zhou 3 Industria & Systems Engineering, Department of Mechanica Engineering,

More information

Math 124B January 31, 2012

Math 124B January 31, 2012 Math 124B January 31, 212 Viktor Grigoryan 7 Inhomogeneous boundary vaue probems Having studied the theory of Fourier series, with which we successfuy soved boundary vaue probems for the homogeneous heat

More information

Separation of Variables and a Spherical Shell with Surface Charge

Separation of Variables and a Spherical Shell with Surface Charge Separation of Variabes and a Spherica She with Surface Charge In cass we worked out the eectrostatic potentia due to a spherica she of radius R with a surface charge density σθ = σ cos θ. This cacuation

More information

STABLE GRAPHS BENJAMIN OYE

STABLE GRAPHS BENJAMIN OYE STABLE GRAPHS BENJAMIN OYE Abstract. In Reguarity Lemmas for Stabe Graphs [1] Maiaris and Sheah appy toos from mode theory to obtain stronger forms of Ramsey's theorem and Szemeredi's reguarity emma for

More information

Copyright information to be inserted by the Publishers. Unsplitting BGK-type Schemes for the Shallow. Water Equations KUN XU

Copyright information to be inserted by the Publishers. Unsplitting BGK-type Schemes for the Shallow. Water Equations KUN XU Copyright information to be inserted by the Pubishers Unspitting BGK-type Schemes for the Shaow Water Equations KUN XU Mathematics Department, Hong Kong University of Science and Technoogy, Cear Water

More information

Integrating Factor Methods as Exponential Integrators

Integrating Factor Methods as Exponential Integrators Integrating Factor Methods as Exponentia Integrators Borisav V. Minchev Department of Mathematica Science, NTNU, 7491 Trondheim, Norway Borko.Minchev@ii.uib.no Abstract. Recenty a ot of effort has been

More information

Asynchronous Control for Coupled Markov Decision Systems

Asynchronous Control for Coupled Markov Decision Systems INFORMATION THEORY WORKSHOP (ITW) 22 Asynchronous Contro for Couped Marov Decision Systems Michae J. Neey University of Southern Caifornia Abstract This paper considers optima contro for a coection of

More information

Throughput Optimal Scheduling for Wireless Downlinks with Reconfiguration Delay

Throughput Optimal Scheduling for Wireless Downlinks with Reconfiguration Delay Throughput Optima Scheduing for Wireess Downinks with Reconfiguration Deay Vineeth Baa Sukumaran vineethbs@gmai.com Department of Avionics Indian Institute of Space Science and Technoogy. Abstract We consider

More information

1D Heat Propagation Problems

1D Heat Propagation Problems Chapter 1 1D Heat Propagation Probems If the ambient space of the heat conduction has ony one dimension, the Fourier equation reduces to the foowing for an homogeneous body cρ T t = T λ 2 + Q, 1.1) x2

More information

Cryptanalysis of PKP: A New Approach

Cryptanalysis of PKP: A New Approach Cryptanaysis of PKP: A New Approach Éiane Jaumes and Antoine Joux DCSSI 18, rue du Dr. Zamenhoff F-92131 Issy-es-Mx Cedex France eiane.jaumes@wanadoo.fr Antoine.Joux@ens.fr Abstract. Quite recenty, in

More information

Bayesian Learning. You hear a which which could equally be Thanks or Tanks, which would you go with?

Bayesian Learning. You hear a which which could equally be Thanks or Tanks, which would you go with? Bayesian Learning A powerfu and growing approach in machine earning We use it in our own decision making a the time You hear a which which coud equay be Thanks or Tanks, which woud you go with? Combine

More information

$, (2.1) n="# #. (2.2)

$, (2.1) n=# #. (2.2) Chapter. Eectrostatic II Notes: Most of the materia presented in this chapter is taken from Jackson, Chap.,, and 4, and Di Bartoo, Chap... Mathematica Considerations.. The Fourier series and the Fourier

More information

Mat 1501 lecture notes, penultimate installment

Mat 1501 lecture notes, penultimate installment Mat 1501 ecture notes, penutimate instament 1. bounded variation: functions of a singe variabe optiona) I beieve that we wi not actuay use the materia in this section the point is mainy to motivate the

More information

MATH 172: MOTIVATION FOR FOURIER SERIES: SEPARATION OF VARIABLES

MATH 172: MOTIVATION FOR FOURIER SERIES: SEPARATION OF VARIABLES MATH 172: MOTIVATION FOR FOURIER SERIES: SEPARATION OF VARIABLES Separation of variabes is a method to sove certain PDEs which have a warped product structure. First, on R n, a inear PDE of order m is

More information

C. Fourier Sine Series Overview

C. Fourier Sine Series Overview 12 PHILIP D. LOEWEN C. Fourier Sine Series Overview Let some constant > be given. The symboic form of the FSS Eigenvaue probem combines an ordinary differentia equation (ODE) on the interva (, ) with a

More information

Age of Information: The Gamma Awakening

Age of Information: The Gamma Awakening Age of Information: The Gamma Awakening Eie Najm and Rajai Nasser LTHI, EPFL, Lausanne, Switzerand Emai: {eie.najm, rajai.nasser}@epf.ch arxiv:604.086v [cs.it] 5 Apr 06 Abstract Status update systems is

More information

Explicit overall risk minimization transductive bound

Explicit overall risk minimization transductive bound 1 Expicit overa risk minimization transductive bound Sergio Decherchi, Paoo Gastado, Sandro Ridea, Rodofo Zunino Dept. of Biophysica and Eectronic Engineering (DIBE), Genoa University Via Opera Pia 11a,

More information

Reliability: Theory & Applications No.3, September 2006

Reliability: Theory & Applications No.3, September 2006 REDUNDANCY AND RENEWAL OF SERVERS IN OPENED QUEUING NETWORKS G. Sh. Tsitsiashvii M.A. Osipova Vadivosto, Russia 1 An opened queuing networ with a redundancy and a renewa of servers is considered. To cacuate

More information

BP neural network-based sports performance prediction model applied research

BP neural network-based sports performance prediction model applied research Avaiabe onine www.jocpr.com Journa of Chemica and Pharmaceutica Research, 204, 6(7:93-936 Research Artice ISSN : 0975-7384 CODEN(USA : JCPRC5 BP neura networ-based sports performance prediction mode appied

More information

A Comparison Study of the Test for Right Censored and Grouped Data

A Comparison Study of the Test for Right Censored and Grouped Data Communications for Statistica Appications and Methods 2015, Vo. 22, No. 4, 313 320 DOI: http://dx.doi.org/10.5351/csam.2015.22.4.313 Print ISSN 2287-7843 / Onine ISSN 2383-4757 A Comparison Study of the

More information

NIKOS FRANTZIKINAKIS. N n N where (Φ N) N N is any Følner sequence

NIKOS FRANTZIKINAKIS. N n N where (Φ N) N N is any Følner sequence SOME OPE PROBLEMS O MULTIPLE ERGODIC AVERAGES IKOS FRATZIKIAKIS. Probems reated to poynomia sequences In this section we give a ist of probems reated to the study of mutipe ergodic averages invoving iterates

More information

(f) is called a nearly holomorphic modular form of weight k + 2r as in [5].

(f) is called a nearly holomorphic modular form of weight k + 2r as in [5]. PRODUCTS OF NEARLY HOLOMORPHIC EIGENFORMS JEFFREY BEYERL, KEVIN JAMES, CATHERINE TRENTACOSTE, AND HUI XUE Abstract. We prove that the product of two neary hoomorphic Hece eigenforms is again a Hece eigenform

More information

Some Measures for Asymmetry of Distributions

Some Measures for Asymmetry of Distributions Some Measures for Asymmetry of Distributions Georgi N. Boshnakov First version: 31 January 2006 Research Report No. 5, 2006, Probabiity and Statistics Group Schoo of Mathematics, The University of Manchester

More information

FRST Multivariate Statistics. Multivariate Discriminant Analysis (MDA)

FRST Multivariate Statistics. Multivariate Discriminant Analysis (MDA) 1 FRST 531 -- Mutivariate Statistics Mutivariate Discriminant Anaysis (MDA) Purpose: 1. To predict which group (Y) an observation beongs to based on the characteristics of p predictor (X) variabes, using

More information

A Statistical Framework for Real-time Event Detection in Power Systems

A Statistical Framework for Real-time Event Detection in Power Systems 1 A Statistica Framework for Rea-time Event Detection in Power Systems Noan Uhrich, Tim Christman, Phiip Swisher, and Xichen Jiang Abstract A quickest change detection (QCD) agorithm is appied to the probem

More information

6 Wave Equation on an Interval: Separation of Variables

6 Wave Equation on an Interval: Separation of Variables 6 Wave Equation on an Interva: Separation of Variabes 6.1 Dirichet Boundary Conditions Ref: Strauss, Chapter 4 We now use the separation of variabes technique to study the wave equation on a finite interva.

More information

Improving the Reliability of a Series-Parallel System Using Modified Weibull Distribution

Improving the Reliability of a Series-Parallel System Using Modified Weibull Distribution Internationa Mathematica Forum, Vo. 12, 217, no. 6, 257-269 HIKARI Ltd, www.m-hikari.com https://doi.org/1.12988/imf.217.611155 Improving the Reiabiity of a Series-Parae System Using Modified Weibu Distribution

More information

Reichenbachian Common Cause Systems

Reichenbachian Common Cause Systems Reichenbachian Common Cause Systems G. Hofer-Szabó Department of Phiosophy Technica University of Budapest e-mai: gszabo@hps.ete.hu Mikós Rédei Department of History and Phiosophy of Science Eötvös University,

More information

Lecture Note 3: Stationary Iterative Methods

Lecture Note 3: Stationary Iterative Methods MATH 5330: Computationa Methods of Linear Agebra Lecture Note 3: Stationary Iterative Methods Xianyi Zeng Department of Mathematica Sciences, UTEP Stationary Iterative Methods The Gaussian eimination (or

More information

On the construction of bivariate exponential distributions with an arbitrary correlation coefficient

On the construction of bivariate exponential distributions with an arbitrary correlation coefficient Downoaded from orbit.dtu.d on: Sep 23, 208 On the construction of bivariate exponentia distributions with an arbitrary correation coefficient Badt, Mogens; Niesen, Bo Friis Pubication date: 2008 Document

More information

Online Load Balancing on Related Machines

Online Load Balancing on Related Machines Onine Load Baancing on Reated Machines ABSTRACT Sungjin Im University of Caifornia at Merced Merced, CA, USA sim3@ucmerced.edu Debmaya Panigrahi Duke University Durham, NC, USA debmaya@cs.duke.edu We give

More information

14 Separation of Variables Method

14 Separation of Variables Method 14 Separation of Variabes Method Consider, for exampe, the Dirichet probem u t = Du xx < x u(x, ) = f(x) < x < u(, t) = = u(, t) t > Let u(x, t) = T (t)φ(x); now substitute into the equation: dt

More information

Appendix of the Paper The Role of No-Arbitrage on Forecasting: Lessons from a Parametric Term Structure Model

Appendix of the Paper The Role of No-Arbitrage on Forecasting: Lessons from a Parametric Term Structure Model Appendix of the Paper The Roe of No-Arbitrage on Forecasting: Lessons from a Parametric Term Structure Mode Caio Ameida cameida@fgv.br José Vicente jose.vaentim@bcb.gov.br June 008 1 Introduction In this

More information

Smoothness equivalence properties of univariate subdivision schemes and their projection analogues

Smoothness equivalence properties of univariate subdivision schemes and their projection analogues Numerische Mathematik manuscript No. (wi be inserted by the editor) Smoothness equivaence properties of univariate subdivision schemes and their projection anaogues Phiipp Grohs TU Graz Institute of Geometry

More information

First-Order Corrections to Gutzwiller s Trace Formula for Systems with Discrete Symmetries

First-Order Corrections to Gutzwiller s Trace Formula for Systems with Discrete Symmetries c 26 Noninear Phenomena in Compex Systems First-Order Corrections to Gutzwier s Trace Formua for Systems with Discrete Symmetries Hoger Cartarius, Jörg Main, and Günter Wunner Institut für Theoretische

More information

Chemical Kinetics Part 2

Chemical Kinetics Part 2 Integrated Rate Laws Chemica Kinetics Part 2 The rate aw we have discussed thus far is the differentia rate aw. Let us consider the very simpe reaction: a A à products The differentia rate reates the rate

More information

Week 6 Lectures, Math 6451, Tanveer

Week 6 Lectures, Math 6451, Tanveer Fourier Series Week 6 Lectures, Math 645, Tanveer In the context of separation of variabe to find soutions of PDEs, we encountered or and in other cases f(x = f(x = a 0 + f(x = a 0 + b n sin nπx { a n

More information

FRIEZE GROUPS IN R 2

FRIEZE GROUPS IN R 2 FRIEZE GROUPS IN R 2 MAXWELL STOLARSKI Abstract. Focusing on the Eucidean pane under the Pythagorean Metric, our goa is to cassify the frieze groups, discrete subgroups of the set of isometries of the

More information

On the Goal Value of a Boolean Function

On the Goal Value of a Boolean Function On the Goa Vaue of a Booean Function Eric Bach Dept. of CS University of Wisconsin 1210 W. Dayton St. Madison, WI 53706 Lisa Heerstein Dept of CSE NYU Schoo of Engineering 2 Metrotech Center, 10th Foor

More information

4 1-D Boundary Value Problems Heat Equation

4 1-D Boundary Value Problems Heat Equation 4 -D Boundary Vaue Probems Heat Equation The main purpose of this chapter is to study boundary vaue probems for the heat equation on a finite rod a x b. u t (x, t = ku xx (x, t, a < x < b, t > u(x, = ϕ(x

More information

A simple reliability block diagram method for safety integrity verification

A simple reliability block diagram method for safety integrity verification Reiabiity Engineering and System Safety 92 (2007) 1267 1273 www.esevier.com/ocate/ress A simpe reiabiity bock diagram method for safety integrity verification Haitao Guo, Xianhui Yang epartment of Automation,

More information

Stochastic Automata Networks: Product Forms and Iterative Solutions.

Stochastic Automata Networks: Product Forms and Iterative Solutions. INSTITUT NATIONAL DE REHERHE EN INFORMATIQUE ET EN AUTOMATIQUE Stochastic Automata Networks: Product Forms and Iterative Soutions. rigitte Pateau and Wiiam J. Stewart N 2939 2 juiet 996 THÈME apport de

More information

Chemical Kinetics Part 2. Chapter 16

Chemical Kinetics Part 2. Chapter 16 Chemica Kinetics Part 2 Chapter 16 Integrated Rate Laws The rate aw we have discussed thus far is the differentia rate aw. Let us consider the very simpe reaction: a A à products The differentia rate reates

More information

A SIMPLIFIED DESIGN OF MULTIDIMENSIONAL TRANSFER FUNCTION MODELS

A SIMPLIFIED DESIGN OF MULTIDIMENSIONAL TRANSFER FUNCTION MODELS A SIPLIFIED DESIGN OF ULTIDIENSIONAL TRANSFER FUNCTION ODELS Stefan Petrausch, Rudof Rabenstein utimedia Communications and Signa Procesg, University of Erangen-Nuremberg, Cauerstr. 7, 958 Erangen, GERANY

More information

Efficiently Generating Random Bits from Finite State Markov Chains

Efficiently Generating Random Bits from Finite State Markov Chains 1 Efficienty Generating Random Bits from Finite State Markov Chains Hongchao Zhou and Jehoshua Bruck, Feow, IEEE Abstract The probem of random number generation from an uncorreated random source (of unknown

More information

6.434J/16.391J Statistics for Engineers and Scientists May 4 MIT, Spring 2006 Handout #17. Solution 7

6.434J/16.391J Statistics for Engineers and Scientists May 4 MIT, Spring 2006 Handout #17. Solution 7 6.434J/16.391J Statistics for Engineers and Scientists May 4 MIT, Spring 2006 Handout #17 Soution 7 Probem 1: Generating Random Variabes Each part of this probem requires impementation in MATLAB. For the

More information

Approximate Bandwidth Allocation for Fixed-Priority-Scheduled Periodic Resources (WSU-CS Technical Report Version)

Approximate Bandwidth Allocation for Fixed-Priority-Scheduled Periodic Resources (WSU-CS Technical Report Version) Approximate Bandwidth Aocation for Fixed-Priority-Schedued Periodic Resources WSU-CS Technica Report Version) Farhana Dewan Nathan Fisher Abstract Recent research in compositiona rea-time systems has focused

More information

The EM Algorithm applied to determining new limit points of Mahler measures

The EM Algorithm applied to determining new limit points of Mahler measures Contro and Cybernetics vo. 39 (2010) No. 4 The EM Agorithm appied to determining new imit points of Maher measures by Souad E Otmani, Georges Rhin and Jean-Marc Sac-Épée Université Pau Veraine-Metz, LMAM,

More information

Approximation and Fast Calculation of Non-local Boundary Conditions for the Time-dependent Schrödinger Equation

Approximation and Fast Calculation of Non-local Boundary Conditions for the Time-dependent Schrödinger Equation Approximation and Fast Cacuation of Non-oca Boundary Conditions for the Time-dependent Schrödinger Equation Anton Arnod, Matthias Ehrhardt 2, and Ivan Sofronov 3 Universität Münster, Institut für Numerische

More information

arxiv:hep-ph/ v1 15 Jan 2001

arxiv:hep-ph/ v1 15 Jan 2001 BOSE-EINSTEIN CORRELATIONS IN CASCADE PROCESSES AND NON-EXTENSIVE STATISTICS O.V.UTYUZH AND G.WILK The Andrzej So tan Institute for Nucear Studies; Hoża 69; 00-689 Warsaw, Poand E-mai: utyuzh@fuw.edu.p

More information

Strauss PDEs 2e: Section Exercise 2 Page 1 of 12. For problem (1), complete the calculation of the series in case j(t) = 0 and h(t) = e t.

Strauss PDEs 2e: Section Exercise 2 Page 1 of 12. For problem (1), complete the calculation of the series in case j(t) = 0 and h(t) = e t. Strauss PDEs e: Section 5.6 - Exercise Page 1 of 1 Exercise For probem (1, compete the cacuation of the series in case j(t = and h(t = e t. Soution With j(t = and h(t = e t, probem (1 on page 147 becomes

More information

Efficient Generation of Random Bits from Finite State Markov Chains

Efficient Generation of Random Bits from Finite State Markov Chains Efficient Generation of Random Bits from Finite State Markov Chains Hongchao Zhou and Jehoshua Bruck, Feow, IEEE Abstract The probem of random number generation from an uncorreated random source (of unknown

More information

Asymptotic Properties of a Generalized Cross Entropy Optimization Algorithm

Asymptotic Properties of a Generalized Cross Entropy Optimization Algorithm 1 Asymptotic Properties of a Generaized Cross Entropy Optimization Agorithm Zijun Wu, Michae Koonko, Institute for Appied Stochastics and Operations Research, Caustha Technica University Abstract The discrete

More information

4 Separation of Variables

4 Separation of Variables 4 Separation of Variabes In this chapter we describe a cassica technique for constructing forma soutions to inear boundary vaue probems. The soution of three cassica (paraboic, hyperboic and eiptic) PDE

More information

Maintenance activities planning and grouping for complex structure systems

Maintenance activities planning and grouping for complex structure systems Maintenance activities panning and grouping for compex structure systems Hai Canh u, Phuc Do an, Anne Barros, Christophe Bérenguer To cite this version: Hai Canh u, Phuc Do an, Anne Barros, Christophe

More information

Partial permutation decoding for MacDonald codes

Partial permutation decoding for MacDonald codes Partia permutation decoding for MacDonad codes J.D. Key Department of Mathematics and Appied Mathematics University of the Western Cape 7535 Bevie, South Africa P. Seneviratne Department of Mathematics

More information

T.C. Banwell, S. Galli. {bct, Telcordia Technologies, Inc., 445 South Street, Morristown, NJ 07960, USA

T.C. Banwell, S. Galli. {bct, Telcordia Technologies, Inc., 445 South Street, Morristown, NJ 07960, USA ON THE SYMMETRY OF THE POWER INE CHANNE T.C. Banwe, S. Gai {bct, sgai}@research.tecordia.com Tecordia Technoogies, Inc., 445 South Street, Morristown, NJ 07960, USA Abstract The indoor power ine network

More information

THE REACHABILITY CONES OF ESSENTIALLY NONNEGATIVE MATRICES

THE REACHABILITY CONES OF ESSENTIALLY NONNEGATIVE MATRICES THE REACHABILITY CONES OF ESSENTIALLY NONNEGATIVE MATRICES by Michae Neumann Department of Mathematics, University of Connecticut, Storrs, CT 06269 3009 and Ronad J. Stern Department of Mathematics, Concordia

More information

Theory of Generalized k-difference Operator and Its Application in Number Theory

Theory of Generalized k-difference Operator and Its Application in Number Theory Internationa Journa of Mathematica Anaysis Vo. 9, 2015, no. 19, 955-964 HIKARI Ltd, www.m-hiari.com http://dx.doi.org/10.12988/ijma.2015.5389 Theory of Generaized -Difference Operator and Its Appication

More information

Determining The Degree of Generalization Using An Incremental Learning Algorithm

Determining The Degree of Generalization Using An Incremental Learning Algorithm Determining The Degree of Generaization Using An Incrementa Learning Agorithm Pabo Zegers Facutad de Ingeniería, Universidad de os Andes San Caros de Apoquindo 22, Las Condes, Santiago, Chie pzegers@uandes.c

More information

SUPPLEMENTARY MATERIAL TO INNOVATED SCALABLE EFFICIENT ESTIMATION IN ULTRA-LARGE GAUSSIAN GRAPHICAL MODELS

SUPPLEMENTARY MATERIAL TO INNOVATED SCALABLE EFFICIENT ESTIMATION IN ULTRA-LARGE GAUSSIAN GRAPHICAL MODELS ISEE 1 SUPPLEMENTARY MATERIAL TO INNOVATED SCALABLE EFFICIENT ESTIMATION IN ULTRA-LARGE GAUSSIAN GRAPHICAL MODELS By Yingying Fan and Jinchi Lv University of Southern Caifornia This Suppementary Materia

More information

Rate-Distortion Theory of Finite Point Processes

Rate-Distortion Theory of Finite Point Processes Rate-Distortion Theory of Finite Point Processes Günther Koiander, Dominic Schuhmacher, and Franz Hawatsch, Feow, IEEE Abstract We study the compression of data in the case where the usefu information

More information

with a unit root is discussed. We propose a modication of the Block Bootstrap which

with a unit root is discussed. We propose a modication of the Block Bootstrap which The Continuous-Path Bock-Bootstrap E. PAPARODITIS and D. N. POLITIS University of Cyprus and University of Caifornia, San Diego Abstract - The situation where the avaiabe data arise from a genera inear

More information

A. Distribution of the test statistic

A. Distribution of the test statistic A. Distribution of the test statistic In the sequentia test, we first compute the test statistic from a mini-batch of size m. If a decision cannot be made with this statistic, we keep increasing the mini-batch

More information

THE THREE POINT STEINER PROBLEM ON THE FLAT TORUS: THE MINIMAL LUNE CASE

THE THREE POINT STEINER PROBLEM ON THE FLAT TORUS: THE MINIMAL LUNE CASE THE THREE POINT STEINER PROBLEM ON THE FLAT TORUS: THE MINIMAL LUNE CASE KATIE L. MAY AND MELISSA A. MITCHELL Abstract. We show how to identify the minima path network connecting three fixed points on

More information

FOURIER SERIES ON ANY INTERVAL

FOURIER SERIES ON ANY INTERVAL FOURIER SERIES ON ANY INTERVAL Overview We have spent considerabe time earning how to compute Fourier series for functions that have a period of 2p on the interva (-p,p). We have aso seen how Fourier series

More information

Maximum likelihood decoding of trellis codes in fading channels with no receiver CSI is a polynomial-complexity problem

Maximum likelihood decoding of trellis codes in fading channels with no receiver CSI is a polynomial-complexity problem 1 Maximum ikeihood decoding of treis codes in fading channes with no receiver CSI is a poynomia-compexity probem Chun-Hao Hsu and Achieas Anastasopouos Eectrica Engineering and Computer Science Department

More information

Data Mining Technology for Failure Prognostic of Avionics

Data Mining Technology for Failure Prognostic of Avionics IEEE Transactions on Aerospace and Eectronic Systems. Voume 38, #, pp.388-403, 00. Data Mining Technoogy for Faiure Prognostic of Avionics V.A. Skormin, Binghamton University, Binghamton, NY, 1390, USA

More information

Two-Stage Least Squares as Minimum Distance

Two-Stage Least Squares as Minimum Distance Two-Stage Least Squares as Minimum Distance Frank Windmeijer Discussion Paper 17 / 683 7 June 2017 Department of Economics University of Bristo Priory Road Compex Bristo BS8 1TU United Kingdom Two-Stage

More information

( ) is just a function of x, with

( ) is just a function of x, with II. MULTIVARIATE CALCULUS The first ecture covered functions where a singe input goes in, and a singe output comes out. Most economic appications aren t so simpe. In most cases, a number of variabes infuence

More information

An Algorithm for Pruning Redundant Modules in Min-Max Modular Network

An Algorithm for Pruning Redundant Modules in Min-Max Modular Network An Agorithm for Pruning Redundant Modues in Min-Max Moduar Network Hui-Cheng Lian and Bao-Liang Lu Department of Computer Science and Engineering, Shanghai Jiao Tong University 1954 Hua Shan Rd., Shanghai

More information

An explicit Jordan Decomposition of Companion matrices

An explicit Jordan Decomposition of Companion matrices An expicit Jordan Decomposition of Companion matrices Fermín S V Bazán Departamento de Matemática CFM UFSC 88040-900 Forianópois SC E-mai: fermin@mtmufscbr S Gratton CERFACS 42 Av Gaspard Coriois 31057

More information

Multi-server queueing systems with multiple priority classes

Multi-server queueing systems with multiple priority classes Muti-server queueing systems with mutipe priority casses Mor Harcho-Bater Taayui Osogami Aan Scheer-Wof Adam Wierman Abstract We present the first near-exact anaysis of an M/PH/ queue with m > 2 preemptive-resume

More information

Componentwise Determination of the Interval Hull Solution for Linear Interval Parameter Systems

Componentwise Determination of the Interval Hull Solution for Linear Interval Parameter Systems Componentwise Determination of the Interva Hu Soution for Linear Interva Parameter Systems L. V. Koev Dept. of Theoretica Eectrotechnics, Facuty of Automatics, Technica University of Sofia, 1000 Sofia,

More information

MA 201: Partial Differential Equations Lecture - 10

MA 201: Partial Differential Equations Lecture - 10 MA 201: Partia Differentia Equations Lecture - 10 Separation of Variabes, One dimensiona Wave Equation Initia Boundary Vaue Probem (IBVP) Reca: A physica probem governed by a PDE may contain both boundary

More information

NOISE-INDUCED STABILIZATION OF STOCHASTIC DIFFERENTIAL EQUATIONS

NOISE-INDUCED STABILIZATION OF STOCHASTIC DIFFERENTIAL EQUATIONS NOISE-INDUCED STABILIZATION OF STOCHASTIC DIFFERENTIAL EQUATIONS TONY ALLEN, EMILY GEBHARDT, AND ADAM KLUBALL 3 ADVISOR: DR. TIFFANY KOLBA 4 Abstract. The phenomenon of noise-induced stabiization occurs

More information

Higher dimensional PDEs and multidimensional eigenvalue problems

Higher dimensional PDEs and multidimensional eigenvalue problems Higher dimensiona PEs and mutidimensiona eigenvaue probems 1 Probems with three independent variabes Consider the prototypica equations u t = u (iffusion) u tt = u (W ave) u zz = u (Lapace) where u = u

More information

Global Optimality Principles for Polynomial Optimization Problems over Box or Bivalent Constraints by Separable Polynomial Approximations

Global Optimality Principles for Polynomial Optimization Problems over Box or Bivalent Constraints by Separable Polynomial Approximations Goba Optimaity Principes for Poynomia Optimization Probems over Box or Bivaent Constraints by Separabe Poynomia Approximations V. Jeyakumar, G. Li and S. Srisatkunarajah Revised Version II: December 23,

More information

AALBORG UNIVERSITY. The distribution of communication cost for a mobile service scenario. Jesper Møller and Man Lung Yiu. R June 2009

AALBORG UNIVERSITY. The distribution of communication cost for a mobile service scenario. Jesper Møller and Man Lung Yiu. R June 2009 AALBORG UNIVERSITY The distribution of communication cost for a mobie service scenario by Jesper Møer and Man Lung Yiu R-29-11 June 29 Department of Mathematica Sciences Aaborg University Fredrik Bajers

More information

Wave Equation Dirichlet Boundary Conditions

Wave Equation Dirichlet Boundary Conditions Wave Equation Dirichet Boundary Conditions u tt x, t = c u xx x, t, < x 1 u, t =, u, t = ux, = fx u t x, = gx Look for simpe soutions in the form ux, t = XxT t Substituting into 13 and dividing

More information

Expectation-Maximization for Estimating Parameters for a Mixture of Poissons

Expectation-Maximization for Estimating Parameters for a Mixture of Poissons Expectation-Maximization for Estimating Parameters for a Mixture of Poissons Brandon Maone Department of Computer Science University of Hesini February 18, 2014 Abstract This document derives, in excrutiating

More information

A natural differential calculus on Lie bialgebras with dual of triangular type

A natural differential calculus on Lie bialgebras with dual of triangular type Centrum voor Wiskunde en Informatica REPORTRAPPORT A natura differentia cacuus on Lie biagebras with dua of trianguar type N. van den Hijigenberg and R. Martini Department of Anaysis, Agebra and Geometry

More information

LECTURE NOTES 9 TRACELESS SYMMETRIC TENSOR APPROACH TO LEGENDRE POLYNOMIALS AND SPHERICAL HARMONICS

LECTURE NOTES 9 TRACELESS SYMMETRIC TENSOR APPROACH TO LEGENDRE POLYNOMIALS AND SPHERICAL HARMONICS MASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Department Physics 8.07: Eectromagnetism II October 7, 202 Prof. Aan Guth LECTURE NOTES 9 TRACELESS SYMMETRIC TENSOR APPROACH TO LEGENDRE POLYNOMIALS AND SPHERICAL

More information

Biometrics Unit, 337 Warren Hall Cornell University, Ithaca, NY and. B. L. Raktoe

Biometrics Unit, 337 Warren Hall Cornell University, Ithaca, NY and. B. L. Raktoe NONISCMORPHIC CCMPLETE SETS OF ORTHOGONAL F-SQ.UARES, HADAMARD MATRICES, AND DECCMPOSITIONS OF A 2 4 DESIGN S. J. Schwager and w. T. Federer Biometrics Unit, 337 Warren Ha Corne University, Ithaca, NY

More information

Stochastic Variational Inference with Gradient Linearization

Stochastic Variational Inference with Gradient Linearization Stochastic Variationa Inference with Gradient Linearization Suppementa Materia Tobias Pötz * Anne S Wannenwetsch Stefan Roth Department of Computer Science, TU Darmstadt Preface In this suppementa materia,

More information

A CLUSTERING LAW FOR SOME DISCRETE ORDER STATISTICS

A CLUSTERING LAW FOR SOME DISCRETE ORDER STATISTICS J App Prob 40, 226 241 (2003) Printed in Israe Appied Probabiity Trust 2003 A CLUSTERING LAW FOR SOME DISCRETE ORDER STATISTICS SUNDER SETHURAMAN, Iowa State University Abstract Let X 1,X 2,,X n be a sequence

More information

Research Article On the Lower Bound for the Number of Real Roots of a Random Algebraic Equation

Research Article On the Lower Bound for the Number of Real Roots of a Random Algebraic Equation Appied Mathematics and Stochastic Anaysis Voume 007, Artice ID 74191, 8 pages doi:10.1155/007/74191 Research Artice On the Lower Bound for the Number of Rea Roots of a Random Agebraic Equation Takashi

More information

STA 216 Project: Spline Approach to Discrete Survival Analysis

STA 216 Project: Spline Approach to Discrete Survival Analysis : Spine Approach to Discrete Surviva Anaysis November 4, 005 1 Introduction Athough continuous surviva anaysis differs much from the discrete surviva anaysis, there is certain ink between the two modeing

More information