Stochastic Perishable Inventory Systems: Dual-Balancing and Look-Ahead Approaches

Size: px
Start display at page:

Download "Stochastic Perishable Inventory Systems: Dual-Balancing and Look-Ahead Approaches"

Transcription

1 Sochasic Perishable Invenory Sysems: Dual-Balancing and Look-Ahead Approaches by Yuhe Diao A hesis presened o he Universiy Of Waerloo in fulfilmen of he hesis requiremen for he degree of Maser of Applied Science in Managemen Sciences Waerloo, Onario, Canada, 2016 Yuhe Diao 2016

2 Auhor's Declaraion I hereby declare ha I am he sole auhor of his hesis. This is a rue copy of he hesis, including any required final revisions, as acceped by my examiners. I undersand ha my hesis may be made elecronically available o he public. ii

3 Absrac We sudy a single-iem, muli-period, sochasic perishable invenory problem under boh backlogging and los-sales circumsances, wih and wihou an order capaciy consrain in each period. We firs model he problem as a dynamic program and hen develop wo heurisics namely, Dual-Balancing (DB) and Look-Ahead (LA) policies, o approximae he opimal invenory level a he beginning of each period. To characerize he holding and backlog cos funcions under he proposed polices, we inroduce a runcaed marginal holding cos for he marginal cos accouning scheme. Our numerical examples demonsrae ha boh DB and LA policies have a possible wors-case performance guaranee of wo in perishable invenory sysems under differen assumpions, and he LA policy significanly ouperforms he DB policy in mos siuaions. We also analyze he arge invenory level in each period (he invenory level a he beginning of each period) under differen policies. We observe ha he arge invenory level under he LA policy is no larger han he opimal one in each period in sysems wihou an order capaciy consrain. iii

4 Acknowledgemens Firs and foremos, I would like o express my sincere graiude o my supervisor Dr. Hossein Abouee Mehrizi for his consan encouragemen, professional insrucions and paience. Wihou his help, he compleion of his hesis would no have been possible. I feel very lucky o have been a suden of such an insighful professor. I also hank Dr. Bissan Ghaddar and Dr. Fama Gzara for heir ime and paience in reviewing my hesis. Secondly, I wan o hank all my friends a Waerloo for he grea ime. Yangzi and Haoran have also helped me a lo. Finally, my hanks would go o my beloved Mum, Dad and Xiaohan, who have always been helping and supporing me ou of difficulies. iv

5 Table of Conens Lis of Figures... vii Lis of Tables... viii 1 Inroducion Lieraure Review Models and Mehodologies Perishable Invenory Sysems Sysem Evoluion in Each Period Objecive Cos Transformaion Dynamic Programming Approach Heurisic Policies Marginal Cos Accouning Dual-Balancing and Look-Ahead Approaches Algorihms Dynamic Programming Algorihm Dual-Balancing Algorihm Look-Ahead Algorihm Numerical Sudy Sysems wih Backlogs and No Ordering Capaciy Consrain Performance of DB and LA Policies Opimal and Heurisic Policies v

6 4.2 Sysems wih Los-Sales and No Ordering Capaciy Consrain Performance of DB and LA Policies Opimal and Heurisic Policies Sysems wih Backlogs and an Ordering Capaciy Consrain Performance of DB and LA Policies Opimal and heurisic Policies Sysems wih Los-Sales and an Ordering Capaciy Consrain Performance of DB and LA Policies Opimal and Heurisic Policies Conclusion References vi

7 Lis of Figures Figure 1. ρ DB for Uniform Demand Figure 2. ρ DB for Binomial Demand Figure 3. ρ DB for Uniform Demand Figure 4. ρ DB for Binomial Demand Figure 5. Cos for h=10, o= Figure 6. Cos for h=0.1, o= Figure 7. ρ DB for Uniform Demand Figure 8. ρ DB for Binomial Demand Figure 9. ρ DB for Uniform Demand Figure 10. ρ DB for Binomial Demand Figure 11. Cos for h=10, o=5 wih an Order Capaciy Consrain vii

8 Lis of Tables Table 1. Per-Uni Cos Combinaions Table 2. ρ LA and ρ DB wih Backlogs and No Capaciy Consrain (%, T=6) Table 3. ρ LA and ρ DB wih Backlogs and No Capaciy Consrain (%, T=8) Table 4. Average Opimal Targe Invenory Level (T=6) Table 5. Average Opimal Targe Invenory Level (T=8) Table 6. Targe Invenory Level under LA and Opimal Policies (T=6) Table 7. Targe Invenory Level under LA and Opimal Policies (T=8) Table 8. ρ LA and ρ DB wih Los-Sales and No Capaciy Consrain (%, T=6) Table 9. ρ LA and ρ DB wih Los-Sales and No Capaciy Consrain (%, T=8) Table 10. Average Opimal Targe Invenory Level (T=6) Table 11. Average Opimal Targe Invenory Level (T=8) Table 12. Targe Invenory Level under LA and Opimal Policies (T=6) Table 13. Targe Invenory Level under LA and Opimal Policies (T=8) Table 14. ρ LA and ρ DB wih Backlogs and a Capaciy Consrain (%) Table 15. Targe Invenory Level under he Opimal Policy Table 16. Comparison of Opimal Expeced Toal Cos Table 17. Targe Invenory Level under he LA Policy Table 18. Targe Invenory Level under LA and Opimal Policies Table 19. ρ LA and ρ DB wih Los-Sales and a Capaciy Consrain (%) Table 20. Targe Invenory Level under he Opimal Policy Table 21. Comparison of Opimal Expeced Toal Cos Table 22. Targe Invenory Level under LA and Opimal Policies viii

9 1 Inroducion Perishable producs are very common in he real life. The main difference beween perishable and non-perishable producs is ha he former mus be consumed wihin a relaively shor period of ime before hey become oudaed, while he laer one can be kep in invenory for a long ime, unil hey are used o saisfied demand. For example, food producs, such as fresh frui and dairy produc, usually have a shor qualiy guaranee period. Besides, blood producs are also perishable and need o be used in a shor ime. For example, blood plaeles are very precious and mus be used wihin six-days [1], and blood plasma also deerioraes wih ime [2]. If perishable producs are no consumed in ime, oudaing can lead o a significan financial loss. For example, abou 321,000 unis of apheresis plaeles oudaed in US in 2011, accouning for 12.8% of he oal processed/produced amoun [3]. Anoher repor indicaes ha he average oudaing rae of perishable producs for reailers and disribuors was around 1.21% in 2007 [4]. On he oher hand, approximaely 50% of he causes of oudaing across a value chain can be miigaed effecively hrough improved planning pracices [4]. These facs indicae he imporance of developing efficien invenory managemen policies for perishable producs. In his sudy, we sudy a sochasic periodic-review invenory sysem over a finie planning horizon wih a perishable produc, which has a fixed lifeime. Inuiively, he lifeime ha we consider needs o be no longer han he lengh of he planning horizon. Also, if he lifeime is 1 period, here will be no excess invenory lef a he end of each period, no maer wheher all demand has been saisfied or no, and we can solve he problem for one period a a ime, separaely. Therefore, in our sudy, we assume ha he lifeime of he producs is longer han 1 period. We assume ha producs in he invenory are used based on he firs-ordered firsconsumed (FOFC) policy o simplify he modeling of he problem. This is a common 1

10 assumpion in he vas majoriy of he lieraure ha producs in he invenory are used based on he FOFC policy [5], which is quie reasonable in pracice. For example, grocery managers usually pu he oldes producs in he mos convenien place on shelves for cusomers, in order o reduce he oudaing rae. The FOFC policy is also he foundaion of many algorihms and policies for perishable invenory problems. For example, i is a criical assumpion for he marginal cos accouning scheme ha maches he available producs and demand in perishable sysems [6], [7]. In general, here are wo differen ways o rea unme demand a he end of each period, eiher backlog hem (sysems wih backlogs) or lose hem (sysems wih los-sales). In backlogging problems, unme demand is backlogged and should be saisfied in he nex period while in los-sales problems, unme demand is los. This difference impacs he period sae ransiion and he amoun of demand a he beginning of each period, which is zero in los-sales problems and he number of backlogs in backlogging problems. Order capaciy is anoher feaure in he invenory problems, which can limi he order quaniy in each period and, hus, affec he ordering policy and he oal cos of he sysem. In his research, we consider boh backlogging and los-sales problems, wih and wihou an order capaciy consrain in each period. The opimal invenory policy for he muli-period invenory sysems wih a perishable produc can be characerized using a dynamic program. However, he dimension of he dynamic program depends on he lengh of he lifeime, and he opimal order quaniy depends on boh he age disribuion of he on-hand invenory and he ime of he curren period. Therefore, he dynamic problem usually has a very large number of sae variables, acion spaces and oucome spaces [1], and he compuaion of he program is usually inracable due o he curse-ofdimensionaliy. Thus, many effors were pu ino developing and examining effecive and efficien heurisics o approximae he opimal policy. For example, he Look-Ahead (LA) 2

11 policy has been proposed for sysems wihou perishable producs [8] and he Dual-Balancing (DB) policy has been applied o boh problems wih and wihou perishable producs [6], [7]. In his sudy, we firs propose algorihms based on boh DB and LA policies for he invenory conrol of he perishable producs wih limied order capaciy. Then, we invesigae he performance of he proposed algorihms using numerical examples. The main conribuions of his hesis are summarized as follows. Algorihms. We use he marginal cos accouning scheme and nesed srucure o exend DB and LA policies for perishable invenory problems wih differen modeling feaures. When here is no order capaciy consrain, he marginal cos accouning was presened in he lieraure [7], [9], and a nesed srucure was used o model problems wih perishable producs [7]. We exend he marginal cos accouning scheme and nesed srucure o perishable problems wih an order capaciy consrain in each period. Wors-case Performance Based on Numerical Examples. For a perishable invenory sysem wihou an order capaciy consrain, he DB policy was proved o have a wors-case performance guaranee of wo in he lieraure [10]. In his sudy, based on our numerical examples, we observe ha he expeced oal coss under he LA policy in perishable sysems wihou an order capaciy consrain, as well as boh DB and LA policies in sysems wih an order capaciy consrain, are less han wice he opimal expeced oal cos. Performance of LA and DB Policies. We compare performance of boh DB and LA policies in perishable invenory sysems under differen assumpions. We observe ha in perishable invenory sysems wihou an order capaciy consrain, he LA policy has much beer performance, on average, han he DB policy, under boh backlogging and los-sales assumpions. In problems wih an order capaciy consrain, while he LA policy ouperforms he DB policy under he backlogging assumpion, i does no always ouperform he DB policy under he los-sales assumpion. 3

12 Truncaed marginal holding cos. We inroduce a new runcaed marginal holding cos for he marginal cos accouning scheme. We show ha wih he runcaed marginal holding cos, we do no need o consider he impac of one decision on he holding cos in all following periods. We also examine he effeciveness of his runcaed marginal holding cos. The res of his hesis is organized as follows. In Chaper 2, we review he lieraure relaed o he sochasic periodic-review perishable invenory sysems and he DB and LA policies. In Chaper 3, we model a perishable invenory problem wih boh backlogs and los-sales, wih and wihou an order capaciy consrain. In Chaper 4, we provide some insighs ino he perishable invenory sysems. In Chaper 5, we conclude our resuls and presen some possible direcions for he fuure sudy. 4

13 2 Lieraure Review Sochasic periodic-review perishable invenory sysems have araced many researchers since 1960 s. A early sages, he focus of he sudies was on developing opimal policies using dynamic programs. The fundamenal characerizaion of he opimal ordering policy wih a woperiod lifeime is provided by Nahmias and Pierskalla [11]. Then, Nahmias [12] and Fries [13] independenly sudy he opimal policy for he general lifeime problem, focusing on backlogging and los-sales problems, respecively. Nahmias [12] uses a per-uni oudaing cos, ordering cos, holding cos, and per-uni per-period shorage cos, which is a sandard cos srucure [14]. They show ha, due o perishabiliy, he srucure of he opimal policy is quie complex. Also, he compuaion of he opimal policy using a dynamic problem is racable only if he lifeime of he produc is shor. This complexiy of he opimal policy is reinforced by Cohen [15], who derives an explici closed-form soluion for a wo-period problem, and discusses procedures o obain he soluion for he m-period case. Many researchers pu effor ino developing heurisic policies for boh backlogging and los-sales models, which would be easier o define and implemen while remaining close o he opimal policy [14]. Brodheim e al. [16] evaluae a class of invenory policies and obain some measures for he invenory sysem, such as he probabiliy of shorage, as well as he easily compuable bounds for he policies. Nahmias [17] proposes hree heurisics, including he criical number and linear approximaions, and compares heir performance using a simulaion model. Nahmias [18] consrucs a bound on he oudaing cos which is only a funcion of he oal on-hand invenory, and develops a myopic policy by subsiuing he oudaing cos wih he bound. Then, Nahmias [19] compares he wo models developed by Nahmias [12] and Fries [13] wih he myopic policy, and shows ha when he remaining periods in he horizon is more han 5

14 he lifeime, he opimal policy in Fries [13] is he same as a discouned version of he one characerized in Nahmias [12]. Nahmias [20] uses a bounded expeced oudaing funcion and a refined ransfer funcion o approximae he problem, focusing on a wo-period problem, which yields o a one-dimensional sae variable. For problems wih los-sales, Nandakumar and Moron [21] derive upper and lower bounds for he order quaniy in each period, and examine Nahmias approximaions wih a weighed average of he upper and lower bounds. The resuls show ha all heurisics have good performance. Cooper [22] considers a fixedcriical number ordering policy and derives new families of bounds on he expeced number of oudaes per period. The performance of he criical-number policies is close o he opimal policy based on his numerical sudy. Recenly, Chen e al. [23] propose wo heurisic policies for join invenory and pricing conrol problems wih boh coninuous and discree demand, and boh backlogging and los-sales cases. They indicae some monooniciy properies of he opimal policy and idenify bounds for he opimal order-up-o invenory level. Li e al. [24] analyze he opimal soluion-srucure of a wo-period lifeime problem and develop a basesock/lis-price heurisic policy. In addiion o he above lieraure focusing on periodic-review perishable invenory sysems, many oher aspecs of perishable invenory problems were considered as well. We refer readers o Kempf e al. [14] for a comprehensive lieraure review. Recenly, Levi e al. [6] uilize a marginal cos accouning scheme and cos balancing idea o develop an approximae policy for sochasic periodic-review invenory problems, and propose a Dual-Balancing as well as a randomized Dual-Balancing algorihms. They show ha boh he Dual-Balancing and randomized Dual-Balancing policies have a wors-case performance guaranee of wo. Then, Levi e al. [25] apply he marginal cos accouning scheme and cos balancing idea o a sochasic invenory problem wih an order capaciy consrain. For los-sales problems, Levi e al. [10] prove ha boh he Dual-Balancing policy and randomized 6

15 Dual-Balancing policy have a wors-case performance guaranee of wo. They also consider a capaciaed model using forced los-sales cos and prove ha he Dual-Balancing policy has a wors-case performance guaranee of wo, even when here is a capaciy consrain on ordering in each period. The marginal cos accouning scheme and Dual-Balancing policy have been also used in perishable invenory sysems. Chao e al. [7] develop he firs approximaion algorihm for a periodic-review perishable invenory sysem wih a perishable produc for boh backlogging and los-sales models. They propose a proporional-balancing policy and a dualbalancing policy wih discouned marginal holding cos, and indicae ha boh policies have a compuaional complexiy of O(mT), which is very efficien compared o he opimal policy. Zhang e al. [9] develop a marginal-cos Dual-Balancing policy for a periodic-review, fixedlifeime perishable invenory conrol problem, and prove ha i has a wors-case performance guaranee of wo. This resul works for boh backlogging and los-sales problems due o a zero lead ime. Then, hey prove ha a myopic policy under he marginal cos accouning scheme provides a lower bound on he opimal ordering quaniy, based on a given on-hand invenory. Anoher approach relaed o our sudy is he Look-Ahead opimizaion approach. Truong [8] sudies a single-iem, muli-period sochasic invenory problem and shows ha he Look- Ahead policy has a wors-case performance guaranee of wo. Also, he Look-Ahead policy, on average, performs well wihin 3.9% of he opimal policy, compared o 19.7% for he Dual- Balancing policy, based on he numerical sudy. 7

16 3 Models and Mehodologies 3.1 Perishable Invenory Sysems We consider a periodic-review sochasic perishable invenory problem, in a finie planning horizon wih T periods, where each period is denoed by, =1,, T. We disinguish a random variable and is realizaion by using capial leers and lower case leers, respecively. The demand in period is denoed by D, =1,, T, which is a random variable. A he beginning of each period, an informaion se f can be observed, conaining he informaion we can have a ha ime. The informaion se has wo pars: realized demands d 1,, d -1, and all oher informaion w 1,, w relaed o he problem in each period. Thus, he informaion se f is one specific realizaion of he random vecor F =(D 1,, D -1 ; W 1,, W ). Furhermore, when f is observed, he fuure demands will have a known condiional join disribuion, including he upcoming curren demand D s, which means E[ D f ] is welldefined and finie for each period s, s=,, T. When we consider differen policies, he informaion se f and he condiional join disribuion is independen of a specific policy [6]. Le lead ime be 0, so he unis ordered in one period can arrive and become available o mee demand in he same period before demand occurs. We apply he firs-ordered firsconsumed (FOFC) policy o calculae he on-hand invenory a he end of each period. This is a common assumpion in perishable problems, because unis in perishable problem have a lifeime, and hose having a shorer remaining lifeime are usually used o saisfy cusomer demand prior o hose arriving laer. The lifeime of he produc is denoed by m, m=2,, T. This means ha afer one uni eners he sysem, i can say for a mos m periods if i has no been used o mee cusomer demand by he end of he m-h periods. For example, suppose here is no on-hand invenory a he beginning of period, and q unis are ordered. These q unis 8

17 arrive in he curren period and are used o mee he demand d. Le m equal o 3 and assume ha he invenory level a he beginning of period is zero. Then, if q >d +d +1 +d +2, here will be sill some excess unis of q a he end of period +2. Bu, hey canno be reained in he sysem because hey have been saying for 3 periods,, +1 and +2. As a resul, he remaining q -(d +d +1 +d +2 ) unis will perish and leave he sysem a he end of period +2. On he oher hand, if q <d +d +1 +d +2, all q unis have been used before hey reach heir lifeime and perish. Some new unis may arrive in period +1 and +2, bu hey need o be kep in he sysem before all q unis arriving in period have been consumed, due o he FOFC policy. A he beginning of period, he on-hand invenory consiss of m-1 ypes of unis, whose remaining lifeime is from m-1 periods o 1 period. Le x,i represen he number of unis having a remaining lifeime of i periods, i=1,, m-1, in he invenory a he beginning of period. Then, he oal invenory a he beginning of period is X m1, i i1 x. Le non-negaive ineger q denoe he number of unis ordered in period. Because lead ime is 0, he q unis will arrive and become available in curren period before demand occurs, wih a lifeime of m periods. So in period, he oal invenory used o mee demand is Y =X +q, which is he basesock level in period. As we discussed in Secion 1, wo general feaures of invenory problems in he lieraure are wheher excess demand is backlogged or los, and wheher here is a limi on he ordering level in each period. In problems wihou an ordering capaciy consrain, he number of ordered unis can be any non-negaive ineger, bu in problems wih an ordering capaciy consrain u, he order quaniy canno surpass he capaciy limi. In a periodic-review sochasic perishable invenory conrol problem, several ypes of cos are incurred. A per-uni ordering cos c ˆ occurs when one uni is ordered in period. A per- ˆ uni holding cos is denoed by h, which occurs a he end of period for each excess uni ha 9

18 is sill in invenory, represening a cos of keeping his uni from period o +1. These excess unis can be held in invenory for a long ime if he demands in he following periods are no high enough o consume all available unis. Also, because one uni can say in he sysem for a mos m periods, he holding cos incurred by one order occurs in a mos m-1 following periods. If he demand in period canno be saisfied immediaely, here will be a penaly, which is ˆ b per uni of unsaisfied demand. This cos has differen meanings in he backlogging and lossales problems. In a backlogging problem, he unsaisfied demand is backlogged in he sysem and can be eliminaed by ordering more in he following period. Therefore, his penaly means he compensaion o he unsaisfied demand and will exis in only one period, if here is no ordering capaciy consrain. Sysems wih an ordering capaciy consrain in each period will be discussed in more deail laer. In a los-sales problem, he unfulfilled demand will leave he sysem a he end of he period. Thus, he penaly is he loss of value caused by he unsaisfied demand and will occur once. For simpliciy, he per-uni penaly in boh ypes of problems are denoed by ˆ b. The fourh cos in a perishable problem is a per-uni oudaing cos o ˆ. I occurs when a uni has been in he sysem for m periods and perishes, before been used o mee demand Sysem Evoluion in Each Period Sysems wih Backlogs and No Ordering Capaciy Consrain In his secion, we discuss a complee process of a perishable invenory sysem wih backlogs, ineger demand and order quaniy, and no ordering capaciy consrain. 1. A he beginning of period, =1,, T, informaion se f and on-hand invenory x,i, i=1,, m-1, are observed, where i, x i 1,, m 1 0, 1, m1, i x 0, x 0, i1,, m 2 {2,, T}, means backlogs in period 1 i, {2,, }, means x 0, i1,, m1 T no backlogs in period 1 Therefore, he oal invenory level is X 10 m1, i i1 x. The condiional join disribuion of.

19 upcoming demand D is known based on he informaion se f. 2. The sysem orders q unis, and hey arrive and become available, having a lifeime of m periods. In sysems wihou an ordering capaciy consrain, he feasible range of q is {0,, }. The oal cos of ordering is cq ˆ. Also, he on-hand invenory level increases from X o Y =X +q. 3. Demand d occurs, and he on-hand invenory is used o mee his demand. The x,1 unis having only 1-period lifeime, are consumed firs. If x,1 >d, he remaining x,1 -d unis will perish and an oudaing cos of,1 oˆ ( x d ) occurs. Oherwise, he x,2 unis having 2- period lifeime are used o mee he excess d -x,1 unis of demand, and so on, unil all invenory is consumed. If d >Y, he demand canno be saisfied compleely by he curren invenory. The remaining d -Y unis of demand becomes backlogs ha need o be me in period +1, wih a backlogging penaly of bˆ ( d Y ). If d <Y, he remaining Y -d unis resul in a holding cos of hˆ ( Y d ). 4. A he end of period, he remaining lifeime of he excess unis of he on-hand invenory 1 decreases by 1 period. Noe ha if, i i, j x d x j1, where x max{ x,0}, all x,i 1 unis will be used o mee demand. Oherwise,, i i, j x d x j1 unis remain for period +1, becoming x +1,i-1, i=2,, m-1. In period, he q unis have a lifeime of m periods, and hey will mee he excess demand afer all X have been used. Therefore, in a m1, j backlogging problem, if q d x j1 m1, j backlogs is q d x j1, backlogs occur and he number of. These backlogs will reain in he sysem for period +1 as a negaive invenory. So he new invenory level a he beginning of period +1 can be calculaed using (similar o Eq. (1) in [7]) 11

20 i 1, i, i1, j x x d x, i {1,, m 2}, j1 x q d x j1 m1 1, m1, j. (1) Then, d will be added ino f +1, which will be he given informaion for period +1, as well as x +1,i, i=1,, m Sysems wih Los-Sales and No Ordering Capaciy Consrain The process of a los-sales problem wihou an order capaciy consrain can be considered as follows. 1. A he beginning of period, =1,, T, informaion se f and on-hand invenory x,i are observed, where i, x i1,, m1 x i, 0, 1. 0, i 1,, m 1 {2,, T} The oal invenory level is X m1, i i1 x. The condiional join disribuion of he upcoming demand D is known based on he informaion se f. 2. The sysem orders q unis, and hey arrive and become available, having a lifeime of m periods. In sysems wihou an order capaciy consrain, he feasible range of q is {0,, }. The oal cos of ordering is cq ˆ. Also, he on-hand invenory level increases from X o Y =X +q. 3. Demand d occurs, and he on-hand invenory is used o mee his demand. The x,1 unis having only 1-period lifeime are consumed firs. If x,1 >d, he remaining x,1 -d unis will perish and an oudaing cos of,1 oˆ ( x d ) occurs. Oherwise, he x,2 unis having 2- period lifeime are used o mee he excess d -x,1 unis of demand, and so on, unil all invenory is consumed. If d >Y, he demand canno be saisfied compleely by he curren invenory. The remaining d -Y unis of demand will be los and leave he sysem, wih a 12

21 los-sales penaly of bˆ ( d Y ). If d <Y, he remaining Y -d unis resul in a holding cos of hˆ ( Y d ). 4. A he end of period, he remaining lifeime of he excess unis of he on-hand invenory 1 decreases by 1 period. Noe ha if, i i, j x d x j1 1 demand. Oherwise,, i i, j x d x j1, all x,i unis will be used o mee unis remain for period +1, becoming x +1,i-1, i=2,, m-1. In period, he q unis have a lifeime of m periods, and hey will mee he excess demand afer all X have been used. Therefore, in a los-sales problem, if 1, m j1, he excess demand will leave he sysem. If 1 m, j q d x j1 q d x j m1, j he q d x j1, remaining unis will become x +1,m-1. So he new invenory level a he beginning of period +1 can be calculaed using (similar o Eq. (18) in [7]) i 1, i, i1, j x x d x, i {1,, m 2}, j1 x q d x j1 m1 1, m1, j. (2) Sysems wih Backlogs and an Ordering Capaciy Consrain In a backlogging problem wih an order capaciy consrain, he order quaniy is limied and he complee process in period can be lised as follows. 1. A he beginning of period, =1,, T, informaion se f and on-hand invenory x,i, i=1,, m-1, are observed, where i, x i 1,, m 1 0, 1, m1, i x 0, x 0, i1,, m 2 {2,, T}, means backlogs in period 1 i, {2,, }, means x 0, i1,, m1 T no backlogs in period 1 The oal invenory level is X m1, i i1 x and he condiional join disribuion of. upcoming demand D is known based on he informaion se f. 2. The sysem orders q unis, and hey arrive and become available, having a lifeime of m 13

22 periods. In a sysem wih an order capaciy consrain u, he feasible range of q is {0,, u }. The oal cos of ordering is cq ˆ. Also, he on-hand invenory level increases from X o Y =X +q. 3. Demand d occurs, and he on-hand invenory is used o mee his demand. The x,1 unis having only 1-period lifeime are consumed firs. If x,1 >d, he remaining x,1 -d unis will perish and an oudaing cos of,1 oˆ ( x d ) occurs. Oherwise, he x,2 unis having 2- period lifeime are used o mee he excess d -x,1 unis of demand, and so on, unil all invenory is consumed. If d >Y, he demand canno be saisfied compleely using he curren invenory. Therefore, he remaining d -Y unis of demand becomes backlogs ha need o be me in period +1, and a backlogging penaly of bˆ ( d Y ) occurs. If d <Y, he remaining Y -d unis resul in a holding cos of hˆ ( Y d ). 4. A he end of period, he remaining lifeime of he excess unis of he on-hand invenory 1 decreases by 1 period. Noe ha if, i i, j x d x j1 1 demand. Oherwise,, i i, j x d x j1, all x,i unis will be used o mee unis remain for period +1, becoming x +1,i-1, i=2,, m-1. In period, he q unis have a lifeime of m periods, and hey will mee he excess demand afer all X have been used. Therefore, in a backlogging problem, if j1 m1, j q d x, backlogs occur and he number of backlogs is m1, j q d x j1. These backlogs will reain in he sysem for period +1 as a negaive invenory. So he new invenory level a he beginning of period +1 can be calculaed using (similar o Eq. (1)) i 1, i, i1, j x x d x, i {1,, m 2}, j1 x q d x j1 m1 1, m1, j 14.

23 Sysems wih Los-Sales and an Ordering Capaciy Consrain The process of a los-sales problem wih an order capaciy consrain can be considered as follows. 1. A he beginning of period, =1,, T, informaion se f and on-hand invenory x,i are observed, where i, x i1,, m1 x i, 0, 1. 0, i 1,, m 1 {2,, T} Then oal invenory level is X m1, i x and he condiional join disribuion of he i1 upcoming demand D is known based on he informaion se f. 2. The sysem orders q unis, and hey arrive and become available, having a lifeime of m periods. In a sysem wih an order capaciy consrain u, he feasible range of q is {0,, u }. The oal cos of ordering is cq ˆ. Also, he on-hand invenory level increases from X o Y =X +q. 3. Demand d occurs, and he on-hand invenory is used o mee his demand. The x,1 unis having only 1-period lifeime are consumed firs. If x,1 >d, he remaining x,1 -d unis will perish and an oudaing cos of,1 oˆ ( x d ) occurs. Oherwise, he x,2 unis having 2- period lifeime are used o mee he excess d -x,1 unis of demand, and so on, unil all invenory is consumed. If d >Y, he demand canno be saisfied compleely using he curren invenory. Therefore, he remaining d -Y unis of demand will be los and leave he sysem, wih a los-sales penaly of bˆ ( d Y ). If d <Y, he remaining Y -d unis resul in a holding cos of hˆ ( Y d ). 4. A he end of period, he remaining lifeime of he excess unis of he on-hand invenory 1 decreases by 1 period. Noe ha if, i i, j x d x j1 1 demand. Oherwise,, i i, j x d x j1, all x,i unis will be used o mee unis remain for period +1, becoming x +1,i-1, 15

24 i=2,, m-1. In period, he q unis have a lifeime of m periods, and hey will mee he excess demand afer all X have been used. Therefore, in a los-sales problem, if 1, m j1, he excess demand will leave he sysem. If 1 m, j q d x j1 q d x j m1, j he q d x j1, remaining unis of q will become x +1,m-1. Therefore, he new invenory level a he beginning of period +1 can be calculaed using (similar o Eq. (2)) i 1, i, i1, j x x d x, i {1,, m 2}, j1 x q d x j1 m1 1, m1, j Objecive We consider a finie horizon problem wih T periods. The oal cos occurring in each period includes holding, backlogging/los-sales penaly, and oudaing coss. A he beginning of period, he informaion se f and on-hand invenory x,i, i=1,, m-1 are observed, and order quaniy q is deermined, which causes an ordering cos of cq ˆ. The holding cos occurs when here are some unis lef afer fulfilling he demand, which means ha Y -d >0. Therefore, he holding cos in period is ˆ ( h Y ) D. If Y -d <0, no all demands are saisfied, and a backlogging/los-sales penaly of ˆ ( b D Y ) occurs. Also, in each period, perishable unis cause an oudaing cos of,1 oˆ ( x D ). Le α, 0 1, denoe he discoun facor. Then, because demand in each period is uncerain, he expeced oal cos in period under policy P can be obained using ˆ ˆ ˆ ˆ, 1,1 S (P) E c q h ( Y D ) b ( D Y ) o ( x D ) f where m1, i. i1 Y x q A he end of T periods, some unis may remain in he sysem, which can be sold in a salvage marke. Le ˆv denoe he per-uni salvage value. Then, he oal salvage value a he 16

25 end of T periods is T m1 T 1, i vˆ i1 x. The objecive is o find he opimal policy based on he informaion se f and on-hand invenory x,i, i=1,, m-1, ha minimizes he expeced oal cos over T periods (similar o Eq. (2) in [7]): C(P) T m1 T T 1, i = S (P) E vˆ x 1 i1 T m1 1 ˆ ˆ,1 T T 1, i E cˆ q h ( Y D ) b ( D Y ) oˆ ( x D ) vˆ x. 1 i1 (3) Cos Transformaion Suppose ha all per-uni coss are consan over he T periods and denoe ĉ, ĥ, ˆb and ô as he consan per-uni ordering cos, holding cos, backlogging/los-sales penaly and oudaing cos. Then, ˆ ˆ ˆ ˆ,, ˆ ˆ c c h h b b, oˆ oˆ, {1,, T}. A cos ransformaion can be carried ou o make ordering cos c=0. Le salvage value ˆv be equal o ĉ, and se non-negaive cos parameers o c=0, h hˆ (1 ) cˆ, b bˆ (1 ) cˆ and o oˆ cˆ. Then, i has been proved ha he expeced oal cos described in (3) can be rewrien as (see Proposiion 1 in [7]) C(P) E h( Y D ) b( D Y ) o( x D ) ce D 1 1 T T 1,1 1 ˆ. Noe ha 1 T cˆe D 1 is a consan and does no depend on he invenory policy. To simplify he noaion, we le α=1 and only consider he expeced cos per period and expeced oal cos as: S h Y D b D Y o x D f,1 (P) E ( ) ( ) ( ), T C S h Y D b D Y o x D f T,1 (P) (P)=E ( ) ( ) ( )

26 3.2 Dynamic Programming Approach Dynamic programming approach is usually used o describe a muli-period sochasic invenory conrol problem, which is defined recursively and deermines he order quaniy q in period o minimize he expeced oal cos over he ime horizon. Noe ha he sae a he beginning of period consiss of he informaion se f and he on-hand invenory x,i, i=1,, m-1, and he expeced cos in period is.,1 S (P) E h( Y D ) b( D Y ) o( x D ) f Le V x x f,1, m1 (,,, ) denoe he opimal expeced cos over periods o T, and suppose ha V x x f j,1 j, m1 j(,,, j), j=+1,, T, have been calculaed. The Bellman formulaion for calculaing he opimal expeced cos in period for he perishable invenory problem wihou an order capaciy consrain is V x x f,1, m1 (,,, ) 1,1 1, m1 min S E V 1(,,, 1) x x F f q {0,, },1 E h( Y D ) b( D Y ) o( x D ) f min, q {0,, } 1,1 1, m1 E V 1( x,, x, F 1) f and for he perishable invenory problem wih an order capaciy consrain is V x x f,1, m1 (,,, ) 1,1 1, m1 min S E V 1(,,, 1) x x F f q {0,, u },1 E h( Y D ) b( D Y ) o( x D ) f min. q {0,, u } 1,1 1, m1 E V 1( x,, x, F 1) f Considering Eq. (1) and (2), he sae ransiion equaions for a problem wih backlogs are i 1, i, i1, j x x D x, i {1,, m 2}, j1 x q D x j1 m1 1, m1, j and for a problem wih los-sales are 18,

27 i 1, i, i1, j x x d x, i {1,, m 2}, j1 x q d x j1 m1 1, m1, j. Due o he curse-of-dimensionaliy, solving he dynamic program discussed above o obain he opimal order quaniy in each period is no sraighforward [6]. Therefore, in he nex secion, we propose wo simple efficien policies based on he marginal cos accouning scheme o approximae he opimal invenory level a he beginning of each period. 3.3 Heurisic Policies Marginal Cos Accouning In he marginal cos accouning scheme, he cos incurred by a decision is no considered periodically. Insead, he oal cos incurred by one decision over all he periods from when he decision is made o he end of he planning horizon is calculaed. There are hree differen ypes of coss caused by ordering q unis in period, which are he holding, backlogging/los-sales penaly, and oudaing coss. When hey occur, he number of periods hey will las for depends on he specific model assumpions. Also, noe ha he cos occurring in one period may conain he marginal coss caused by decisions made in previous periods. In his secion, we firs characerize hese marginal coss. Recall ha he FOFC policy is applied here, meaning ha unis ordered earlier should be used o saisfy demand prior o hose ordered laer. One imporan feaure abou using he marginal cos accouning is ha once one decision is made, he coss caused by i from he curren period o he end of he planning horizon are only affeced by he demand, bu no he decisions in he subsequen periods. This feaure ensures ha he expeced oal coss associaed o differen decisions can be obained before he decision is made, so we can compare he 19

28 differen decisions and heir corresponding expeced oal coss, and choose he bes one Sysems wih Backlogs and No Ordering Capaciy Consrain a) Marginal holding cos accouning Holding cos occurs when some unis have no been consumed a he end of one period. Le he order quaniy in period be q. Because hese q unis have a lifeime of m periods in a perishable problem when hey become available, hey can say in he sysem for a mos m periods, from period o +m-1, if hey are no used o mee demand. Therefore, in he marginal holding cos accouning, he holding cos incurred by hese q unis may occurs a mos from period o period +m-1. In order o calculae he marginal holding cos incurred by hese q unis, we need o know how many unis of hese q unis can reain in he sysem a he end of each period. Recall ha a he beginning of period, he on-hand available invenory is X m1, i x. i1 These X unis need o be consumed prior o q due o he FOFC policy. Therefore, wheher q unis will be used in period depends on he sign of D -X. If D -X >0, he on-hand invenory is no enough o mee demand in period, hus, some unis of q are needed. If D -X 0, all q unis will be kep in he sysem for he nex period. Therefore, he remaining unis of q a he end of period is m 1, i q D x i1. In period +1, afer unused unis of x,1, if here are any, have perished, he number of unme D is,1 ( D x ), which will be me by m1 i, i2 x. Then, D +1 occurs and his will also be me by m1 i, x before q is used. Therefore, he remaining unis of q a he end of period +1 is i2 m 1 1,1, i q D D x x i2. The number of remaining unis of q a he end of period +2,, +m-1 can be calculaed 20

29 similarly. For he simpliciy of expression, a nesed srucure is used o describe hese resuls (see e.g., [7]). Le B ( x, i) denoe he ousanding demand in period +i-1 afer he unis having a lifeime of i periods or less have been depleed. Then, (similar o Eq. (5) in [7]) B ( x, i) : 0 i 0. (4) (, 1) {1,, 1} i 1, i D B x i x i m Therefore, he number of remaining unis of q a he end of period j, j=,, +m-1, can be calculaed, using m 1 j, i q D B x, j x. i j1 Considering ha here are T periods in he horizon, he oal marginal holding cos over period o +m-1 incurred by ordering q unis in period can be defined as (similar o Eq. (4) in [7]) where x y min{ x, y}. ( m1) T m 1 j, i H ( q ) : h q D B ( x, j ) x j, i j1 In pracical siuaions, a runcaed marginal holding cos can be calculaed. Insead of adding all holding cos incurred by q from period o ( m 1) T, he runcaed marginal holding cos considers he holding cos from period o +l-1, l=1,, m. Le, l H ( q ), l=1,, m, denoe he runcaed marginal holding cos, which represens he marginal holding cos incurred by q over periods o ( m 1) T. When l=1, only he marginal holding cos in period is considered. Thus, he runcaed marginal holding cos incurred by q in period wih l=1 is m 1,1, i H ( q ) h q D x i1. When l=2 and <T,,2 H ( q ) can be obained by adding he marginal holding cos in period 21

30 and +1, which is H,2 ( q ) h q D B ( x, j ) x j i2 ( 1) T m1 j, i m1 m1, i 1,1, i h q D x h q D D x x i1 i2 If =T, we only need o consider he marginal holding cos in period T, and i is easy o see ha in period T, no maer wha l is, he runcaed marginal holding cos will be. m 1 T, l T T T T, i H ( q ) h q D x i1. For l=3,, m-1, he runcaed marginal holding cos can be obained using a similar mehod, and when l=m, he runcaed marginal holding cos is exacly equal o he marginal holding cos. In general, he runcaed marginal holding cos is ( l1) T m 1, l j, i H ( q ) : h q D B ( x, j ) x, l {1,, m} j. i j1 An inuiive explanaion for he runcaed marginal holding cos is ha because he ordering cos is consan, ordering more in one period for he fuure demand will cos more han ordering in fuure periods. Therefore, here is no moivaion o order more in he curren period for he fuure demand unless he variabiliy of demand is high. Hereafer, we will use, l H ( q ), l=1,, m, insead of H ( q ). b) Marginal backlogging penaly accouning The marginal backlogging penaly caused by ordering q unis only occurs in period in a backlogging problem wihou an order capaciy consrain. Before he demand is observed, he available on-hand invenory ha can be used o mee he demand is m1, i i1 x q. Therefore, he unsaisfied demand is 22

31 m1, i D x q, i1 and he marginal backlogging penaly is (similar o Eq. (7) in [7]) m1, i ( q ) : bd x q. (5) i1 c) Marginal oudaing cos accouning The marginal oudaing cos incurred by ordering q in period occurs in period +m-1 if hese q unis have no been consumed compleely by he end of period +m-1. This oudaing cos can only be incurred by orders placed from period o T-m+1 due o he m-period lifeime. The unis arriving in period T-m+2 o T will have no reached heir lifeime of m periods by he end of period T, so he oudaing cos for hese orders is 0. For he q unis ordered in period, =1,, T-m+1, hey will be used o mee demand only when all m1 i, x have been used up. Noe ha x,i unis have a lifeime of i periods, so hey i1 will perish and leave he sysem a he end of period +i-1 if hey have no been used by hen. Therefore, he nesed srucure in he marginal holding cos is also needed o calculae he oudaing cos. Wih he same definiion of B ( x, i) given in he marginal holding cos, he marginal oudaing cos can be obained using (similar o eq. (6) in [7]) m1 ( ): (, 1), {1,, 1} O q o q D B x m T m. (6) Sysems wih Los-Sales and No Ordering Capaciy Consrain a) Marginal holding cos accouning The difference beween sysems wih backlog and los-sales is he way he excess demand is reaed. However, holding cos only occurs when excess invenory exiss. In period, if he on-hand invenory m1, i x q is larger han d, here will be no excess demand a he end i1 of he period. Therefore, he runcaed marginal holding cos in los-sales problems wihou an 23

32 order capaciy consrain is ( l1) T m 1, l j, i H ( q ) : h q D B ( x, j ) x, l {1,, m} j. i j1 b) Marginal los-sales penaly accouning When he demand in period is no saisfied compleely, he excess par is los a he end of period and leads o a los-sales penaly. This los-sales penaly only occurs once a he end of period. Alhough in a backlogging problem his excess demand reain in he sysem, he calculaion of he amoun of he excess demand is he same for boh sysems wih backlogs and los-sales, which is m1, i D x q. i1 Therefore, he marginal los-sales penaly for ordering q unis in period is (see Eq. (5)) m1, i ( q ) : bd x q. i1 c) Marginal oudaing cos accouning The marginal oudaing cos incurred by ordering q in period occurs in period +m-1 if hese q unis have no been consumed compleely by he end of period +m-1. If an oudaing cos occurs, he calculaion is he same as he one given in Secion , which is (see Eq. (6)) m1 ( ): (, 1), {1,, 1} O q o q D B x m T m Sysems wih Backlogs and an Ordering Capaciy Consrain a) Marginal holding cos accouning For any given order quaniy a he beginning of period, q, he available invenory a he beginning of period is m1, i i1 x q. Therefore, he runcaed marginal holding cos, which is based on he available invenory and observed demand, follows he discussion in Secion 24

33 and can be obained using ( l1) T m 1, l j, i H ( q ) : h q D B ( x, j ) x, l {1,, m} j. i j1 b) Marginal backlogging penaly accouning To obain he marginal backlogging penaly in sysems wih an ordering capaciy consrain, a forced shorage is inroduced in [25]. Noe ha if he order capaciy in period, u, is less han he order quaniy q P under policy P in a problem wihou an ordering capaciy consrain, he order capaciy makes q o be less han q P. Therefore, he number of backlogs, m1, i ha should be considered o calculae he backlogging penaly is no D x q i1 because he order quaniy q canno reach q P due o he capaciy consrain. Also, he penaly of he unsaisfied demand beyond u canno be calculaed as he backlogging penaly caused by he decision of ordering q unis. Le q denoe he slack capaciy in period, where q u q, 0 q u. In period, we use a forced shorage, which is he number of unis ha could have been ordered more o avoid he backlogging penaly, o calculae he cos of he backlogs, insead of comparing he available invenory and demand direcly. Le W,s denoe he forced shorage, which occurs in period s, s=,, T, because of he ordering decision in period. Then, when s=, we ge m 1,, i W min q, D x q i1. (7) If all demand in period is saisfied, or we order u unis in period, W, will be 0. If backlog occurs in period, he capaciy u +1 should be considered o see wheher he sum of he backlogs and new demand d +1 in period +1 can be saisfied by ordering q +1 unis under he order capaciy u +1. If so, no backlog will exis a he end of period +1. If no, backlogs from period should be saisfied wih a higher prioriy due o he FOFC policy. Afer x,1 unis have 25

34 perished a he end of period, he remaining demand is m1, i i2 x increases o D x,1, and he invenory is q. Suppose ha u +1 unis are ordered in period +1. Then, he available invenory m1, i 1 i2 x q u +1 increases o D 1 D x,1 1,1 m1, i 1 D D x x q u. Afer demand D +1 is observed, he oal demand in period. Comparing he available invenory and demand, if i2, ordering q unis in period causes no backlogging penaly in period +1. Oherwise, ordering q unis in period is no enough o mee he cumulaive demand D [, 1], where D [, s] s i i D, s=,, T, even hough he unis of he maximum feasible order quaniy u +1 are ordered in period +1. Noe ha m1 i2 1,1, i 1 D D x x q u is he shorage of invenory in period +1. Similar o period, we compare his value wih he slack capaciy caused by ordering q unis, which is q o deermine he forced shorage m1, 1 1,1, i 1 W min q, D D x x q u i2. Because u +1 is he maximum feasible order quaniy in period +1, he forced shorage W,+1 using u +1 is no larger han ha using oher order quaniies under he ordering capaciy consrain in period +1. The analysis for he subsequen periods o period +1 is similar. We define a new vecor x, which consiss of unis ha are already in he invenory in period and ha will come ino he invenory, as,1, m 1 1 T : x,, x x ; q ; u,, u, where =1,, T-1. If =T, no u unis will be ordered by he end of T periods, so we have T T,1 T, m 1 T x : x,, x ; q. Noe ha x consiss of invenory x,i, q and order capaciy u j, for j=+1,...,t. Therefore, here are T-+m componens in oal. The orders of q and u j, j=+1,...,t, unis are placed from period 26

35 o T, from when hey will have a lifeime of m periods. We can consider he periods from o he period when he q and u j, j=+1,...,t, unis are ordered as an exra lifeime for hese unis. For example, consider period and assume ha u +1 unis will be ordered in period +1. The one period from o +1 can be considered as an exra lifeime for he u +1 unis. Thus, in period, we can hink ha he u +1 unis have a lifeime of m+1 periods. Therefore, generally, in period, he I, x unis, I=1,,T-+m, have a lifeime from 1 period o T-+m periods. Bu, when hey come ino he invenory and become available, hey all have a lifeime of only m periods. Then, we modify he definiion of B ( x, i) in Eq. (4) o B ( x, I), using x and I, as B ( x, I) : 0 I 0 1, I I D B ( x, I 1) x I 1,, T m. Noe ha B ( x, I) denoes he ousanding demand in period +I-1 afer he unis having a lifeime of I periods or less have been depleed, similar o B ( x, i). The difference beween B ( x, I) and B ( x, i) is ha he former considers order quaniies placed in periods +1,, T. Therefore, a general expression for he forced shorage in period s, s=,, T, which is caused by ordering q unis in period is s m, s s, I W min q, D B ( x, s ) x, s {,, T}, I s1 and he corresponding backlogging penaly in he problem wih an order capaciy consrain is T, s ( q ) : b W. s c) Marginal oudaing cos accouning In a backlogging problem wih an order capaciy consrain, he order capaciy does no impac he calculaion of he marginal oudaing cos. If he order quaniy is smaller han ha in a problem wihou an order capaciy consrain, he probabiliy of oudaing decreases. Therefore, he oudaing cos can be calculaed using (see Eq. (6)) 27

36 m1 ( ): (, 1), {1,, 1} O q o q D B x m T m Sysems wih Los-Sales and an Ordering Capaciy Consrain a) Marginal holding cos accouning The derivaion of he runcaed marginal holding cos is based on he available invenory a he beginning of a period and observed demand during ha period. Therefore, following he discussion in Secion , we ge ( l1) T m 1, l j, i H ( q ) : h q D B ( x, j ) x, l {1,, m} j. i j1 b) Marginal los-sales penaly accouning In a los-sales problem wih an order capaciy consrain, similar o Secion , he excess demand in period canno be used o calculae he los-sales penaly incurred by ordering q unis in period direcly. This is because his excess demand may surpass he difference beween q and he order capaciy u. Therefore, he forced shorage is needed o obain he marginal los-sales penaly. Similar o he backlogging problem wih an ordering capaciy consrain, if we always order a he capaciy level in period s, s=+1,, T, and unsaisfied demand occurs in a period, his los-sales penaly is associaed o he shorage in period. 1, In period, he excess demand is m D x i q. We compare his excess demand i1 wih he slack capaciy q u q, 0 q u, o deermine he exac los-sales quaniy ha is caused by ordering q unis in period, which is (see Eq. (7)) m 1,, i W min q, D x q i1. Because he excess demand will be los in a los-sales problem, so in period +1 he only demand ha needs o be me is D +1. Noe ha he unused unis of x,1 will perish a he end of period, and should be deduced from he invenory. Therefore, a he end of period, afer x,1 unis have perished, he excess invenory is 28

T L. t=1. Proof of Lemma 1. Using the marginal cost accounting in Equation(4) and standard arguments. t )+Π RB. t )+K 1(Q RB

T L. t=1. Proof of Lemma 1. Using the marginal cost accounting in Equation(4) and standard arguments. t )+Π RB. t )+K 1(Q RB Elecronic Companion EC.1. Proofs of Technical Lemmas and Theorems LEMMA 1. Le C(RB) be he oal cos incurred by he RB policy. Then we have, T L E[C(RB)] 3 E[Z RB ]. (EC.1) Proof of Lemma 1. Using he marginal

More information

Inventory Control of Perishable Items in a Two-Echelon Supply Chain

Inventory Control of Perishable Items in a Two-Echelon Supply Chain Journal of Indusrial Engineering, Universiy of ehran, Special Issue,, PP. 69-77 69 Invenory Conrol of Perishable Iems in a wo-echelon Supply Chain Fariborz Jolai *, Elmira Gheisariha and Farnaz Nojavan

More information

Vehicle Arrival Models : Headway

Vehicle Arrival Models : Headway Chaper 12 Vehicle Arrival Models : Headway 12.1 Inroducion Modelling arrival of vehicle a secion of road is an imporan sep in raffic flow modelling. I has imporan applicaion in raffic flow simulaion where

More information

3.1.3 INTRODUCTION TO DYNAMIC OPTIMIZATION: DISCRETE TIME PROBLEMS. A. The Hamiltonian and First-Order Conditions in a Finite Time Horizon

3.1.3 INTRODUCTION TO DYNAMIC OPTIMIZATION: DISCRETE TIME PROBLEMS. A. The Hamiltonian and First-Order Conditions in a Finite Time Horizon 3..3 INRODUCION O DYNAMIC OPIMIZAION: DISCREE IME PROBLEMS A. he Hamilonian and Firs-Order Condiions in a Finie ime Horizon Define a new funcion, he Hamilonian funcion, H. H he change in he oal value of

More information

Applying Genetic Algorithms for Inventory Lot-Sizing Problem with Supplier Selection under Storage Capacity Constraints

Applying Genetic Algorithms for Inventory Lot-Sizing Problem with Supplier Selection under Storage Capacity Constraints IJCSI Inernaional Journal of Compuer Science Issues, Vol 9, Issue 1, No 1, January 2012 wwwijcsiorg 18 Applying Geneic Algorihms for Invenory Lo-Sizing Problem wih Supplier Selecion under Sorage Capaciy

More information

Inventory Analysis and Management. Multi-Period Stochastic Models: Optimality of (s, S) Policy for K-Convex Objective Functions

Inventory Analysis and Management. Multi-Period Stochastic Models: Optimality of (s, S) Policy for K-Convex Objective Functions Muli-Period Sochasic Models: Opimali of (s, S) Polic for -Convex Objecive Funcions Consider a seing similar o he N-sage newsvendor problem excep ha now here is a fixed re-ordering cos (> 0) for each (re-)order.

More information

Energy Storage Benchmark Problems

Energy Storage Benchmark Problems Energy Sorage Benchmark Problems Daniel F. Salas 1,3, Warren B. Powell 2,3 1 Deparmen of Chemical & Biological Engineering 2 Deparmen of Operaions Research & Financial Engineering 3 Princeon Laboraory

More information

2-Approximation Algorithms for Perishable Inventory Control When FIFO Is an Optimal Issuing Policy arxiv: v2 [math.

2-Approximation Algorithms for Perishable Inventory Control When FIFO Is an Optimal Issuing Policy arxiv: v2 [math. Vol. 00, No. 0, Xxxxx 0000, pp. 000 000 issn0000-0000 eissn0000-0000 00 0000 0001 INFORMS doi 10.1287/xxxx.0000.0000 c 0000 INFORMS 2-Approximaion Algorihms for Perishable Invenory Conrol When FIFO Is

More information

10. State Space Methods

10. State Space Methods . Sae Space Mehods. Inroducion Sae space modelling was briefly inroduced in chaper. Here more coverage is provided of sae space mehods before some of heir uses in conrol sysem design are covered in he

More information

Excel-Based Solution Method For The Optimal Policy Of The Hadley And Whittin s Exact Model With Arma Demand

Excel-Based Solution Method For The Optimal Policy Of The Hadley And Whittin s Exact Model With Arma Demand Excel-Based Soluion Mehod For The Opimal Policy Of The Hadley And Whiin s Exac Model Wih Arma Demand Kal Nami School of Business and Economics Winson Salem Sae Universiy Winson Salem, NC 27110 Phone: (336)750-2338

More information

An Inventory Model for Time Dependent Weibull Deterioration with Partial Backlogging

An Inventory Model for Time Dependent Weibull Deterioration with Partial Backlogging American Journal of Operaional Research 0, (): -5 OI: 0.593/j.ajor.000.0 An Invenory Model for Time ependen Weibull eerioraion wih Parial Backlogging Umakana Mishra,, Chaianya Kumar Tripahy eparmen of

More information

Physics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle

Physics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle Chaper 2 Newonian Mechanics Single Paricle In his Chaper we will review wha Newon s laws of mechanics ell us abou he moion of a single paricle. Newon s laws are only valid in suiable reference frames,

More information

Competitive and Cooperative Inventory Policies in a Two-Stage Supply-Chain

Competitive and Cooperative Inventory Policies in a Two-Stage Supply-Chain Compeiive and Cooperaive Invenory Policies in a Two-Sage Supply-Chain (G. P. Cachon and P. H. Zipkin) Presened by Shruivandana Sharma IOE 64, Supply Chain Managemen, Winer 2009 Universiy of Michigan, Ann

More information

Final Spring 2007

Final Spring 2007 .615 Final Spring 7 Overview The purpose of he final exam is o calculae he MHD β limi in a high-bea oroidal okamak agains he dangerous n = 1 exernal ballooning-kink mode. Effecively, his corresponds o

More information

CENTRALIZED VERSUS DECENTRALIZED PRODUCTION PLANNING IN SUPPLY CHAINS

CENTRALIZED VERSUS DECENTRALIZED PRODUCTION PLANNING IN SUPPLY CHAINS CENRALIZED VERSUS DECENRALIZED PRODUCION PLANNING IN SUPPLY CHAINS Georges SAHARIDIS* a, Yves DALLERY* a, Fikri KARAESMEN* b * a Ecole Cenrale Paris Deparmen of Indusial Engineering (LGI), +3343388, saharidis,dallery@lgi.ecp.fr

More information

Approximation algorithms for capacitated stochastic inventory systems with setup costs

Approximation algorithms for capacitated stochastic inventory systems with setup costs Approximaion algorihms for capaciaed sochasic invenory sysems wih seup coss The MIT Faculy has made his aricle openly available. Please share how his access benefis you. Your sory maers. Ciaion As Published

More information

Approximation Algorithms for Perishable Inventory Systems

Approximation Algorithms for Perishable Inventory Systems Approximaion Algorihms for Perishable Invenory Sysems Xiuli Chao Deparmen of Indusrial and Operaions Engineering, Universiy of Michigan, Ann Arbor, MI 48109, xchao@umich.edu Xiing Gong Deparmen of Sysems

More information

Diebold, Chapter 7. Francis X. Diebold, Elements of Forecasting, 4th Edition (Mason, Ohio: Cengage Learning, 2006). Chapter 7. Characterizing Cycles

Diebold, Chapter 7. Francis X. Diebold, Elements of Forecasting, 4th Edition (Mason, Ohio: Cengage Learning, 2006). Chapter 7. Characterizing Cycles Diebold, Chaper 7 Francis X. Diebold, Elemens of Forecasing, 4h Ediion (Mason, Ohio: Cengage Learning, 006). Chaper 7. Characerizing Cycles Afer compleing his reading you should be able o: Define covariance

More information

Application of a Stochastic-Fuzzy Approach to Modeling Optimal Discrete Time Dynamical Systems by Using Large Scale Data Processing

Application of a Stochastic-Fuzzy Approach to Modeling Optimal Discrete Time Dynamical Systems by Using Large Scale Data Processing Applicaion of a Sochasic-Fuzzy Approach o Modeling Opimal Discree Time Dynamical Sysems by Using Large Scale Daa Processing AA WALASZE-BABISZEWSA Deparmen of Compuer Engineering Opole Universiy of Technology

More information

Seminar 4: Hotelling 2

Seminar 4: Hotelling 2 Seminar 4: Hoelling 2 November 3, 211 1 Exercise Par 1 Iso-elasic demand A non renewable resource of a known sock S can be exraced a zero cos. Demand for he resource is of he form: D(p ) = p ε ε > A a

More information

20. Applications of the Genetic-Drift Model

20. Applications of the Genetic-Drift Model 0. Applicaions of he Geneic-Drif Model 1) Deermining he probabiliy of forming any paricular combinaion of genoypes in he nex generaion: Example: If he parenal allele frequencies are p 0 = 0.35 and q 0

More information

Macroeconomic Theory Ph.D. Qualifying Examination Fall 2005 ANSWER EACH PART IN A SEPARATE BLUE BOOK. PART ONE: ANSWER IN BOOK 1 WEIGHT 1/3

Macroeconomic Theory Ph.D. Qualifying Examination Fall 2005 ANSWER EACH PART IN A SEPARATE BLUE BOOK. PART ONE: ANSWER IN BOOK 1 WEIGHT 1/3 Macroeconomic Theory Ph.D. Qualifying Examinaion Fall 2005 Comprehensive Examinaion UCLA Dep. of Economics You have 4 hours o complee he exam. There are hree pars o he exam. Answer all pars. Each par has

More information

Deteriorating Inventory Model with Time. Dependent Demand and Partial Backlogging

Deteriorating Inventory Model with Time. Dependent Demand and Partial Backlogging Applied Mahemaical Sciences, Vol. 4, 00, no. 7, 36-369 Deerioraing Invenory Model wih Time Dependen Demand and Parial Backlogging Vinod Kumar Mishra Deparmen of Compuer Science & Engineering Kumaon Engineering

More information

Essential Microeconomics : OPTIMAL CONTROL 1. Consider the following class of optimization problems

Essential Microeconomics : OPTIMAL CONTROL 1. Consider the following class of optimization problems Essenial Microeconomics -- 6.5: OPIMAL CONROL Consider he following class of opimizaion problems Max{ U( k, x) + U+ ( k+ ) k+ k F( k, x)}. { x, k+ } = In he language of conrol heory, he vecor k is he vecor

More information

Lecture 4 Notes (Little s Theorem)

Lecture 4 Notes (Little s Theorem) Lecure 4 Noes (Lile s Theorem) This lecure concerns one of he mos imporan (and simples) heorems in Queuing Theory, Lile s Theorem. More informaion can be found in he course book, Bersekas & Gallagher,

More information

Errata (1 st Edition)

Errata (1 st Edition) P Sandborn, os Analysis of Elecronic Sysems, s Ediion, orld Scienific, Singapore, 03 Erraa ( s Ediion) S K 05D Page 8 Equaion (7) should be, E 05D E Nu e S K he L appearing in he equaion in he book does

More information

STATE-SPACE MODELLING. A mass balance across the tank gives:

STATE-SPACE MODELLING. A mass balance across the tank gives: B. Lennox and N.F. Thornhill, 9, Sae Space Modelling, IChemE Process Managemen and Conrol Subjec Group Newsleer STE-SPACE MODELLING Inroducion: Over he pas decade or so here has been an ever increasing

More information

An introduction to the theory of SDDP algorithm

An introduction to the theory of SDDP algorithm An inroducion o he heory of SDDP algorihm V. Leclère (ENPC) Augus 1, 2014 V. Leclère Inroducion o SDDP Augus 1, 2014 1 / 21 Inroducion Large scale sochasic problem are hard o solve. Two ways of aacking

More information

Online Appendix to Solution Methods for Models with Rare Disasters

Online Appendix to Solution Methods for Models with Rare Disasters Online Appendix o Soluion Mehods for Models wih Rare Disasers Jesús Fernández-Villaverde and Oren Levinal In his Online Appendix, we presen he Euler condiions of he model, we develop he pricing Calvo block,

More information

Chapter 2. First Order Scalar Equations

Chapter 2. First Order Scalar Equations Chaper. Firs Order Scalar Equaions We sar our sudy of differenial equaions in he same way he pioneers in his field did. We show paricular echniques o solve paricular ypes of firs order differenial equaions.

More information

SPH3U: Projectiles. Recorder: Manager: Speaker:

SPH3U: Projectiles. Recorder: Manager: Speaker: SPH3U: Projeciles Now i s ime o use our new skills o analyze he moion of a golf ball ha was ossed hrough he air. Le s find ou wha is special abou he moion of a projecile. Recorder: Manager: Speaker: 0

More information

This document was generated at 1:04 PM, 09/10/13 Copyright 2013 Richard T. Woodward. 4. End points and transversality conditions AGEC

This document was generated at 1:04 PM, 09/10/13 Copyright 2013 Richard T. Woodward. 4. End points and transversality conditions AGEC his documen was generaed a 1:4 PM, 9/1/13 Copyrigh 213 Richard. Woodward 4. End poins and ransversaliy condiions AGEC 637-213 F z d Recall from Lecure 3 ha a ypical opimal conrol problem is o maimize (,,

More information

Lecture 2-1 Kinematics in One Dimension Displacement, Velocity and Acceleration Everything in the world is moving. Nothing stays still.

Lecture 2-1 Kinematics in One Dimension Displacement, Velocity and Acceleration Everything in the world is moving. Nothing stays still. Lecure - Kinemaics in One Dimension Displacemen, Velociy and Acceleraion Everyhing in he world is moving. Nohing says sill. Moion occurs a all scales of he universe, saring from he moion of elecrons in

More information

Biol. 356 Lab 8. Mortality, Recruitment, and Migration Rates

Biol. 356 Lab 8. Mortality, Recruitment, and Migration Rates Biol. 356 Lab 8. Moraliy, Recruimen, and Migraion Raes (modified from Cox, 00, General Ecology Lab Manual, McGraw Hill) Las week we esimaed populaion size hrough several mehods. One assumpion of all hese

More information

On a Discrete-In-Time Order Level Inventory Model for Items with Random Deterioration

On a Discrete-In-Time Order Level Inventory Model for Items with Random Deterioration Journal of Agriculure and Life Sciences Vol., No. ; June 4 On a Discree-In-Time Order Level Invenory Model for Iems wih Random Deerioraion Dr Biswaranjan Mandal Associae Professor of Mahemaics Acharya

More information

5. Stochastic processes (1)

5. Stochastic processes (1) Lec05.pp S-38.45 - Inroducion o Teleraffic Theory Spring 2005 Conens Basic conceps Poisson process 2 Sochasic processes () Consider some quaniy in a eleraffic (or any) sysem I ypically evolves in ime randomly

More information

Comparing Means: t-tests for One Sample & Two Related Samples

Comparing Means: t-tests for One Sample & Two Related Samples Comparing Means: -Tess for One Sample & Two Relaed Samples Using he z-tes: Assumpions -Tess for One Sample & Two Relaed Samples The z-es (of a sample mean agains a populaion mean) is based on he assumpion

More information

A Study of Inventory System with Ramp Type Demand Rate and Shortage in The Light Of Inflation I

A Study of Inventory System with Ramp Type Demand Rate and Shortage in The Light Of Inflation I Inernaional Journal of Mahemaics rends and echnology Volume 7 Number Jan 5 A Sudy of Invenory Sysem wih Ramp ype emand Rae and Shorage in he Ligh Of Inflaion I Sangeea Gupa, R.K. Srivasava, A.K. Singh

More information

Maintenance Models. Prof. Robert C. Leachman IEOR 130, Methods of Manufacturing Improvement Spring, 2011

Maintenance Models. Prof. Robert C. Leachman IEOR 130, Methods of Manufacturing Improvement Spring, 2011 Mainenance Models Prof Rober C Leachman IEOR 3, Mehods of Manufacuring Improvemen Spring, Inroducion The mainenance of complex equipmen ofen accouns for a large porion of he coss associaed wih ha equipmen

More information

23.2. Representing Periodic Functions by Fourier Series. Introduction. Prerequisites. Learning Outcomes

23.2. Representing Periodic Functions by Fourier Series. Introduction. Prerequisites. Learning Outcomes Represening Periodic Funcions by Fourier Series 3. Inroducion In his Secion we show how a periodic funcion can be expressed as a series of sines and cosines. We begin by obaining some sandard inegrals

More information

SZG Macro 2011 Lecture 3: Dynamic Programming. SZG macro 2011 lecture 3 1

SZG Macro 2011 Lecture 3: Dynamic Programming. SZG macro 2011 lecture 3 1 SZG Macro 2011 Lecure 3: Dynamic Programming SZG macro 2011 lecure 3 1 Background Our previous discussion of opimal consumpion over ime and of opimal capial accumulaion sugges sudying he general decision

More information

Information Relaxations and Duality in Stochastic Dynamic Programs

Information Relaxations and Duality in Stochastic Dynamic Programs OPERATIONS RESEARCH Vol. 58, No. 4, Par 1 of 2, July Augus 2010, pp. 785 801 issn 0030-364X eissn 1526-5463 10 5804 0785 informs doi 10.1287/opre.1090.0796 2010 INFORMS Informaion Relaxaions and Dualiy

More information

Appendix 14.1 The optimal control problem and its solution using

Appendix 14.1 The optimal control problem and its solution using 1 Appendix 14.1 he opimal conrol problem and is soluion using he maximum principle NOE: Many occurrences of f, x, u, and in his file (in equaions or as whole words in ex) are purposefully in bold in order

More information

E β t log (C t ) + M t M t 1. = Y t + B t 1 P t. B t 0 (3) v t = P tc t M t Question 1. Find the FOC s for an optimum in the agent s problem.

E β t log (C t ) + M t M t 1. = Y t + B t 1 P t. B t 0 (3) v t = P tc t M t Question 1. Find the FOC s for an optimum in the agent s problem. Noes, M. Krause.. Problem Se 9: Exercise on FTPL Same model as in paper and lecure, only ha one-period govenmen bonds are replaced by consols, which are bonds ha pay one dollar forever. I has curren marke

More information

Robust estimation based on the first- and third-moment restrictions of the power transformation model

Robust estimation based on the first- and third-moment restrictions of the power transformation model h Inernaional Congress on Modelling and Simulaion, Adelaide, Ausralia, 6 December 3 www.mssanz.org.au/modsim3 Robus esimaion based on he firs- and hird-momen resricions of he power ransformaion Nawaa,

More information

Let us start with a two dimensional case. We consider a vector ( x,

Let us start with a two dimensional case. We consider a vector ( x, Roaion marices We consider now roaion marices in wo and hree dimensions. We sar wih wo dimensions since wo dimensions are easier han hree o undersand, and one dimension is a lile oo simple. However, our

More information

Two Popular Bayesian Estimators: Particle and Kalman Filters. McGill COMP 765 Sept 14 th, 2017

Two Popular Bayesian Estimators: Particle and Kalman Filters. McGill COMP 765 Sept 14 th, 2017 Two Popular Bayesian Esimaors: Paricle and Kalman Filers McGill COMP 765 Sep 14 h, 2017 1 1 1, dx x Bel x u x P x z P Recall: Bayes Filers,,,,,,, 1 1 1 1 u z u x P u z u x z P Bayes z = observaion u =

More information

On the Optimal Policy Structure in Serial Inventory Systems with Lost Sales

On the Optimal Policy Structure in Serial Inventory Systems with Lost Sales On he Opimal Policy Srucure in Serial Invenory Sysems wih Los Sales Woonghee Tim Huh, Columbia Universiy Ganesh Janakiraman, New York Universiy May 21, 2008 Revised: July 30, 2008; December 23, 2008 Absrac

More information

Problem 1 / 25 Problem 2 / 20 Problem 3 / 10 Problem 4 / 15 Problem 5 / 30 TOTAL / 100

Problem 1 / 25 Problem 2 / 20 Problem 3 / 10 Problem 4 / 15 Problem 5 / 30 TOTAL / 100 eparmen of Applied Economics Johns Hopkins Universiy Economics 602 Macroeconomic Theory and Policy Miderm Exam Suggesed Soluions Professor Sanjay hugh Fall 2008 NAME: The Exam has a oal of five (5) problems

More information

d 1 = c 1 b 2 - b 1 c 2 d 2 = c 1 b 3 - b 1 c 3

d 1 = c 1 b 2 - b 1 c 2 d 2 = c 1 b 3 - b 1 c 3 and d = c b - b c c d = c b - b c c This process is coninued unil he nh row has been compleed. The complee array of coefficiens is riangular. Noe ha in developing he array an enire row may be divided or

More information

Linear Response Theory: The connection between QFT and experiments

Linear Response Theory: The connection between QFT and experiments Phys540.nb 39 3 Linear Response Theory: The connecion beween QFT and experimens 3.1. Basic conceps and ideas Q: How do we measure he conduciviy of a meal? A: we firs inroduce a weak elecric field E, and

More information

t is a basis for the solution space to this system, then the matrix having these solutions as columns, t x 1 t, x 2 t,... x n t x 2 t...

t is a basis for the solution space to this system, then the matrix having these solutions as columns, t x 1 t, x 2 t,... x n t x 2 t... Mah 228- Fri Mar 24 5.6 Marix exponenials and linear sysems: The analogy beween firs order sysems of linear differenial equaions (Chaper 5) and scalar linear differenial equaions (Chaper ) is much sronger

More information

Lecture 20: Riccati Equations and Least Squares Feedback Control

Lecture 20: Riccati Equations and Least Squares Feedback Control 34-5 LINEAR SYSTEMS Lecure : Riccai Equaions and Leas Squares Feedback Conrol 5.6.4 Sae Feedback via Riccai Equaions A recursive approach in generaing he marix-valued funcion W ( ) equaion for i for he

More information

Solutions Problem Set 3 Macro II (14.452)

Solutions Problem Set 3 Macro II (14.452) Soluions Problem Se 3 Macro II (14.452) Francisco A. Gallego 04/27/2005 1 Q heory of invesmen in coninuous ime and no uncerainy Consider he in nie horizon model of a rm facing adjusmen coss o invesmen.

More information

A Dynamic Model of Economic Fluctuations

A Dynamic Model of Economic Fluctuations CHAPTER 15 A Dynamic Model of Economic Flucuaions Modified for ECON 2204 by Bob Murphy 2016 Worh Publishers, all righs reserved IN THIS CHAPTER, OU WILL LEARN: how o incorporae dynamics ino he AD-AS model

More information

Asymptotic Optimality of Order-up-to Policies in Lost Sales Inventory Systems

Asymptotic Optimality of Order-up-to Policies in Lost Sales Inventory Systems Asympoic Opimaliy of Order-up-o Policies in Los Sales Invenory Sysems Woonghee Tim Huh 1, Ganesh Janakiraman 2, John A. Mucksad 3, Paa Rusmevichienong 4 December 4, 2006 Subjec Classificaion: Invenory/Producion:

More information

not to be republished NCERT MATHEMATICAL MODELLING Appendix 2 A.2.1 Introduction A.2.2 Why Mathematical Modelling?

not to be republished NCERT MATHEMATICAL MODELLING Appendix 2 A.2.1 Introduction A.2.2 Why Mathematical Modelling? 256 MATHEMATICS A.2.1 Inroducion In class XI, we have learn abou mahemaical modelling as an aemp o sudy some par (or form) of some real-life problems in mahemaical erms, i.e., he conversion of a physical

More information

Two Coupled Oscillators / Normal Modes

Two Coupled Oscillators / Normal Modes Lecure 3 Phys 3750 Two Coupled Oscillaors / Normal Modes Overview and Moivaion: Today we ake a small, bu significan, sep owards wave moion. We will no ye observe waves, bu his sep is imporan in is own

More information

Problem Set 5. Graduate Macro II, Spring 2017 The University of Notre Dame Professor Sims

Problem Set 5. Graduate Macro II, Spring 2017 The University of Notre Dame Professor Sims Problem Se 5 Graduae Macro II, Spring 2017 The Universiy of Nore Dame Professor Sims Insrucions: You may consul wih oher members of he class, bu please make sure o urn in your own work. Where applicable,

More information

Robotics I. April 11, The kinematics of a 3R spatial robot is specified by the Denavit-Hartenberg parameters in Tab. 1.

Robotics I. April 11, The kinematics of a 3R spatial robot is specified by the Denavit-Hartenberg parameters in Tab. 1. Roboics I April 11, 017 Exercise 1 he kinemaics of a 3R spaial robo is specified by he Denavi-Harenberg parameers in ab 1 i α i d i a i θ i 1 π/ L 1 0 1 0 0 L 3 0 0 L 3 3 able 1: able of DH parameers of

More information

1. An introduction to dynamic optimization -- Optimal Control and Dynamic Programming AGEC

1. An introduction to dynamic optimization -- Optimal Control and Dynamic Programming AGEC This documen was generaed a :45 PM 8/8/04 Copyrigh 04 Richard T. Woodward. An inroducion o dynamic opimizaion -- Opimal Conrol and Dynamic Programming AGEC 637-04 I. Overview of opimizaion Opimizaion is

More information

International Journal of Scientific & Engineering Research, Volume 4, Issue 10, October ISSN

International Journal of Scientific & Engineering Research, Volume 4, Issue 10, October ISSN Inernaional Journal of Scienific & Engineering Research, Volume 4, Issue 10, Ocober-2013 900 FUZZY MEAN RESIDUAL LIFE ORDERING OF FUZZY RANDOM VARIABLES J. EARNEST LAZARUS PIRIYAKUMAR 1, A. YAMUNA 2 1.

More information

Class Meeting # 10: Introduction to the Wave Equation

Class Meeting # 10: Introduction to the Wave Equation MATH 8.5 COURSE NOTES - CLASS MEETING # 0 8.5 Inroducion o PDEs, Fall 0 Professor: Jared Speck Class Meeing # 0: Inroducion o he Wave Equaion. Wha is he wave equaion? The sandard wave equaion for a funcion

More information

Final Exam Advanced Macroeconomics I

Final Exam Advanced Macroeconomics I Advanced Macroeconomics I WS 00/ Final Exam Advanced Macroeconomics I February 8, 0 Quesion (5%) An economy produces oupu according o α α Y = K (AL) of which a fracion s is invesed. echnology A is exogenous

More information

Matlab and Python programming: how to get started

Matlab and Python programming: how to get started Malab and Pyhon programming: how o ge sared Equipping readers he skills o wrie programs o explore complex sysems and discover ineresing paerns from big daa is one of he main goals of his book. In his chaper,

More information

Stability and Bifurcation in a Neural Network Model with Two Delays

Stability and Bifurcation in a Neural Network Model with Two Delays Inernaional Mahemaical Forum, Vol. 6, 11, no. 35, 175-1731 Sabiliy and Bifurcaion in a Neural Nework Model wih Two Delays GuangPing Hu and XiaoLing Li School of Mahemaics and Physics, Nanjing Universiy

More information

Presentation Overview

Presentation Overview Acion Refinemen in Reinforcemen Learning by Probabiliy Smoohing By Thomas G. Dieerich & Didac Busques Speaer: Kai Xu Presenaion Overview Bacground The Probabiliy Smoohing Mehod Experimenal Sudy of Acion

More information

14 Autoregressive Moving Average Models

14 Autoregressive Moving Average Models 14 Auoregressive Moving Average Models In his chaper an imporan parameric family of saionary ime series is inroduced, he family of he auoregressive moving average, or ARMA, processes. For a large class

More information

A Hop Constrained Min-Sum Arborescence with Outage Costs

A Hop Constrained Min-Sum Arborescence with Outage Costs A Hop Consrained Min-Sum Arborescence wih Ouage Coss Rakesh Kawara Minnesoa Sae Universiy, Mankao, MN 56001 Email: Kawara@mnsu.edu Absrac The hop consrained min-sum arborescence wih ouage coss problem

More information

Module 2 F c i k c s la l w a s o s f dif di fusi s o i n

Module 2 F c i k c s la l w a s o s f dif di fusi s o i n Module Fick s laws of diffusion Fick s laws of diffusion and hin film soluion Adolf Fick (1855) proposed: d J α d d d J (mole/m s) flu (m /s) diffusion coefficien and (mole/m 3 ) concenraion of ions, aoms

More information

On Measuring Pro-Poor Growth. 1. On Various Ways of Measuring Pro-Poor Growth: A Short Review of the Literature

On Measuring Pro-Poor Growth. 1. On Various Ways of Measuring Pro-Poor Growth: A Short Review of the Literature On Measuring Pro-Poor Growh 1. On Various Ways of Measuring Pro-Poor Growh: A Shor eview of he Lieraure During he pas en years or so here have been various suggesions concerning he way one should check

More information

BU Macro BU Macro Fall 2008, Lecture 4

BU Macro BU Macro Fall 2008, Lecture 4 Dynamic Programming BU Macro 2008 Lecure 4 1 Ouline 1. Cerainy opimizaion problem used o illusrae: a. Resricions on exogenous variables b. Value funcion c. Policy funcion d. The Bellman equaion and an

More information

Christos Papadimitriou & Luca Trevisan November 22, 2016

Christos Papadimitriou & Luca Trevisan November 22, 2016 U.C. Bereley CS170: Algorihms Handou LN-11-22 Chrisos Papadimiriou & Luca Trevisan November 22, 2016 Sreaming algorihms In his lecure and he nex one we sudy memory-efficien algorihms ha process a sream

More information

Two New Uncertainty Programming Models of Inventory with Uncertain Costs

Two New Uncertainty Programming Models of Inventory with Uncertain Costs Journal of Informaion & Compuaional Science 8: 2 (211) 28 288 Available a hp://www.joics.com Two New Uncerainy Programming Models of Invenory wih Uncerain Coss Lixia Rong Compuer Science and Technology

More information

Approximation Algorithms for Perishable Inventory Systems with Setup Costs

Approximation Algorithms for Perishable Inventory Systems with Setup Costs Approximaion Algorihms for Perishable Invenory Sysems wih Seup Coss Huanan Zhang, Cong Shi, Xiuli Chao Indusrial and Operaions Engineering, Universiy of Michigan, MI 48105, {zhanghn, shicong, xchao}@umich.edu

More information

Simulation-Solving Dynamic Models ABE 5646 Week 2, Spring 2010

Simulation-Solving Dynamic Models ABE 5646 Week 2, Spring 2010 Simulaion-Solving Dynamic Models ABE 5646 Week 2, Spring 2010 Week Descripion Reading Maerial 2 Compuer Simulaion of Dynamic Models Finie Difference, coninuous saes, discree ime Simple Mehods Euler Trapezoid

More information

Numerical Dispersion

Numerical Dispersion eview of Linear Numerical Sabiliy Numerical Dispersion n he previous lecure, we considered he linear numerical sabiliy of boh advecion and diffusion erms when approimaed wih several spaial and emporal

More information

Logic in computer science

Logic in computer science Logic in compuer science Logic plays an imporan role in compuer science Logic is ofen called he calculus of compuer science Logic plays a similar role in compuer science o ha played by calculus in he physical

More information

ACE 562 Fall Lecture 5: The Simple Linear Regression Model: Sampling Properties of the Least Squares Estimators. by Professor Scott H.

ACE 562 Fall Lecture 5: The Simple Linear Regression Model: Sampling Properties of the Least Squares Estimators. by Professor Scott H. ACE 56 Fall 005 Lecure 5: he Simple Linear Regression Model: Sampling Properies of he Leas Squares Esimaors by Professor Sco H. Irwin Required Reading: Griffihs, Hill and Judge. "Inference in he Simple

More information

Basic definitions and relations

Basic definitions and relations Basic definiions and relaions Lecurer: Dmiri A. Molchanov E-mail: molchan@cs.u.fi hp://www.cs.u.fi/kurssi/tlt-2716/ Kendall s noaion for queuing sysems: Arrival processes; Service ime disribuions; Examples.

More information

SUPPLEMENTARY INFORMATION

SUPPLEMENTARY INFORMATION SUPPLEMENTARY INFORMATION DOI: 0.038/NCLIMATE893 Temporal resoluion and DICE * Supplemenal Informaion Alex L. Maren and Sephen C. Newbold Naional Cener for Environmenal Economics, US Environmenal Proecion

More information

Single-Pass-Based Heuristic Algorithms for Group Flexible Flow-shop Scheduling Problems

Single-Pass-Based Heuristic Algorithms for Group Flexible Flow-shop Scheduling Problems Single-Pass-Based Heurisic Algorihms for Group Flexible Flow-shop Scheduling Problems PEI-YING HUANG, TZUNG-PEI HONG 2 and CHENG-YAN KAO, 3 Deparmen of Compuer Science and Informaion Engineering Naional

More information

Modal identification of structures from roving input data by means of maximum likelihood estimation of the state space model

Modal identification of structures from roving input data by means of maximum likelihood estimation of the state space model Modal idenificaion of srucures from roving inpu daa by means of maximum likelihood esimaion of he sae space model J. Cara, J. Juan, E. Alarcón Absrac The usual way o perform a forced vibraion es is o fix

More information

A New Perturbative Approach in Nonlinear Singularity Analysis

A New Perturbative Approach in Nonlinear Singularity Analysis Journal of Mahemaics and Saisics 7 (: 49-54, ISSN 549-644 Science Publicaions A New Perurbaive Approach in Nonlinear Singulariy Analysis Ta-Leung Yee Deparmen of Mahemaics and Informaion Technology, The

More information

Predator - Prey Model Trajectories and the nonlinear conservation law

Predator - Prey Model Trajectories and the nonlinear conservation law Predaor - Prey Model Trajecories and he nonlinear conservaion law James K. Peerson Deparmen of Biological Sciences and Deparmen of Mahemaical Sciences Clemson Universiy Ocober 28, 213 Ouline Drawing Trajecories

More information

Notes on Kalman Filtering

Notes on Kalman Filtering Noes on Kalman Filering Brian Borchers and Rick Aser November 7, Inroducion Daa Assimilaion is he problem of merging model predicions wih acual measuremens of a sysem o produce an opimal esimae of he curren

More information

ODEs II, Lecture 1: Homogeneous Linear Systems - I. Mike Raugh 1. March 8, 2004

ODEs II, Lecture 1: Homogeneous Linear Systems - I. Mike Raugh 1. March 8, 2004 ODEs II, Lecure : Homogeneous Linear Sysems - I Mike Raugh March 8, 4 Inroducion. In he firs lecure we discussed a sysem of linear ODEs for modeling he excreion of lead from he human body, saw how o ransform

More information

Random Walk with Anti-Correlated Steps

Random Walk with Anti-Correlated Steps Random Walk wih Ani-Correlaed Seps John Noga Dirk Wagner 2 Absrac We conjecure he expeced value of random walks wih ani-correlaed seps o be exacly. We suppor his conjecure wih 2 plausibiliy argumens and

More information

Guest Lectures for Dr. MacFarlane s EE3350 Part Deux

Guest Lectures for Dr. MacFarlane s EE3350 Part Deux Gues Lecures for Dr. MacFarlane s EE3350 Par Deux Michael Plane Mon., 08-30-2010 Wrie name in corner. Poin ou his is a review, so I will go faser. Remind hem o go lisen o online lecure abou geing an A

More information

Sliding Mode Extremum Seeking Control for Linear Quadratic Dynamic Game

Sliding Mode Extremum Seeking Control for Linear Quadratic Dynamic Game Sliding Mode Exremum Seeking Conrol for Linear Quadraic Dynamic Game Yaodong Pan and Ümi Özgüner ITS Research Group, AIST Tsukuba Eas Namiki --, Tsukuba-shi,Ibaraki-ken 5-856, Japan e-mail: pan.yaodong@ais.go.jp

More information

Lecture Notes 2. The Hilbert Space Approach to Time Series

Lecture Notes 2. The Hilbert Space Approach to Time Series Time Series Seven N. Durlauf Universiy of Wisconsin. Basic ideas Lecure Noes. The Hilber Space Approach o Time Series The Hilber space framework provides a very powerful language for discussing he relaionship

More information

KINEMATICS IN ONE DIMENSION

KINEMATICS IN ONE DIMENSION KINEMATICS IN ONE DIMENSION PREVIEW Kinemaics is he sudy of how hings move how far (disance and displacemen), how fas (speed and velociy), and how fas ha how fas changes (acceleraion). We say ha an objec

More information

12: AUTOREGRESSIVE AND MOVING AVERAGE PROCESSES IN DISCRETE TIME. Σ j =

12: AUTOREGRESSIVE AND MOVING AVERAGE PROCESSES IN DISCRETE TIME. Σ j = 1: AUTOREGRESSIVE AND MOVING AVERAGE PROCESSES IN DISCRETE TIME Moving Averages Recall ha a whie noise process is a series { } = having variance σ. The whie noise process has specral densiy f (λ) = of

More information

Chapter Floating Point Representation

Chapter Floating Point Representation Chaper 01.05 Floaing Poin Represenaion Afer reading his chaper, you should be able o: 1. conver a base- number o a binary floaing poin represenaion,. conver a binary floaing poin number o is equivalen

More information

Learning a Class from Examples. Training set X. Class C 1. Class C of a family car. Output: Input representation: x 1 : price, x 2 : engine power

Learning a Class from Examples. Training set X. Class C 1. Class C of a family car. Output: Input representation: x 1 : price, x 2 : engine power Alpaydin Chaper, Michell Chaper 7 Alpaydin slides are in urquoise. Ehem Alpaydin, copyrigh: The MIT Press, 010. alpaydin@boun.edu.r hp://www.cmpe.boun.edu.r/ ehem/imle All oher slides are based on Michell.

More information

5.1 - Logarithms and Their Properties

5.1 - Logarithms and Their Properties Chaper 5 Logarihmic Funcions 5.1 - Logarihms and Their Properies Suppose ha a populaion grows according o he formula P 10, where P is he colony size a ime, in hours. When will he populaion be 2500? We

More information

On Multicomponent System Reliability with Microshocks - Microdamages Type of Components Interaction

On Multicomponent System Reliability with Microshocks - Microdamages Type of Components Interaction On Mulicomponen Sysem Reliabiliy wih Microshocks - Microdamages Type of Componens Ineracion Jerzy K. Filus, and Lidia Z. Filus Absrac Consider a wo componen parallel sysem. The defined new sochasic dependences

More information

2. Nonlinear Conservation Law Equations

2. Nonlinear Conservation Law Equations . Nonlinear Conservaion Law Equaions One of he clear lessons learned over recen years in sudying nonlinear parial differenial equaions is ha i is generally no wise o ry o aack a general class of nonlinear

More information

ACE 564 Spring Lecture 7. Extensions of The Multiple Regression Model: Dummy Independent Variables. by Professor Scott H.

ACE 564 Spring Lecture 7. Extensions of The Multiple Regression Model: Dummy Independent Variables. by Professor Scott H. ACE 564 Spring 2006 Lecure 7 Exensions of The Muliple Regression Model: Dumm Independen Variables b Professor Sco H. Irwin Readings: Griffihs, Hill and Judge. "Dumm Variables and Varing Coefficien Models

More information

1 Review of Zero-Sum Games

1 Review of Zero-Sum Games COS 5: heoreical Machine Learning Lecurer: Rob Schapire Lecure #23 Scribe: Eugene Brevdo April 30, 2008 Review of Zero-Sum Games Las ime we inroduced a mahemaical model for wo player zero-sum games. Any

More information