An Extension of the Multi-machine Multi-level Proportional Lot Sizing and Scheduling Model for Product-Life-Cycle Demand

Size: px
Start display at page:

Download "An Extension of the Multi-machine Multi-level Proportional Lot Sizing and Scheduling Model for Product-Life-Cycle Demand"

Transcription

1 International Journal of Operations Research Vol., No., (00) An Extension of the Multi-machine Multi-level Proportional Lot Sizing an Scheuling Moel for Prouct-Life-Cycle Ping-eng Chang *, Ching-Hsiang Chang, an Chia-sung Chen Department of Inustrial Engineering an Enterprise Information, Box 98, unghai University, aichung 07, aiwan, ROC Abstract Often lot sizing an scheuling moels have neglecte the consieration of influence of prouct life cycles (PLCs). But, accoringly as esirable in a firm s prouction strategy, ecisions of lot sizing scheuling may have to respon to the stage changes in a prouct s life cycle. In this paper, we propose an extension of the lot sizing scheuling, the multi-item, multi-level prouct-structure an multi-machine proportional lot sizing scheuling moel by Kimms (999), by the concept of PLC. he avantages have been realize an inclue both the setup an holing cost improvements an an improve lot sizing scheuling that better matches the proucts emans in the PLCs. In aition, the combine effects of the setup (cost) learning effect an cash flow are examine in the moel an the genetic algorithm is aopte as the solution ai tool. A numerical example is provie an emonstrates the avantages of this moel an the effects of the factors consiere. Keywors Proportional lot sizing an scheuling problem; Prouct life cycle; Cash flow, Learning curve; Genetic algorithms. INRODUCION Lot sizing an scheuling or etermining which quantities of what items have to be prouce at what time to meet the emans of the proucts is two of the most important problems in prouction planning. Many researchers have evelope moels an approaches to these problems (e.g., see Drexl an Kimms, 997 for a goo review). Four types of problems have been reporte an investigate an inclue the capacitate lot sizing problems (e.g. Barany et al., 98; Chen an hizy, 990; Haase, 99), iscrete lot sizing an scheuling problems (e.g. Fleschman, 99), continuous setup lot sizing problems (e.g. Bitran an Matsuo, 98;Karnarkar an Schrage, 98), an proportional lot sizing an scheuling problems (PLSP) (e.g. Drexl an Haase, 99; Kimms, 99). In this paper, we inten to propose an extension of the PLSP with multi-item multi-level-prouct-structure an multi-machines (PLSP-MM) (Kimms, 999) for the concern of prouct-life-cycle (PLC) emans. As the PLC concern can be an important factor, but it has been neglecte in the existing moel. In this PLSP-MM problem, we make the same assumptions as those one by Kimms (999) besies the PLC emans. We assume a finite horizon an which is ivie into iscrete perios with ynamic an eterministic external emans for the items. No shortage or backlog is allowe. Items may consist of other items that are to be prouce before the former items can be prouce. he intermeiate items may in turn cause internal emans for other items an so on. Also, we assume that the capacity per perio of each machine is constraine. he prouction of one item consumes an item-specific amount of machine capacity. o prouce an item the machine has to be set up for the item. Every setup causes item-specific setup costs. Setup time an sequence epenent setup cost are ignore. If a remaining capacity occurs in a perio it can be use to scheule another item. An ile perios between lots of a same item o not cause aitional setup costs. Items that are prouce within a perio without a eman in that perio are store in inventory to meet future emans, while inventory is assume to be uncapacitate but have holing costs. herefore, the obective is to fin the prouction plan feasible an cheap, i.e., the sum of setup an holing costs must be as low as possible. Furthermore, uner prouct life cycles (Polli an Cook, 99; Rink an Swan, 979), a firm often has to change its market strategies several times (Lambkin an Day, 989), (Levitt, 9), which is not only ue to the competitors challenges but also that the chance of the prouct s market economic state an its role as acting in the market may have change (Aaker, 99). herefore, a firm often has to prepare a series of strategies for ifferent stages of PLCs (Kurawarwala an Matsou, 998) an the strategies for prouction to meet the ifferent stages of PLCs. Existing lot sizing an scheuling moels often have focuse on the problems without the consieration of PLC stages an changes of emans in the PLCs. hough Hill (99) has propose a moel with the PLC concern, but his assumption was mae on a single item, single-prouct-level moel an cannot solve the multi-item multi-level prouct structure problems. Especially it is not consiere that the lot sizes of the proucts therefore may have to respon to * Corresponing author, ptchang@ie.thu.eu.tw.

2 IJOR Vol., No., (00) the stage changes of the PLC. It is only a pure inventory moel instea of multi-item scheuling an inventory moel. herefore, in this paper the PLSP moel will be consiere with the PLC concept. Besies that, the effects of the cash flow an learning curve (Argote an Epple, 990) will also be examine. he technique of genetic algorithms (GAs) will be aopte as the solution tool for the current moel. he remainer of this paper is organize as: Section provies the moel evelopment. Section provies the solution proceure by incorporating the genetic algorithm. Section provies a numerical example for emonstration. In section we raw conclusions.. MODEL DEVELOPMEN. he multi-level PLSP with multiple machines (PLSP-MM) In this subsection the PLSP-MM moel by Kimms (999) will be reviewe. Aitional assumptions in aition to those given in the last section will be provie as follows. In a multi-level prouct structure, the preceence relationships among the items can be efine as an acyclic gozinto-structure of a general type. Several resources (or machines) can be available for manufacturing the items an each item is prouce on item-specific machines (or resources). Also, positive lea-times of the items can be efine ue to the technological restrictions such as cooling or transportation for instance. Moreover, items may share common machines (or resources), some of them may be scarce, an capacities of these resources may vary over time. For this problem, able efines the ecision variables. able provies the problem parameters neee. Using these symbols, the PLSP-MM may be moelle as a mixe integer program (999) as follows: able. Decision variables Symbol Definition I t Inventory for item at the en of perio t. q t Prouction quantity of item in perio t. x t y t J = i= Binary variable for inicating whether a setup for item occurs in perio t x = ) or not (x t = 0). Binary variable which inicates whether machine m is set up for item at the en of perio t (y t = ) or not (y t = 0). Min ( sx + hi ), () t t t I able. Parameters for the PLSP-MM Symbol Definition J Number of items. M Number of machines. Number of perios. a i Gozinto -factor. Its value is zero if item i is not an immeiate successor of item. Otherwise, it is the quantity of item to prouce one item i. C mt Available capacity of machine m in perio t. t External eman for item in perio t. b Non-negative holing cost for having one unit of item one perio in inventory. m Machines on which items is prouce. p Capacity nees for proucing one unit of item. s Non-negative setup cost for item. v Positive an integral lea time of item. I 0 Initial inventory for item. y 0 Unique initial setup state. Set of all items that share the machine m, i.e., α m β { {,,..., } } α m = J m = m. Set of immeiate successors of item, i.e., { {,,..., } i 0} β = i J a >. { t+ v } a q min, t i iτ i β τ= t+,0 =,,, J, t =,,,, J α m y t m =,,, M, t =,,,, x y y t t ( t ) =,,, J, t =,,,, ( t m t ( t ) t ) pq C y y =,,, J, t =,,,, α m pq C t mt m =,,, M, t =,,, { }, y t 0,, i =,,, I, =,,, J, t =,,,, () () () () (7) (8) subect to I = I + q a q, t ( t ) t t i it i β =,,, J, t =,,,, () I, q, x t t t 0, =,,, J, t =,,,, i =,,, I, =,,, J, he obective () is to minimize the sum of setup an holing costs. Eqs. () are the inventory balances. At the (9)

3 IJOR Vol., No., (00) en of a perio t there shoul be in inventory what was in there at the en of perio t plus what is prouce minus external an internal emans. Constraints () guarantee a lower boun of inventory of item at the en of a perio for fulfillment of the internal emans for the amount of items require to prouce other items but within the next v perios of its positive lea time. Constraints () make sure that the setup state of each machine is uniquely efine at the en of each perio. hose perios in which setups of a machine m for item happen are spotte by Eqs. () that if the machine is not set up for the particular item at the beginning of a perio, y (t-) = 0, but it is at the en of that perio (y t = ), then a setup must have occurre in perio t, i.e., x t 0 =. Ile perios may occur in orer to save setup costs. Due to Eqs. (), prouction can take place only if there is a proper setup state either at the beginning or at the en of a particular perio. Note that if the machine is not set up for an item at the beginning of a perio but it is at the en, the setup takes place in that perio, which means that prouction may take place in that perio. As the consequence of Eqs. (), at most two items can be manufacture on each machine per perio. Capacity constraints are formulate in Eqs. (7). Eqs. (8) efine the binary-value setup state variables, while Eqs. (9) are simply non-negativity conitions. One shoul easily convince himself that ue to Eqs. () in combination with Eq. () setup variables x t are inee zero-one value. Hence non-negativity conitions are sufficient for these variables. By letting inventory variables I t be nonnegative, backlogs cannot occur. he following will introuce the present extene moel.. he propose PLC-base PLSP moel his subsection introuces the propose PLC-base PLSP-MM moel. he aitional parameters ue to this consieration will be first introuce. hen, the extension of the PLSP-MM moel will be given... Notations for the PLC concerns In applications of PLCs, important steps relevant to this research inclue etermining the transitional points between the consecutive stages, such as the introuction, growth, maturity, an ecline stages an forecasting the emans of the PLCs. Useful technique has been foun in (Chang an Chang, 000), (Chang an Chang, 00) an the reviews therein. able provies the aitional parameters to the PLC-concerne moelling ue to the concern of emans of the PLCs. Here we shoul exten the original PLSP-MM moel by extening it to allow for ifferent aitional capacity an lot sizing constraints for the various stages of PLCs. he moelling resultant may have a better scheuling for the PLC emans... Aitional capacity-allocation constraints propose In lot sizing scheuling moels, often, the capacity constraints have been treate equally for all items competing for it. No consierations have been place on the proucts eman-stage changes. As a result when the PLC stages of the proucts exist strongly, the allocation of capacities for ifferent stages of proucts may occur inaequately for fulfillment of the emans. Higher inventories have to be use in some perios. However, a reality can be observe as well that in practice a company usually can have ifferent strategies for capacity allocation an ifferent proucts in ifferent stages for any preictable market changes. herefore, if PLC is not consiere, the solutions of the conventional lot sizing scheuling moels may become unsuitable or suboptimal. o overcome this problem, here a terme stage capacity-allocation weight for each stage of the proucts base on their PLCs may be introuce. It may be efine by Stage capacity-allocation weight Stage total eman PLC total time length = PLC total eman Stage total time length, (0) able. Aitional parameters for the PLC base PLSP moelling Symbol Definition B Starting perio of the PLC of item. ransitional perio between the introuction an In growth stages of item. ransitional perio between the growth an Gr maturity stages of item. ransitional perio between the maturity an Ma ecline stages of item. E En perio of PLC of item. In Number of perios in the introuction stage of item. Gr Number of perios in the growth stage of item. Ma Number of perios in the maturity stage of item. De Number of perios in the ecline stage of item. Number of perios of the PLC of item. n Number of immeiate successors of item. an taking into consieration the PLC stage emans an lengths of the proucts in capacity allocation of the lot-sizing scheuling. herefore, by the stage capacity-allocation weights, the aitional constraints of capacity allocations for each item in each perio can be efine (will be enote as C mt) an also replace their original C m t in Eqs. (). Also, for consieration of the original constraints C mt, C mt can be efine by multiplying the original C mt by the respective stage capacity-allocation weight of the items at each stage. Hence, for introuction stage,

4 IJOR Vol., No., (00) In t B t Cmt Cmt E In t B for =,,, J, t =, +,, B B In () number of immeiate successors of an item into consieration too, to avoi shortages for the emans. he aitional require lot-sizing constraints may be efine by Aitional lot sizing constraint Stage total eman = No. of immeiate successors Stage total time length () For growth stage, Gr t In + t Cmt Cmt E Gr t B for =,,, J, t = +,, In Gr () Hence For introuction stage: In B 0 q t n, In for =,,, J, t =, +,, For growth stage: t B B In () For maturity stage, Ma t Gr + t Cmt Cmt E Ma t B for =,,, J, t = +,, Gr Ma () Gr In + 0 q t n, Gr for =,,..., J, t = +,...,, For maturity stage: t In Gr (7) For ecline stage, E t Ma + t Cmt Cmt E De t B for =,,, J, t = +,, Ma E () In Eqs. ()-(), the minimum is taken to the purpose of maintaining the original capacity constraints un-exceee. An the results fulfill the strategic reassignment of capacity for each item in each perio at each stage. herefore, the stage capacity-allocation weight will reflect the capacity assignment requirement (constraints) of the items in ifferent stages... Lot sizing constraints propose Analogous to the arguments in capacity allocation constraints, PLC emans vary with stages. he lot-sizing ecisions of the proucts shoul also be given with the consieration of the stage variations too not only the capacity availability. o supply this iea into the lot sizing ecision of ifferent stage prouction requirements of the items, aitional lot-sizing constraints may be efine. Note that an item may be require for manufacturing other item(s). For these lot-sizing constraints to be feasible, the lot sizing an thus resultant inventories shoul take the Ma Gr + 0 q t n, Ma for =,,..., J, t = +,...,, For ecline stage: E Ma + t Gr Ma 0 q t n, De for =,,..., J, t = +,...,. t Ma.. Cash flow an learning curve consierations E (8) (9) Practically time value of money is consiere necessitate in any investment, revenue an expenses. Some lot-sizing moels have investigate a cash-flow-oriente moel to replace the non-cash-flow one (e.g., see Hofmann, 998; Ram an homas, 99). he moel propose here will also examine the cash flow effects for a further investigation. Moreover, as an item s prouction usually exhibits the learning curve, which is important in PLCs, the effects of learning curve cost will also be examine in the current moel. he setup costs may be consiere with the learning curve effect. Following these results (Argote an Epple, 990; Ram an homas, 99; Yelle, 979), the setup cost

5 IJOR Vol., No., (00) for an item may be moelle in a power function: (Setup cost) S = S t t b, { t min} for =,,..., J, t S S (0) where the learning rate b = 0.98 is use in the numerical example. For the effects of cash flows, iscounting of the costs at the various perios with a rate of return i per perio to the equivalent present worth or proecting them into the future worth or equivalent annual worth performs the task (William et al., 997) he equivalent present worth is aopte here. hey inclue: (Setup costs at each perio) α m pq C t mt, m =,,, M, t =,,,, y t { 0, }, i =,,, I, =,,, J, t =,,,, 0.98 S t = S t, for =,,..., J, In t B t Cmt Cmt E In t B for =,,, J, t =, +,,, B B In J = S t x t t, for t =,,..., ( + i ) (Holing costs at each perio) J = hi t t, for t =,,..., ( + i ) () () he entire moel propose, PLC-base PLSP-MM, can be efine as, by collecting the above formulations: Min J = i= subect to s t x h I t + t t + t + ( i) ( i) I I q a q t = ( t ) + t t i it, i β =,,, J, t =,,,, I t i β { t+ v } a q min, τ =+ t =,,, J, t =,,,, J α m y t, i iτ m =,,, M, t =,,,, x y y, t t ( t ) =,,, J, t =,,,, ( t m t ( t ) t ) pq C y y =,,, J, t =,,,,,, () Gr t In + t Cmt Cmt E Gr t B for =,,, J, t = +,,, In Gr Ma t Gr + mt = mt mt E Ma t B C min C, C, for =,,, J, t = +,,, Gr Ma E t Ma + t Cmt Cmt E De t B for =,,, J, t = +,,, Gr In + Ma E t 0 q t n, Gr for =,,..., J, t = +,...,, E Ma + In Gr 0 q t n, De for =,,..., J, t = +,...,, t Ma E I, q, x 0, =,,, J, t =,,,. t t t his moel will also be enote as PLC-PLSP-MM-CL, where CL stans for the cash flow ( C ) an learning curve ( L ) effect. Also a similar enotation will be use for

6 IJOR Vol., No., (00) the original moel, PLSP-MM, an therefore, PLSP-MM-CL enotes the original PLSP-MM moel with the cash flow an learning curve effects. Also, for the comparison purpose, PLC-PLSP-MM will enote the current moel without the cash flow an learning curve effects. he following introuces the solution proceure of the moel.. HE SOLUION PROCEDURE For optimization problems, genetic algorithms have been prove powerful an general for any type of problem an particularly for those of complex nature an moelling. With the vast literature an application cases, here the reaer is referre to the references (Disney et al., 000), (Hyun et al., 998), (Ip et al., 000), (Khoua et al., 998), (Kim an Kim, 99), (Kimms, 999), (Chang an Lo, 00), (S ik )provie for the basic iea. he esign of the GA for aiing in the solution of the moel above is provie as the solution proceure as follows. Step. GA representation/coing: Let S M enote a chromosome consisting of M binary genes s mt, m =,,, M, t =,,. hese genes are arrange in the orer of machines an time perios as Figure. Each gene represents whether a set up of the corresponing machine occurs or not in the particular perio t. Step. Initial population of chromosomes: An initial population of chromosomes S mt of a esirable size (PSIZE) is generate ranomly an each S M represents a feasible solution. PSIZE = 0 has been use in the numerical example. Step. Selection rules for setting up for items: For the above moel, we have use a backwar proceure for eciing the lot size scheuling, i.e., starting from the perio an then working backwars towars perio. Meanwhile, for each s mt = (t =,,,, m =,,, M) of each SM, selecting the proper item that requires the machine to be set up has been the central of this step. Here five rules are use an applie in orer of rule (δ ) to rule (δ ). If with a selection rule there are more-than-one items that can be selecte, the next rule is applie. able gives the parameters to be use in these selection rules. Rule (δ ) Base on minimizing the holing cost: Step (δ -): Fin the set α mt of items, which require the machine m to prouce an have unfulfille eman at time t. αmt = { αm CDt ( + ) > 0} () αm Nr q τ > 0. τ =+ t able. Parameters for the selection rules Symbol Definition CD t Cumulative eman for item that has not been met yet in perio t, i.e., = ( ) CD t q an CD = 0. τ= t τ τ = t+ τ + Nr Net requirement of item, i.e., = τ = τ Nr. pat Path of an item (with the meaning to be explaine in rule ). ep Depth of an item (with the meaning to be explaine in rule ). IS k Set of items selecte by selection rule k. Step (δ -): Evaluate the holing cost that can result from each of the items in α mt. he iea is that in orer to let all immeiate successors i of an item can be prouce in the next perio (t = t + ), the prouction of item at machine m must be complete in perio t, or otherwise it will cause holing cost of the item. herefore, the holing cost that an item α mt coul cause may be evaluate as h CD (t+). his rule is to select the item that can result in the possible maximal holing cost an to scheule it first in orer to avoi high holing costs. Namely, IS δ an ( CDl( t+ ) ) max{ αmt h ( CD ( t+ ) )} h l = l αmt, (a) =, if ISδ, S mt = (b) 0, otherwise, where S mt inicates the prouction an setup status of item in perio t whether it is scheule on machine m or not. If IS δ, apply the next rule δ. Rule (δ ) Base on minimizing the setup cost: Step (δ -): he iea for this rule is that for the same machine if an item can be prouce in consecutive perios aitional setup can be avoie. For this rule, one woul like to scheule an item first if it has alreay been scheule for prouction on machine m in the next perio (t + ). Namely, this rule will select the item(s) as: resetting α mt for δ, Machine m M Perio t Gene s mt Figure. Chromosome representation.

7 IJOR Vol., No., (00) α mt = { δ } m( t ) IS S = α, (a) + mt with the four items. herefore, this rule realizes that 7 an { αm ( t+ ) t 0} αmt = CD + > αm Nr q τ > 0. τ =+ t (b) he purpose of (b) is to ensure that the item selecte has an unfulfille eman incluing that for perio t. If in (a) { ISδ S + = } = m( t ), let ISδ = IS δ, omit the next step an go to the next rule. Otherwise, continue to the next step. Step (δ -): Select the item from α mt with the minimal setup cost. if ISδ =, S mt = { l αmt sl = min { αmt s }}, (7) 0, otherwise. Again if IS δ, apply the next rule δ. Rule (δ ) Base on the maximal epth of the items in the prouct structure: Assume a prouct manufacturing structure is as shown in Figure. he epths of the items can be etermine as ep =, ep =, ep =, ep =. In orer to reuce the generating of inferior or even infeasible solutions, we woul look for the eeper-epth or eeper-associate item to scheule first. his rule realizes that if ISδ =, S mt = { l ISδ ep = max { δ } l IS ep }, (8) 0, otherwise. Again if IS δ, apply the next rule δ. Figure. An example of prouct structure. Rule (δ ) Base on the longest path of the items: By a similar iea to rule, rule focuses on the entire path of the items in the prouct structure for avoiing inferior solutions ue to long influence of an item to the inventory an setup costs if it is not scheule first. Again by referring to Figure, the path lengths of the items as shown for instance can be etermine as pat = pat = pat = pat = as they are all in a path if ISδ =, S mt = { l ISδ pat = max{ δ } l IS pat } (9) 0, otherwise. If again IS δ, apply the final rule δ. Rule (δ ) Ranom selection: If applying the above rules still cannot settle the choice, a ranom selection may be applie. Step. Carry out the program () an etermine all other ecision variables values an the obective function value. Step. GA reprouction/selection: he GA reprouction is esigne an performe here with the roulette wheel selection base on the fitness function values. In GA operation, this requires the evaluation of the fitness (obective) function value of each chromosome an which is for maximization. he fitness function is efine as the reciprocal of the obective function above for each chromosome. Step. GA crossover: he crossover operation is esigne with a crossover probability P c for actually performing the crossover operation or otherwise the two ranomly selecte chromosomes enter irectly the next generation. It is performe in two ranomly selecte chromosomes with two ranomly selecte gene-points. hen exchange of the genes between the crossover points of the two chromosomes generates two offspring chromosomes. Step 7. GA mutation: Mutation operation is carrie out with a mutation probability P m on a ranomly selecte gene of a chromosome ranomly selecte here an mutates it from 0 to or vice versa. After completing this operation, a generation of the GA is reache. hen, if the termination conition is met, the GA process terminates. Or otherwise it returns to step. In aition, in the above proceure we also esigne it with a floating crossover an mutation probability mechanism for escaping from entrapment of local optimum possible. his is with that, when the process continues to fin the same value for a number of generations n, the GA increases graually the probabilities P c an P m with (P cu P cd )/GEN an (P mu P md )/GEN, respectively, where P cu an P cd enote the upper limit an efault, respectively, of P c an similarly for P m, an GEN enotes the number of generations. ill the solution changes, they are then reset to the efault value. If the number of generations GEN is reache an still the same solution is foun, the process terminates. After a number of empirical trial runs with all moels, these probability parameters, (P cu P cd )/GEN = (0.9 0.)/00, (P mu P md )/GEN = (0. 0.)/00, GEN = 00, an n = 0 were use.. NUMERICAL EXAMPLE Assume a firm is planning for prouction of two

8 IJOR Vol., No., (00) 8 proucts X or an Y or for the next ten perios for emans as shown in able. Figure. shows the manufacturing structure of the two proucts accoring to the bills of materials, where the numbers along the arrows inicate the quantities of the items require to make a unit of another item. able shows the manufacturing parameters of all these items. here are three machines available an each machine has 80 units of capacity per perio of time. able 7 shows the manufacturing relationships between these machines an items. It shows that which item can be prouce on what machine(s). For example, item can be prouce on machine or. able 8 shows the optimal chromosomes of machine setups S M with each moel introuce after 0 runs of the solution proceures. Other results obtaine by these moels are provie an analyze in the next sections. able. s for proucts X an Y Prouct Perio X Y Figure. he gozinto-structure for manufacturing. able. Manufacturing parameters Parameter Item S h P V I y able 7. Relationships between machines an items Relationship between the machine an item Machine Item 0 a 0 0 a a inicates that an item can be prouce on a machine an 0 inicates that it cannot.. he optimal lot sizing an scheuling he outputs are from the PLC-PLSP-MM-CL moel an the original PLSP-MM, PLC-PLSP-MM, an PLSP-MM-CL moels too. able 9 shows the optimal setup statuses of the machines for the items at the ens of perios, y t, with the applie selection rule by the PLC-PLSP-MM-CL moel. able 8. Optimal chromosomes of machine setups S mt with each moel S mt Machine Original PLSP-MM PLC-PLSP-MM PLSP-MM-CL PLC-PLSP-MM-CL able 9. Optimal setup statuses y t by the PLC-PLSP-MM-CL moel y t Machine Item t = (δ )(δ )(δ )(δ ) (δ ) (δ ) (δ ) (δ ) (δ )(δ ) (δ ) (δ ) 0 (δ ) (δ ) (δ ) (δ ) (δ ) (δ ) (δ ) (δ ) (δ ) (δ ) (δ )(δ ) 0 0 (δ ) 0 (δ ) (δ ) (δ ) (δ ) (δ )(δ ) (δ ) 0 (δ ) (δ ) 0 (δ )(δ ) (δ ) (δ ) (δ ) (δ )(δ )(δ ) (δ ) (δ ) 0 (δ ) (δ )(δ )(δ ) herefore, able 0 shows the optimal machine setups that are to be occurring in the ifferent perios for items, x t of the PLC-PLSP-MM-CL. Furthermore, the optimal lot sizing q t an optimal inventory ecisions I t of the various items by the PLC-PLSP-MM-CL moel are shown in ables an. he results of the original PLSP-MM, PLC-PLSP-MM, an PLSP-MM-CL moels are also obtaine; but for space limitations, they will not be provie here in etail. But, Figures -7 provie the final lot sizing an scheuling of the four moels in respective Gantt charts. In the next subsections, comparisons of these moels shall be mae.

9 IJOR Vol., No., (00) 9 able 0. Optimal machine setups to occur in ifferent perios for the items, x t, of the PLC-PLSP-MM-CL x t Machine Item t = able. Optimal lot-sizing ecisions q t of the items by the PLC-PLSP-MM-CL moel q t Item t = able. Optimal inventory ecisions I t of the items by the PLC-PLSP-MM-CL I t Item t = Machine Perio : Unuse Capacity : Unuse Perio Figure. Optimal lot sizing an scheuling by the PLC-PLSP-MM-CL.

10 IJOR Vol., No., (00) 0 Machine : Unuse Capacity Unuse Perio Figure. Lot sizing an scheuling by the original PLSP-MM. Perio Machine Perio : Unuse Capacity : Unuse Perio Figure. Optimal lot sizing an scheuling by the PLC-PLSP-MM. Machine Perio : Unuse Capacity : Unuse Perio Figure 7. Optimal lot sizing an scheuling by the PLSP-MM-CL.

11 IJOR Vol., No., (00). Comparison of the lot sizing scheuling costs of the moels able shows the setup cost, holing cost, an their total resulting with each moel. he following important remarks can be mae. () he effect of PLC concern: Examining the two pairs of moel s costs, original PLSP-MM moel with the PLC- PLSP-MM an PLSP-MM-CL with the PLC-PLSP-MM- CL moel (able ) inicates that both pairs iffer in whether the PLC is consiere or not an that not only the holing cost but also the setup cost are improve by the PLC consierations. Particularly a greater saving can be even given in the holing cost. he pairs of comparison of moels costs are summarize in able. () he combine effect of cash flow an setup-cost learning: Although the two moels, original PLSP-MM moel an PLSP-MM-CL (an similarly PLC-PLSP-MM moel an PLC-PLSP-MM-CL) may not be irectly compare ue to the reason of ifferent cost bases, one shoul still observe that the improvement of the total cost ue to the PLC consieration combine with the cash flow an setup cost learning (i.e., PLSP-MM-CL PLC-PLSP- MM-CL) is greater than that with only the PLC factor (i.e., original PLSP-MM PLC-PLSP-MM) in able. his inicates that the CL has also an effect on the moel an therefore shoul be taken into consieration in the moel. able. Optimal cost of the moels Moel Setup Holing Cost Cost otal Cost Original PLSP-MM PLC-PLSP-MM 00 8 PLSP-MM-CL PLC-PLSP-MM-CL 7 7 able. Comparison of the optimal cost improvement of the moels Reuction Setup Cost HolingCost otal Cost Original PLSP-MM PLC-PLSP-MM (-.%) (-.%) (-7.0%) PLSP-MM-CL PLC-PLSP-MM-CL 8 7 (-.%) 8 7 (-0.%) 8 (-.%). Comparison of eviations between scheule outputs of the moels an emans his subsection further compares the eviations between the scheule output an items emans by the moels. his task can be one either in a single-perio-by-single-perio manner or by a cumulate manner. For space limitations, only the cumulate results are provie here. In able, first the emans of the items in each perio are summarize. herefore, the cumulate emans an cumulate scheule outputs for all items (in units of capacity) up to each perio by the four moels are shown in able. he sum of square ifferences, SSD, between the cumulate emans an outputs of each moel can be shown in able 7. Similar remarks to that in Section. can be mae. he SSDs by the PLC-PLSP-MM-CL moel an PLC-PLSP-MM moel are much smaller than that by the PLSP-MM-CL an PLSP-MM, respectively. Furthermore, although the SSD by the PLSP-MM-CL is higher than that by the PLSP-MM, the SSD by the PLC-PLSP-MM-CL is smaller than that by PLC-PLSP-MM. So, again we observe that the PLC-PLSP- MM moel prouces a lot-sizing scheule that is better matching the proucts emans in the PLCs than that by the PLSP-MM moel. An the combine effect of the cash flow an setup cost learning curve (i.e., the PLC- PLSP-MM-CL moel) again reinforces this result. he extension of the PLSP-MM by the PLC-PLSP-MM-CL moel provies benefits to the PLSP problem with the prouct-life-cycle consierations. able. s of the items for each perio Item t = able. Cumulate emans an scheule outputs of the moels Cumulate s Cumulate Outputs: Original PLSP-MM PLC-PLSP-MM PLSP-MM-CL PLC-PLSP-MM-CL able 7. Sum of square ifferences (SSD) between cumulate emans an scheule outputs by each moel Original PLC-PLSP-MM PLSP-MM-CL PLC-PLSP-MM-CL PLSP-MM SSD CONCLUSIONS his paper has presente an extene PLSP-MM moel base on the prouct life cycle concept. It has introuce the proper PLC prouction capacity-allocation constraints an lot sizing constraints for the proucts for ifferent stages as the aitional constraints. he results of the example showe that the benefits of the extene moel provie not only reuce both the setup an holing costs but also that the lot sizing scheuling was significantly much better matching the proucts emans in the PLCs than the PLSP-MM moel. Also, in the moel

12 IJOR Vol., No., (00) the setup cost learning an cash flow were also consiere. he results also showe that these factors shoul be inclue in the moel. Further research from this research may apply a similar or re-evelope approach to the moels of other lot sizing an scheuling problems such as that reviewe in this paper or (Drexl an Kimms, 997). ACKNOWLEDGEMENS his research was partially supporte by the National Science Council of aiwan, ROC uner NSC 90-8-E REFERENCES. Aaker, D.A. (99). Strategic Market Management. Wiley, New York.. Argote, L. an Epple, D. (990). Learning curve in manufacturing. Science, 7: Barany, I., Van Roy,. J., an Wolsey, L. A. (98). Strong formulations for multi-item capacitate lot sizing. Management Science, 0(0):-.. Bitran, G.R. an Matsuo, H. (98). Approximation formulations for the single-prouct capacitate lot sizing problem. Operations Research, (): -7.. Chang, P.-. an Chang, C.-H.. (Dec., 000). A prouct life cycle moelling with stage characteristic-preserving capability, th Annual International Conference of Inustrial Engineering heory, Applications & Practice, aiwan, pp.-.. Chang, P.-. an Chang, C.-H. (00). A stage characteristic-preserving prouct life cycle moeling. Mathematical an Computer Moelling, 7: Chang, P.-. an Lo, Y.-. (00).Moelling of ob-shop scheuling with multiple quantitative an qualitative obectives an a GA/S mixture approach. International Journal of Computer Integrate Manufacturing, : Chen, W.H. an hizy, J.M.. (990). Analysis of relaxations for the multi-item capacitate lot-sizing problem. Annals of Operations Research, : Disney, S.M., Naim, M. N., an owill, D. R.. (000). Genetic algorithm optimisation of a class of inventory control systems. International Journal of Prouction Economics, 8: Drexl, A. an Haase, K. (99). Proportional Lotsizing an Scheuling. International Journal of Prouction Economics, 0(): Drexl, A. an Kimms, A. (997). Lot sizing an scheuling Survey an extensions. European Journal of Operational Research, 99: -.. Fleischman, B. (99). he iscrete lot sizing an scheuling problem with sequence-epenent setup costs. European Journal of Operational Research, 7: Haase, K. (99). Capacitate lot-sizing with sequence epenent setup costs. OR Spektrum, B. 8: -9.. Hill, R. M. (99). Batching policies for a prouct life cycle. International Journal of Prouction Economics, : -7.. Hofmann, C. (998). Investments in moern prouction technology an the cash flow-oriente EPQ moel. International Journal of Prouction Economics, (): Hyun, C.J., Kim, Y., an Kim, Y. K. (998). A genetic algorithm for multiple obective sequencing problems in mixe moel assembly lines. Computers an Operations Research, : W. H., Li, Y., Man, K. F., an ang, K. S. (000). Multi-prouct planning an scheuling using genetic algorithm approach. Computers an Inustrial Engineering, 8: Karnarkar, U.S. an Schrage, L. (98). he eterministic ynamic prouct cycling problem. Operations Research, : Khoua, M., Michalewicz, Z., an Wilmot, M.. (998). he use of genetic algorithms to solve the economic lot size scheuling problem. European Journal of Operational Research, 0: Kim, J. an Kim, Y. (99). Simulate annealing an genetic algorithms for scheuling proucts with multi-level prouct structure. Computers an Operations Research, : Kimms, A. (99). Multi-level, single-machine lot sizing an scheuling (with initial inventory). European Journal of Operational Research, 89: Kimms, A. (999). A genetic algorithm for multi-level, multi-machine lot sizing an scheuling. Computers an Operations Research, : Kurawarwala, A.A. an Matsou, H. (998). Prouct growth moels for meium-term forecasting of short-life-cycle proucts. echnological Forecasting an Social Changes, 7: Lambkin, M. an Day, G.S. (989). Evolutionary process in competitive markets: Beyon the prouct life cycle. Journal of Marketing, : -0.. Levitt,. (9). Exploit the prouct life cycle. Harvar Business Review, : Polli, R. an Cook, V. (99). Valiity of the prouct life cycle. Journal of Business, : Ram, R. an homas, J. S.. (99). Optimal an heuristic policies for lot sizing with learning in setups, Journal of Operations Management, : Rink, D.R. an Swan, J.E. (979). Prouct life cycle research: A literature review, Journal of Business Research, 7: Sikora, R. (99). A genetic algorithm for integrating lot-sizing an sequencing in scheuling a capacitate flow line. Computers an Inustrial Engineering, 0(): William, G.S., James, A.B., an Elin, M.W. (997). Engineering Economy. Prentice Hall.. Yelle, L. E. (979). he learning curve: Historical review an comprehensive survey. Decision Sciences, 0: 0-.

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 CONTROL OF A PRODUCTION SYSTEM WITH INVENTORY-LEVEL-DEPENDENT DEMAND

OPTIMAL CONTROL OF A PRODUCTION SYSTEM WITH INVENTORY-LEVEL-DEPENDENT DEMAND Applie Mathematics E-Notes, 5(005), 36-43 c ISSN 1607-510 Available free at mirror sites of http://www.math.nthu.eu.tw/ amen/ OPTIMAL CONTROL OF A PRODUCTION SYSTEM WITH INVENTORY-LEVEL-DEPENDENT DEMAND

More information

Scott E. Grasman 1, Zaki Sari 2 and Tewfik Sari 3

Scott E. Grasman 1, Zaki Sari 2 and Tewfik Sari 3 RAIRO Operations Research RAIRO Oper. Res. 41 (27) 455 464 DOI: 1.151/ro:2731 NEWSVENDOR SOLUTIONS WITH GENERAL RANDOM YIELD DISTRIBUTIONS Scott E. Grasman 1, Zaki Sari 2 an Tewfik Sari 3 Abstract. Most

More information

CONTROL CHARTS FOR VARIABLES

CONTROL CHARTS FOR VARIABLES UNIT CONTOL CHATS FO VAIABLES Structure.1 Introuction Objectives. Control Chart Technique.3 Control Charts for Variables.4 Control Chart for Mean(-Chart).5 ange Chart (-Chart).6 Stanar Deviation Chart

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

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

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

Designing of Acceptance Double Sampling Plan for Life Test Based on Percentiles of Exponentiated Rayleigh Distribution

Designing of Acceptance Double Sampling Plan for Life Test Based on Percentiles of Exponentiated Rayleigh Distribution International Journal of Statistics an Systems ISSN 973-675 Volume, Number 3 (7), pp. 475-484 Research Inia Publications http://www.ripublication.com Designing of Acceptance Double Sampling Plan for Life

More information

A note on asymptotic formulae for one-dimensional network flow problems Carlos F. Daganzo and Karen R. Smilowitz

A note on asymptotic formulae for one-dimensional network flow problems Carlos F. Daganzo and Karen R. Smilowitz A note on asymptotic formulae for one-imensional network flow problems Carlos F. Daganzo an Karen R. Smilowitz (to appear in Annals of Operations Research) Abstract This note evelops asymptotic formulae

More information

Monopoly Part III: Multiple Local Equilibria and Adaptive Search

Monopoly Part III: Multiple Local Equilibria and Adaptive Search FH-Kiel University of Applie Sciences Prof Dr Anreas Thiemer, 00 e-mail: anreasthiemer@fh-kiele Monopoly Part III: Multiple Local Equilibria an Aaptive Search Summary: The information about market eman

More information

Expected Value of Partial Perfect Information

Expected Value of Partial Perfect Information Expecte Value of Partial Perfect Information Mike Giles 1, Takashi Goa 2, Howar Thom 3 Wei Fang 1, Zhenru Wang 1 1 Mathematical Institute, University of Oxfor 2 School of Engineering, University of Tokyo

More information

Situation awareness of power system based on static voltage security region

Situation awareness of power system based on static voltage security region The 6th International Conference on Renewable Power Generation (RPG) 19 20 October 2017 Situation awareness of power system base on static voltage security region Fei Xiao, Zi-Qing Jiang, Qian Ai, Ran

More information

LATTICE-BASED D-OPTIMUM DESIGN FOR FOURIER REGRESSION

LATTICE-BASED D-OPTIMUM DESIGN FOR FOURIER REGRESSION The Annals of Statistics 1997, Vol. 25, No. 6, 2313 2327 LATTICE-BASED D-OPTIMUM DESIGN FOR FOURIER REGRESSION By Eva Riccomagno, 1 Rainer Schwabe 2 an Henry P. Wynn 1 University of Warwick, Technische

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

A simple model for the small-strain behaviour of soils

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

More information

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

Cascaded redundancy reduction

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

More information

Balancing Expected and Worst-Case Utility in Contracting Models with Asymmetric Information and Pooling

Balancing Expected and Worst-Case Utility in Contracting Models with Asymmetric Information and Pooling Balancing Expecte an Worst-Case Utility in Contracting Moels with Asymmetric Information an Pooling R.B.O. erkkamp & W. van en Heuvel & A.P.M. Wagelmans Econometric Institute Report EI2018-01 9th January

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

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

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

Developing a Method for Increasing Accuracy and Precision in Measurement System Analysis: A Fuzzy Approach

Developing a Method for Increasing Accuracy and Precision in Measurement System Analysis: A Fuzzy Approach Journal of Inustrial Engineering 6(00 5-3 Developing a Metho for Increasing Accuracy an Precision in Measurement System Analysis: A Fuzzy Approach Abolfazl Kazemi a, Hassan Haleh a, Vahi Hajipour b, Seye

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

Lyapunov Functions. V. J. Venkataramanan and Xiaojun Lin. Center for Wireless Systems and Applications. School of Electrical and Computer Engineering,

Lyapunov Functions. V. J. Venkataramanan and Xiaojun Lin. Center for Wireless Systems and Applications. School of Electrical and Computer Engineering, On the Queue-Overflow Probability of Wireless Systems : A New Approach Combining Large Deviations with Lyapunov Functions V. J. Venkataramanan an Xiaojun Lin Center for Wireless Systems an Applications

More information

This module is part of the. Memobust Handbook. on Methodology of Modern Business Statistics

This module is part of the. Memobust Handbook. on Methodology of Modern Business Statistics This moule is part of the Memobust Hanbook on Methoology of Moern Business Statistics 26 March 2014 Metho: Balance Sampling for Multi-Way Stratification Contents General section... 3 1. Summary... 3 2.

More information

1 The Derivative of ln(x)

1 The Derivative of ln(x) Monay, December 3, 2007 The Derivative of ln() 1 The Derivative of ln() The first term or semester of most calculus courses will inclue the it efinition of the erivative an will work out, long han, a number

More information

A Review of Multiple Try MCMC algorithms for Signal Processing

A Review of Multiple Try MCMC algorithms for Signal Processing A Review of Multiple Try MCMC algorithms for Signal Processing Luca Martino Image Processing Lab., Universitat e València (Spain) Universia Carlos III e Mari, Leganes (Spain) Abstract Many applications

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

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

Online Appendix for Trade Policy under Monopolistic Competition with Firm Selection

Online Appendix for Trade Policy under Monopolistic Competition with Firm Selection Online Appenix for Trae Policy uner Monopolistic Competition with Firm Selection Kyle Bagwell Stanfor University an NBER Seung Hoon Lee Georgia Institute of Technology September 6, 2018 In this Online

More information

How to Minimize Maximum Regret in Repeated Decision-Making

How to Minimize Maximum Regret in Repeated Decision-Making How to Minimize Maximum Regret in Repeate Decision-Making Karl H. Schlag July 3 2003 Economics Department, European University Institute, Via ella Piazzuola 43, 033 Florence, Italy, Tel: 0039-0-4689, email:

More information

Short Intro to Coordinate Transformation

Short Intro to Coordinate Transformation Short Intro to Coorinate Transformation 1 A Vector A vector can basically be seen as an arrow in space pointing in a specific irection with a specific length. The following problem arises: How o we represent

More information

A simplified macroscopic urban traffic network model for model-based predictive control

A simplified macroscopic urban traffic network model for model-based predictive control Delft University of Technology Delft Center for Systems an Control Technical report 9-28 A simplifie macroscopic urban traffic network moel for moel-base preictive control S. Lin, B. De Schutter, Y. Xi,

More information

Agmon Kolmogorov Inequalities on l 2 (Z d )

Agmon Kolmogorov Inequalities on l 2 (Z d ) Journal of Mathematics Research; Vol. 6, No. ; 04 ISSN 96-9795 E-ISSN 96-9809 Publishe by Canaian Center of Science an Eucation Agmon Kolmogorov Inequalities on l (Z ) Arman Sahovic Mathematics Department,

More information

CMSC 313 Preview Slides

CMSC 313 Preview Slides CMSC 33 Preview Slies These are raft slies. The actual slies presente in lecture may be ifferent ue to last minute changes, scheule slippage,... UMBC, CMSC33, Richar Chang CMSC 33 Lecture

More information

Lecture 6: Calculus. In Song Kim. September 7, 2011

Lecture 6: Calculus. In Song Kim. September 7, 2011 Lecture 6: Calculus In Song Kim September 7, 20 Introuction to Differential Calculus In our previous lecture we came up with several ways to analyze functions. We saw previously that the slope of a linear

More information

Web Appendix to Firm Heterogeneity and Aggregate Welfare (Not for Publication)

Web Appendix to Firm Heterogeneity and Aggregate Welfare (Not for Publication) Web ppeni to Firm Heterogeneity an ggregate Welfare Not for Publication Marc J. Melitz Harvar University, NBER, an CEPR Stephen J. Reing Princeton University, NBER, an CEPR March 6, 203 Introuction his

More information

Text S1: Simulation models and detailed method for early warning signal calculation

Text S1: Simulation models and detailed method for early warning signal calculation 1 Text S1: Simulation moels an etaile metho for early warning signal calculation Steven J. Lae, Thilo Gross Max Planck Institute for the Physics of Complex Systems, Nöthnitzer Str. 38, 01187 Dresen, Germany

More information

Capacity Analysis of MIMO Systems with Unknown Channel State Information

Capacity Analysis of MIMO Systems with Unknown Channel State Information Capacity Analysis of MIMO Systems with Unknown Channel State Information Jun Zheng an Bhaskar D. Rao Dept. of Electrical an Computer Engineering University of California at San Diego e-mail: juzheng@ucs.eu,

More information

Optimal Variable-Structure Control Tracking of Spacecraft Maneuvers

Optimal Variable-Structure Control Tracking of Spacecraft Maneuvers Optimal Variable-Structure Control racking of Spacecraft Maneuvers John L. Crassiis 1 Srinivas R. Vaali F. Lanis Markley 3 Introuction In recent years, much effort has been evote to the close-loop esign

More information

Predictive Control of a Laboratory Time Delay Process Experiment

Predictive Control of a Laboratory Time Delay Process Experiment Print ISSN:3 6; Online ISSN: 367-5357 DOI:0478/itc-03-0005 Preictive Control of a aboratory ime Delay Process Experiment S Enev Key Wors: Moel preictive control; time elay process; experimental results

More information

3.2 Shot peening - modeling 3 PROCEEDINGS

3.2 Shot peening - modeling 3 PROCEEDINGS 3.2 Shot peening - moeling 3 PROCEEDINGS Computer assiste coverage simulation François-Xavier Abaie a, b a FROHN, Germany, fx.abaie@frohn.com. b PEENING ACCESSORIES, Switzerlan, info@peening.ch Keywors:

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

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

Perfect Matchings in Õ(n1.5 ) Time in Regular Bipartite Graphs

Perfect Matchings in Õ(n1.5 ) Time in Regular Bipartite Graphs Perfect Matchings in Õ(n1.5 ) Time in Regular Bipartite Graphs Ashish Goel Michael Kapralov Sanjeev Khanna Abstract We consier the well-stuie problem of fining a perfect matching in -regular bipartite

More information

Power Generation and Distribution via Distributed Coordination Control

Power Generation and Distribution via Distributed Coordination Control Power Generation an Distribution via Distribute Coorination Control Byeong-Yeon Kim, Kwang-Kyo Oh, an Hyo-Sung Ahn arxiv:407.4870v [math.oc] 8 Jul 204 Abstract This paper presents power coorination, power

More information

Cost-based Heuristics and Node Re-Expansions Across the Phase Transition

Cost-based Heuristics and Node Re-Expansions Across the Phase Transition Cost-base Heuristics an Noe Re-Expansions Across the Phase Transition Elan Cohen an J. Christopher Beck Department of Mechanical & Inustrial Engineering University of Toronto Toronto, Canaa {ecohen, jcb}@mie.utoronto.ca

More information

Lecture 22. Lecturer: Michel X. Goemans Scribe: Alantha Newman (2004), Ankur Moitra (2009)

Lecture 22. Lecturer: Michel X. Goemans Scribe: Alantha Newman (2004), Ankur Moitra (2009) 8.438 Avance Combinatorial Optimization Lecture Lecturer: Michel X. Goemans Scribe: Alantha Newman (004), Ankur Moitra (009) MultiFlows an Disjoint Paths Here we will survey a number of variants of isjoint

More information

TEMPORAL AND TIME-FREQUENCY CORRELATION-BASED BLIND SOURCE SEPARATION METHODS. Yannick DEVILLE

TEMPORAL AND TIME-FREQUENCY CORRELATION-BASED BLIND SOURCE SEPARATION METHODS. Yannick DEVILLE TEMPORAL AND TIME-FREQUENCY CORRELATION-BASED BLIND SOURCE SEPARATION METHODS Yannick DEVILLE Université Paul Sabatier Laboratoire Acoustique, Métrologie, Instrumentation Bât. 3RB2, 8 Route e Narbonne,

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

Acute sets in Euclidean spaces

Acute sets in Euclidean spaces Acute sets in Eucliean spaces Viktor Harangi April, 011 Abstract A finite set H in R is calle an acute set if any angle etermine by three points of H is acute. We examine the maximal carinality α() of

More information

Lecture Introduction. 2 Examples of Measure Concentration. 3 The Johnson-Lindenstrauss Lemma. CS-621 Theory Gems November 28, 2012

Lecture Introduction. 2 Examples of Measure Concentration. 3 The Johnson-Lindenstrauss Lemma. CS-621 Theory Gems November 28, 2012 CS-6 Theory Gems November 8, 0 Lecture Lecturer: Alesaner Mąry Scribes: Alhussein Fawzi, Dorina Thanou Introuction Toay, we will briefly iscuss an important technique in probability theory measure concentration

More information

Leaving Randomness to Nature: d-dimensional Product Codes through the lens of Generalized-LDPC codes

Leaving Randomness to Nature: d-dimensional Product Codes through the lens of Generalized-LDPC codes Leaving Ranomness to Nature: -Dimensional Prouct Coes through the lens of Generalize-LDPC coes Tavor Baharav, Kannan Ramchanran Dept. of Electrical Engineering an Computer Sciences, U.C. Berkeley {tavorb,

More information

Constraint Reformulation and a Lagrangian Relaxation based Solution Algorithm for a Least Expected Time Path Problem Abstract 1.

Constraint Reformulation and a Lagrangian Relaxation based Solution Algorithm for a Least Expected Time Path Problem Abstract 1. Constraint Reformulation an a Lagrangian Relaation base Solution Algorithm for a Least Epecte Time Path Problem Liing Yang State Key Laboratory of Rail Traffic Control an Safety, Being Jiaotong University,

More information

Slovak University of Technology in Bratislava Institute of Information Engineering, Automation, and Mathematics PROCEEDINGS

Slovak University of Technology in Bratislava Institute of Information Engineering, Automation, and Mathematics PROCEEDINGS lovak University of echnology in Bratislava Institute of Information Engineering, Automation, an athematics PROCEEDING 7 th International Conference on Process Control 009 Hotel Baník, Štrbské Pleso, lovakia,

More information

arxiv: v1 [math.co] 29 May 2009

arxiv: v1 [math.co] 29 May 2009 arxiv:0905.4913v1 [math.co] 29 May 2009 simple Havel-Hakimi type algorithm to realize graphical egree sequences of irecte graphs Péter L. Erős an István Miklós. Rényi Institute of Mathematics, Hungarian

More information

We G Model Reduction Approaches for Solution of Wave Equations for Multiple Frequencies

We G Model Reduction Approaches for Solution of Wave Equations for Multiple Frequencies We G15 5 Moel Reuction Approaches for Solution of Wave Equations for Multiple Frequencies M.Y. Zaslavsky (Schlumberger-Doll Research Center), R.F. Remis* (Delft University) & V.L. Druskin (Schlumberger-Doll

More information

Total Energy Shaping of a Class of Underactuated Port-Hamiltonian Systems using a New Set of Closed-Loop Potential Shape Variables*

Total Energy Shaping of a Class of Underactuated Port-Hamiltonian Systems using a New Set of Closed-Loop Potential Shape Variables* 51st IEEE Conference on Decision an Control December 1-13 212. Maui Hawaii USA Total Energy Shaping of a Class of Uneractuate Port-Hamiltonian Systems using a New Set of Close-Loop Potential Shape Variables*

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

arxiv:hep-th/ v1 3 Feb 1993

arxiv:hep-th/ v1 3 Feb 1993 NBI-HE-9-89 PAR LPTHE 9-49 FTUAM 9-44 November 99 Matrix moel calculations beyon the spherical limit arxiv:hep-th/93004v 3 Feb 993 J. Ambjørn The Niels Bohr Institute Blegamsvej 7, DK-00 Copenhagen Ø,

More information

Modeling the effects of polydispersity on the viscosity of noncolloidal hard sphere suspensions. Paul M. Mwasame, Norman J. Wagner, Antony N.

Modeling the effects of polydispersity on the viscosity of noncolloidal hard sphere suspensions. Paul M. Mwasame, Norman J. Wagner, Antony N. Submitte to the Journal of Rheology Moeling the effects of polyispersity on the viscosity of noncolloial har sphere suspensions Paul M. Mwasame, Norman J. Wagner, Antony N. Beris a) epartment of Chemical

More information

Optimistic Planning with Long Sequences of Identical Actions: An Extended Theoretical and Experimental Study

Optimistic Planning with Long Sequences of Identical Actions: An Extended Theoretical and Experimental Study Volume 56, Number 1-2, 2015 27 Optimistic Planning with Long Sequences of Ientical Actions: An Extene Theoretical an Experimental Stuy Koppány Máthé, Lucian Bușoniu, Liviu Miclea epartment of Automation,

More information

A Model of Discovery

A Model of Discovery Moel of Discovery Michele olrin an Davi K. Levine 1 his version: December 26, 2008 First version: December 13, 2008 Prepare for the 2008 merican Economic ssociation Meetings. Session itle: Intellectual

More information

A Novel Decoupled Iterative Method for Deep-Submicron MOSFET RF Circuit Simulation

A Novel Decoupled Iterative Method for Deep-Submicron MOSFET RF Circuit Simulation A Novel ecouple Iterative Metho for eep-submicron MOSFET RF Circuit Simulation CHUAN-SHENG WANG an YIMING LI epartment of Mathematics, National Tsing Hua University, National Nano evice Laboratories, an

More information

An Anisotropic Hardening Model for Springback Prediction

An Anisotropic Hardening Model for Springback Prediction An Anisotropic Harening Moel for Springback Preiction Danielle Zeng an Z. Ceric Xia Scientific Research Laboratories For Motor Company Dearborn, MI 48 Abstract. As more Avance High-Strength Steels (AHSS

More information

Experiment 2, Physics 2BL

Experiment 2, Physics 2BL Experiment 2, Physics 2BL Deuction of Mass Distributions. Last Upate: 2009-05-03 Preparation Before this experiment, we recommen you review or familiarize yourself with the following: Chapters 4-6 in Taylor

More information

The canonical controllers and regular interconnection

The canonical controllers and regular interconnection Systems & Control Letters ( www.elsevier.com/locate/sysconle The canonical controllers an regular interconnection A.A. Julius a,, J.C. Willems b, M.N. Belur c, H.L. Trentelman a Department of Applie Mathematics,

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

Non-deterministic Social Laws

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

More information

Markov Chains in Continuous Time

Markov Chains in Continuous Time Chapter 23 Markov Chains in Continuous Time Previously we looke at Markov chains, where the transitions betweenstatesoccurreatspecifietime- steps. That it, we mae time (a continuous variable) avance in

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

State-Space Model for a Multi-Machine System

State-Space Model for a Multi-Machine System State-Space Moel for a Multi-Machine System These notes parallel section.4 in the text. We are ealing with classically moele machines (IEEE Type.), constant impeance loas, an a network reuce to its internal

More information

TIME-DELAY ESTIMATION USING FARROW-BASED FRACTIONAL-DELAY FIR FILTERS: FILTER APPROXIMATION VS. ESTIMATION ERRORS

TIME-DELAY ESTIMATION USING FARROW-BASED FRACTIONAL-DELAY FIR FILTERS: FILTER APPROXIMATION VS. ESTIMATION ERRORS TIME-DEAY ESTIMATION USING FARROW-BASED FRACTIONA-DEAY FIR FITERS: FITER APPROXIMATION VS. ESTIMATION ERRORS Mattias Olsson, Håkan Johansson, an Per öwenborg Div. of Electronic Systems, Dept. of Electrical

More information

Logarithmic spurious regressions

Logarithmic spurious regressions Logarithmic spurious regressions Robert M. e Jong Michigan State University February 5, 22 Abstract Spurious regressions, i.e. regressions in which an integrate process is regresse on another integrate

More information

The Uniqueness of the Overall Assurance Interval for Epsilon in DEA Models by the Direction Method

The Uniqueness of the Overall Assurance Interval for Epsilon in DEA Models by the Direction Method Available online at http://nrm.srbiau.ac.ir Vol., No., Summer 5 Journal of New Researches in Mathematics Science an Research Branch (IAU) he Uniqueness of the Overall Assurance Interval for Epsilon in

More information

Improving Estimation Accuracy in Nonrandomized Response Questioning Methods by Multiple Answers

Improving Estimation Accuracy in Nonrandomized Response Questioning Methods by Multiple Answers International Journal of Statistics an Probability; Vol 6, No 5; September 207 ISSN 927-7032 E-ISSN 927-7040 Publishe by Canaian Center of Science an Eucation Improving Estimation Accuracy in Nonranomize

More information

Introduction to Markov Processes

Introduction to Markov Processes Introuction to Markov Processes Connexions moule m44014 Zzis law Gustav) Meglicki, Jr Office of the VP for Information Technology Iniana University RCS: Section-2.tex,v 1.24 2012/12/21 18:03:08 gustav

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

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

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

Multi-robot Formation Control Using Reinforcement Learning Method

Multi-robot Formation Control Using Reinforcement Learning Method Multi-robot Formation Control Using Reinforcement Learning Metho Guoyu Zuo, Jiatong Han, an Guansheng Han School of Electronic Information & Control Engineering, Beijing University of Technology, Beijing

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

A Course in Machine Learning

A Course in Machine Learning A Course in Machine Learning Hal Daumé III 12 EFFICIENT LEARNING So far, our focus has been on moels of learning an basic algorithms for those moels. We have not place much emphasis on how to learn quickly.

More information

Comparative Approaches of Calculation of the Back Water Curves in a Trapezoidal Channel with Weak Slope

Comparative Approaches of Calculation of the Back Water Curves in a Trapezoidal Channel with Weak Slope Proceeings of the Worl Congress on Engineering Vol WCE, July 6-8,, Lonon, U.K. Comparative Approaches of Calculation of the Back Water Curves in a Trapezoial Channel with Weak Slope Fourar Ali, Chiremsel

More information

The new concepts of measurement error s regularities and effect characteristics

The new concepts of measurement error s regularities and effect characteristics The new concepts of measurement error s regularities an effect characteristics Ye Xiaoming[1,] Liu Haibo [3,,] Ling Mo[3] Xiao Xuebin [5] [1] School of Geoesy an Geomatics, Wuhan University, Wuhan, Hubei,

More information

Placement and tuning of resonance dampers on footbridges

Placement and tuning of resonance dampers on footbridges Downloae from orbit.tu.k on: Jan 17, 19 Placement an tuning of resonance ampers on footbriges Krenk, Steen; Brønen, Aners; Kristensen, Aners Publishe in: Footbrige 5 Publication ate: 5 Document Version

More information

Including the Consumer Function. 1.0 Constant demand as consumer utility function

Including the Consumer Function. 1.0 Constant demand as consumer utility function Incluin the Consumer Function What we i in the previous notes was solve the cost minimization problem. In these notes, we want to (a) see what such a solution means in the context of solvin a social surplus

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

Summary: Differentiation

Summary: Differentiation Techniques of Differentiation. Inverse Trigonometric functions The basic formulas (available in MF5 are: Summary: Differentiation ( sin ( cos The basic formula can be generalize as follows: Note: ( sin

More information

UNIFYING PCA AND MULTISCALE APPROACHES TO FAULT DETECTION AND ISOLATION

UNIFYING PCA AND MULTISCALE APPROACHES TO FAULT DETECTION AND ISOLATION UNIFYING AND MULISCALE APPROACHES O FAUL DEECION AND ISOLAION Seongkyu Yoon an John F. MacGregor Dept. Chemical Engineering, McMaster University, Hamilton Ontario Canaa L8S 4L7 yoons@mcmaster.ca macgreg@mcmaster.ca

More information

Examining Geometric Integration for Propagating Orbit Trajectories with Non-Conservative Forcing

Examining Geometric Integration for Propagating Orbit Trajectories with Non-Conservative Forcing Examining Geometric Integration for Propagating Orbit Trajectories with Non-Conservative Forcing Course Project for CDS 05 - Geometric Mechanics John M. Carson III California Institute of Technology June

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

Non-Linear Bayesian CBRN Source Term Estimation

Non-Linear Bayesian CBRN Source Term Estimation Non-Linear Bayesian CBRN Source Term Estimation Peter Robins Hazar Assessment, Simulation an Preiction Group Dstl Porton Down, UK. probins@stl.gov.uk Paul Thomas Hazar Assessment, Simulation an Preiction

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

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

Hybrid Fusion for Biometrics: Combining Score-level and Decision-level Fusion

Hybrid Fusion for Biometrics: Combining Score-level and Decision-level Fusion Hybri Fusion for Biometrics: Combining Score-level an Decision-level Fusion Qian Tao Raymon Velhuis Signals an Systems Group, University of Twente Postbus 217, 7500AE Enschee, the Netherlans {q.tao,r.n.j.velhuis}@ewi.utwente.nl

More information

UNDERSTANDING INTEGRATION

UNDERSTANDING INTEGRATION UNDERSTANDING INTEGRATION Dear Reaer The concept of Integration, mathematically speaking, is the "Inverse" of the concept of result, the integration of, woul give us back the function f(). This, in a way,

More information

NON-LINEAR CONTROL STRATEGY FOR STRIP CASTING PROCESSES

NON-LINEAR CONTROL STRATEGY FOR STRIP CASTING PROCESSES NON-LINER CONROL SREGY FOR SRIP CSING PROCESSES V.. Oliveira, R. R. Nascimento, N. S. D. rrifano E. Gesualo, J. P. V. osetti Departamento e Engenharia Elétrica, USP São Carlos CP 359, CEP 3560-970 São

More information

Research Article When Inflation Causes No Increase in Claim Amounts

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

More information

A Novel Unknown-Input Estimator for Disturbance Estimation and Compensation

A Novel Unknown-Input Estimator for Disturbance Estimation and Compensation A Novel Unknown-Input Estimator for Disturbance Estimation an Compensation Difan ang Lei Chen Eric Hu School of Mechanical Engineering he University of Aelaie Aelaie South Australia 5005 Australia leichen@aelaieeuau

More information