A Fuzzy Model for the Multiobjective Emergency Facility Location Problem with A-Distance

Size: px
Start display at page:

Download "A Fuzzy Model for the Multiobjective Emergency Facility Location Problem with A-Distance"

Transcription

1 The Open Cybernecs and Sysemcs Journal, 007, 1, A Fuzzy Model for he Mulobecve Emergency Facly Locaon Problem wh A-Dsance T. Uno *, H. Kaagr and K. Kao Deparmen of Arfcal Complex Sysems Engneerng, Graduae School of Engneerng, Hroshma Unversy, Hgashhroshma, Hroshma, Japan Absrac: Ths paper proposes a new mulobecve locaon problem for emergency facles, e.g. ambulance servce saons and fre saons. In he problem, one of he obecves s o mnmze he maxmal dsance of he pahs from he locaed emergency facles o hospals va accdens. The oher s o maxmze he frequency of accdens ha he emergency facles can respond quckly. An neracve fuzzy sasfcng mehod wh parcle swarm opmzaon (PSO) mehod s proposed for fndng a sasfcng locaon for he problem. Compuaonal resuls llusrae he mehod wh numercal examples of he mulobecve locaon problem. INTRODUCTION In hs paper, we consder a new emergency facly locaon problem (EFLP), such ha ambulance servce saons [1-3], fre saons [1, 4, 5], ec. Masuom and Ish [] consdered EFLPs wh he suaon ha f an accden occurs, he neares emergency facly sends ambulances o and nured people are conveyed o he neares hospal. We propose a new mulobecve EFLP by exendng Masuom and Ish s EFLP. EFLPs have he followng wo mporan facors. One s dsance (or norm); for deals of he relaon beween facly locaon and norm, he reader can refer o he sudy of Marn [6]. There are wo norms wdely used n sudes abou he EFLPs. One s he Eucldean norm [7, 8], whch s assumed ha can be raveled o any orenaons a any pons. However, hs assumpon does no usually hold for facly locaon n urban areas. The oher s he block norm [9-11], whch s assumed ha can be raveled o gven several allowable orenaons of movemen wh weghs a any pons. Reclnear dsance [8] s regarded as one of he block norms such ha here are wo allowable orenaons whch cross a rgh angles wh he same weghs. The EFLPs wh he reclnear dsance are ofen suded [7, 1]. The A-dsance defned by Wdmayer e al. [13] s also regarded as one of block norms such ha here are several allowable orenaons of movemen wh he same weghs. Masuom and Ish [] consder an EFLP wh he A- dsance. In hs paper, we propose a new EFLP based on he EFLP wh he A-dsance. The oher s creron of opmaly for facly locaon. In general EFLPs [, 4, 7, 14, 15], an obecve for facly locaon s o mnmze he maxmal dsance beween emergency facles and he scenes of accdens. In hs paper, we nroduce anoher new obecve, whch s o maxmze *Address correspondence o hs auhor a he Deparmen of Arfcal Complex Sysems Engneerng, Graduae School of Engneerng, Hroshma Unversy, Hgashhroshma, Hroshma, Japan; E-mal: unoake@hroshma-u.ac.p frequency of accdens ha emergency facles can respond quckly. Then, a new mulobecve EFLP wh he wo obecves s formulaed. Mos mulobecve EFLPs do no have complee opmal soluons. For fndng a sasfcng soluon of he mulobecve EFLP for he decson maker (DM), we apply neracve fuzzy sasfcng mehod proposed by Sakawa and Yano [16]. Kaagr e al. [17] recenly proposed neracve fuzzy sasfcng mehod for mulobecve fuzzy random lnear programmng problems. In hs mehod, we need o fnd an opmal soluon for each of he mnmax problems wh he correspondng reference membershp values. Parcle swarm opmzaon (PSO), whch s proposed by Kennedy and Eberhar [18] and mproved by Masu e al. [19], s one of he effcen soluon mehods for nonlnear programmng problem. We propose o apply he PSO mproved by ulzng some characerscs of he EFLPs. The organzaon of he paper s as follows. In Secon, we gve he defnon of he A-dsance and s properes. In Secon 3, we formulae a mulobecve EFLP wh he A- dsance. For he formulaed EFLP, frs we propose he mehod o compue he obecve values for each locaon n Secon 4. In order o fnd a sasfcng soluon for he DM, we nroduce he neracve fuzzy sasfcng mehod proposed by Sakawa and Yano [16] n Secon 5. In order o solve he mnmax problems n hs mehod, we proposed a PSO mehod consderng characerscs of he EFLP n Secon 6. In secon 7, we show he resuls for applyng he mehod o numercal examples of he mulobecve EFLPs. Fnally, we make menon of conclusons and fuure remarks n Secon 8. A-DISTANCE In hs secon, we descrbe he defnon of A-dsance and s properes. We consder he suaon ha here are a orenaons whch can only move n he plane R. The orenaons are represened as he angles beween he correspondng sragh lnes o orenaons and he Caresan x-axs; for example, orenaon 0 s he x-axs and orenaon / s X/ Benham Scence Publshers Ld.

2 The Open Cybernecs and Sysemcs Journal, 007, Volume 1 Uno e al. he y-axs. Le A = { 1,..., a } be a se of orenaons such ha 0 1 a <. A lne, a half lne, or a lne segmen s called A-orened f s orenaon s one of hose n A. Then, he A-dsance beween wo pons p 1 and s represened as follows: p R d, p ), f p 1 and p are n an A-orened lne, d A, p ):= (1) mn d p 3 { R A, p 3 ) + d A (p 3, p )}, oherwse where d (, ) means he Eucldean dsance. Fg. (1) shows an example of he A-dsance beween p 1 and p R. The reclnear dsance s represened as he A-dsance wh A = { 0, / }. Fg. (). An A-crcle ( A = { / 3, /,3 / 4} ). The sdes and verces of he Vorono polygons are called Vorono edges and Vorono pons, respecvely. The se of all Vorono polygons, whch can be regarded as a paron of R, s called Vorono dagram wh he A-dsance. Vorono dagram for Q s denoed by VD A (Q). Fg. (1). An example of he A-dsance. Wdmayer e al. [13] show ha f pons p 1 and p are no n any A-orened lnes, here exss a leas one pon p 3 R such ha d A, p ) = d, p 3 ) + d (p 3, p ) () Theorem 1: A Vorono dagram n he A-dsance for Q can be consruced n O ( n log n) mes usng O (n) space, whch s asympocally opmal n he wors case. Proof: Correcness follows from he consderaon n Secon 3 of he reference [13]. Opmaly s due, for nsance, o a reducon of sorng [0]. Fg. (3) shows an example of he Vorono dagram wh A = { 0, / }. For a pon p and a dsance d > 0, he locus of all pons p ' wh d A ( p, p' ) = d s called he A-crcle wh cener p and radus d. As shown n Fg. (), A-crcle has s boundary of a-gon whose corner pons are he nersecons of he crcle wh cener p and radus d and he A-orened lnes hrough p. For wo pons p 1 and p, he bsecor of p 1 and p wh he A-dsance s defned as follows: B p, p ) = { p d ( p, p ) = d ( p, )}. (3) A( 1 A 1 A p Le Q = { p 1,..., p n } be a se of n pons p 1,..., p n R. Then, he Vorono polygon V A ( p ), = 1,..., n, wh he A- dsance s defned as follows: V A (p ) = p d A (p, p ) d A (p, p ) (4) { } Fg. (3). Vorono dagram wh he A-dsance. FORMULATION OF MULTI-OBJECTIVE EFLP In hs secon, we formulae a mulobecve EFLP wh he A-dsance. Le S R be a closed convex polygon n whch accdens occur and he DM needs o locae emergency facles. We consder he suaon ha f an accden occurs a a pon, he neares emergency facly o he pon sends ambulances o he pon and hen nured people n he accden are conveyed from he pon o he neares hospal.

3 Mulobecve Emergency Facly Locaon Problem wh A-Dsance The Open Cybernecs and Sysemcs Journal, 007, Volume 1 3 Frs, we show he mnmax creron abou pah from he emergency facles o he hospals va he pons of accdens. Le h 1,..., h m S be ses of m hospals, le y 1,..., y n S be ses of n emergency facles, and le n Y = ( y 1,..., yn ) S. Then, f an accden occurs a a pon p S, he A-dsance for he above pah s represened as follows: u(y, = mn d A (y, + mn d A (p, h ) (5) =1,...,n =1,...,m Because he DM does no know where accdens occur n S beforehand, one of our obecves s nerpreed as copng wh any accden pons n S quckly. Then, he frs obecve funcon s represened as follows: f 1 (Y ):= max u(y, (6) ps Secondly, we show a new creron abou frequency of accdens. We assume ha he DM knows pons where accdens frequenly occur n S, called accden pons. There are k accden pons whose ses are denoed by a 1,...,a k S, and each of her accden pons has a wegh abou frequency of accdens, denoed by w 1,..., w k > 0, respecvely. Le > 0 be an upper lm of he dsance from he emergency facles o he hospals such ha a medcal reamen for nured people can be n me. The oher of our obecves s nerpreed as maxmzng he sum of he weghs of frequency of accden pons ha he emergency facles can cover for a gven. Then, he second obecve funcon s represened as follows: f (Y ):= (7) w I (Y ) where I (Y ):= a u(y, { } (8) Therefore, a mulobecve EFLP s formulaed as follows: mnmze f 1 (Y ) maxmze f (Y ) subec o Y = (y 1,..., y n ) S n COMPUTATION OF THE OBJECTIVE VALUES OF MULTIOBJECTIVE EFLP In order o fnd an opmal soluon of (9), we need o compue he wo obecve values for each soluon. For he second obecve funcon, we can compue s obecve value by measurng he A-dsances from he emergency facles o he hospals va he k accdens pons. In he followng par of hs secon, we propose he mehod o compue he frs obecve value for each locaon. For he case n = 1, Masuom and Ish [] showed he followng heorem. Theorem : If n = 1, p S whch maxmzes u ( Y, s one of he followng pons: (9) Verces of he boundary of S, Inersecons of Vorono edges of each V h ),..., ( h ) and he boundary of S. If he DM locaes wo or more emergency facles, we need o consder whch of he emergency facles s used a any pons n S. From he defnon of Vorono polygon, V A ( y ) s he se of pons whch uses emergency facly. Whle p S whch maxmzes u ( Y, s only on he boundary of S f n = 1, p S whch maxmzes u ( Y, may be n he neror of S and on he Vorono edges of V A ( h 1 ),..., V A ( h m ) and V A ( y 1 ),..., V A ( y m ) f n. Then, Theorem 1 can be exended o he followng corollary: Corollary 1: If n, p S whch maxmzes u ( Y, s one of he followng pons: Verces of he boundary of S, Inersecons of Vorono edges of each V h ),..., ( h ) and he boundary of S, Vorono pons of each V h ),..., ( h ), Vorono pons of each V y ),..., ( y ), Inersecons of Vorono edges of each V y ),..., ( y ) and he boundary of S, Inersecons of Vorono edges of each V h ),..., ( h ) and Vorono edges of each V A ( 1 ),..., V A ( y m y ). The above pons can be found by drawng Vorono dagram for hospals and Vorono dagram for each locaon of emergency facles. Then, we can fnd he frs obecve value by compung he maxmal dsance for pahs from emergency facles o hospals va hese pons. INTERACTIVE FUZZY SATISFICING APPROACH In hs secon, we nroduce he neracve fuzzy sasfcng mehod proposed by Sakawa and Yano [16] n order o fnd a sasfcng soluon of (9) for he DM. For decson makng n real world, he DM usually prefer o make an obecve funcon value more/less han a ceran value raher han o maxmze/mnmze s obecve funcon value. Such an obecve, called a fuzzy obecve, ncludes vagueness based upon udgmen of he DM. In hs paper, we represen he wo obecves of (9) as fuzzy obecves provded by membershp funcons, denoed by μ 1 and μ. Now we nroduce an example of membershp funcons for each obecve funcon. Le d e denoe he dsance such ha he DM s que sasfed f he frs obecve value s less han d e, and d l denoe he dsance such ha she/he s sasfed o a ceran degree f s obecve value s more han d bu less han e

4 4 The Open Cybernecs and Sysemcs Journal, 007, Volume 1 Uno e al. d l. Then, we use he followng lnear membershp funcon for he former obecve: 1, f f 1 (Y ) < d e, f μ 1 ( f 1 (Y )) := 1 (Y ) d e, f d e f 1 (Y ) < d l, d l d e 0, f d l f 1 (Y ) (10) Nex, one of he smples ways o provde membershp funcon for he laer obecve s as follows: μ ( f (Y )) := f (Y ) (11) k w =1 Then, (9) s ransformed as he followng mulobecve fuzzy programmng problem: mnmze μ 1 ( f 1 (Y )) maxmze μ (f (Y )) subec o Y = (y 1,..., y n ) S n (1) Snce here s generally no complee opmal soluon for mulobecve programmng problem ncludng (1), he concep of he M-Pareo opmal soluon s usually used for mulobecve fuzzy programmng problems. Defnon 1: Soluon Y * s an M-Pareo opmal soluon of (1) f and only f here does no exs any soluons n S Y such ha μ ( f ( Y )) μ ( f ( Y*)) for all = 1, and μ ( f ( Y )) > μ ( f ( Y*)) for a leas one {1, }. The neracve fuzzy sasfcng mehod [16] s o fnd a sasfcng M-Pareo opmal soluon hrough neracon o he DM. Le ( μ 1, μ) be a par of nal reference membershp levels of membershp funcon μ 1 and μ, respecvely. The neracve fuzzy sasfcng mehod for (1) can be descrbed as follows: Algorhm 1: Ineracve Fuzzy Sasfcng Mehod Sep 1. Provde wo membershp funcons μ 1 and μ accordng o (10) and (11). Sep. Se he nal reference membershp levels μ, μ ) (1,1). Sep 3. ( 1 = For he gven par of reference membershp levels ( μ 1, μ), solve he followng correspondng mnmax problem: mnmze max{μ μ ( f (Y )) =1, + (μ μ ( f (Y )))} =1 subec o Y = (y 1,..., y n ) S n where, s a suffcenly small posve number. Sep 1. (13) If he DM s sasfed wh he curren levels of he M-Pareo opmal soluon, STOP. Then he curren M-Pareo opmal soluon s a sasfcng soluon for he DM. Sep. Updae he par of curren reference membershp levels ( μ 1, μ) based on nformaon of preference of he DM, he curren values of he membershp funcons, ec. Reurn o Sep 3. In he neracve fuzzy sasfcng mehod, we need o solve he mnmax problems n Sep 3 effcenly. In he nex secon, we propose an effcen soluon mehod for he mnmax problem. A SOLUTION ALGORITHM FOR MINIMAX PROB- LEM A PSO mehod proposed by Kennedy and Eberhar [18] s based on he socal behavor ha a populaon of ndvduals adaps o s envronmen by reurnng o promsng regons ha were prevously dscovered [1]. Ths adapaon o he envronmen s a sochasc process ha depends upon boh he memory of each ndvdual, called parcle, and he knowledge ganed by he populaon, called swarm. In he numercal mplemenaon of hs smplfed socal model, each parcle has he followng hree arbues: he poson vecor n he search space, he velocy vecor and he bes poson n s rack, and he bes poson of he swarm. The process can be oulned as follows. Algorhm : Oulne of he PSO Mehod Sep 1. Generae he nal swarm nvolvng N parcles a random. Sep. Calculae he new velocy vecor for each parcle, based on s arbues. Sep 3. Calculae he new poson of each parcle from he curren poson and s new velocy vecor. Sep 4. If he ermnal condon s sasfed, STOP. The bes soluon gven n he searchng hsory s an approxmae opmal soluon. Oherwse, go o Sep. To be more specfc, for he poson and he velocy vecor of he -h parcle a me, denoed by x and v, respecvely, he new velocy vecor of he -h parcle a me + 1 s calculaed by he followng scheme nroduced by Sh and Eberhar []. v + 1 : 1 1 g where me, and = v + c R ( p x ) + c R ( p x ) (14) R 1 and R are random numbers beween 0 and 1 a p s he bes poson of he -h parcle n s rack p g s he bes poson of he swarm. There are hree problem-dependen parameers, he nera of he parcle, and wo rus parameers c 1 and c. Then, he new poson of he -h parcle a me + 1 s calculaed from he followng equaon: + 1 : + 1 x = x + v (15)

5 Mulobecve Emergency Facly Locaon Problem wh A-Dsance The Open Cybernecs and Sysemcs Journal, 007, Volume 1 5 The -h parcle calculaes he nex search drecon vecor v by (14) n consderaon of he curren search drec- +1 on vecor v, he drecon vecor gong from he curren search poson x o he bes poson n s rack p and he drecon vecor gong from he curren search poson x o he bes poson of he swarm p g, moves from he +1 curren poson x o he nex search poson x calculaed by (14). The parameer conrols he amoun of he move o search globally n early sage and o search locally by decreasng gradually. I s defned by follows: poson of he swarm. Moreover, hey proposed he mulplex srechng mehod, whch s he exended verson of he srechng mehod proposed by Parsopoulos and Varahas [3]. In order o sele he second problem, Masu e al. [19] proposed o generae nal parcles n he feasble se by ulzng he homomorphsm proposed by Kozel and Mchalewcz [4]. Moreover, here are ofen cases ha a parcle afer move s no always nfeasble f we use he updang equaon of search poson menoned above. To deal wh such a suaon, Masu e al. [19] dvded he swarm no wo subswarms; one s he move of a parcle o he nfeasble regon s acceped and he oher s no. 0 T max 0 ( ) : = (16) 0.75T max 0 where T max s he number of maxmum eraon mes, Tmax s an nal value a he me eraon, and s he las value. The searchng procedure of PSO s shown n Fg. (4). Fg. (4). Movemen model for PSO. Comparng he evaluaon value of a parcle afer +1 movemen, denoed by f ( x ), wh ha of he bes poson n s rack, denoed by f ( p ), f f ( x ) s beer +1 han ) f ( x f ( p ), hen he bes poson n s rack s updaed as f ( p ). Oherwse, hen he bes poson n he + 1 swarm s updaed as f ( p ) f ( p ). In hs way, a parcle ges nformaon of he bes poson of new oneself and he swarm, and moves accordng o (14) and (15) agan, and searches newly. The summary of he PSO mehod s shown n Fg. (5). Such a PSO echnque ncludes wo drawbacks. One s ha parcles concenrae on he bes search poson of he swarm and hey canno easly escape from he local opmal +1 soluon snce he move drecon vecor v calculaed by (14) always ncludes he drecon vecor o he bes search poson of he swarm. Anoher s ha a parcle afer move s no always feasble for problems wh consrans. In order o sele he frs ssue, Masu e al. [19] proposed he leavng acs for parcles whch are on he bes Fg. (5). Summary of PSO algorhm. IMPROVEMENT OF PSO METHOD In hs sudy, based upon he PSO mehod nroduced n he prevous secon, we proposed o add a new fourh movemen n (14). Le q S be maxmzer of u ( x,, such a pon can be found by usng Theorem 1. In (9) approachng he neares emergency facly o q, obecve value of f 1( Y ) may be mproved. Therefore we nroduce such a movemen o our soluon mehod. Le p a be a poson such ha for x, a se of he neares facly s changed o q and ses of he oher facles are fxed. Then, our proposng new velocy vecor of he -h parcle a me + 1 s represened as follows: v +1 := v + c 1 R 1 (p x ) + c R (p g x ) +c 3 R 3 (p a x ) where c 3 s a rus parameer and beween 0 and 1. NUMERICAL EXPERIMENTS (17) R 3 s a random number In hs secon, we apply neracve fuzzy sasfcng approach wh he PSO mehod o an example of our proposng mulobecve EFLPs. In hs example, we consder an EFLP

6 6 The Open Cybernecs and Sysemcs Journal, 007, Volume 1 Uno e al. for wo emergency facles, ha s n =. We represen S as a convex hull ncludng 0 pons gven n [ 0,100] [0,100] randomly. For he A-dsance, we se A = { 0, / 4, /,3 / 4}. For hospals, we se m = 3 and her ses are gven n S randomly. For he frequency of accdens, we se = 15, and for each of 100 accden pons, s se and wegh are randomly gven n S and (0,1], respecvely. We llusrae he neracve fuzzy sasfcng approach for he above example of he mulobecve EFLP. For parameers of PSO, we use he same parameer as he sudy of Masu e al. [19], ha s, s populaon sze s 40, s generaon s 500, and c 1 = c =. Moreover, we se he new parameer c 3 = A Sep 1, n order o represen fuzzness abou wo obecves, we use membershp funcons (10) and (11) n Secon 5 wh seng d = 5 and d = 10. e A Sep 3, we solve mnmax problem (13) for each gven μ, ) by solvng he PSO mehod n secon 6. We se ( 1 μ = In order o verfy effcency of PSO mehod, we apply he soluon mehod GENOCOP [4], whch s a genec algorhm for numercal opmzaon for consraned problem proposed by Kozel and Mchalewcz [4], o mnmax problems. For parameers of GENOCOP, we se he same populaon sze and generaon. Moreover, we se he oher parameers smlarly o he sudy of Kozel and Mchalewcz [4]. Compuaonal resuls for each gven ( μ 1, μ) a 0 mes by PSO and GENOCOP are shown n Tables 1 and, respecvely, where hese resuls are gven by usng DELL Opplex GX60 (CPU: Penum(R) 4.33GHz, RAM: 51MB). Table 1. Compuaonal Resuls by PSO l Table. Compuaonal Resuls by GENOCOP Mnmax Problem 1 3 μ μ Bes Mean Wors Mean CPU Tme (Sec) curren reference membershp levels ( μ 1, μ) by consderng he curren values of he membershp funcons, and resolve mnmax problem o Sep 3. In hs example of EFLP, we assume ha he DM hnks ha he frs obecve s more mporan han he second obecve. Then he DM hopes o mprove he value of μ 1 even f he value of μ s changed for he worse. However, he DM does no hope o go he value of μ oo bad. Then, an example of he neracve fuzzy sasfcng mehods s shown n Table 3. Table 3. Resuls of Ineracve Fuzzy Sasfcng Approach Ieraon 1 3 μ μ μ 1 ( f 1 (Y*)) μ ( f (Y*)) Mnmax Problem 1 3 μ μ Bes Mean Wors Mean CPU Tme (Sec) From Tables 1 and, PSO can fnd beer soluons han GENOCOP by meanngs of boh mean and sably. Ths means ha effcency of PSO for such mnmax problems. A Sep 4, he DM evaluaes wheher he M-Pareo opmal soluon gven by solvng he mnmax problem a Sep 3 s sasfed or no. If s soluon sasfes he DM, hs algorhm s ermnaed. Oherwse, ask he DM o updae he In Table 3, he DM consders as follows: A Ieraon 1, he DM s no sasfed M-Pareo opmal soluon because he frs obecve value s bad. A Ieraon, for mprovng he frs obecve value, he DM decreases μ, whch s reference membershp levels for he second obecve. Then, he M- Pareo opmal soluon gven a Ieraon s mproved for he frs obecve. However, she/he s no sasfed because he second obecve value s oo bad. A Ieraon 3, for mprovng he second obecve value a lle, he DM decreases μ 1, whch s reference membershp levels for he frs obecve. Then, she/he obans a sasfcng soluon for boh wo obecves, so he algorhm s ermnaed. CONCLUSIONS AND FUTURE STUDIES In hs paper, we have proposed a new mulobecve EFLP wh he A-dsance. In order o oban a sasfcng

7 Mulobecve Emergency Facly Locaon Problem wh A-Dsance The Open Cybernecs and Sysemcs Journal, 007, Volume 1 7 soluon for he DM, we have proposed an neracve fuzzy sasfcng mehod whch nvolves he procedure of solvng mnmax problems by he PSO mehod. By applyng an example of mulobecve EFLPs, we showed effcency of PSO and llusraed he neracve fuzzy sasfcng mehod. In he mulobecve EFLPs, we assume ha he regon of facly locaon S s a convex polygon because Theorem and Corollary 1 use he assumpon. However, n order o apply he EFLPs o more general cases, we need o consder varous shapes of S whch are non-convex, non-conneced, ec. To consruc our soluon mehod for general shapes of S s a fuure sudy. Moreover, f S ncludes many hospals and he DM locaes many emergency facles n S, we need o fnd an opmal soluon more effcenly for he mnmax problems n he neracve fuzzy sasfcng mehods. To consder an effcen soluon mehod for large-scale mulobecve EFLPs s also a fuure sudy. REFERENCES [1] C. Araz, H. Selm, and I. Ozkarahan, A fuzzy mul-obecve coverng-based vehcle locaon model for emergency servces, Compu. Opera. Res., vol. 34, no. 3, pp , March 007. [] T. Masuom and H. Ish, Mnmax locaon problem wh A- dsance, J. Opl. Res. Soc., vol. 41, pp , June [3] H. K. Raagopalan, C. Saydam, and J. Xao, A mulperod se coverng locaon model for dynamc redeploymen of ambulances, Compu. Opera. Res., vol. 35, no. 3, pp , March 008. [4] D. R. Plane and T. E. Hendrc, Mahemacal programmng and he locaon of fre companes for he Denver fre deparmen, Opns. Res., vol. 5, pp , July-Augus [5] L. Yang, B. F. Jones, and S. H. Yang, A fuzzy mul-obecve programmng for opmzaon of fre saon locaons hrough genec algorhms, Eur. J. Oper. Res., vol. 181, no., pp , Sepember 007. [6] H. Marn and A. Schöbel, Medan hyperplanes n normed spaces - a survey, Dscree Appl. Mah., vol. 89, pp , December [7] J. Elznga and D. W. Hearn, Geomercal soluons for some mnmax locaon problems, Trans. Sc., vol. 6, pp , November 197. [8] S. Özen and S. Güne, Effec of feaure-ype n selecng dsance measure for an arfcal mmune sysem as a paern recognzer Dgal Sgnal Processng (In Press.) [9] B. Pelegrn and F. R. Fernandez, Deermnaon of effcen pons n mulple-obecve locaon problems, Navel Res. Log., vol. 35, pp , December [10] J. F. Thsse, J. E. Hendrc, and R. E. Wendell, Some properes of locaon problems wh block and round norm, Opns. Res., vol. 3, pp , November-December [11] J. E. Ward and R. E. Wendell, Usng block norm for locaon modelng, Opns. Res., vol. 33, pp , Sepember-Ocober [1] G. O. Wesolowsky, Recangular dsance locaon under he mnmax opmaly creron, Trans. Sc., vol. 6, pp , February 197. [13] P. Wdmayer, Y. F. Wu, and C. K. Wong, On some dsance problems n fxed orenaons, SIAM J. Compu., vol. 16, pp , Augus [14] R. L. Francs, A geomercal soluon procedure for a reclnear mnmax locaon problem, AIIE Trans., vol. 4, pp , December 197. [15] D. K. Kulshresha, A mn-max locaon problem wh demand pons arbrarly dsrbued n a compac conneced space, J. Opl. Res. Soc., vol. 38, pp , May [16] M. Sakawa and H. Yano, An neracve fuzzy sasfcng mehod usng augmened mnmax problems and s applcaon o envronmenal sysems, IEEE Trans. Sys. Man Cybern., vol. SMC-15, pp , November-December [17] H. Kaagr, M. Sakawa, K. Kao, and I. Nshzak, Ineracve mulobecve fuzzy random lnear programmng: Maxmzaon of possbly and probably, Eur. J. Oper. Res., (In Press.) [18] J. Kennedy and R. C. Eberhar, Parcle swarm opmzaon, n Proc. of IEEE In. Conf. Neural Neworks, Pscaaway, NJ, pp , [19] T. Masu, M. Sakawa, T. Uno, K. Kao, M. Hgashmor, and M. Kaneko, Parcle swarm opmzaon for ump hegh maxmzaon of a seral lnk robo, Journal of Advanced Compuaonal Inellgence and Inellgen Informacs, vol. 11, no. 8, pp , Ocorber 007. [0] F. Hwang, An O(n log n) algorhm for reclnear mnmal spannng rees, J. Assoc. Compu. Mach., vol. 6, pp , Aprl [1] J. Kennedy, W. M. Spears, Machng algorhms o problems: an expermenal es of he parcle swarm and some genec algorhms on he mulmodal problem generaor, n Proc. of IEEE In. Conf. Evoluonary Compuaon, Anchorage, Alaska, [] Y. Sh and R. C. Eberhar, A modfed parcle swarm opmzer, n Proc. of IEEE In. Conf. Evoluonary Compuaon, Anchorage, Alaska, [3] K. E. Parsopoulos and M. N. Varahas, Recen approaches o global opmzaon problems hrough Parcle Swarm Opmzaon, Na. Compu., vol. 1, pp , June 00. [4] S. Kozel and Z. Mchalewcz, Evoluonary algorhms, homomorphous mappngs, and consraned parameer opmzaon, Evol. Compu., vol. 7, no. 1, pp , Sprng Receved: Ocober 5, 007 Revsed: Ocober 9, 007 Acceped: November 8, 007

PARTICLE SWARM OPTIMIZATION FOR INTERACTIVE FUZZY MULTIOBJECTIVE NONLINEAR PROGRAMMING. T. Matsui, M. Sakawa, K. Kato, T. Uno and K.

PARTICLE SWARM OPTIMIZATION FOR INTERACTIVE FUZZY MULTIOBJECTIVE NONLINEAR PROGRAMMING. T. Matsui, M. Sakawa, K. Kato, T. Uno and K. Scenae Mahemacae Japoncae Onlne, e-2008, 1 13 1 PARTICLE SWARM OPTIMIZATION FOR INTERACTIVE FUZZY MULTIOBJECTIVE NONLINEAR PROGRAMMING T. Masu, M. Sakawa, K. Kao, T. Uno and K. Tamada Receved February

More information

/07/$ IEEE

/07/$ IEEE An Application of Interactive Fuzzy Satisficing Approach with Particle Swarm Optimization for Multiobjective Emergency Facility Location Problem with A-distance Takeshi Uno, Kosuke Kato, Hideki Katagiri

More information

GENERATING CERTAIN QUINTIC IRREDUCIBLE POLYNOMIALS OVER FINITE FIELDS. Youngwoo Ahn and Kitae Kim

GENERATING CERTAIN QUINTIC IRREDUCIBLE POLYNOMIALS OVER FINITE FIELDS. Youngwoo Ahn and Kitae Kim Korean J. Mah. 19 (2011), No. 3, pp. 263 272 GENERATING CERTAIN QUINTIC IRREDUCIBLE POLYNOMIALS OVER FINITE FIELDS Youngwoo Ahn and Kae Km Absrac. In he paper [1], an explc correspondence beween ceran

More information

ISSN MIT Publications

ISSN MIT Publications MIT Inernaonal Journal of Elecrcal and Insrumenaon Engneerng Vol. 1, No. 2, Aug 2011, pp 93-98 93 ISSN 2230-7656 MIT Publcaons A New Approach for Solvng Economc Load Dspach Problem Ansh Ahmad Dep. of Elecrcal

More information

Approximate Analytic Solution of (2+1) - Dimensional Zakharov-Kuznetsov(Zk) Equations Using Homotopy

Approximate Analytic Solution of (2+1) - Dimensional Zakharov-Kuznetsov(Zk) Equations Using Homotopy Arcle Inernaonal Journal of Modern Mahemacal Scences, 4, (): - Inernaonal Journal of Modern Mahemacal Scences Journal homepage: www.modernscenfcpress.com/journals/jmms.aspx ISSN: 66-86X Florda, USA Approxmae

More information

HEAT CONDUCTION PROBLEM IN A TWO-LAYERED HOLLOW CYLINDER BY USING THE GREEN S FUNCTION METHOD

HEAT CONDUCTION PROBLEM IN A TWO-LAYERED HOLLOW CYLINDER BY USING THE GREEN S FUNCTION METHOD Journal of Appled Mahemacs and Compuaonal Mechancs 3, (), 45-5 HEAT CONDUCTION PROBLEM IN A TWO-LAYERED HOLLOW CYLINDER BY USING THE GREEN S FUNCTION METHOD Sansław Kukla, Urszula Sedlecka Insue of Mahemacs,

More information

Dynamic Team Decision Theory. EECS 558 Project Shrutivandana Sharma and David Shuman December 10, 2005

Dynamic Team Decision Theory. EECS 558 Project Shrutivandana Sharma and David Shuman December 10, 2005 Dynamc Team Decson Theory EECS 558 Proec Shruvandana Sharma and Davd Shuman December 0, 005 Oulne Inroducon o Team Decson Theory Decomposon of he Dynamc Team Decson Problem Equvalence of Sac and Dynamc

More information

Existence and Uniqueness Results for Random Impulsive Integro-Differential Equation

Existence and Uniqueness Results for Random Impulsive Integro-Differential Equation Global Journal of Pure and Appled Mahemacs. ISSN 973-768 Volume 4, Number 6 (8), pp. 89-87 Research Inda Publcaons hp://www.rpublcaon.com Exsence and Unqueness Resuls for Random Impulsve Inegro-Dfferenal

More information

Outline. Probabilistic Model Learning. Probabilistic Model Learning. Probabilistic Model for Time-series Data: Hidden Markov Model

Outline. Probabilistic Model Learning. Probabilistic Model Learning. Probabilistic Model for Time-series Data: Hidden Markov Model Probablsc Model for Tme-seres Daa: Hdden Markov Model Hrosh Mamsuka Bonformacs Cener Kyoo Unversy Oulne Three Problems for probablsc models n machne learnng. Compung lkelhood 2. Learnng 3. Parsng (predcon

More information

Solution in semi infinite diffusion couples (error function analysis)

Solution in semi infinite diffusion couples (error function analysis) Soluon n sem nfne dffuson couples (error funcon analyss) Le us consder now he sem nfne dffuson couple of wo blocks wh concenraon of and I means ha, n a A- bnary sysem, s bondng beween wo blocks made of

More information

Cubic Bezier Homotopy Function for Solving Exponential Equations

Cubic Bezier Homotopy Function for Solving Exponential Equations Penerb Journal of Advanced Research n Compung and Applcaons ISSN (onlne: 46-97 Vol. 4, No.. Pages -8, 6 omoopy Funcon for Solvng Eponenal Equaons S. S. Raml *,,. Mohamad Nor,a, N. S. Saharzan,b and M.

More information

On One Analytic Method of. Constructing Program Controls

On One Analytic Method of. Constructing Program Controls Appled Mahemacal Scences, Vol. 9, 05, no. 8, 409-407 HIKARI Ld, www.m-hkar.com hp://dx.do.org/0.988/ams.05.54349 On One Analyc Mehod of Consrucng Program Conrols A. N. Kvko, S. V. Chsyakov and Yu. E. Balyna

More information

In the complete model, these slopes are ANALYSIS OF VARIANCE FOR THE COMPLETE TWO-WAY MODEL. (! i+1 -! i ) + [(!") i+1,q - [(!

In the complete model, these slopes are ANALYSIS OF VARIANCE FOR THE COMPLETE TWO-WAY MODEL. (! i+1 -! i ) + [(!) i+1,q - [(! ANALYSIS OF VARIANCE FOR THE COMPLETE TWO-WAY MODEL The frs hng o es n wo-way ANOVA: Is here neracon? "No neracon" means: The man effecs model would f. Ths n urn means: In he neracon plo (wh A on he horzonal

More information

Variants of Pegasos. December 11, 2009

Variants of Pegasos. December 11, 2009 Inroducon Varans of Pegasos SooWoong Ryu bshboy@sanford.edu December, 009 Youngsoo Cho yc344@sanford.edu Developng a new SVM algorhm s ongong research opc. Among many exng SVM algorhms, we wll focus on

More information

Lecture 11 SVM cont

Lecture 11 SVM cont Lecure SVM con. 0 008 Wha we have done so far We have esalshed ha we wan o fnd a lnear decson oundary whose margn s he larges We know how o measure he margn of a lnear decson oundary Tha s: he mnmum geomerc

More information

V.Abramov - FURTHER ANALYSIS OF CONFIDENCE INTERVALS FOR LARGE CLIENT/SERVER COMPUTER NETWORKS

V.Abramov - FURTHER ANALYSIS OF CONFIDENCE INTERVALS FOR LARGE CLIENT/SERVER COMPUTER NETWORKS R&RATA # Vol.) 8, March FURTHER AALYSIS OF COFIDECE ITERVALS FOR LARGE CLIET/SERVER COMPUTER ETWORKS Vyacheslav Abramov School of Mahemacal Scences, Monash Unversy, Buldng 8, Level 4, Clayon Campus, Wellngon

More information

MANY real-world applications (e.g. production

MANY real-world applications (e.g. production Barebones Parcle Swarm for Ineger Programmng Problems Mahamed G. H. Omran, Andres Engelbrech and Ayed Salman Absrac The performance of wo recen varans of Parcle Swarm Opmzaon (PSO) when appled o Ineger

More information

CS434a/541a: Pattern Recognition Prof. Olga Veksler. Lecture 4

CS434a/541a: Pattern Recognition Prof. Olga Veksler. Lecture 4 CS434a/54a: Paern Recognon Prof. Olga Veksler Lecure 4 Oulne Normal Random Varable Properes Dscrmnan funcons Why Normal Random Varables? Analycally racable Works well when observaon comes form a corruped

More information

Robustness Experiments with Two Variance Components

Robustness Experiments with Two Variance Components Naonal Insue of Sandards and Technology (NIST) Informaon Technology Laboraory (ITL) Sascal Engneerng Dvson (SED) Robusness Expermens wh Two Varance Componens by Ana Ivelsse Avlés avles@ns.gov Conference

More information

THE PREDICTION OF COMPETITIVE ENVIRONMENT IN BUSINESS

THE PREDICTION OF COMPETITIVE ENVIRONMENT IN BUSINESS THE PREICTION OF COMPETITIVE ENVIRONMENT IN BUSINESS INTROUCTION The wo dmensonal paral dfferenal equaons of second order can be used for he smulaon of compeve envronmen n busness The arcle presens he

More information

M. Y. Adamu Mathematical Sciences Programme, AbubakarTafawaBalewa University, Bauchi, Nigeria

M. Y. Adamu Mathematical Sciences Programme, AbubakarTafawaBalewa University, Bauchi, Nigeria IOSR Journal of Mahemacs (IOSR-JM e-issn: 78-578, p-issn: 9-765X. Volume 0, Issue 4 Ver. IV (Jul-Aug. 04, PP 40-44 Mulple SolonSoluons for a (+-dmensonalhroa-sasuma shallow waer wave equaon UsngPanlevé-Bӓclund

More information

WiH Wei He

WiH Wei He Sysem Idenfcaon of onlnear Sae-Space Space Baery odels WH We He wehe@calce.umd.edu Advsor: Dr. Chaochao Chen Deparmen of echancal Engneerng Unversy of aryland, College Par 1 Unversy of aryland Bacground

More information

Clustering (Bishop ch 9)

Clustering (Bishop ch 9) Cluserng (Bshop ch 9) Reference: Daa Mnng by Margare Dunham (a slde source) 1 Cluserng Cluserng s unsupervsed learnng, here are no class labels Wan o fnd groups of smlar nsances Ofen use a dsance measure

More information

An introduction to Support Vector Machine

An introduction to Support Vector Machine An nroducon o Suppor Vecor Machne 報告者 : 黃立德 References: Smon Haykn, "Neural Neworks: a comprehensve foundaon, second edon, 999, Chaper 2,6 Nello Chrsann, John Shawe-Tayer, An Inroducon o Suppor Vecor Machnes,

More information

J i-1 i. J i i+1. Numerical integration of the diffusion equation (I) Finite difference method. Spatial Discretization. Internal nodes.

J i-1 i. J i i+1. Numerical integration of the diffusion equation (I) Finite difference method. Spatial Discretization. Internal nodes. umercal negraon of he dffuson equaon (I) Fne dfference mehod. Spaal screaon. Inernal nodes. R L V For hermal conducon le s dscree he spaal doman no small fne spans, =,,: Balance of parcles for an nernal

More information

Genetic Algorithm in Parameter Estimation of Nonlinear Dynamic Systems

Genetic Algorithm in Parameter Estimation of Nonlinear Dynamic Systems Genec Algorhm n Parameer Esmaon of Nonlnear Dynamc Sysems E. Paeraks manos@egnaa.ee.auh.gr V. Perds perds@vergna.eng.auh.gr Ah. ehagas kehagas@egnaa.ee.auh.gr hp://skron.conrol.ee.auh.gr/kehagas/ndex.hm

More information

SOME NOISELESS CODING THEOREMS OF INACCURACY MEASURE OF ORDER α AND TYPE β

SOME NOISELESS CODING THEOREMS OF INACCURACY MEASURE OF ORDER α AND TYPE β SARAJEVO JOURNAL OF MATHEMATICS Vol.3 (15) (2007), 137 143 SOME NOISELESS CODING THEOREMS OF INACCURACY MEASURE OF ORDER α AND TYPE β M. A. K. BAIG AND RAYEES AHMAD DAR Absrac. In hs paper, we propose

More information

EEL 6266 Power System Operation and Control. Chapter 5 Unit Commitment

EEL 6266 Power System Operation and Control. Chapter 5 Unit Commitment EEL 6266 Power Sysem Operaon and Conrol Chaper 5 Un Commmen Dynamc programmng chef advanage over enumeraon schemes s he reducon n he dmensonaly of he problem n a src prory order scheme, here are only N

More information

Planar truss bridge optimization by dynamic programming and linear programming

Planar truss bridge optimization by dynamic programming and linear programming IABSE-JSCE Jon Conference on Advances n Brdge Engneerng-III, Augus 1-, 015, Dhaka, Bangladesh. ISBN: 978-984-33-9313-5 Amn, Oku, Bhuyan, Ueda (eds.) www.abse-bd.org Planar russ brdge opmzaon by dynamc

More information

Relative controllability of nonlinear systems with delays in control

Relative controllability of nonlinear systems with delays in control Relave conrollably o nonlnear sysems wh delays n conrol Jerzy Klamka Insue o Conrol Engneerng, Slesan Techncal Unversy, 44- Glwce, Poland. phone/ax : 48 32 37227, {jklamka}@a.polsl.glwce.pl Keywor: Conrollably.

More information

Math 128b Project. Jude Yuen

Math 128b Project. Jude Yuen Mah 8b Proec Jude Yuen . Inroducon Le { Z } be a sequence of observed ndependen vecor varables. If he elemens of Z have a on normal dsrbuon hen { Z } has a mean vecor Z and a varancecovarance marx z. Geomercally

More information

Attribute Reduction Algorithm Based on Discernibility Matrix with Algebraic Method GAO Jing1,a, Ma Hui1, Han Zhidong2,b

Attribute Reduction Algorithm Based on Discernibility Matrix with Algebraic Method GAO Jing1,a, Ma Hui1, Han Zhidong2,b Inernaonal Indusral Informacs and Compuer Engneerng Conference (IIICEC 05) Arbue educon Algorhm Based on Dscernbly Marx wh Algebrac Mehod GAO Jng,a, Ma Hu, Han Zhdong,b Informaon School, Capal Unversy

More information

Robust and Accurate Cancer Classification with Gene Expression Profiling

Robust and Accurate Cancer Classification with Gene Expression Profiling Robus and Accurae Cancer Classfcaon wh Gene Expresson Proflng (Compuaonal ysems Bology, 2005) Auhor: Hafeng L, Keshu Zhang, ao Jang Oulne Background LDA (lnear dscrmnan analyss) and small sample sze problem

More information

Comb Filters. Comb Filters

Comb Filters. Comb Filters The smple flers dscussed so far are characered eher by a sngle passband and/or a sngle sopband There are applcaons where flers wh mulple passbands and sopbands are requred Thecomb fler s an example of

More information

Online Supplement for Dynamic Multi-Technology. Production-Inventory Problem with Emissions Trading

Online Supplement for Dynamic Multi-Technology. Production-Inventory Problem with Emissions Trading Onlne Supplemen for Dynamc Mul-Technology Producon-Invenory Problem wh Emssons Tradng by We Zhang Zhongsheng Hua Yu Xa and Baofeng Huo Proof of Lemma For any ( qr ) Θ s easy o verfy ha he lnear programmng

More information

FTCS Solution to the Heat Equation

FTCS Solution to the Heat Equation FTCS Soluon o he Hea Equaon ME 448/548 Noes Gerald Reckenwald Porland Sae Unversy Deparmen of Mechancal Engneerng gerry@pdxedu ME 448/548: FTCS Soluon o he Hea Equaon Overvew Use he forward fne d erence

More information

CHAPTER 10: LINEAR DISCRIMINATION

CHAPTER 10: LINEAR DISCRIMINATION CHAPER : LINEAR DISCRIMINAION Dscrmnan-based Classfcaon 3 In classfcaon h K classes (C,C,, C k ) We defned dscrmnan funcon g j (), j=,,,k hen gven an es eample, e chose (predced) s class label as C f g

More information

Tight results for Next Fit and Worst Fit with resource augmentation

Tight results for Next Fit and Worst Fit with resource augmentation Tgh resuls for Nex F and Wors F wh resource augmenaon Joan Boyar Leah Epsen Asaf Levn Asrac I s well known ha he wo smple algorhms for he classc n packng prolem, NF and WF oh have an approxmaon rao of

More information

Chapter Lagrangian Interpolation

Chapter Lagrangian Interpolation Chaper 5.4 agrangan Inerpolaon Afer readng hs chaper you should be able o:. dere agrangan mehod of nerpolaon. sole problems usng agrangan mehod of nerpolaon and. use agrangan nerpolans o fnd deraes and

More information

Advanced Macroeconomics II: Exchange economy

Advanced Macroeconomics II: Exchange economy Advanced Macroeconomcs II: Exchange economy Krzyszof Makarsk 1 Smple deermnsc dynamc model. 1.1 Inroducon Inroducon Smple deermnsc dynamc model. Defnons of equlbrum: Arrow-Debreu Sequenal Recursve Equvalence

More information

Single-loop System Reliability-Based Design & Topology Optimization (SRBDO/SRBTO): A Matrix-based System Reliability (MSR) Method

Single-loop System Reliability-Based Design & Topology Optimization (SRBDO/SRBTO): A Matrix-based System Reliability (MSR) Method 10 h US Naonal Congress on Compuaonal Mechancs Columbus, Oho 16-19, 2009 Sngle-loop Sysem Relably-Based Desgn & Topology Opmzaon (SRBDO/SRBTO): A Marx-based Sysem Relably (MSR) Mehod Tam Nguyen, Junho

More information

Notes on the stability of dynamic systems and the use of Eigen Values.

Notes on the stability of dynamic systems and the use of Eigen Values. Noes on he sabl of dnamc ssems and he use of Egen Values. Source: Macro II course noes, Dr. Davd Bessler s Tme Seres course noes, zarads (999) Ineremporal Macroeconomcs chaper 4 & Techncal ppend, and Hamlon

More information

Short-Term Load Forecasting Using PSO-Based Phase Space Neural Networks

Short-Term Load Forecasting Using PSO-Based Phase Space Neural Networks Proceedngs of he 5h WSEAS In. Conf. on SIMULATION, MODELING AND OPTIMIZATION, Corfu, Greece, Augus 7-9, 005 (pp78-83) Shor-Term Load Forecasng Usng PSO-Based Phase Space Neural Neworks Jang Chuanwen, Fang

More information

Lecture 18: The Laplace Transform (See Sections and 14.7 in Boas)

Lecture 18: The Laplace Transform (See Sections and 14.7 in Boas) Lecure 8: The Lalace Transform (See Secons 88- and 47 n Boas) Recall ha our bg-cure goal s he analyss of he dfferenal equaon, ax bx cx F, where we emloy varous exansons for he drvng funcon F deendng on

More information

Reactive Methods to Solve the Berth AllocationProblem with Stochastic Arrival and Handling Times

Reactive Methods to Solve the Berth AllocationProblem with Stochastic Arrival and Handling Times Reacve Mehods o Solve he Berh AllocaonProblem wh Sochasc Arrval and Handlng Tmes Nsh Umang* Mchel Berlare* * TRANSP-OR, Ecole Polyechnque Fédérale de Lausanne Frs Workshop on Large Scale Opmzaon November

More information

Performance Analysis for a Network having Standby Redundant Unit with Waiting in Repair

Performance Analysis for a Network having Standby Redundant Unit with Waiting in Repair TECHNI Inernaonal Journal of Compung Scence Communcaon Technologes VOL.5 NO. July 22 (ISSN 974-3375 erformance nalyss for a Nework havng Sby edundan Un wh ang n epar Jendra Sngh 2 abns orwal 2 Deparmen

More information

TSS = SST + SSE An orthogonal partition of the total SS

TSS = SST + SSE An orthogonal partition of the total SS ANOVA: Topc 4. Orhogonal conrass [ST&D p. 183] H 0 : µ 1 = µ =... = µ H 1 : The mean of a leas one reamen group s dfferen To es hs hypohess, a basc ANOVA allocaes he varaon among reamen means (SST) equally

More information

Volatility Interpolation

Volatility Interpolation Volaly Inerpolaon Prelmnary Verson March 00 Jesper Andreasen and Bran Huge Danse Mares, Copenhagen wan.daddy@danseban.com brno@danseban.com Elecronc copy avalable a: hp://ssrn.com/absrac=69497 Inro Local

More information

[ ] 2. [ ]3 + (Δx i + Δx i 1 ) / 2. Δx i-1 Δx i Δx i+1. TPG4160 Reservoir Simulation 2018 Lecture note 3. page 1 of 5

[ ] 2. [ ]3 + (Δx i + Δx i 1 ) / 2. Δx i-1 Δx i Δx i+1. TPG4160 Reservoir Simulation 2018 Lecture note 3. page 1 of 5 TPG460 Reservor Smulaon 08 page of 5 DISCRETIZATIO OF THE FOW EQUATIOS As we already have seen, fne dfference appromaons of he paral dervaves appearng n he flow equaons may be obaned from Taylor seres

More information

Bayes rule for a classification problem INF Discriminant functions for the normal density. Euclidean distance. Mahalanobis distance

Bayes rule for a classification problem INF Discriminant functions for the normal density. Euclidean distance. Mahalanobis distance INF 43 3.. Repeon Anne Solberg (anne@f.uo.no Bayes rule for a classfcaon problem Suppose we have J, =,...J classes. s he class label for a pxel, and x s he observed feaure vecor. We can use Bayes rule

More information

@FMI c Kyung Moon Sa Co.

@FMI c Kyung Moon Sa Co. Annals of Fuzzy Mahemacs and Informacs Volume 8, No. 2, (Augus 2014), pp. 245 257 ISSN: 2093 9310 (prn verson) ISSN: 2287 6235 (elecronc verson) hp://www.afm.or.kr @FMI c Kyung Moon Sa Co. hp://www.kyungmoon.com

More information

This document is downloaded from DR-NTU, Nanyang Technological University Library, Singapore.

This document is downloaded from DR-NTU, Nanyang Technological University Library, Singapore. Ths documen s downloaded from DR-NTU, Nanyang Technologcal Unversy Lbrary, Sngapore. Tle A smplfed verb machng algorhm for word paron n vsual speech processng( Acceped verson ) Auhor(s) Foo, Say We; Yong,

More information

Method of upper lower solutions for nonlinear system of fractional differential equations and applications

Method of upper lower solutions for nonlinear system of fractional differential equations and applications Malaya Journal of Maemak, Vol. 6, No. 3, 467-472, 218 hps://do.org/1.26637/mjm63/1 Mehod of upper lower soluons for nonlnear sysem of fraconal dfferenal equaons and applcaons D.B. Dhagude1 *, N.B. Jadhav2

More information

HEURISTICS FOR ONE CLASS OF MINIMAL COVERING PROBLEM IN CASE OF LOCATING UNDESIRABLE GOODS

HEURISTICS FOR ONE CLASS OF MINIMAL COVERING PROBLEM IN CASE OF LOCATING UNDESIRABLE GOODS HEURSTCS FOR ONE CLASS OF MNMAL COVERNG PROBLEM N CASE OF LOCATNG UNDESRABLE GOODS BRANKA DMTRJEVC MLORAD VDOVC The Faculy of Transpor and Traffc Engneerng, The Unversy of Belgrade, Belgrade, Serba and

More information

New M-Estimator Objective Function. in Simultaneous Equations Model. (A Comparative Study)

New M-Estimator Objective Function. in Simultaneous Equations Model. (A Comparative Study) Inernaonal Mahemacal Forum, Vol. 8, 3, no., 7 - HIKARI Ld, www.m-hkar.com hp://dx.do.org/.988/mf.3.3488 New M-Esmaor Objecve Funcon n Smulaneous Equaons Model (A Comparave Sudy) Ahmed H. Youssef Professor

More information

Particle Swarm Optimization Algorithm with Reverse-Learning and Local-Learning Behavior

Particle Swarm Optimization Algorithm with Reverse-Learning and Local-Learning Behavior 35 JOURNAL OF SOFTWARE, VOL. 9, NO. 2, FEBRUARY 214 Parcle Swarm Opmzaon Algorhm wh Reverse-Learnng and Local-Learnng Behavor Xuewen Xa Naonal Engneerng Research Cener for Saelle Posonng Sysem, Wuhan Unversy,

More information

( ) () we define the interaction representation by the unitary transformation () = ()

( ) () we define the interaction representation by the unitary transformation () = () Hgher Order Perurbaon Theory Mchael Fowler 3/7/6 The neracon Represenaon Recall ha n he frs par of hs course sequence, we dscussed he chrödnger and Hesenberg represenaons of quanum mechancs here n he chrödnger

More information

First-order piecewise-linear dynamic circuits

First-order piecewise-linear dynamic circuits Frs-order pecewse-lnear dynamc crcus. Fndng he soluon We wll sudy rs-order dynamc crcus composed o a nonlnear resse one-por, ermnaed eher by a lnear capacor or a lnear nducor (see Fg.. Nonlnear resse one-por

More information

John Geweke a and Gianni Amisano b a Departments of Economics and Statistics, University of Iowa, USA b European Central Bank, Frankfurt, Germany

John Geweke a and Gianni Amisano b a Departments of Economics and Statistics, University of Iowa, USA b European Central Bank, Frankfurt, Germany Herarchcal Markov Normal Mxure models wh Applcaons o Fnancal Asse Reurns Appendx: Proofs of Theorems and Condonal Poseror Dsrbuons John Geweke a and Gann Amsano b a Deparmens of Economcs and Sascs, Unversy

More information

UNIVERSITAT AUTÒNOMA DE BARCELONA MARCH 2017 EXAMINATION

UNIVERSITAT AUTÒNOMA DE BARCELONA MARCH 2017 EXAMINATION INTERNATIONAL TRADE T. J. KEHOE UNIVERSITAT AUTÒNOMA DE BARCELONA MARCH 27 EXAMINATION Please answer wo of he hree quesons. You can consul class noes, workng papers, and arcles whle you are workng on he

More information

Effective Task Scheduling and Dynamic Resource Optimization based on Heuristic Algorithms in Cloud Computing Environment

Effective Task Scheduling and Dynamic Resource Optimization based on Heuristic Algorithms in Cloud Computing Environment KSII TRANSACTIONS ON INTERNET AND INFORMATION SYSTEMS VOL 11, NO 12, Dec 2017 5780 Copyrgh c2017 KSII Effecve Task Schedulng and Dynamc Resource Opmzaon based on Heursc Algorhms n Cloud Compung Envronmen

More information

Graduate Macroeconomics 2 Problem set 5. - Solutions

Graduate Macroeconomics 2 Problem set 5. - Solutions Graduae Macroeconomcs 2 Problem se. - Soluons Queson 1 To answer hs queson we need he frms frs order condons and he equaon ha deermnes he number of frms n equlbrum. The frms frs order condons are: F K

More information

Comparison of Differences between Power Means 1

Comparison of Differences between Power Means 1 In. Journal of Mah. Analyss, Vol. 7, 203, no., 5-55 Comparson of Dfferences beween Power Means Chang-An Tan, Guanghua Sh and Fe Zuo College of Mahemacs and Informaon Scence Henan Normal Unversy, 453007,

More information

On computing differential transform of nonlinear non-autonomous functions and its applications

On computing differential transform of nonlinear non-autonomous functions and its applications On compung dfferenal ransform of nonlnear non-auonomous funcons and s applcaons Essam. R. El-Zahar, and Abdelhalm Ebad Deparmen of Mahemacs, Faculy of Scences and Humanes, Prnce Saam Bn Abdulazz Unversy,

More information

Mechanics Physics 151

Mechanics Physics 151 Mechancs Physcs 5 Lecure 9 Hamlonan Equaons of Moon (Chaper 8) Wha We Dd Las Tme Consruced Hamlonan formalsm H ( q, p, ) = q p L( q, q, ) H p = q H q = p H = L Equvalen o Lagrangan formalsm Smpler, bu

More information

Hongyuan Gao* and Ming Diao

Hongyuan Gao* and Ming Diao In. J. odellng, Idenfcaon and Conrol, Vol. X, No. Y, 200X Culural frework algorhm and s applcaon for dgal flers desgn Hongyuan Gao* and ng Dao College of Informaon and Communcaon Engneerng, Harbn Engneerng

More information

Solving the multi-period fixed cost transportation problem using LINGO solver

Solving the multi-period fixed cost transportation problem using LINGO solver Inernaonal Journal of Pure and Appled Mahemacs Volume 119 No. 12 2018, 2151-2157 ISSN: 1314-3395 (on-lne verson) url: hp://www.pam.eu Specal Issue pam.eu Solvng he mul-perod fxed cos ransporaon problem

More information

NATIONAL UNIVERSITY OF SINGAPORE PC5202 ADVANCED STATISTICAL MECHANICS. (Semester II: AY ) Time Allowed: 2 Hours

NATIONAL UNIVERSITY OF SINGAPORE PC5202 ADVANCED STATISTICAL MECHANICS. (Semester II: AY ) Time Allowed: 2 Hours NATONAL UNVERSTY OF SNGAPORE PC5 ADVANCED STATSTCAL MECHANCS (Semeser : AY 1-13) Tme Allowed: Hours NSTRUCTONS TO CANDDATES 1. Ths examnaon paper conans 5 quesons and comprses 4 prned pages.. Answer all

More information

Lecture 6: Learning for Control (Generalised Linear Regression)

Lecture 6: Learning for Control (Generalised Linear Regression) Lecure 6: Learnng for Conrol (Generalsed Lnear Regresson) Conens: Lnear Mehods for Regresson Leas Squares, Gauss Markov heorem Recursve Leas Squares Lecure 6: RLSC - Prof. Sehu Vjayakumar Lnear Regresson

More information

( t) Outline of program: BGC1: Survival and event history analysis Oslo, March-May Recapitulation. The additive regression model

( t) Outline of program: BGC1: Survival and event history analysis Oslo, March-May Recapitulation. The additive regression model BGC1: Survval and even hsory analyss Oslo, March-May 212 Monday May 7h and Tuesday May 8h The addve regresson model Ørnulf Borgan Deparmen of Mahemacs Unversy of Oslo Oulne of program: Recapulaon Counng

More information

Mechanics Physics 151

Mechanics Physics 151 Mechancs Physcs 5 Lecure 9 Hamlonan Equaons of Moon (Chaper 8) Wha We Dd Las Tme Consruced Hamlonan formalsm Hqp (,,) = qp Lqq (,,) H p = q H q = p H L = Equvalen o Lagrangan formalsm Smpler, bu wce as

More information

CS286.2 Lecture 14: Quantum de Finetti Theorems II

CS286.2 Lecture 14: Quantum de Finetti Theorems II CS286.2 Lecure 14: Quanum de Fne Theorems II Scrbe: Mara Okounkova 1 Saemen of he heorem Recall he las saemen of he quanum de Fne heorem from he prevous lecure. Theorem 1 Quanum de Fne). Le ρ Dens C 2

More information

II. Light is a Ray (Geometrical Optics)

II. Light is a Ray (Geometrical Optics) II Lgh s a Ray (Geomercal Opcs) IIB Reflecon and Refracon Hero s Prncple of Leas Dsance Law of Reflecon Hero of Aleandra, who lved n he 2 nd cenury BC, posulaed he followng prncple: Prncple of Leas Dsance:

More information

Including the ordinary differential of distance with time as velocity makes a system of ordinary differential equations.

Including the ordinary differential of distance with time as velocity makes a system of ordinary differential equations. Soluons o Ordnary Derenal Equaons An ordnary derenal equaon has only one ndependen varable. A sysem o ordnary derenal equaons consss o several derenal equaons each wh he same ndependen varable. An eample

More information

ON THE WEAK LIMITS OF SMOOTH MAPS FOR THE DIRICHLET ENERGY BETWEEN MANIFOLDS

ON THE WEAK LIMITS OF SMOOTH MAPS FOR THE DIRICHLET ENERGY BETWEEN MANIFOLDS ON THE WEA LIMITS OF SMOOTH MAPS FOR THE DIRICHLET ENERGY BETWEEN MANIFOLDS FENGBO HANG Absrac. We denfy all he weak sequenal lms of smooh maps n W (M N). In parcular, hs mples a necessary su cen opologcal

More information

Way nding under Uncertainty in Continuous Time and Space by Dynamic Programming

Way nding under Uncertainty in Continuous Time and Space by Dynamic Programming Way ndng under Uncerany n Connuous Tme and Space by Dynamc Programmng Dr. Serge P. Hoogendoorn (s.hoogendoorn@c.udelf.nl) Transporaon and Tra c Engneerng Secon Faculy of Cvl Engneerng and Geoscences -

More information

Introduction ( Week 1-2) Course introduction A brief introduction to molecular biology A brief introduction to sequence comparison Part I: Algorithms

Introduction ( Week 1-2) Course introduction A brief introduction to molecular biology A brief introduction to sequence comparison Part I: Algorithms Course organzaon Inroducon Wee -2) Course nroducon A bref nroducon o molecular bology A bref nroducon o sequence comparson Par I: Algorhms for Sequence Analyss Wee 3-8) Chaper -3, Models and heores» Probably

More information

e-journal Reliability: Theory& Applications No 2 (Vol.2) Vyacheslav Abramov

e-journal Reliability: Theory& Applications No 2 (Vol.2) Vyacheslav Abramov June 7 e-ournal Relably: Theory& Applcaons No (Vol. CONFIDENCE INTERVALS ASSOCIATED WITH PERFORMANCE ANALYSIS OF SYMMETRIC LARGE CLOSED CLIENT/SERVER COMPUTER NETWORKS Absrac Vyacheslav Abramov School

More information

Ordinary Differential Equations in Neuroscience with Matlab examples. Aim 1- Gain understanding of how to set up and solve ODE s

Ordinary Differential Equations in Neuroscience with Matlab examples. Aim 1- Gain understanding of how to set up and solve ODE s Ordnary Dfferenal Equaons n Neuroscence wh Malab eamples. Am - Gan undersandng of how o se up and solve ODE s Am Undersand how o se up an solve a smple eample of he Hebb rule n D Our goal a end of class

More information

APOC #232 Capacity Planning for Fault-Tolerant All-Optical Network

APOC #232 Capacity Planning for Fault-Tolerant All-Optical Network APOC #232 Capacy Plannng for Faul-Toleran All-Opcal Nework Mchael Kwok-Shng Ho and Kwok-wa Cheung Deparmen of Informaon ngneerng The Chnese Unversy of Hong Kong Shan, N.T., Hong Kong SAR, Chna -mal: kwcheung@e.cuhk.edu.hk

More information

ABSTRACT KEYWORDS. Bonus-malus systems, frequency component, severity component. 1. INTRODUCTION

ABSTRACT KEYWORDS. Bonus-malus systems, frequency component, severity component. 1. INTRODUCTION EERAIED BU-MAU YTEM ITH A FREQUECY AD A EVERITY CMET A IDIVIDUA BAI I AUTMBIE IURACE* BY RAHIM MAHMUDVAD AD HEI HAAI ABTRACT Frangos and Vronos (2001) proposed an opmal bonus-malus sysems wh a frequency

More information

PSO Algorithm Particle Filters for Improving the Performance of Lane Detection and Tracking Systems in Difficult Roads

PSO Algorithm Particle Filters for Improving the Performance of Lane Detection and Tracking Systems in Difficult Roads Sensors 2012, 12, 17168-17185; do:10.3390/s121217168 Arcle OPEN ACCESS sensors ISSN 1424-8220 www.mdp.com/journal/sensors PSO Algorhm Parcle Flers for Improvng he Performance of Lane Deecon and Trackng

More information

Chapter 6: AC Circuits

Chapter 6: AC Circuits Chaper 6: AC Crcus Chaper 6: Oulne Phasors and he AC Seady Sae AC Crcus A sable, lnear crcu operang n he seady sae wh snusodal excaon (.e., snusodal seady sae. Complee response forced response naural response.

More information

Testing a new idea to solve the P = NP problem with mathematical induction

Testing a new idea to solve the P = NP problem with mathematical induction Tesng a new dea o solve he P = NP problem wh mahemacal nducon Bacground P and NP are wo classes (ses) of languages n Compuer Scence An open problem s wheher P = NP Ths paper ess a new dea o compare he

More information

A Principled Approach to MILP Modeling

A Principled Approach to MILP Modeling A Prncpled Approach o MILP Modelng John Hooer Carnege Mellon Unvers Augus 008 Slde Proposal MILP modelng s an ar, bu need no be unprncpled. Slde Proposal MILP modelng s an ar, bu need no be unprncpled.

More information

3. OVERVIEW OF NUMERICAL METHODS

3. OVERVIEW OF NUMERICAL METHODS 3 OVERVIEW OF NUMERICAL METHODS 3 Inroducory remarks Ths chaper summarzes hose numercal echnques whose knowledge s ndspensable for he undersandng of he dfferen dscree elemen mehods: he Newon-Raphson-mehod,

More information

Lecture VI Regression

Lecture VI Regression Lecure VI Regresson (Lnear Mehods for Regresson) Conens: Lnear Mehods for Regresson Leas Squares, Gauss Markov heorem Recursve Leas Squares Lecure VI: MLSC - Dr. Sehu Vjayakumar Lnear Regresson Model M

More information

Modeling and Solving of Multi-Product Inventory Lot-Sizing with Supplier Selection under Quantity Discounts

Modeling and Solving of Multi-Product Inventory Lot-Sizing with Supplier Selection under Quantity Discounts nernaonal ournal of Appled Engneerng Research SSN 0973-4562 Volume 13, Number 10 (2018) pp. 8708-8713 Modelng and Solvng of Mul-Produc nvenory Lo-Szng wh Suppler Selecon under Quany Dscouns Naapa anchanaruangrong

More information

An Optimal Control Approach to the Multi-agent Persistent Monitoring Problem

An Optimal Control Approach to the Multi-agent Persistent Monitoring Problem An Opmal Conrol Approach o he Mul-agen Perssen Monorng Problem Chrsos.G. Cassandras, Xuchao Ln and Xu Chu Dng Dvson of Sysems Engneerng and Cener for Informaon and Sysems Engneerng Boson Unversy, cgc@bu.edu,

More information

Appendix to Online Clustering with Experts

Appendix to Online Clustering with Experts A Appendx o Onlne Cluserng wh Expers Furher dscusson of expermens. Here we furher dscuss expermenal resuls repored n he paper. Ineresngly, we observe ha OCE (and n parcular Learn- ) racks he bes exper

More information

Fall 2010 Graduate Course on Dynamic Learning

Fall 2010 Graduate Course on Dynamic Learning Fall 200 Graduae Course on Dynamc Learnng Chaper 4: Parcle Flers Sepember 27, 200 Byoung-Tak Zhang School of Compuer Scence and Engneerng & Cognve Scence and Bran Scence Programs Seoul aonal Unversy hp://b.snu.ac.kr/~bzhang/

More information

P R = P 0. The system is shown on the next figure:

P R = P 0. The system is shown on the next figure: TPG460 Reservor Smulaon 08 page of INTRODUCTION TO RESERVOIR SIMULATION Analycal and numercal soluons of smple one-dmensonal, one-phase flow equaons As an nroducon o reservor smulaon, we wll revew he smples

More information

PARTICLE SWARM OPTIMIZATION BASED ON BOTTLENECK MACHINE FOR JOBSHOP SCHEDULING

PARTICLE SWARM OPTIMIZATION BASED ON BOTTLENECK MACHINE FOR JOBSHOP SCHEDULING Proceedng 7 h Inernaonal Semnar on Indusral Engneerng and Managemen PARTICLE SWARM OPTIMIZATION BASED ON BOTTLENECK MACHINE FOR JOBSHOP SCHEDULING Rahm Mauldya Indusral Engneerng Deparmen, Indusral Engneerng

More information

Hongbin Dong Computer Science and Technology College Harbin Engineering University Harbin, China

Hongbin Dong Computer Science and Technology College Harbin Engineering University Harbin, China Hongbn Dong Compuer Scence and Harbn Engneerng Unversy Harbn, Chna donghongbn@hrbeu.edu.cn Xue Yang Compuer Scence and Harbn Engneerng Unversy Harbn, Chna yangxue@hrbeu.edu.cn Xuyang Teng Compuer Scence

More information

Mechanics Physics 151

Mechanics Physics 151 Mechancs Physcs 5 Lecure 0 Canoncal Transformaons (Chaper 9) Wha We Dd Las Tme Hamlon s Prncple n he Hamlonan formalsm Dervaon was smple δi δ Addonal end-pon consrans pq H( q, p, ) d 0 δ q ( ) δq ( ) δ

More information

Particle Swarm Procedure for the Capacitated Open Pit Mining Problem

Particle Swarm Procedure for the Capacitated Open Pit Mining Problem Parcle Swarm Procedure for he Capacaed Open P Mnng Problem Jacques A Ferland, Jorge Amaya 2, Melody Suzy Dumo Déparemen d nformaque e de recherche opéraonnelle, Unversé de Monréal, Monréal, Québec, Canada

More information

FI 3103 Quantum Physics

FI 3103 Quantum Physics /9/4 FI 33 Quanum Physcs Aleander A. Iskandar Physcs of Magnesm and Phooncs Research Grou Insu Teknolog Bandung Basc Conces n Quanum Physcs Probably and Eecaon Value Hesenberg Uncerany Prncle Wave Funcon

More information

(,,, ) (,,, ). In addition, there are three other consumers, -2, -1, and 0. Consumer -2 has the utility function

(,,, ) (,,, ). In addition, there are three other consumers, -2, -1, and 0. Consumer -2 has the utility function MACROECONOMIC THEORY T J KEHOE ECON 87 SPRING 5 PROBLEM SET # Conder an overlappng generaon economy le ha n queon 5 on problem e n whch conumer lve for perod The uly funcon of he conumer born n perod,

More information

A Deterministic Algorithm for Summarizing Asynchronous Streams over a Sliding Window

A Deterministic Algorithm for Summarizing Asynchronous Streams over a Sliding Window A Deermnsc Algorhm for Summarzng Asynchronous Sreams over a Sldng ndow Cosas Busch Rensselaer Polyechnc Insue Srkana Trhapura Iowa Sae Unversy Oulne of Talk Inroducon Algorhm Analyss Tme C Daa sream: 3

More information

Linear Response Theory: The connection between QFT and experiments

Linear Response Theory: The connection between QFT and experiments Phys540.nb 39 3 Lnear Response Theory: The connecon beween QFT and expermens 3.1. Basc conceps and deas Q: ow do we measure he conducvy of a meal? A: we frs nroduce a weak elecrc feld E, and hen measure

More information