MULTI-ITEM PRODUCTION PLANNING AND MANAGEMENT SYSTEM BASED ON UNFULFILLED ORDER RATE IN SUPPLY CHAIN

Size: px
Start display at page:

Download "MULTI-ITEM PRODUCTION PLANNING AND MANAGEMENT SYSTEM BASED ON UNFULFILLED ORDER RATE IN SUPPLY CHAIN"

Transcription

1 Journal of the Operatons Research Socety of Japan 2007, Vol. 50, No. 3, MULTI-ITEM PRODUCTION PLANNING AND MANAGEMENT SYSTEM BASED ON UNFULFILLED ORDER RATE IN SUPPLY CHAIN Nobuyuk Ueno Koj Okuhara Hroak Ish Hroak Shbuk Toshak Kuramoto Prefectural Unversty Osaka Unversty Mazda Motor YNS INC. of Hroshma Corporaton (Receved May 26, 2005; Revsed September 19, 2006) Abstract In the automoble ndustry, a usual busness model has a problem to realze mass customzaton, because t s dffcult to satsfy the dversfed customer needs. Ths paper proposes a mult-tem producton and nventory plannng method of the mass customzaton wth the consderaton of the restrcton of daly manufacturng capacty and so on. Ths model s formulated as a stochastc programmng problem, and then the sub problem as a lnear programmng problem. An effcent and practcal algorthm for the mult-tem model s developed. Keywords: Inventory, producton plannng, mass-customzaton, mathematcal programmng, stochastc model 1. Introducton In recent years, the automoble ndustry has pushed forward wth reducton n cost, nducton of foregn captal, competton between supplers, and so on. It becomes very mportant to remove waste of the producton actvtes and to meet the demands of an ndvdual customer and a changng market. In partcular, we must satsfy the varety of customer specfcaton n product and servce, wthout droppng the productve effcency n mass customzaton [1,5,10]. For example, n case of a certan car manufacturer, t s sad that there s possblty of 300,000,000 ways of specfcaton for a customzed product n one model of a car. On manufacturng of ths customzed product, levelng of producton load becomes very dffcult because manufacturers have to meet orders of customers based on each specfcaton. Then the productvty decreases, fnshed goods n stock ncrease, and t becomes dffcult to deal wth customer specfcaton. There have already been some manufacturng strateges for mass customzaton. However, the desgn method of producton plannng and management system for t has not yet been establshed. Under the precondton that delvery lead tme whch s expected by customers s longer than producton lead tme whch s necessary for manufacturers, Make-to-Order management system n whch manufacturng starts after recevng an order has been appled for a varety of customer specfcatons as producton plannng and management system for mass customzaton. However, t s often the case that delvery lead tme becomes shorter than producton lead tme. Therefore manufacturers have to start manufacturng before they receve an order from customers. For a varety of customer specfcatons, Make-to-Order management system and Parts Orented Producton System (POPS) type module manufacturng system are proposed [9]. 201

2 202 N. Ueno, K. Okuhara, H. Ish, H. Shbuk & T. Kuramoto The Materal Requrement Plannng (MRP) [11] and Advanced Plannng & Schedulng (APS) [7] are presented to plan and manage such systems. Ths MRP satsfes varous demands from customers by promotng modularzaton n producton. By makng use of modularzaton, the MRP makes producton lead tme shorter, and t manages Make-to- Order management system on precondton that delvery lead tme s longer than producton lead tme. When delvery lead tme s shorter than producton lead tme, we consder the mnmum stock as necessary beforehand. But n these management systems, the mechansm between substtutablty of unfulflled order and stocks on mass customzaton has not been dscussed quanttatvely. On precondton that delvery lead tme s shorter than producton lead tme, manufacture seat system for both prospectve stocks and order stocks has been proposed [13]. A study on the analyss of manufacture seat system coverage has also been done [6]. However, t has been ponted out that further theoretcal work about the decson method of proper manufacture seat should be done [13]. In ths paper, we take up producton plannng on the suppler sde [4,14] among manufacturer and assembly supplers n automoble ndustry. On precondton that delvery lead tme s shorter than producton lead tme, our targets are both prospectve stocks and order stocks. Our producton plannng and management system does nether dscrmnate between them nor depend on producton seat. We regard mass customzaton as fluctuaton of order quanttes from manufacturer to suppler, and we deal wth order quanttes from manufacturer, that s, the demand for suppler, stochastcally. At frst, we defne the producton plannng and management system as the stochastc programmng problem [2] to fnd producton plan mnmzng manufacture and nventory cost under both unfulflled order constrant and producton constrant. Next we show how to get practcal and effectve algorthm to the problem. Fnally, we descrbe a desgn procedure of the whole producton plannng and management system for mplementng mass customzaton. 2. Mass Customzaton Envronment The conceptual fgure of collaboraton model between manufacturer and supplers for mplementng mass customzaton for each tem (part) s shown n Fgure 1. The example n Fgure 1 shows that there are 3 types of order nformaton (forecast order) from manufacturer to suppler. The monthly forecast order gves the prospected order 3 months before, the weekly forecast order gves the predcton value 2 weeks before, and the delvery nstructon to suppler s gven as the frm order 1-3 days before delvery due date. In the mass customzaton envronment, the ncrease of customer specfcatons means the ncrease of producton specfcatons n manufacturer. Ths mples that the quanttes of frm order reflectng customer needs have large fluctuaton lke the bottom mage of Fgure 1. For suppler, the standard producton lead tme for knds of products (or parts) usually takes about 1 week. Therefore producton lead tme becomes longer than delvery lead tme. Suppler must start producton n advance based on producton plan by usng MRP, for example, accordng to weekly forecast order. Furthermore the suppler must perform the producton plannng whle at the same tme avodng unfulflled order to the frm order whch s gven 1-3 days before delvery due date. We regard mass customzaton as fluctuaton of order from customer to manufacturer and from manufacturer to suppler. Then t can be grasped from the suppler s pont of vew as the fluctuaton of order quanttes of target part from manufacturer to suppler. Ths

3 Plannng Based on Unfulflled Order Rate 203 Fgure 1: Future collaboraton model between manufacturer and suppler for mplementng mass customzaton Fgure 2: Dfference between MTO and ATO means that the fluctuaton of order quanttes can be represented by standard devaton around average that forecast order gves. Fgure 2 shows the dfference between make-toorder (MTO) and assemble-to-order (ATO) n envronment under mass customzaton of producton management system. We focus on ths case, and suppler should manufacture accordng to weekly forecast order. As s shown n Fgure 3, assumng that σ 0 denotes the order fluctuaton of concerned model of car from customer to manufacturer and σ denotes the order fluctuaton about concerned producton (or parts) from manufacturer to supplers. Then the producton plannng and management system mplementng mass customzaton n supplers sde can be consdered by the problem such as how to product and to stock ndvdual parts n advance n order not to run short of supples for σ whch gves fluctuaton of frm order from manufacturer. In ths paper, we assume that σ s gven and demand s defned as a normal dstrbuton wth tme varant average and tme varant devaton of order. We here focus on the suppler s producton plannng system, so the effect of σ only are dscussed n ths model. However t s mportant to nvestgate the mpact of σ 0 on the model, whch should be ncluded n the future research.

4 204 N. Ueno, K. Okuhara, H. Ish, H. Shbuk & T. Kuramoto Fgure 3: Flow of order among customer, manufacturer and suppler 3. Unfulflled Order Rate for Suppler s Producton Plannng and Management System 3.1. Dervaton of unfulflled order rate In mass customzaton envronment, we formulate problem whch determnes proper producton quanttes among n perods to optmze nventory change at fnal producton stage of suppler, correspondng to fluctuaton of demand from manufacturer. In ths secton, we dscuss about one certan tem. [Notaton] : Perod ( n). d : Frm order of manufacturer at perod. x : Producton quantty at perod of suppler. S : Inventory quantty at perod of suppler. p : Manufacturng cost per module at perod of suppler. h : Inventory holdng cost per module at perod of suppler. r : Total product quantty untl perod n of suppler. R : Set of lnear producton constrants of suppler. SO : Unfulflled order rate untl perod of suppler. β : Unfulflled order rate of planned target of suppler. It s assumed that d obeys an average d and standard devaton ω where d, d j ( j) are ndependent each other, and ω = σ d, where σ s the devaton of order. Let ntal nventory be S 0. R s denoted by (x 1, x 2,, x n ) R where R denotes lnear producton constrant, and s convex. SO s probablty functon whch s n short of delvery to the frm order at least untl perod. The nventory quantty S at perod s gven by S = S 0 + x t d t (3.1) t=1 t=1 where d s random varable so that S becomes random varable. It obeys normal dstrbuton whch has the followng tme varant average and varance. S can be replaced by Average m = S 0 + x t d t, (3.2) t=1 t=1 V arance σ 2 = ωt 2. (3.3) t=1 y = S m σ. (3.4)

5 Plannng Based on Unfulflled Order Rate 205 The probablty functon M whch satsfes nventory quantty S derved by 0 at perod s M = 0 1 e (S m ) 2σ 2 ds = 2πσ 2 1 m σ e y 2 2π 2 dy. (3.5) Thus () M = () M = 0.5 These equatons can be represented by m σ 0 m σ m σ M = sgn(m ) e y 2 2 dy (m 0), (3.6) 2π 1 e y 2 2 dy (m < 0). (3.7) 2π 1 e y 2 2 dy. (3.8) 2π We consder an unfulflled order rate SO n whch represents probablty that S 1 0, S 2 0, S 3 0,, S n 0 s not satsfed. Then, the upper bound of the unfulflled order rate SO n s gven by 1 n M t. Snce t s dffcult to obtan the true value of the unfulflled t=1 order rate, we defne the unfulflled order rate SO n by n SO n = 1 M t. (3.9) t=1 By consderng such unfulflled order rate, we can construct the producton plannng on safety sde. In addton, SO n can be rewrtten by applyng the ntegraton by parts to M t as follows; SO n = 1 A (3.10) where A = n t=1 1 2 ( ) 2 m t σ t sgn(m t ) e 2π k=1 m t σ t 2k+1. (3.11) (2k + 1) 3.2. Property of unfulflled order rate We clarfy the property of unfulflled order rate to understand effect of some varables. Lemma 1 SO n s monotonous decrease functon of x. Lemma 2 SO n s monotonous decrease functon of S 0. Lemma 3 SO n s monotonous ncrease functon of d. Lemma 4 SO n s monotonous ncrease functon of ω n m > 0 ( n). Proof We apply partal dervatve SO n wth respect to x, S 0, d and ω. Let α be arbtrary varable. SO n α = M 1 α M 2 M n M 1 M 2 α M 3 M n M 1 M 2 M n 1 M n α. (3.12)

6 206 N. Ueno, K. Okuhara, H. Ish, H. Shbuk & T. Kuramoto Snce M > 0 and f(y )dy = f(g(α)) g(α) α g(α) α, (3.13) by substtutng α nto x, S 0, d and ω, we can derve the followng propertes. M α = 1 2πσ e 1 2 M α = 1 2πσ e 1 2 M α = 1 2π e 1 2 ( m σ ) 2 > 0 (α x, S 0 ), (3.14) ( m σ ) 2 < 0 (α d ), (3.15) ( ) 2 m σ m ω < 0 (α ω σ 2, m > 0). (3.16) σ Equatons (3.14) (3.16) and M > 0 ( n) lead lemmas Producton Plannng for Implementng Mass Customzaton wth Mult-tem 4.1. Problem formulaton of mass customzaton wth mult-tem In ths secton, we defne suppler s producton plannng and management system as the stochastc programmng problem to fnd producton plan mnmzng manufacturng and nventory holdng cost under both unfulflled order rate constrant and producton constrant. In order to consder the mult-tem case, prevous notatons are extended by addng nformaton about jth part of product (j m). For example, d j, x j and S j denote the frm order about jth part of manufacturer at perod wth tme-varant average d j, the producton quantty of jth part at perod and nventory quantty of jth part at perod respectvely. Also p j and h j denote manufacturng cost per module about jth part of suppler at perod and nventory holdng cost per module about jth part of suppler at perod, respectvely. Those parameters, p j and h j, wll be regarded as the same value about each part wthout losng generalty and practcalty n order to smplfy a transformaton of objectve functon. And S j 0 denotes ntal nventory quantty of jth part, β j s targeted unfulflled order rate of jth part, and Q denotes upper bound of total product quantty at perod. [Formula 1] mnmze m n n E[ { p j x j + h j S j }] (4.1) j=1 =1 =1 s.t. S j 0 + x j t d j t 0 (, j) (4.2) t=1 t=1 n x j = r j ( j) (4.3) =1 SOn j β j ( j) (4.4) m x j Q ( ) (4.5) j=1 (x j 1, x j 2,, x j n) R j ( j) (4.6) x j 0 (, j) (4.7) The evaluaton functon (Equaton (4.1)) expresses that we must fnd optmal solutons (x j 1, x j 2,, x j n), (j = 1, 2,, m) that mnmze the expected value (expectaton) of the

7 Plannng Based on Unfulflled Order Rate 207 sum of producton cost and nventory cost. Equaton (4.2) s non negatve condton about nventory quantty that the demand s not over the forecast order. Equaton (4.3) s constrant about target producton quantty. Equaton (4.4) gves constrant about unfulflled order rate, Equaton (4.5) denotes lmt of producton capacty at each perod and Equaton (4.6) represents lnear constrant about producton. The above problem s a manufacture/nventory problem where demand changes stochastcally [2]. Although the (s, S) strategy [3,12] s known as an effectve tool for such a problem, t s not applcable when the problem has many condtons about producton. Thus we propose practcal and effectve algorthm by solvng lnear problem as partal problems repeatedly Procedure of algorthm We frst substtute E[S j ] = S j 0 + t=1 x j t t=1 dj t nto total cost MC of Equaton (4.1), MC = = ( m n ) ( m n ) p j x j + h j S j 0 j=1 =1 j=1 =1 { ( m n )} { ( m n )} + h j x j t h j d j t j=1 =1 t=1 j=1 =1 t=1 ( m n ) { ( m n )} p j x j + h j x j t + MC 1 (4.8) j=1 =1 j=1 =1 t=1 where { ( m n m n MC 1 = h j S j 0 h j j=1 =1 j=1 =1 t=1 d j t )} (4.9) Assumng that h j = h j, p j = p j, = 1,, n wthout losng generalty and practcalty on smplfyng a transformaton, frst term of Equaton (4.8) can be rewrtten by m j=1 p j r j because of m j=1 p j n =1 x j. Then by usng MC s calculated by m MC 2 = h j n (n t + 1) x j t, (4.10) j=1 t=1 m MC = MC 1 + MC 2 + p j r j (4.11) j=1 where MC 1 and m j=1 p j r j are constants. Here we formulate Mult-tem P problem applyng relaxaton strategy by exceptng Equaton (4.4) n Formula 1. [Mult-tem P problem] m mnmze h j n (n t + 1)x j t (4.12) j=1 t=1 s.t. S j 0 + x j t d j t 0 (, j) (4.13) t=1 t=1

8 208 N. Ueno, K. Okuhara, H. Ish, H. Shbuk & T. Kuramoto Fgure 4: Basc dea of proposed algorthm for a mult-tem producton plannng model n x j = r j ( j) (4.14) =1 m x j Q ( ) (4.15) j=1 (x j 1, x j 2,, x j n) R j ( j) (4.16) x j 0 (, j) (4.17) When we obtan (x j 1, x j 2,, x j n ), (j = 1, 2,, m) as optmal soluton of mult-tem P problem, t corresponds to global optmal soluton f t satsfes Equaton (4.4). If t does not satsfy Equaton (4.4) then we consder the followng mult-tem MP problem where the unfulflled order rate holds smaller value than certan target β j about jth part at perod. [Mult-tem MP problem] mnmze m h j n (n t + 1)x j t (4.18) j=1 t=1 s.t. S j 0 + x j t d j t 0 (, j) (4.19) t=1 t=1 S j 0 + x j t d j t K j ( B j, j) (4.20) t=1 t=1 n x j = r j ( j) (4.21) =1 m x j Q ( ) (4.22) j=1 (x j 1, x j 2,, x j n) R j ( j) (4.23) x j 0 (, j) (4.24) Compared wth the Mult-tem P problem, the Mult-tem MP problem nvolves Equaton (4.20) newly. From Fgure 4, K j and B j n Equaton (4.20) can be understood as follows. We defne the target unfulflled order rate β j for each perod n plan from unfulflled

9 order rate of jth part β j as follows; Plannng Based on Unfulflled Order Rate β j = n 1 β j (4.25) Let unfulflled order rate whch s calculated at the frst teraton on the mult-tem MP problem be β j(1). By usng β j(1), we derve β j(2) = β j(1) sgn(β j(1) β j ) βj(1) β j, (4.26) INC (2) j = arg max n (βj(2) ) (4.27) and B j = { (2) j }. In order to hold generally that the probablty, nventory quantty S j s less than y j at perod, s smaller than the unfulflled order rate β j n Fgure 4, we obtan y j β j = 1 2σ 2πσ j e 2 ds j, (4.28) K j = m j y j, (4.29) (Sj mj )2 B j denote the set of perod ndces, whch the unfulflled order rate of certan perod becomes smaller value than β j. Form Equaton (4.26), (4.27), (4.28) and (4.29), we fnd K j(2) satsfyng β j(2) = y j(2) 1 2πσ j e (S j mj(1) ) 2 2σ j 2 ds j, (4.30) K j(2) = m j(1) y j(2). (4.31) Thus we can reconstruct new mult-tem MP problem that the probablty, nventory quantty S j whch s less than 0 of jth part at perod (2) j, takes smaller value than β j(2) by applyng K j(2) and (2) j B j to Equaton (4.20). We wll show the reason why by solvng mult-tem MP problem repeatedly, we can fnd the approxmate soluton. Whenever Equaton (4.20) s added sequentally, we can get the mproved producton quantty x j(2) about jth part at perod that the probablty, whch nventory quantty S j(2) of the jth part at perod (2) j s less than 0, takes smaller value than β j(2) by K j(2) and (2) j β j to Equaton (4.20). Then mproved producton quantty x j(2) at one teraton makes both unfulflled order rate β j(2) of jth part at perod and the unfulflled order rate of jth part β j(2) monotonously decrease by Lemma 1 at secton 3.2. These mechansm s repeatedly by changng next step targeted unfulflled order rate β j(l) (l > 2) of jth part at perod decreasng, whch s defned lke Equaton (4.26) when l = 2. Now that, we ntroduce update parameter INC to next step unfulflled order rate n calculaton (Equaton (4.26)) of the target unfulflled order rate β j at each perod n plan. Fgure 5 shows nput factors between manufacturer and suppler. Purpose of research s to develop a suppler s producton plannng system whch s mnmzng total cost subject to unfulflled order rate and producton constrants when weekly forecast order and devaton of order σ are gven.

10 210 N. Ueno, K. Okuhara, H. Ish, H. Shbuk & T. Kuramoto Fgure 5: Input factors between manufacturer and suppler 4.3. Algorthm In ths secton, we descrbe algorthm for fndng the approxmate soluton of Formula 1. Step 1 Let l = 1, T j = {1, 2,, n} ( j), B j = {} ( j) and ϵ = T j denotes set of perod ndex. β j = 1 n 1 β j ( n, j m) (4.32) Step 2 Fndng a soluton of the mult-tem P problem, then we obtan the soluton (x j 1, x j 2,, x j n ), cost C (l) as the objectve value of Equaton (4.18), unfulflled order rate β j(l) and unfulflled order rate at each perod β j(l). Step 3 () If β j(l) β j, ( j) holds, then (x j 1, x j 2,, x j n ) s global approxmate soluton (END). () Otherwse, for certan j holdng β j(l) > β j, t s updated by then β j(l+1) = β j(l) sgn(β j(l) β j ) βj(l) β j INC, (4.33) (l+1) j = arg max n (βj(l+1) ) (4.34) {B j } = {B j } { (l+1) j }, {T j } = {T j } { (l+1) j } (4.35) and reconstruct the mult-tem MP problem by usng β j(l+1) = y j(l+1) 1 2πσ j e (S j mj(l) ) 2 2σ j 2 ds j, (4.36) K j(l+1) = m j(l) y j(l+1). (4.37) Step 4 If we get a soluton (x j 1, x j 2,, x j n ), cost C (l+1), and the unfulflled order rate β j(l+1) for the mult-tem MP problem, and the followng relaton holds then procedure s end. β j(l+1) β j(l) < ε j ( j) (4.38)

11 Plannng Based on Unfulflled Order Rate 211 Table 1 d, ω, σ used n numercal example d ω σ Fgure 6: Result of proposed algorthm If we can not get a soluton (x j 1, x j 2,, x j n ), cost C (l+1), and the unfulflled order rate β j(l+1) for the mult-tem MP problem, then procedure s end. Step 5 () If β j(l+1) β j, ( j) holds, then (x j 1, x j 2,, x j n ) s global approxmate soluton (END). () Otherwse, let l = l + 1, then go to Step 3 (). 5. Numercal Experments 5.1. Effect of proposed unfulflled order rate as crtera In ths secton, n order to confrm the effectveness of crtera of unfulflled order rate we propose n ths paper, we apply our algorthm to actual operaton data n auto parts suppler. Table 1 shows d, ω and σ about forecast order. It s assumed that S 0 = 24, p = 1, h = 1, r = 82, ω = d 0.35 for each and β = We set some parameters based on the actual data such as manufacturng and nventory holdng costs to show the effcency of ths approach. Fgure 6 shows comparson of nventory change between the result of proposal method n case of m = 1, INC = 2 and the actual producton plan. We call the result of real operaton planned by producton manager n the suppler base case. Here base case s consequence of actual operaton data n suppler. From the results of Fgure 6, we can fnd that our method provde lower nventory change because of consderng nventory cost n frst half, and provde hgher nventory change for respondng to fluctuaton of demand n last half. Fgure 7 shows the transton of both unfulflled order rate and the cost at each teraton by usng proposed algorthm. The unfulflled order rate SO n decreases monotonously and soluton s obtaned after 21 teratons. The cost s plotted by relatve comparson for settng the cost of ntal state to 100%. We can fnd that the unfulflled order rate s mproved at each teraton. We can also fnd that the exstence of trade-off between the unfulflled order rate and the cost s shown n Fgure 7. Fgure 8 shows mprovement of unfulflled order rate by proposed technque when ntal nventory s changed. When ntal nventory s gven by 24, then unfulflled order rate takes

12 212 N. Ueno, K. Okuhara, H. Ish, H. Shbuk & T. Kuramoto Fgure 7: Result of unfulflled order rate and cost for each teraton ( : base case, : proposed algorthm, : proposed algorthm wth same cost of base case) Fgure 8: Performance of proposal algorthm 11.7% by base case, but t takes 5.8% by proposal algorthm. Partcularly, from the result about proposed method wth same cost of base case, t takes 8.5% then we can fnd our method mprovng both unfulflled order rate and cost Behavor of unfulflled order rate n two-tem case Here we consder the case that suppler s product conssts of two tems. Table 2 shows d j used n an example system of 2 tems, where S j 0 = 24, p j = 1, h j = 1, r 1 = 82, r 2 = 135, Q = 50, ω j = d j 0.35 and β j = Fgure 9 shows results of nventory by proposed algorthm at 1st and 7th teratons, and fnal soluton. Compared wth obtaned result at 1st teraton, we can fnd the tendency of becomng hgher about nventory level when teraton ncreases. Fgure 9 shows typcal result. Fgure 10 shows result of producton quantty by proposed algorthm. Compared wth base case, we can fnd that proposed algorthm plans necessarly the precedence producton at 4th, 5th, 6th and 7th for tem number 1, also t plans necessarly the precedence producton at 2nd for tem number 2. Table 2 d j used n an example system of 2 tems d d

13 Plannng Based on Unfulflled Order Rate 213 < tem number 1 > < tem number 2 > Fgure 9: Results of nventory by proposed algorthm for each tem Fgure 10: Result of producton quantty by proposed algorthm

14 214 N. Ueno, K. Okuhara, H. Ish, H. Shbuk & T. Kuramoto Fgure 11: Result of cost and unfulflled order rate for each teraton Fgure 12: Results for each teraton by changng a step sze for update rule Fgure 11 shows result of cost and unfulflled order rate for each teraton n the case that INC =2. The unfulflled order rate decreases gradually, and soluton s obtaned after 48th teraton. Cost s represented by rate compared wth the cost obtaned at 1st teraton. From the result, we can fnd that proposed algorthm can mprove the unfulflled order rate, and we can understand exstence of alternatve between cost and the unfulflled order rate. Fgure 12 shows results for each teraton by changng a step sze for update rule INC where Q = 50. It can be found that algorthm, where INC=2, s converged faster than another. Table 3 shows results by changng Q and INC. Result when INC = 2 s qute nearly converged same as one when INC = 10 for each Q. Fgure 13 shows performance of proposed algorthm for each tem. For the tem number 1, the unfulflled order rate of base case takes 11.7%, but the unfulflled order rate obtaned by proposed method becomes 5.8%. Especally, when we compare t under the same cost of base case, proposed method realzes 8.3% about the unfulflled order rate. Smlarly, for tem number 2, the unfulflled order rate s mproved from 52.1% n base case to 14.5% n Table 3 Results by changng Q and INC SOn 1 SOn 2 cost Q INC = 2 INC = 10 INC = 2 INC = 10 INC = 2 INC =

15 Plannng Based on Unfulflled Order Rate 215 Fgure 13: Performance of proposed algorthm for each tem Table 4 d j used n an example system of 5 tems d d d d d proposed method. Furthermore, even f we compare t under the same cost of base case, proposed method realzes 36.1% about the unfulflled order rate. Proposed algorthm s better than base case from the pont of both unfulflled order rate and cost. From these results, compared wth base case, we can fnd that the proposed method can mprove both cost and unfulflled order rate under keepng manufacturng capacty restrctons and so on Behavor of unfulflled order rate n mult-tem case In ths secton, the behavor of unfulflled order rate n mult-tem case s nvestgated. Table 4 shows d j used n ths numercal example wth 5 tems. The other parameters are assumed to be S 1 0 = 18, S 2 0 = 24, S 3 0 = 14, S 4 0 = 24, S 5 0 = 18, p j = 1, h j = 1, r 1 = 82, r 2 = 135, r 3 = 112, r 4 = 132, r 5 = 140, Q = 100, ω j = d j 0.20 and β j = Fgure 14 shows results when system apples INC = 2 and INC = 10 as the step sze for update rule. From ts results, we can understand that the unfulflled order rate changes decrease monotonously and converges also n case of 5 tems Proposal desgn method n actual case In ths secton, we show proposal desgn method wth 2 steps n actual case. As the desgn method s llustrated n Fgure 15, we frst solve the producton plannng problem by the proposed algorthm n the prevous secton, and we call ths algorthm MCPS (Masscustomzaton Producton Plannng and Management System), and by ths result we modfy orgnal due date for delvery. Almost all the case t means precedence producton. Next we execute schedulng based on modfed due date for delvery. So, we frst determne the producton plannng of fnal stage for each tem n order to obtan mnmum nventory level whle keepng target unfulflled order rate. For the due date of delvery that s gven by each order, t can be realzed by mprovng a schedule that fnshes the fnal stage. For example, fnal stage s set by the fnal process about suppler s assembly, or the delvery to manufacturer s warehouse from suppler. That s, such schedule

16 216 N. Ueno, K. Okuhara, H. Ish, H. Shbuk & T. Kuramoto (a) INC=2 (b) INC=10 Fgure 14: Result of unfulflled order rate and cost for each teraton ( : date that order s assgned and to be manufactured) Fgure 15: Proposal desgn method n actual case

17 Plannng Based on Unfulflled Order Rate 217 ams both the proper necessary precedence producton for each tem and the holdng of mnmum nventory. Next we schedule upstream process from fnal state by usng commercal software that supports system based on Make-to-Order producton system. We wll select the best software that can take nto account about the condton of suppler s producton restrcton n detal. In order to realze mass customzaton, n the case that we apply the exstng technque, for example, MRP II [11], APS [7] drectly, t s necessary to set up safety stock n advance. However, safety stock method s not realstc from the fact that forecast order changes n the planned perod usually happen n the automoble ndustry. By applyng the proposed algorthm, t s not necessary to gve the safety stock n advance. When the planned unfulflled order rate s gven n advance, then the nearly optmal producton quantty and smultaneously the mnmum nventory for planned term, nstead of safety stock, are decded. In sum, both MCPS and schedulng software solve sequentally, and realze mass customzaton. 6. Concludng Remarks In ths paper, we took up producton plannng of parts assembly suppler n auto ndustry, and defned suppler s mult-tem producton plannng and management system for mplementng mass customzaton. The problem about proper nventory quantty of mult-tem producton was formulated by stochastc programmng problem. We proposed the effectve algorthm to obtan nearly optmal nventory quantty for the producton plannng and management system mplementng mass customzaton. Further work s an applcaton of the decomposton prncple [8] to fnd the soluton of the mult-tem MP problem whch was derved, because f the number of tems becomes large then the sze of sub problem for teraton become huge. Acknowledgments The authors would lke to thank referees for ther valuable comments on the orgnal verson whch mproved ths paper very much. In addton, the authors wsh to express apprecaton to Mr. Kyozo Furuta for hs assstance. Ths research was supported by Grant-n-Ad for Scentfc Research (c) References [1] S. Bller, E.K. Bsh, and A. Murel: Impact of manufacturng flexblty on supply chan performance n the automotve ndustry. In J.S. Song and D.D. Yao (eds.): Supply Chan Structures Coordnaton, Informaton and Optmzaton (Kluwer s Internatonal Seres, 2001). [2] J.R. Brge and F. Louveaux: Introducton to Stochastc Programmng (Sprnger-Verlag, 1997). [3] J. Bramel and D. Smch-Lev: The Logc of Logstcs (Sprnger-Verlag, 1997). [4] S. Chopra and P. Mendl: Supply Chan Management Strategy, Plannng and Operaton (Prentce-Hall, 2000). [5] J.H. Glmore and B.J. Pne: Markets of One : Creatng Customer-Unque Value through Mass Customzaton (Harvard Busness School Press, 2000). [6] Y. Kobayash and H. Tsubone: Producton seat bookng system for the combnaton of make-to-order and make-to-stock as producton envronment. Journal of Japan Industral Management Assocaton, 52 (2001),

18 218 N. Ueno, K. Okuhara, H. Ish, H. Shbuk & T. Kuramoto [7] M. Kuroda: MRP to APS - nnovaton of producton management and the new role of schedulng. Proceedngs of the Schedulng Symposum 2002 (2002), 2 13 (n Japanese). [8] L.S. Lasdon: Optmzaton Theory for Large Systems (Dover Publcatons, 2002). [9] J. Nakane (ed.): BTO Producton System (Nkkan Kogyo Shmbun, 2000). [10] B.J. Pne: Mass Customzaton (Harvard Busness School Press, 1993). [11] K. Shekh: Manufacturng Resource Plannng (MRP II) (McGraw-Hll, 2003). [12] E.A. Slver, D.F. Pyke, and R. Peterson: Inventory Management and Producton Plannng and Schedulng, 3rd edton (Jon Wlley & Sons, 1998). [13] T. Tamura and S. Fujta: Producton seat system for producton management. Journal of Japan Industral Management Assocaton, 4 (1994), [14] S. Tayur, R. Ganeshan, and M. Magazne: Quanttatve Models for Supply Chan Management (Kluwer s Internatonal Seres, 1998). Nobuyuk Ueno Faculty of Management and Informaton Systems Prefectural Unversty of Hroshma Ujnahgash Mnam-ku Hroshma , Japan E-mal: ueno@pu-hroshma.ac.jp

Module 3 LOSSY IMAGE COMPRESSION SYSTEMS. Version 2 ECE IIT, Kharagpur

Module 3 LOSSY IMAGE COMPRESSION SYSTEMS. Version 2 ECE IIT, Kharagpur Module 3 LOSSY IMAGE COMPRESSION SYSTEMS Verson ECE IIT, Kharagpur Lesson 6 Theory of Quantzaton Verson ECE IIT, Kharagpur Instructonal Objectves At the end of ths lesson, the students should be able to:

More information

COMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS

COMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS Avalable onlne at http://sck.org J. Math. Comput. Sc. 3 (3), No., 6-3 ISSN: 97-537 COMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS

More information

On the Multicriteria Integer Network Flow Problem

On the Multicriteria Integer Network Flow Problem BULGARIAN ACADEMY OF SCIENCES CYBERNETICS AND INFORMATION TECHNOLOGIES Volume 5, No 2 Sofa 2005 On the Multcrtera Integer Network Flow Problem Vassl Vasslev, Marana Nkolova, Maryana Vassleva Insttute of

More information

The Minimum Universal Cost Flow in an Infeasible Flow Network

The Minimum Universal Cost Flow in an Infeasible Flow Network Journal of Scences, Islamc Republc of Iran 17(2): 175-180 (2006) Unversty of Tehran, ISSN 1016-1104 http://jscencesutacr The Mnmum Unversal Cost Flow n an Infeasble Flow Network H Saleh Fathabad * M Bagheran

More information

Winter 2008 CS567 Stochastic Linear/Integer Programming Guest Lecturer: Xu, Huan

Winter 2008 CS567 Stochastic Linear/Integer Programming Guest Lecturer: Xu, Huan Wnter 2008 CS567 Stochastc Lnear/Integer Programmng Guest Lecturer: Xu, Huan Class 2: More Modelng Examples 1 Capacty Expanson Capacty expanson models optmal choces of the tmng and levels of nvestments

More information

Assortment Optimization under MNL

Assortment Optimization under MNL Assortment Optmzaton under MNL Haotan Song Aprl 30, 2017 1 Introducton The assortment optmzaton problem ams to fnd the revenue-maxmzng assortment of products to offer when the prces of products are fxed.

More information

Lecture Notes on Linear Regression

Lecture Notes on Linear Regression Lecture Notes on Lnear Regresson Feng L fl@sdueducn Shandong Unversty, Chna Lnear Regresson Problem In regresson problem, we am at predct a contnuous target value gven an nput feature vector We assume

More information

CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE

CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE Analytcal soluton s usually not possble when exctaton vares arbtrarly wth tme or f the system s nonlnear. Such problems can be solved by numercal tmesteppng

More information

Online Appendix. t=1 (p t w)q t. Then the first order condition shows that

Online Appendix. t=1 (p t w)q t. Then the first order condition shows that Artcle forthcomng to ; manuscrpt no (Please, provde the manuscrpt number!) 1 Onlne Appendx Appendx E: Proofs Proof of Proposton 1 Frst we derve the equlbrum when the manufacturer does not vertcally ntegrate

More information

ECE559VV Project Report

ECE559VV Project Report ECE559VV Project Report (Supplementary Notes Loc Xuan Bu I. MAX SUM-RATE SCHEDULING: THE UPLINK CASE We have seen (n the presentaton that, for downlnk (broadcast channels, the strategy maxmzng the sum-rate

More information

Interactive Bi-Level Multi-Objective Integer. Non-linear Programming Problem

Interactive Bi-Level Multi-Objective Integer. Non-linear Programming Problem Appled Mathematcal Scences Vol 5 0 no 65 3 33 Interactve B-Level Mult-Objectve Integer Non-lnear Programmng Problem O E Emam Department of Informaton Systems aculty of Computer Scence and nformaton Helwan

More information

EEE 241: Linear Systems

EEE 241: Linear Systems EEE : Lnear Systems Summary #: Backpropagaton BACKPROPAGATION The perceptron rule as well as the Wdrow Hoff learnng were desgned to tran sngle layer networks. They suffer from the same dsadvantage: they

More information

Simultaneous Optimization of Berth Allocation, Quay Crane Assignment and Quay Crane Scheduling Problems in Container Terminals

Simultaneous Optimization of Berth Allocation, Quay Crane Assignment and Quay Crane Scheduling Problems in Container Terminals Smultaneous Optmzaton of Berth Allocaton, Quay Crane Assgnment and Quay Crane Schedulng Problems n Contaner Termnals Necat Aras, Yavuz Türkoğulları, Z. Caner Taşkın, Kuban Altınel Abstract In ths work,

More information

Chapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems

Chapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems Numercal Analyss by Dr. Anta Pal Assstant Professor Department of Mathematcs Natonal Insttute of Technology Durgapur Durgapur-713209 emal: anta.bue@gmal.com 1 . Chapter 5 Soluton of System of Lnear Equatons

More information

Transfer Functions. Convenient representation of a linear, dynamic model. A transfer function (TF) relates one input and one output: ( ) system

Transfer Functions. Convenient representation of a linear, dynamic model. A transfer function (TF) relates one input and one output: ( ) system Transfer Functons Convenent representaton of a lnear, dynamc model. A transfer functon (TF) relates one nput and one output: x t X s y t system Y s The followng termnology s used: x y nput output forcng

More information

Feature Selection: Part 1

Feature Selection: Part 1 CSE 546: Machne Learnng Lecture 5 Feature Selecton: Part 1 Instructor: Sham Kakade 1 Regresson n the hgh dmensonal settng How do we learn when the number of features d s greater than the sample sze n?

More information

MMA and GCMMA two methods for nonlinear optimization

MMA and GCMMA two methods for nonlinear optimization MMA and GCMMA two methods for nonlnear optmzaton Krster Svanberg Optmzaton and Systems Theory, KTH, Stockholm, Sweden. krlle@math.kth.se Ths note descrbes the algorthms used n the author s 2007 mplementatons

More information

Generalized Linear Methods

Generalized Linear Methods Generalzed Lnear Methods 1 Introducton In the Ensemble Methods the general dea s that usng a combnaton of several weak learner one could make a better learner. More formally, assume that we have a set

More information

Structure and Drive Paul A. Jensen Copyright July 20, 2003

Structure and Drive Paul A. Jensen Copyright July 20, 2003 Structure and Drve Paul A. Jensen Copyrght July 20, 2003 A system s made up of several operatons wth flow passng between them. The structure of the system descrbes the flow paths from nputs to outputs.

More information

The Value of Demand Postponement under Demand Uncertainty

The Value of Demand Postponement under Demand Uncertainty Recent Researches n Appled Mathematcs, Smulaton and Modellng The Value of emand Postponement under emand Uncertanty Rawee Suwandechocha Abstract Resource or capacty nvestment has a hgh mpact on the frm

More information

ON A DETERMINATION OF THE INITIAL FUNCTIONS FROM THE OBSERVED VALUES OF THE BOUNDARY FUNCTIONS FOR THE SECOND-ORDER HYPERBOLIC EQUATION

ON A DETERMINATION OF THE INITIAL FUNCTIONS FROM THE OBSERVED VALUES OF THE BOUNDARY FUNCTIONS FOR THE SECOND-ORDER HYPERBOLIC EQUATION Advanced Mathematcal Models & Applcatons Vol.3, No.3, 2018, pp.215-222 ON A DETERMINATION OF THE INITIAL FUNCTIONS FROM THE OBSERVED VALUES OF THE BOUNDARY FUNCTIONS FOR THE SECOND-ORDER HYPERBOLIC EUATION

More information

A PROBABILITY-DRIVEN SEARCH ALGORITHM FOR SOLVING MULTI-OBJECTIVE OPTIMIZATION PROBLEMS

A PROBABILITY-DRIVEN SEARCH ALGORITHM FOR SOLVING MULTI-OBJECTIVE OPTIMIZATION PROBLEMS HCMC Unversty of Pedagogy Thong Nguyen Huu et al. A PROBABILITY-DRIVEN SEARCH ALGORITHM FOR SOLVING MULTI-OBJECTIVE OPTIMIZATION PROBLEMS Thong Nguyen Huu and Hao Tran Van Department of mathematcs-nformaton,

More information

Some modelling aspects for the Matlab implementation of MMA

Some modelling aspects for the Matlab implementation of MMA Some modellng aspects for the Matlab mplementaton of MMA Krster Svanberg krlle@math.kth.se Optmzaton and Systems Theory Department of Mathematcs KTH, SE 10044 Stockholm September 2004 1. Consdered optmzaton

More information

A new Approach for Solving Linear Ordinary Differential Equations

A new Approach for Solving Linear Ordinary Differential Equations , ISSN 974-57X (Onlne), ISSN 974-5718 (Prnt), Vol. ; Issue No. 1; Year 14, Copyrght 13-14 by CESER PUBLICATIONS A new Approach for Solvng Lnear Ordnary Dfferental Equatons Fawz Abdelwahd Department of

More information

Chapter - 2. Distribution System Power Flow Analysis

Chapter - 2. Distribution System Power Flow Analysis Chapter - 2 Dstrbuton System Power Flow Analyss CHAPTER - 2 Radal Dstrbuton System Load Flow 2.1 Introducton Load flow s an mportant tool [66] for analyzng electrcal power system network performance. Load

More information

A Hybrid Variational Iteration Method for Blasius Equation

A Hybrid Variational Iteration Method for Blasius Equation Avalable at http://pvamu.edu/aam Appl. Appl. Math. ISSN: 1932-9466 Vol. 10, Issue 1 (June 2015), pp. 223-229 Applcatons and Appled Mathematcs: An Internatonal Journal (AAM) A Hybrd Varatonal Iteraton Method

More information

LOW BIAS INTEGRATED PATH ESTIMATORS. James M. Calvin

LOW BIAS INTEGRATED PATH ESTIMATORS. James M. Calvin Proceedngs of the 007 Wnter Smulaton Conference S G Henderson, B Bller, M-H Hseh, J Shortle, J D Tew, and R R Barton, eds LOW BIAS INTEGRATED PATH ESTIMATORS James M Calvn Department of Computer Scence

More information

Appendix for Causal Interaction in Factorial Experiments: Application to Conjoint Analysis

Appendix for Causal Interaction in Factorial Experiments: Application to Conjoint Analysis A Appendx for Causal Interacton n Factoral Experments: Applcaton to Conjont Analyss Mathematcal Appendx: Proofs of Theorems A. Lemmas Below, we descrbe all the lemmas, whch are used to prove the man theorems

More information

Stanford University CS359G: Graph Partitioning and Expanders Handout 4 Luca Trevisan January 13, 2011

Stanford University CS359G: Graph Partitioning and Expanders Handout 4 Luca Trevisan January 13, 2011 Stanford Unversty CS359G: Graph Parttonng and Expanders Handout 4 Luca Trevsan January 3, 0 Lecture 4 In whch we prove the dffcult drecton of Cheeger s nequalty. As n the past lectures, consder an undrected

More information

Econ107 Applied Econometrics Topic 3: Classical Model (Studenmund, Chapter 4)

Econ107 Applied Econometrics Topic 3: Classical Model (Studenmund, Chapter 4) I. Classcal Assumptons Econ7 Appled Econometrcs Topc 3: Classcal Model (Studenmund, Chapter 4) We have defned OLS and studed some algebrac propertes of OLS. In ths topc we wll study statstcal propertes

More information

ANSWERS. Problem 1. and the moment generating function (mgf) by. defined for any real t. Use this to show that E( U) var( U)

ANSWERS. Problem 1. and the moment generating function (mgf) by. defined for any real t. Use this to show that E( U) var( U) Econ 413 Exam 13 H ANSWERS Settet er nndelt 9 deloppgaver, A,B,C, som alle anbefales å telle lkt for å gøre det ltt lettere å stå. Svar er gtt . Unfortunately, there s a prntng error n the hnt of

More information

CSci 6974 and ECSE 6966 Math. Tech. for Vision, Graphics and Robotics Lecture 21, April 17, 2006 Estimating A Plane Homography

CSci 6974 and ECSE 6966 Math. Tech. for Vision, Graphics and Robotics Lecture 21, April 17, 2006 Estimating A Plane Homography CSc 6974 and ECSE 6966 Math. Tech. for Vson, Graphcs and Robotcs Lecture 21, Aprl 17, 2006 Estmatng A Plane Homography Overvew We contnue wth a dscusson of the major ssues, usng estmaton of plane projectve

More information

A new construction of 3-separable matrices via an improved decoding of Macula s construction

A new construction of 3-separable matrices via an improved decoding of Macula s construction Dscrete Optmzaton 5 008 700 704 Contents lsts avalable at ScenceDrect Dscrete Optmzaton journal homepage: wwwelsevercom/locate/dsopt A new constructon of 3-separable matrces va an mproved decodng of Macula

More information

A New Algorithm for Finding a Fuzzy Optimal. Solution for Fuzzy Transportation Problems

A New Algorithm for Finding a Fuzzy Optimal. Solution for Fuzzy Transportation Problems Appled Mathematcal Scences, Vol. 4, 200, no. 2, 79-90 A New Algorthm for Fndng a Fuzzy Optmal Soluton for Fuzzy Transportaton Problems P. Pandan and G. Nataraan Department of Mathematcs, School of Scence

More information

Amiri s Supply Chain Model. System Engineering b Department of Mathematics and Statistics c Odette School of Business

Amiri s Supply Chain Model. System Engineering b Department of Mathematics and Statistics c Odette School of Business Amr s Supply Chan Model by S. Ashtab a,, R.J. Caron b E. Selvarajah c a Department of Industral Manufacturng System Engneerng b Department of Mathematcs Statstcs c Odette School of Busness Unversty of

More information

College of Computer & Information Science Fall 2009 Northeastern University 20 October 2009

College of Computer & Information Science Fall 2009 Northeastern University 20 October 2009 College of Computer & Informaton Scence Fall 2009 Northeastern Unversty 20 October 2009 CS7880: Algorthmc Power Tools Scrbe: Jan Wen and Laura Poplawsk Lecture Outlne: Prmal-dual schema Network Desgn:

More information

Linear Approximation with Regularization and Moving Least Squares

Linear Approximation with Regularization and Moving Least Squares Lnear Approxmaton wth Regularzaton and Movng Least Squares Igor Grešovn May 007 Revson 4.6 (Revson : March 004). 5 4 3 0.5 3 3.5 4 Contents: Lnear Fttng...4. Weghted Least Squares n Functon Approxmaton...

More information

The Order Relation and Trace Inequalities for. Hermitian Operators

The Order Relation and Trace Inequalities for. Hermitian Operators Internatonal Mathematcal Forum, Vol 3, 08, no, 507-57 HIKARI Ltd, wwwm-hkarcom https://doorg/0988/mf088055 The Order Relaton and Trace Inequaltes for Hermtan Operators Y Huang School of Informaton Scence

More information

Numerical Heat and Mass Transfer

Numerical Heat and Mass Transfer Master degree n Mechancal Engneerng Numercal Heat and Mass Transfer 06-Fnte-Dfference Method (One-dmensonal, steady state heat conducton) Fausto Arpno f.arpno@uncas.t Introducton Why we use models and

More information

Difference Equations

Difference Equations Dfference Equatons c Jan Vrbk 1 Bascs Suppose a sequence of numbers, say a 0,a 1,a,a 3,... s defned by a certan general relatonshp between, say, three consecutve values of the sequence, e.g. a + +3a +1

More information

Computing Correlated Equilibria in Multi-Player Games

Computing Correlated Equilibria in Multi-Player Games Computng Correlated Equlbra n Mult-Player Games Chrstos H. Papadmtrou Presented by Zhanxang Huang December 7th, 2005 1 The Author Dr. Chrstos H. Papadmtrou CS professor at UC Berkley (taught at Harvard,

More information

Additional Codes using Finite Difference Method. 1 HJB Equation for Consumption-Saving Problem Without Uncertainty

Additional Codes using Finite Difference Method. 1 HJB Equation for Consumption-Saving Problem Without Uncertainty Addtonal Codes usng Fnte Dfference Method Benamn Moll 1 HJB Equaton for Consumpton-Savng Problem Wthout Uncertanty Before consderng the case wth stochastc ncome n http://www.prnceton.edu/~moll/ HACTproect/HACT_Numercal_Appendx.pdf,

More information

Uncertainty in measurements of power and energy on power networks

Uncertainty in measurements of power and energy on power networks Uncertanty n measurements of power and energy on power networks E. Manov, N. Kolev Department of Measurement and Instrumentaton, Techncal Unversty Sofa, bul. Klment Ohrdsk No8, bl., 000 Sofa, Bulgara Tel./fax:

More information

Introduction to Vapor/Liquid Equilibrium, part 2. Raoult s Law:

Introduction to Vapor/Liquid Equilibrium, part 2. Raoult s Law: CE304, Sprng 2004 Lecture 4 Introducton to Vapor/Lqud Equlbrum, part 2 Raoult s Law: The smplest model that allows us do VLE calculatons s obtaned when we assume that the vapor phase s an deal gas, and

More information

Chapter 13: Multiple Regression

Chapter 13: Multiple Regression Chapter 13: Multple Regresson 13.1 Developng the multple-regresson Model The general model can be descrbed as: It smplfes for two ndependent varables: The sample ft parameter b 0, b 1, and b are used to

More information

On a direct solver for linear least squares problems

On a direct solver for linear least squares problems ISSN 2066-6594 Ann. Acad. Rom. Sc. Ser. Math. Appl. Vol. 8, No. 2/2016 On a drect solver for lnear least squares problems Constantn Popa Abstract The Null Space (NS) algorthm s a drect solver for lnear

More information

Games of Threats. Elon Kohlberg Abraham Neyman. Working Paper

Games of Threats. Elon Kohlberg Abraham Neyman. Working Paper Games of Threats Elon Kohlberg Abraham Neyman Workng Paper 18-023 Games of Threats Elon Kohlberg Harvard Busness School Abraham Neyman The Hebrew Unversty of Jerusalem Workng Paper 18-023 Copyrght 2017

More information

On the correction of the h-index for career length

On the correction of the h-index for career length 1 On the correcton of the h-ndex for career length by L. Egghe Unverstet Hasselt (UHasselt), Campus Depenbeek, Agoralaan, B-3590 Depenbeek, Belgum 1 and Unverstet Antwerpen (UA), IBW, Stadscampus, Venusstraat

More information

Single-Facility Scheduling over Long Time Horizons by Logic-based Benders Decomposition

Single-Facility Scheduling over Long Time Horizons by Logic-based Benders Decomposition Sngle-Faclty Schedulng over Long Tme Horzons by Logc-based Benders Decomposton Elvn Coban and J. N. Hooker Tepper School of Busness, Carnege Mellon Unversty ecoban@andrew.cmu.edu, john@hooker.tepper.cmu.edu

More information

Erratum: A Generalized Path Integral Control Approach to Reinforcement Learning

Erratum: A Generalized Path Integral Control Approach to Reinforcement Learning Journal of Machne Learnng Research 00-9 Submtted /0; Publshed 7/ Erratum: A Generalzed Path Integral Control Approach to Renforcement Learnng Evangelos ATheodorou Jonas Buchl Stefan Schaal Department of

More information

Problem Set 9 Solutions

Problem Set 9 Solutions Desgn and Analyss of Algorthms May 4, 2015 Massachusetts Insttute of Technology 6.046J/18.410J Profs. Erk Demane, Srn Devadas, and Nancy Lynch Problem Set 9 Solutons Problem Set 9 Solutons Ths problem

More information

Module 9. Lecture 6. Duality in Assignment Problems

Module 9. Lecture 6. Duality in Assignment Problems Module 9 1 Lecture 6 Dualty n Assgnment Problems In ths lecture we attempt to answer few other mportant questons posed n earler lecture for (AP) and see how some of them can be explaned through the concept

More information

Perfect Competition and the Nash Bargaining Solution

Perfect Competition and the Nash Bargaining Solution Perfect Competton and the Nash Barganng Soluton Renhard John Department of Economcs Unversty of Bonn Adenauerallee 24-42 53113 Bonn, Germany emal: rohn@un-bonn.de May 2005 Abstract For a lnear exchange

More information

Annexes. EC.1. Cycle-base move illustration. EC.2. Problem Instances

Annexes. EC.1. Cycle-base move illustration. EC.2. Problem Instances ec Annexes Ths Annex frst llustrates a cycle-based move n the dynamc-block generaton tabu search. It then dsplays the characterstcs of the nstance sets, followed by detaled results of the parametercalbraton

More information

The Exact Formulation of the Inverse of the Tridiagonal Matrix for Solving the 1D Poisson Equation with the Finite Difference Method

The Exact Formulation of the Inverse of the Tridiagonal Matrix for Solving the 1D Poisson Equation with the Finite Difference Method Journal of Electromagnetc Analyss and Applcatons, 04, 6, 0-08 Publshed Onlne September 04 n ScRes. http://www.scrp.org/journal/jemaa http://dx.do.org/0.46/jemaa.04.6000 The Exact Formulaton of the Inverse

More information

An Admission Control Algorithm in Cloud Computing Systems

An Admission Control Algorithm in Cloud Computing Systems An Admsson Control Algorthm n Cloud Computng Systems Authors: Frank Yeong-Sung Ln Department of Informaton Management Natonal Tawan Unversty Tape, Tawan, R.O.C. ysln@m.ntu.edu.tw Yngje Lan Management Scence

More information

Lecture 10 Support Vector Machines II

Lecture 10 Support Vector Machines II Lecture 10 Support Vector Machnes II 22 February 2016 Taylor B. Arnold Yale Statstcs STAT 365/665 1/28 Notes: Problem 3 s posted and due ths upcomng Frday There was an early bug n the fake-test data; fxed

More information

Simulated Power of the Discrete Cramér-von Mises Goodness-of-Fit Tests

Simulated Power of the Discrete Cramér-von Mises Goodness-of-Fit Tests Smulated of the Cramér-von Mses Goodness-of-Ft Tests Steele, M., Chaselng, J. and 3 Hurst, C. School of Mathematcal and Physcal Scences, James Cook Unversty, Australan School of Envronmental Studes, Grffth

More information

For now, let us focus on a specific model of neurons. These are simplified from reality but can achieve remarkable results.

For now, let us focus on a specific model of neurons. These are simplified from reality but can achieve remarkable results. Neural Networks : Dervaton compled by Alvn Wan from Professor Jtendra Malk s lecture Ths type of computaton s called deep learnng and s the most popular method for many problems, such as computer vson

More information

1 Convex Optimization

1 Convex Optimization Convex Optmzaton We wll consder convex optmzaton problems. Namely, mnmzaton problems where the objectve s convex (we assume no constrants for now). Such problems often arse n machne learnng. For example,

More information

Inventory Model with Backorder Price Discount

Inventory Model with Backorder Price Discount Vol. No. 7-7 ead Tme and Orderng Cost Reductons are Interdependent n Inventory Model wth Bacorder Prce scount Yu-Jen n Receved: Mar. 7 Frst Revson: May. 7 7 ccepted: May. 7 bstract The stochastc nventory

More information

3.1 Expectation of Functions of Several Random Variables. )' be a k-dimensional discrete or continuous random vector, with joint PMF p (, E X E X1 E X

3.1 Expectation of Functions of Several Random Variables. )' be a k-dimensional discrete or continuous random vector, with joint PMF p (, E X E X1 E X Statstcs 1: Probablty Theory II 37 3 EPECTATION OF SEVERAL RANDOM VARIABLES As n Probablty Theory I, the nterest n most stuatons les not on the actual dstrbuton of a random vector, but rather on a number

More information

Credit Card Pricing and Impact of Adverse Selection

Credit Card Pricing and Impact of Adverse Selection Credt Card Prcng and Impact of Adverse Selecton Bo Huang and Lyn C. Thomas Unversty of Southampton Contents Background Aucton model of credt card solctaton - Errors n probablty of beng Good - Errors n

More information

LINEAR REGRESSION ANALYSIS. MODULE IX Lecture Multicollinearity

LINEAR REGRESSION ANALYSIS. MODULE IX Lecture Multicollinearity LINEAR REGRESSION ANALYSIS MODULE IX Lecture - 30 Multcollnearty Dr. Shalabh Department of Mathematcs and Statstcs Indan Insttute of Technology Kanpur 2 Remedes for multcollnearty Varous technques have

More information

Time-Varying Systems and Computations Lecture 6

Time-Varying Systems and Computations Lecture 6 Tme-Varyng Systems and Computatons Lecture 6 Klaus Depold 14. Januar 2014 The Kalman Flter The Kalman estmaton flter attempts to estmate the actual state of an unknown dscrete dynamcal system, gven nosy

More information

princeton univ. F 17 cos 521: Advanced Algorithm Design Lecture 7: LP Duality Lecturer: Matt Weinberg

princeton univ. F 17 cos 521: Advanced Algorithm Design Lecture 7: LP Duality Lecturer: Matt Weinberg prnceton unv. F 17 cos 521: Advanced Algorthm Desgn Lecture 7: LP Dualty Lecturer: Matt Wenberg Scrbe: LP Dualty s an extremely useful tool for analyzng structural propertes of lnear programs. Whle there

More information

Kernel Methods and SVMs Extension

Kernel Methods and SVMs Extension Kernel Methods and SVMs Extenson The purpose of ths document s to revew materal covered n Machne Learnng 1 Supervsed Learnng regardng support vector machnes (SVMs). Ths document also provdes a general

More information

arxiv:cs.cv/ Jun 2000

arxiv:cs.cv/ Jun 2000 Correlaton over Decomposed Sgnals: A Non-Lnear Approach to Fast and Effectve Sequences Comparson Lucano da Fontoura Costa arxv:cs.cv/0006040 28 Jun 2000 Cybernetc Vson Research Group IFSC Unversty of São

More information

Chapter Newton s Method

Chapter Newton s Method Chapter 9. Newton s Method After readng ths chapter, you should be able to:. Understand how Newton s method s dfferent from the Golden Secton Search method. Understand how Newton s method works 3. Solve

More information

Uncertainty and auto-correlation in. Measurement

Uncertainty and auto-correlation in. Measurement Uncertanty and auto-correlaton n arxv:1707.03276v2 [physcs.data-an] 30 Dec 2017 Measurement Markus Schebl Federal Offce of Metrology and Surveyng (BEV), 1160 Venna, Austra E-mal: markus.schebl@bev.gv.at

More information

The Study of Teaching-learning-based Optimization Algorithm

The Study of Teaching-learning-based Optimization Algorithm Advanced Scence and Technology Letters Vol. (AST 06), pp.05- http://dx.do.org/0.57/astl.06. The Study of Teachng-learnng-based Optmzaton Algorthm u Sun, Yan fu, Lele Kong, Haolang Q,, Helongang Insttute

More information

U.C. Berkeley CS294: Beyond Worst-Case Analysis Luca Trevisan September 5, 2017

U.C. Berkeley CS294: Beyond Worst-Case Analysis Luca Trevisan September 5, 2017 U.C. Berkeley CS94: Beyond Worst-Case Analyss Handout 4s Luca Trevsan September 5, 07 Summary of Lecture 4 In whch we ntroduce semdefnte programmng and apply t to Max Cut. Semdefnte Programmng Recall that

More information

A Simple Inventory System

A Simple Inventory System A Smple Inventory System Lawrence M. Leems and Stephen K. Park, Dscrete-Event Smulaton: A Frst Course, Prentce Hall, 2006 Hu Chen Computer Scence Vrgna State Unversty Petersburg, Vrgna February 8, 2017

More information

Hongyi Miao, College of Science, Nanjing Forestry University, Nanjing ,China. (Received 20 June 2013, accepted 11 March 2014) I)ϕ (k)

Hongyi Miao, College of Science, Nanjing Forestry University, Nanjing ,China. (Received 20 June 2013, accepted 11 March 2014) I)ϕ (k) ISSN 1749-3889 (prnt), 1749-3897 (onlne) Internatonal Journal of Nonlnear Scence Vol.17(2014) No.2,pp.188-192 Modfed Block Jacob-Davdson Method for Solvng Large Sparse Egenproblems Hongy Mao, College of

More information

Solving Nonlinear Differential Equations by a Neural Network Method

Solving Nonlinear Differential Equations by a Neural Network Method Solvng Nonlnear Dfferental Equatons by a Neural Network Method Luce P. Aarts and Peter Van der Veer Delft Unversty of Technology, Faculty of Cvlengneerng and Geoscences, Secton of Cvlengneerng Informatcs,

More information

Week 5: Neural Networks

Week 5: Neural Networks Week 5: Neural Networks Instructor: Sergey Levne Neural Networks Summary In the prevous lecture, we saw how we can construct neural networks by extendng logstc regresson. Neural networks consst of multple

More information

A Local Variational Problem of Second Order for a Class of Optimal Control Problems with Nonsmooth Objective Function

A Local Variational Problem of Second Order for a Class of Optimal Control Problems with Nonsmooth Objective Function A Local Varatonal Problem of Second Order for a Class of Optmal Control Problems wth Nonsmooth Objectve Functon Alexander P. Afanasev Insttute for Informaton Transmsson Problems, Russan Academy of Scences,

More information

Supplement: Proofs and Technical Details for The Solution Path of the Generalized Lasso

Supplement: Proofs and Technical Details for The Solution Path of the Generalized Lasso Supplement: Proofs and Techncal Detals for The Soluton Path of the Generalzed Lasso Ryan J. Tbshran Jonathan Taylor In ths document we gve supplementary detals to the paper The Soluton Path of the Generalzed

More information

e - c o m p a n i o n

e - c o m p a n i o n OPERATIONS RESEARCH http://dxdoorg/0287/opre007ec e - c o m p a n o n ONLY AVAILABLE IN ELECTRONIC FORM 202 INFORMS Electronc Companon Generalzed Quantty Competton for Multple Products and Loss of Effcency

More information

829. An adaptive method for inertia force identification in cantilever under moving mass

829. An adaptive method for inertia force identification in cantilever under moving mass 89. An adaptve method for nerta force dentfcaton n cantlever under movng mass Qang Chen 1, Mnzhuo Wang, Hao Yan 3, Haonan Ye 4, Guola Yang 5 1,, 3, 4 Department of Control and System Engneerng, Nanng Unversty,

More information

Note 10. Modeling and Simulation of Dynamic Systems

Note 10. Modeling and Simulation of Dynamic Systems Lecture Notes of ME 475: Introducton to Mechatroncs Note 0 Modelng and Smulaton of Dynamc Systems Department of Mechancal Engneerng, Unversty Of Saskatchewan, 57 Campus Drve, Saskatoon, SK S7N 5A9, Canada

More information

How Strong Are Weak Patents? Joseph Farrell and Carl Shapiro. Supplementary Material Licensing Probabilistic Patents to Cournot Oligopolists *

How Strong Are Weak Patents? Joseph Farrell and Carl Shapiro. Supplementary Material Licensing Probabilistic Patents to Cournot Oligopolists * How Strong Are Weak Patents? Joseph Farrell and Carl Shapro Supplementary Materal Lcensng Probablstc Patents to Cournot Olgopolsts * September 007 We study here the specal case n whch downstream competton

More information

Lecture 12: Discrete Laplacian

Lecture 12: Discrete Laplacian Lecture 12: Dscrete Laplacan Scrbe: Tanye Lu Our goal s to come up wth a dscrete verson of Laplacan operator for trangulated surfaces, so that we can use t n practce to solve related problems We are mostly

More information

STOCHASTIC INVENTORY MODELS INVOLVING VARIABLE LEAD TIME WITH A SERVICE LEVEL CONSTRAINT * Liang-Yuh OUYANG, Bor-Ren CHUANG 1.

STOCHASTIC INVENTORY MODELS INVOLVING VARIABLE LEAD TIME WITH A SERVICE LEVEL CONSTRAINT * Liang-Yuh OUYANG, Bor-Ren CHUANG 1. Yugoslav Journal of Operatons Research 10 (000), Number 1, 81-98 STOCHASTIC INVENTORY MODELS INVOLVING VARIABLE LEAD TIME WITH A SERVICE LEVEL CONSTRAINT Lang-Yuh OUYANG, Bor-Ren CHUANG Department of Management

More information

NON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS

NON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS IJRRAS 8 (3 September 011 www.arpapress.com/volumes/vol8issue3/ijrras_8_3_08.pdf NON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS H.O. Bakodah Dept. of Mathematc

More information

ALGORITHM FOR THE CALCULATION OF THE TWO VARIABLES CUBIC SPLINE FUNCTION

ALGORITHM FOR THE CALCULATION OF THE TWO VARIABLES CUBIC SPLINE FUNCTION ANALELE ŞTIINŢIFICE ALE UNIVERSITĂŢII AL.I. CUZA DIN IAŞI (S.N.) MATEMATICĂ, Tomul LIX, 013, f.1 DOI: 10.478/v10157-01-00-y ALGORITHM FOR THE CALCULATION OF THE TWO VARIABLES CUBIC SPLINE FUNCTION BY ION

More information

Mean Field / Variational Approximations

Mean Field / Variational Approximations Mean Feld / Varatonal Appromatons resented by Jose Nuñez 0/24/05 Outlne Introducton Mean Feld Appromaton Structured Mean Feld Weghted Mean Feld Varatonal Methods Introducton roblem: We have dstrbuton but

More information

2.3 Nilpotent endomorphisms

2.3 Nilpotent endomorphisms s a block dagonal matrx, wth A Mat dm U (C) In fact, we can assume that B = B 1 B k, wth B an ordered bass of U, and that A = [f U ] B, where f U : U U s the restrcton of f to U 40 23 Nlpotent endomorphsms

More information

Second Order Analysis

Second Order Analysis Second Order Analyss In the prevous classes we looked at a method that determnes the load correspondng to a state of bfurcaton equlbrum of a perfect frame by egenvalye analyss The system was assumed to

More information

Hidden Markov Models

Hidden Markov Models Hdden Markov Models Namrata Vaswan, Iowa State Unversty Aprl 24, 204 Hdden Markov Model Defntons and Examples Defntons:. A hdden Markov model (HMM) refers to a set of hdden states X 0, X,..., X t,...,

More information

Comparison of the Population Variance Estimators. of 2-Parameter Exponential Distribution Based on. Multiple Criteria Decision Making Method

Comparison of the Population Variance Estimators. of 2-Parameter Exponential Distribution Based on. Multiple Criteria Decision Making Method Appled Mathematcal Scences, Vol. 7, 0, no. 47, 07-0 HIARI Ltd, www.m-hkar.com Comparson of the Populaton Varance Estmators of -Parameter Exponental Dstrbuton Based on Multple Crtera Decson Makng Method

More information

k t+1 + c t A t k t, t=0

k t+1 + c t A t k t, t=0 Macro II (UC3M, MA/PhD Econ) Professor: Matthas Kredler Fnal Exam 6 May 208 You have 50 mnutes to complete the exam There are 80 ponts n total The exam has 4 pages If somethng n the queston s unclear,

More information

A Bayes Algorithm for the Multitask Pattern Recognition Problem Direct Approach

A Bayes Algorithm for the Multitask Pattern Recognition Problem Direct Approach A Bayes Algorthm for the Multtask Pattern Recognton Problem Drect Approach Edward Puchala Wroclaw Unversty of Technology, Char of Systems and Computer etworks, Wybrzeze Wyspanskego 7, 50-370 Wroclaw, Poland

More information

A MODIFIED METHOD FOR SOLVING SYSTEM OF NONLINEAR EQUATIONS

A MODIFIED METHOD FOR SOLVING SYSTEM OF NONLINEAR EQUATIONS Journal of Mathematcs and Statstcs 9 (1): 4-8, 1 ISSN 1549-644 1 Scence Publcatons do:1.844/jmssp.1.4.8 Publshed Onlne 9 (1) 1 (http://www.thescpub.com/jmss.toc) A MODIFIED METHOD FOR SOLVING SYSTEM OF

More information

Lossy Compression. Compromise accuracy of reconstruction for increased compression.

Lossy Compression. Compromise accuracy of reconstruction for increased compression. Lossy Compresson Compromse accuracy of reconstructon for ncreased compresson. The reconstructon s usually vsbly ndstngushable from the orgnal mage. Typcally, one can get up to 0:1 compresson wth almost

More information

Basic Statistical Analysis and Yield Calculations

Basic Statistical Analysis and Yield Calculations October 17, 007 Basc Statstcal Analyss and Yeld Calculatons Dr. José Ernesto Rayas Sánchez 1 Outlne Sources of desgn-performance uncertanty Desgn and development processes Desgn for manufacturablty A general

More information

DETERMINATION OF TEMPERATURE DISTRIBUTION FOR ANNULAR FINS WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY BY HPM

DETERMINATION OF TEMPERATURE DISTRIBUTION FOR ANNULAR FINS WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY BY HPM Ganj, Z. Z., et al.: Determnaton of Temperature Dstrbuton for S111 DETERMINATION OF TEMPERATURE DISTRIBUTION FOR ANNULAR FINS WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY BY HPM by Davood Domr GANJI

More information

Support Vector Machines. Vibhav Gogate The University of Texas at dallas

Support Vector Machines. Vibhav Gogate The University of Texas at dallas Support Vector Machnes Vbhav Gogate he Unversty of exas at dallas What We have Learned So Far? 1. Decson rees. Naïve Bayes 3. Lnear Regresson 4. Logstc Regresson 5. Perceptron 6. Neural networks 7. K-Nearest

More information

Fundamental loop-current method using virtual voltage sources technique for special cases

Fundamental loop-current method using virtual voltage sources technique for special cases Fundamental loop-current method usng vrtual voltage sources technque for specal cases George E. Chatzaraks, 1 Marna D. Tortorel 1 and Anastasos D. Tzolas 1 Electrcal and Electroncs Engneerng Departments,

More information

A Robust Method for Calculating the Correlation Coefficient

A Robust Method for Calculating the Correlation Coefficient A Robust Method for Calculatng the Correlaton Coeffcent E.B. Nven and C. V. Deutsch Relatonshps between prmary and secondary data are frequently quantfed usng the correlaton coeffcent; however, the tradtonal

More information