Computers and Mathematics with Applications. Computing a fuzzy shortest path in a network with mixed fuzzy arc lengths using α-cuts

Size: px
Start display at page:

Download "Computers and Mathematics with Applications. Computing a fuzzy shortest path in a network with mixed fuzzy arc lengths using α-cuts"

Transcription

1 Computers and Mathematcs wth Applcatons 60 (00) Contents lsts avalable at ScenceDrect Computers and Mathematcs wth Applcatons journal homepage: Computng a fuzzy shortest path n a network wth mxed fuzzy arc lengths usng α-cuts Al Tajdn a, Iraj Mahdav a,, Nezam Mahdav-Amr b, Bahram Sadeghpour-Gldeh c a Department of Industral Engneerng, Mazandaran Unversty of Scence and Technology, Babol, Iran b Faculty of Mathematcal Scences, Sharf Unversty of Technology, Tehran, Iran c Department of Mathematcs and Statstcs, Mazandaran Unversty, Babolsar, Iran a r t c l e n f o a b s t r a c t Keywords: Fuzzy numbers α-cut Dstance functon Shortest path Lnear least squares Regresson Dynamc programmng We are concerned wth the desgn of a model and an algorthm for computng a shortest path n a network havng varous types of fuzzy arc lengths. Frst, we develop a new technque for the addton of varous fuzzy numbers n a path usng α-cuts by proposng a lnear least squares model to obtan membershp functons for the consdered addtons. Then, usng a recently proposed dstance functon for comparson of fuzzy numbers, we present a dynamc programmng method for fndng a shortest path n the network. Examples are worked out to llustrate the applcablty of the proposed model. 00 Elsever Ltd. All rghts reserved.. Introducton The problem of fndng a shortest path from a specfed source node to any other node s fundamental n graph theory, and s of contnung nterest [,]. Ths problem arses from many applcatons ncludng transportaton, routng, communcatons, supply chan management or models nvolvng agents. Let G = (V, E) be a graph, where V s the set of vertces (nodes) and E s the set of edges (arcs). A path between two nodes s an alternatng sequence of vertces and edges begnnng wth a startng node and endng wth an endng node. The dstance (cost) of a path s the sum of the weghts (arc lengths) of the edges on the path. However, snce there can be more than one path between any two vertces, the problem of fndng a path wth a mnmal cost between two specfed vertces of nterest s the so-called shortest path problem (SPP). Although n conventonal graph theory, the weghts of the edges n an SPP are assumed to be precse real numbers, for most practcal applcatons, these parameters (.e., costs, capactes, demands, tme, etc.) are naturally mprecse. In such cases, an approprate modelng approach may justfably make use of fuzzy numbers, and so does the name fuzzy shortest path problem (FSPP) appear n the lterature [,3,4]. The FSPP, nvolvng addton and comparson of fuzzy numbers, s qute dfferent from the conventonal SPP, whch only nvolves crsp numbers. In an FSPP, the costs beng fuzzy numbers, the task of fndng a path beng smaller than all the others s not straghtforward, as the comparson of fuzzy numbers as an operaton can be defned n a wde varety of ways. Recently, several results have been publshed on the FSPP [5]. The work of Dubos and Prade [6] s one of the frst on ths subject and consders extensons of the classcal Floyd and Ford Moore Bellman (FMB) algorthms. However, t was verfed that the algorthm would compute the shortest path dstance wthout dentfyng an exstng path (see [7]) as was outlned by Klen [8] wth a fuzzy domnance set. Ln and Chern [3] defned the denomnaton of vtal arcs as beng those whose removal from the path resulted n an ncrease of the cost. Another algorthm for ths problem, presented by Okada and Gen [9,0], s a generalzaton of the Djkstra algorthm. In ths algorthm, the weghts of the arcs are consdered to be nterval numbers, and are defned usng a partal order between nterval numbers. Okada and Soper [5] characterzed a Correspondng author. Tel.: ; fax: E-mal address: rajarash@redffmal.com (I. Mahdav) /$ see front matter 00 Elsever Ltd. All rghts reserved. do:0.06/j.camwa

2 990 A. Tajdn et al. / Computers and Mathematcs wth Applcatons 60 (00) soluton not as the shortest path, but as a fuzzy set soluton, where each element of the set s a no domnated path or a Pareto optmal path wth fuzzy edge weghts. However, ths algorthm does not provde decson-makers wth any gudelnes for choosng a best path accordng to ther own vewponts (optmstc/pessmstc, rsky/conservatve) []. Blue et al. [] presented an algorthm whch would fnd a cut value to lmt the number of analyzed paths, and then appled a modfed verson of the k-shortest path (crsp) algorthm proposed by Eppsten [7]. Followng the dea of fndng a fuzzy set soluton, Okada [3] ntroduced the concept of the degree of possblty of an arc beng on a shortest path. Among the most recent work s the one by Nayeem and Pal [4] that proposes an algorthm based on the acceptance ndex ntroduced by Sengupta and Pal [4] and whch gves a sngle fuzzy shortest path or a gudelne for choosng a best fuzzy shortest path accordng to the decson-maker vewpont [5]. Here, we propose a new approach and an algorthm to fnd a shortest path n a network wth varous fuzzy arc lengths. The remander of the paper s organzed as follows. In Secton, basc concepts and defntons are gven. Secton 3 explans ways of computng α-cuts for fuzzy numbers. We present our fuzzy sum operator by use of a lnear least squares model n Secton 4. In Secton 5, usng a recently proposed dstance functon, we present a dynamc programmng algorthm for fndng fuzzy shortest path n a mxed fuzzy network. We conclude n Secton 6.. Defntons We start wth basc defntons of some well-known fuzzy numbers. Defnton. An LR fuzzy number s represented by ã = (m, a, b) LR, wth the membershp functon, µã(x), defned by [ ] m x L x m a µã(x) = [ ] x m R x m, b where L and R are non-ncreasng functons from R + to [0, ], L(0) = R(0) =, m s the center, a s the left spread and b s the rght spread. Note that f L(x) = R(x) = x wth 0 < x <, then x s a trangular fuzzy number and s represented by the trplet ã = (m, a, b), wth the membershp functon, µã(x), defned by ( ) m x x m a µã(x) = ( ) x m x m. b Defnton. A trapezodal fuzzy number ã s shown by ã = (a, a, a 3, a 4 ), wth the membershp functon as follows: 0 x a x a a x a a a µã(x) = a x a 3 a 4 x a 3 x a 4 a 4 a 3 0 a 4 x. A general trapezodal fuzzy number, along wth a cut (to be explaned later n Defnton 4), s shown n Fg.. It s apparent that a trangular fuzzy number s a specal trapezodal fuzzy number wth a = a 3. Defnton 3. If L(x) = R(x) = e x, wth x R, then x s a normal fuzzy number that s shown by (m, σ ) and the correspondng membershp functon s defned to be: x m µã(x) = e ( σ ), x R, where m s the mean and σ s the standard devaton. A normal fuzzy number, along wth a cut (to be explaned later n Defnton 4), s shown n Fg.. Defnton 4. The α-cut and strong α-cut for a fuzzy number ã are shown by ã α and ã + α, respectvely, and for α [0, ] are defned to be: ã α = { x µã(x) α, x X }, ã + α = { x µã(x) > α, x X }, where X s the unversal set.

3 A. Tajdn et al. / Computers and Mathematcs wth Applcatons 60 (00) Fg.. A trapezodal fuzzy number wth an α-cut. Fg.. A normal fuzzy number wth an α-cut. Note that the upper and lower bounds for the α-cut set (ã α ) are shown by sup ã α and nf ã α, respectvely. Here, we assume that the upper and lower bounds of α-cuts are fnte values and for smplcty we show sup ã α by ã R α and nf ã α by ã L α (see Fgs. and ). 3. Computng α-cuts for fuzzy numbers For the LR fuzzy numbers wth L and R nvertble functons, the α-cuts are: [ ] m x α = L m x = L (α) ã L α a a = x = m al (α), [ ] x m α = R x m = R (α) ã R α b b = x = m + br (α). For specfc L and R functons, the followng cases are dscussed. 3.. α-cuts for trapezodal fuzzy numbers Let ã = (a, a, a 3, a 4 ) be a trapezodal fuzzy number. An α-cut for ã, ã α, s computed as: α = x a a a ã L α = x = (a a )α + a, α = a 4 x a 4 a 3 ã R α = x = a 4 (a 4 a 3 )α, ã α = {ãl α = (a a )α + a ã R α = a 4 (a 4 a 3 )α, 0 α, () where ã = α [ãl α, ãr α ] s the correspondng α-cut. The α-cuts for trangular fuzzy numbers are obtaned by usng the above equatons consderng a = a α-cuts for normal fuzzy numbers If ã = (m, σ ) s a normal fuzzy number, then ã α s computed as: m x α = e ( σ ) ln(α) = m x ã L α σ = x = m σ ln α

4 99 A. Tajdn et al. / Computers and Mathematcs wth Applcatons 60 (00) x m α = e ( σ ) ln(α) = x m ã R α σ = x = m + σ ln α {ãl ã α = α = m σ ln α ã R α = m + σ, 0 < α. () ln α 4. Fuzzy approxmate sum operators Here, we propose an approach for summng varous fuzzy numbers approxmately usng α-cuts. The approxmaton s based on fttng an approprate model for the sum usng α-cuts of the addton as the ftness data. Let us dvde the α-nterval [0, ] nto n equal subntervals by lettng α 0 = 0, α = α + α, α =, =,..., n. Ths way, we have a set of n + n equdstant ponts. For the normal fuzzy numbers x (, + ), t s mproper to assume α beng equal to zero. Therefore, n ths case we consder α (0, ], and thus use the nonzero α, n. Here, we ntend to show how to add up a trapezodal fuzzy number wth a normal one. We present a numercal approach to approxmate the sum and ts correspondng membershp functon. 4.. α-cut sum Let ã = (a, a, a 3, a 4 ) and b = (m, σ ) be the trapezodal and normal fuzzy numbers, respectvely. Gven α (0, ], n, the α-cut sum of these fuzzy numbers usng Eqs. () and () s obtaned as follows: L { c c α = ã α + bα c α = α = (ã L α + bl α ) c α R = (ãr α +, n, br α ), where, [ c α L = (a a )α + a + m σ ] ln α, [ c α R = a 4 (a 4 a 3 )α + m + σ ] ln α. (3) Usng Eq. (3), correspondng to α, n, n ponts are obtaned for c (n ponts for the c α L and n ponts for the c α R ). Usng these ponts, t s possble to approxmate the sum of the two fuzzy numbers. An approxmate membershp functon of the sum s computed by fttng an approprate functon usng the α-cut ponts. For the addton of normal and trapezodal fuzzy numbers, the case beng consdered n our examples later on, we propose an exponental membershp functon for approxmatng the sum as follows (later, we wll see that ths choce would ndeed provde a good model for the approxmatng sum of trapezodal and normal fuzzy numbers). Let x = c α R and y = µ( c α R ), and usng the n ponts ( ) (x, y ), n, consder the fttng model as y = e x λ β. The unknown parameters λ and β appear nonlnearly. We lnearze the model, by notng that for any x > λ (as s the case here for the rght hand model), we must have: ( ) x λ ln y =. (4) β Snce 0 < y, then ln(y ) 0, and thus we can wrte, ( ) x λ ln y =, β and hence, β ln y + λ = x. Therefore, we defne a lnear of least squares model for the mnmzaton of error as follows: n ( mn E = β ) ln y + λ x, = where x = c R α s gven by (3). To solve (7), we need to have: E β = E λ = [( ) ln y (β )] ln y + λ x = 0, (8) [ (β )] ln y + λ x = 0. (5) (6) (7)

5 A. Tajdn et al. / Computers and Mathematcs wth Applcatons 60 (00) Hence, we need to solve the followng so-called normal equatons for the unknown parameters β and λ: x ln y, (9) β β ln y + λ ln y + nλ = ln y = x. (0) Usng Cramer s rule to solve (9) and (0), β and λ are explctly determned to be: x ln y ln y x n n x ln y ln y x β = β = ln y ln y n ln y ln y, () ln y ln y n ( ) ln y ln y x ln y x ln y x ( ) x ln y ln y λ = λ = ln y ln y n ln y ln y. () ln y ln y n ( Now, smlarly let x = c α L and ȳ = µ( c α L ), and consder the model y = e λ ) x β. Usng the above approach, we have: Or, ( λ ) x ln y =. β β ln y + λ = x. (3) Thus, solvng mn E = n = smlarly leads to: ( β ln y + λ x ), ln y β = n n x ln y + ln y x ln y ln y, (4) ( x ) ln y ln y λ = ln y n x ln y ln y ln y. (5) Thus, the membershp functon s determned to be: ) β x < λ µ c (x) = λ x λ ( ) e x λ β x > λ, e ( λ x (6) wth λ, β, λ and β as defned by (), (), (4) and (5), respectvely.

6 994 A. Tajdn et al. / Computers and Mathematcs wth Applcatons 60 (00) Normal Fuzzy Number Trapezod Fuzzy Number x 8 Fg. 3a. The α-cut ponts of fuzzy numbers. Remarks. An equvalent approach for fttng the least squares model s to consder the bass functons { ln y, }, for the fttng relaton (6). Ths would lead to the mnmzaton model, wth mn θ Aθ x, ln y A =.., ln yn θ = [ β λ ], x = x.. x n. The soluton of (7) s obtaned by solvng the normal equatons, A T Aθ = A T x, yeldng the same results as () and (). Convenently, ths general data fttng approach can be used to consder other membershp functons n cases usng varous types of fuzzy numbers. The key element s thus to decde the approprate basc functons. Now, we provde a numercal llustraton for our proposed model for addton of trapezodal and normal fuzzy numbers. Example. Consder the followng two fuzzy numbers, one beng normal and the other beng trapezodal: ã = (3, 5), b = (4, 0, 7, 6). The dagrams of the numbers at the α-cuts and ther sum, (ã + b), usng Eq. (3), are shown n Fgs. 3a and 3b, respectvely. As ndcated n the dagram, the result s nether trapezodal nor normal. Usng the ponts obtaned from the α-cuts consderng n = 00, the values of (λ, β) and (λ, β ) are obtaned by (9) (5). For ths example, we obtan: ( ) e 3 x x < 3 µ c (x) = 3 x 30 ( ) e x x > 30. (7) 5. An algorthm for fuzzy shortest path n a network 5.. Dstance between fuzzy numbers Knowng that we can obtan a good approxmaton for the addton of varous fuzzy numbers by use of α-cuts, we compute the dstance between two fuzzy numbers usng the resultng ponts from the α-cuts. Assume that ã and b are two fuzzy

7 A. Tajdn et al. / Computers and Mathematcs wth Applcatons 60 (00) Trapezod + Normal < α Output x 50 Fg. 3b. The sum of mxed fuzzy numbers usng α-cuts. numbers. We apply a fuzzy rankng method for fuzzy numbers. We have used ths rankng method effectvely n a recent work [6]. Let us consder fuzzy mn operatons as follows: Mn value (ã, b) = (mn(a, b ), mn(a, b ), mn(a 3, b 3 ), mn(a 4, b 4 )). (8) It s evdent that, for non-comparable fuzzy numbers ã and b, the fuzzy mn operaton results n a fuzzy number dfferent from both of them. For example, for ã = (5, 0, 3, 9) and b = (6, 9, 5, 0), we get from (8) a fuzzy MṼ = Mn value(ã, b) = (5, 9, 3, 9), whch s dfferent from both ã and b. To allevate ths drawback, we use a method based on the dstance between fuzzy numbers. We use the dstance functon ntroduced n [7]. The man advantages of ths dstance functon are the generalty of ts usage on varous fuzzy numbers, and ts relablty n dstngushng unequal fuzzy numbers. Indeed, the usage of the dstance functon below worked out to be qute approprate for our approach. The D p,q -dstance, ndexed by parameters < p < and 0 < q <, between two fuzzy numbers ã and b s a nonnegatve functon gven by: [ ] ( q) a α b p α dα + q a + α b + p p α dα, p <, D p,q (ã, b) = 0 0 (9) ( q) sup a α b α + q nf a + α b + α, p =. 0<α 0<α The analytcal propertes of D p,q depend on the frst parameter p, whle the second parameter q of D p,q characterzes the subjectve weght attrbuted to the end ponts of the support;.e., the a + α and a α of the fuzzy numbers. If there s no reason for dstngushng any sde of the fuzzy numbers, then D p, s recommended. Havng q close to results n consderng the rght sde of the support of the fuzzy numbers more favorably. Snce the sgnfcance of the end ponts of the support of the fuzzy numbers s assumed to be the same, then we consder q =. For two fuzzy numbers ã and b wth correspondng α -cuts, the D p,q dstance s approxmately proportonal to: D p,q (ã, b) = [( q) n a α b p α + q = n a + α b + p α = ] p. (0) If q =, p =, then the above equaton turns nto: D, (ã, b) = n (a α b α ) + n (a + α b + α ). () = =

8 996 A. Tajdn et al. / Computers and Mathematcs wth Applcatons 60 (00) To compare two fuzzy arc lengths ã and b wth α -cuts as ther approxmatons, snce they are supposed to represent postve values, we compare them wth MṼ = (0, 0,..., 0). In fact, we use formula () to compute D, (ã, MṼ) and D, ( b, MṼ) and then use these values for comparson of the two numbers. Here, we consder ã = (6, 9, 5, 0) and b = (5, 0, 3, 9) wth n = 0. Then, the α-cuts for ã are obtaned to be a + α = {9.5, 9, 8.5, 8, 7.5, 7, 6.5, 6, 5.5, 5}. a α = {9, 8.7, 8.4, 8., 7.8, 7.5, 7., 6.9, 6.6, 6.3} and for b we have b + α = {8.4, 7.8, 7., 6.6, 6, 5.4, 4.8, 4., 3.6, 3}. b α = {0, 9.5, 9, 8.5, 8, 7.5, 7, 6.5, 6, 5.5}. As a result, the dstances are D, (ã, MṼ) = D, ( b, MṼ) = Therefore, ã > b. 5.. An algorthm for computng a shortest path The followng dynamc programmng algorthm s for computng the shortest path n a network. The algorthm s based on Floyd s dynamc programmng method to fnd a shortest path, f t exsts, between every par of nodes and j n the network [8]. We make use of the followng optmal value functon f k (, j) and the correspondng labelng functon P k (, j): f k (, j) = length of the shortest path from node to node j when the path s consdered to use only the nodes from the set of nodes {,..., k}, P k (, j) = the last ntermedate node on the shortest path from node to node j usng {,..., k} as ntermedate node, where, s a source node, j s the end node and k refers to an ntermedate node. The dynamc updatng for the optmal path length and ts correspondng labelng are: f k (, j) = mn {f k (, j), f k (, k) + f k (k, j)}, { Pk (, j) f k s not on shortest path from to j usng {,..., k} P k (, j) = P k (k, j) otherwse. We are now ready to gve the steps of the algorthm. Algorthm. A dynamc programmng method for computng a shortest path n a fuzzy network G = (V, A), where V s the set of nodes wth V = N, and A s the set of arcs. The value dj s the fuzzy arc dstance for arc (, j), f t exsts. Below, fk (, j) the shortest path length s set to, when there s no arc. Step : Let k = 0 and fk (, j) = dj, for all (, j) A, fk = (, j) =, for all (, j) A. If an arc exsts from node to node j then let P k (, j) =. Step : Let k = k +. Do the followng steps for =,, 3,..., N, j =,, 3,..., N, j. ]. Compute the value of fk (, j) = mn [ fk (, j), fk (, k) + fk (k, j) (for the addton, our proposed method dscussed n Secton 4. and for comparson of fuzzy numbers the D p,q dstance functon () of Secton 5. are appled).. If node k s not on the shortest path usng the nodes {,,..., k} as ntermedate nodes, then let P k (, j) = P k (, j) else let P k (, j) = P k (k, j). Step 3: If k < N then go to Step. Step 4: Obtan the shortest path usng the P k (, j) values. If f N (, j) =, then there s no path between and j. The shortest path from node to j, f t exsts, s dentfed backwards and read by the nodes: j, P N (, j) = k followed by P N (, k),..., P N (, l) =, where l s the node mmedately after n the path.

9 A. Tajdn et al. / Computers and Mathematcs wth Applcatons 60 (00) (,3,4,5) (4,) (5,4) 4 (4,8,,6) 3 (5,) Fg. 4. A small szed network havng mxed fuzzy arc lengths. Table The f0 (, j) matrx for k = 0. /j 3 4 (, 3, 4, 5) (4, 8,, 6) (4, ) (5, 4) 3 (5, ) Table The P 0 (, j) matrx for k = 0. /j Table 3 The f (, j) matrx for k =. /j 3 4 (, 3, 4, 5) (4, 8,, 6) (4, ) (5, 4) 3 (5, ) Table 4 The P (, j) matrx for k =. /j Termnaton and complexty of the algorthm The proposed algorthm termnates after N outer teratons correspondng to k. A total of N(N ) addtons and comparsons are needed for every k. For each addton, n fuzzy addtons for the α -cuts should be performed resultng n n(n)(n ) addtons. For comparsons, we have (n + )N(N ) addtons and (n + ) N(N ) multplcatons usng (7). Therefore, the total needed operatons are (6n + ) N(N ) addtons and multplcatons, wth N(N ) comparsons. Example. Consder the mxed fuzzy network n Fg. 4 wth four nodes and fve arcs havng two trapezodal and three normal arc lengths as specfed n Table. Step : We set fk (, j) = dj, for k = 0, as specfed n Table. Therefore, wth P k (, j) =, Table s obtaned. Step : Here, we consder k = and compute the value of f k (, j) = mn [f k (, j), f k (, k) + f k (k, j)]. The result s shown n Table 3. Therefore, for P k (, j) =, Table 4 s obtaned. If node k s not on the shortest path usng {,,..., k} as ntermedate nodes, then we consder P k (, j) = P k (, j), otherwse we let P k (, j) = P k (, k). We now report the results obtaned for other values of k n Tables 5 0. Note that, the

10 998 A. Tajdn et al. / Computers and Mathematcs wth Applcatons 60 (00) Table 5 The f (, j) matrx for k =. /j 3 4 (, 3, 4, 5) V W (4, ) (5, 4) 3 (5, ) Table 6 The P (, j) matrx for k =. /j Table 7 The f3 (, j) matrx for k = 3. /j 3 4 (, 3, 4, 5) V W (4, ) (9, ) 3 (5, ) Table 8 The P 3 (, j) matrx for k = 3. /j Table 9 The f4 (, j) matrx for k = 4. /j 3 4 (, 3, 4, 5) V 3 W 3 (4, ) (9, ) 3 (5, ) Table 0 The P 4 (, j) matrx for k = 4. /j sets V and W are the ponts obtaned by α-cut addtons, where the V and W values are obtaned by the α -cuts consderng n = 0. It ncludes 0 ponts for the a α and 0 ponts for the a + α : V = {(4.5857, 0.474), (4.9336, ), (5.074, ), (5.4477, ), ( , ), (5.8858, 9.47), (6.078, 8.897), (6.376, ), (6.5754, ), (7, 8)} W = {(.0303, ), (.55, ), (.9, 4.089), (3.57, 3.489), (4.698,.830), (4.74,.589), (5.3,.6889), (5.905,.0895), (6.606, ), (8, 9)} V = {(4.5857, 0.474), (4.9336, ), (5.074, ), (5.4477, ), ( , ), (5.8858, 9.47), (6.078, 8.897), (6.376, ), (6.5754, ), (7, 8)} W = {(8.0655, ), (8.6673, ), (9.0549, ), ( , 5.545), ( , 5.65), (0.706, 4.894), (0.5056, ), (0.855, 4.448), (.508, 3.749), (, 3)}

11 A. Tajdn et al. / Computers and Mathematcs wth Applcatons 60 (00) x Fg. 5. The membershp functon. Table The analyss of parameter varaton. n Path β β λ λ Table The parameters of non-optmal path wth ther dstance values for Example. n Path β β λ λ D, V 3 = {(4.5857, 0.474), (4.9336, ), (5.074, ), (5.4477, ), ( , ), (5.8858, 9.47), (6.078, 8.897), (6.376, ), (6.5754, ), (7, 8)} W 3 = {(8.0655, ), (8.6673, ), (9.0549, ), ( , 5.545), ( , 5.65), (0.706, 4.894), (0.5056, ), (0.855, 4.448), (.508, 3.749), (, 3)}. Fnally, when k = N, we dentfy the shortest path as follows: Shortest path from to 4: 3 4. Shortest path length from to 4: (8.0655, ), (8.6673, ), (9.0549, ), ( , 5.545), ( , 5.65), (0.706, 4.894), (0.5056, ), (0.855, 4.448), (.508, 3.749), (, 3). Here, we obtan the membershp functon as shown n Fg. 5. To nvestgate the varaton of β, β, λ, and λ wth respect to n, for the membershp functon (6), we solve the least squares problem for dfferent values of n and obtan the left and rght membershp functon parameters. The results are reported n Table. Note that, the path does not change for dfferent szes of n and after n = 00, the varatons n membershp functon parameters are neglgble. Also, for ths optmal path D, = Moreover, the costs of other paths wth n = 00 are obtaned as shown n Table. Next, we consder a problem wth a larger number of nodes to specfy the general nput and output structures usng the proposed algorthm. Example 3. We consder a larger network as shown n Fg. 6. We consder the network havng mxed arc lengths (a combnaton of normal and trapezodal fuzzy numbers) and use our dynamc algorthm to fnd the shortest paths. The arc lengths are specfed n Table 3. Usng the dstance functon D p,q (for q = / and p = ), the shortest path from the source node to the destnaton node 3 s determned to be: To fnd an optmal path, Table 4 s used. For nstance, to fnd an optmal path from node to node 3, accordng to P 3 (, j) matrx, we have P 3 (, 3) and P 3 (, 3) = 8. From P 3 (, 3) = 8, the path 8 3 s obtaned.

12 000 A. Tajdn et al. / Computers and Mathematcs wth Applcatons 60 (00) Fg. 6. A network. Table 3 The arc lengths. Arc Fuzzy number Arc Fuzzy number Arc Fuzzy number (, ) (, 3, 5, 7) (, 3) (40, ) (, 4) (8, 0,, 3) (, 5) (7, 8, 9, 0) (, 6) (35, 0) (, 7) (6,,, 3) (3, 8) (40, ) (4, 7) (7, 0,, 4) (4, ) (6, 0, 3, 4) (5, 8) (9, 9) (5, ) (7, 0, 3, 4) (5, ) (0, 3, 5, 7) (6, 9) (6, 8, 0, ) (6, 0) (35, ) (7, 0) (9, 0,, 3) (7, ) (6, 7, 8, 9) (8, ) (5, 8, 9, 0) (8, 3) (0, 5) (9, 6) (6, 7, 9, 0) (0, 6) (40, 3) (0, 7) (5, 9, 0, ) (, 4) (8, 9,, 3) (, 7) (8, 9) (, 4) (3, 4, 6, 8) (, 5) (, 4, 5, 6) (3, 5) (37, ) (3, 9) (7, 8, 9, 0) (4, ) (,, 3, 4) (5, 8) (8, 9,, 3) (5, 9) (5, 7) (6, 0) (38, ) (7, 0) (7, 0,, ) (7, ) (6, 7, 8, 0) (8, ) (5, 7, 8, 9) (8, ) (6, 5) (8, 3) (5,5) (9, ) (5, 6, 7, 9) (0, 3) (3, 4, 6, 7) (, 3) (, 5, 7, 8) (, 3) (0, 5) Table 4 P 3 (, j) matrx. k/j Also, P 3 (, 8) = 5, P 3 (, 5) =, and P 3 (, ) = 5. Snce P 3 (, 5) =, then we Stop. Therefore, the shortest path s: Usng the least squares model, we regress the estmated ponts and obtan the followng equaton as membershp functon for the shortest path. By addton of varous fuzzy numbers on the correspondng path, the membershp functon s obtaned

13 A. Tajdn et al. / Computers and Mathematcs wth Applcatons 60 (00) α-cut value the regressed model x Fg. 7. The obtaned membershp functon usng the least squares model. Table 5 The results of parameter varaton. n Path β β λ λ Fg. 8. The network. Table 6 The arc lengths of network (both crsp and fuzzy). Arc Fuzzy lengths Crsp lengths Arc Fuzzy lengths Crsp lengths Arc Fuzzy lengths Crsp lengths (, ) (,, 7) (, 3) (, 3, 4) 3 (, 3) (3, 5, 8) 5 (, 4) (,, 6) (, 5) (, 5, 6) 5 (3, 4) (, 3, 4) 3 (4, 5) (, 4, 5) 4 (4, 6) (,, 5) (5, 6) (4, 5, 8) 5 by () (4) as follows: ( ) e x x < 59 µ c (x) = 59 < x < 65 ( ) e x x > 65. The regressed membershp functon s presented n Fg. 7. To analyze the varatons of β, β, λ, and λ wth respect to n, we use dfferent szes of n and obtan the left and rght membershp functon parameters. The results are reported n Table 5. Note that the path does not change for dfferent values of n bgger than 00, and the varatons of membershp functon parameters are neglgble. 6. Dscusson Our approach can easly be used when the dstances are crsp values. Ths can serve as a tool for comparng crsp versus fuzzy arc length for a partcular case. Consder the example shown by Fg. 8. Both the fuzzy and crsp dstances of arcs are gven n Table 6. As an example, we consdered the crsp value of each arc to be the correspondng maxmum membershp value. The obtaned results for the shortest paths correspondng to the

14 00 A. Tajdn et al. / Computers and Mathematcs wth Applcatons 60 (00) two cases show the dfference of the optmal paths. For the fuzzy case, the optmal path s 3 4 6, whle for the crsp case, the optmal path s 4 6. The dstance value of optmal path n the fuzzy network, usng our proposed algorthm s 5.93, whle for the crsp stuaton s 6 (for the path 4 6). Moreover, for the optmal path n the fuzzy network, we wll obtan the value of 8 for the crsp case, showng the dfference of =.07 n favor of the fuzzy case. Ths example shows the possble effectveness of consderng fuzzy arc lengths n a network. 7. Conclusons A novel practcal approach was proposed for computng a shortest path n a fuzzy network havng mxed fuzzy arc lengths. In dong ths, an α-cut method was presented to compute the addton of varous fuzzy numbers as arc lengths. To obtan an approxmaton of the correspondng membershp functon for the addton, we proposed a lnear least squares model. Fnally, usng a recently proposed dstance functon, we showed how to decde dstances for comparson of fuzzy arc lengths to be used n our proposed dynamc programmng algorthm for fndng an optmal (shortest) path. The effectveness of our approach was shown by workng out llustratng examples. The proposed model, whle beng practcally smple, has the flexblty to consder a mxture of varous types of fuzzy arc lengths n a general network. We also gave a comparatve case of fuzzy and crsp fuzzy cases to pont out a possble effectveness of a fuzzy network. References [] T.-N. Chuang, J.-Y. Kung, The fuzzy shortest path length and the correspondng shortest path n a network, Comput. Oper. Res. 3 (005) [] M.T. Takahash, Contrbuções ao estudo de grafos fuzzy: teora e algortmos, Ph.D. Thess, Faculdade de Engenhara Elétrca e de Computação, UNICAMP, 004. [3] C. Ln, M.S. Chern, The fuzzy shortest path problem and ts most vtal arcs, Fuzzy Sets and Systems 58 (993) [4] S.M.A. Nayeem, M. Pal, Shortest path problem on a network wth mprecse edge weght, Fuzzy Optm. Decs. Mak. 4 (005) [5] S. Okada, T. Soper, A shortest path problem on a network wth fuzzy arc lengths, Fuzzy Sets and Systems 09 (000) [6] D. Dubos, H. Prade, Fuzzy Sets and Systems: Theory and Applcatons, Academc Press, New York, 980. [7] D. Eppsten, Fndng the k-shortest paths, n: Proc. IEEE Symp. on Foundatons of Computer Scence, 994, pp [8] C.M. Klen, Fuzzy shortest paths, Fuzzy Sets and Systems 39 (99) 7 4. [9] S. Okada, M. Gen, Order relaton between ntervals and ts applcatons to shortest path problem, n: Proc. 5th Annu. Conf. on Computers and Industral Engneerng, vol. 5, 993, pp [0] S. Okada, M. Gen, Fuzzy shortest path problem, n: Proc. 6th Ann. Conf. on Computers and Industral Engneerng, vol. 7, 994, pp [] D. Dubos, H. Prade, Rankng fuzzy numbers n the settng of possblty theory, Inform. Sc. 30 (983) [] M. Blue, B. Bush, J. Puckett, Unfed approach to fuzzy graph problems, Fuzzy Sets and Systems 5 (00) [3] S. Okada, Fuzzy shortest path problems ncorporatng nteractvty among paths, Fuzzy Sets and Systems 4 (3) (004) [4] A. Sengupta, T.K. Pal, On comparng nterval numbers, European J. Oper. Res. 7 (000) [5] J.A. Moreno, J.M. Moreno, J.L. Verdegay, Fuzzy locaton problems on networks, Fuzzy Sets and Systems 4 (004) [6] I. Mahdav, R. Nourfar, A. Hedarzade, N. Mahdav Amr, A dynamc programmng approach for fndng shortest chans n a fuzzy network, Appl. Soft Comput. 9 (009) [7] B. Sadeghpour Gldeh, D. Gen, La dstance-dp, q et le coffcent de corrélaton entre deux varables aléatores floues, Actes de LFA 00, Monse, Belgum, 00, pp [8] R.W. Floyd, Algorthm 97, shortest path, Commun. ACM 5 (96) 345.

Complement of Type-2 Fuzzy Shortest Path Using Possibility Measure

Complement of Type-2 Fuzzy Shortest Path Using Possibility Measure Intern. J. Fuzzy Mathematcal rchve Vol. 5, No., 04, 9-7 ISSN: 30 34 (P, 30 350 (onlne Publshed on 5 November 04 www.researchmathsc.org Internatonal Journal of Complement of Type- Fuzzy Shortest Path Usng

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 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

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

NUMERICAL DIFFERENTIATION

NUMERICAL DIFFERENTIATION NUMERICAL DIFFERENTIATION 1 Introducton Dfferentaton s a method to compute the rate at whch a dependent output y changes wth respect to the change n the ndependent nput x. Ths rate of change s called the

More information

Errors for Linear Systems

Errors for Linear Systems Errors for Lnear Systems When we solve a lnear system Ax b we often do not know A and b exactly, but have only approxmatons  and ˆb avalable. Then the best thng we can do s to solve ˆx ˆb exactly whch

More information

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

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

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

Foundations of Arithmetic

Foundations of Arithmetic Foundatons of Arthmetc Notaton We shall denote the sum and product of numbers n the usual notaton as a 2 + a 2 + a 3 + + a = a, a 1 a 2 a 3 a = a The notaton a b means a dvdes b,.e. ac = b where c s an

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

Formulas for the Determinant

Formulas for the Determinant page 224 224 CHAPTER 3 Determnants e t te t e 2t 38 A = e t 2te t e 2t e t te t 2e 2t 39 If 123 A = 345, 456 compute the matrx product A adj(a) What can you conclude about det(a)? For Problems 40 43, use

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

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

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 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

FUZZY GOAL PROGRAMMING VS ORDINARY FUZZY PROGRAMMING APPROACH FOR MULTI OBJECTIVE PROGRAMMING PROBLEM

FUZZY GOAL PROGRAMMING VS ORDINARY FUZZY PROGRAMMING APPROACH FOR MULTI OBJECTIVE PROGRAMMING PROBLEM Internatonal Conference on Ceramcs, Bkaner, Inda Internatonal Journal of Modern Physcs: Conference Seres Vol. 22 (2013) 757 761 World Scentfc Publshng Company DOI: 10.1142/S2010194513010982 FUZZY GOAL

More information

Speeding up Computation of Scalar Multiplication in Elliptic Curve Cryptosystem

Speeding up Computation of Scalar Multiplication in Elliptic Curve Cryptosystem H.K. Pathak et. al. / (IJCSE) Internatonal Journal on Computer Scence and Engneerng Speedng up Computaton of Scalar Multplcaton n Ellptc Curve Cryptosystem H. K. Pathak Manju Sangh S.o.S n Computer scence

More information

Determine the Optimal Order Quantity in Multi-items&s EOQ Model with Backorder

Determine the Optimal Order Quantity in Multi-items&s EOQ Model with Backorder Australan Journal of Basc and Appled Scences, 5(7): 863-873, 0 ISSN 99-878 Determne the Optmal Order Quantty n Mult-tems&s EOQ Model wth Backorder Babak Khabr, Had Nasser, 3 Ehsan Ehsan and Nma Kazem Department

More information

Neryškioji dichotominių testo klausimų ir socialinių rodiklių diferencijavimo savybių klasifikacija

Neryškioji dichotominių testo klausimų ir socialinių rodiklių diferencijavimo savybių klasifikacija Neryškoj dchotomnų testo klausmų r socalnų rodklų dferencjavmo savybų klasfkacja Aleksandras KRYLOVAS, Natalja KOSAREVA, Julja KARALIŪNAITĖ Technologcal and Economc Development of Economy Receved 9 May

More information

One-sided finite-difference approximations suitable for use with Richardson extrapolation

One-sided finite-difference approximations suitable for use with Richardson extrapolation Journal of Computatonal Physcs 219 (2006) 13 20 Short note One-sded fnte-dfference approxmatons sutable for use wth Rchardson extrapolaton Kumar Rahul, S.N. Bhattacharyya * Department of Mechancal Engneerng,

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

Dynamic Programming. Preview. Dynamic Programming. Dynamic Programming. Dynamic Programming (Example: Fibonacci Sequence)

Dynamic Programming. Preview. Dynamic Programming. Dynamic Programming. Dynamic Programming (Example: Fibonacci Sequence) /24/27 Prevew Fbonacc Sequence Longest Common Subsequence Dynamc programmng s a method for solvng complex problems by breakng them down nto smpler sub-problems. It s applcable to problems exhbtng the propertes

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

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

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

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

International Journal of Mathematical Archive-3(3), 2012, Page: Available online through ISSN

International Journal of Mathematical Archive-3(3), 2012, Page: Available online through   ISSN Internatonal Journal of Mathematcal Archve-3(3), 2012, Page: 1136-1140 Avalable onlne through www.ma.nfo ISSN 2229 5046 ARITHMETIC OPERATIONS OF FOCAL ELEMENTS AND THEIR CORRESPONDING BASIC PROBABILITY

More information

EEL 6266 Power System Operation and Control. Chapter 3 Economic Dispatch Using Dynamic Programming

EEL 6266 Power System Operation and Control. Chapter 3 Economic Dispatch Using Dynamic Programming EEL 6266 Power System Operaton and Control Chapter 3 Economc Dspatch Usng Dynamc Programmng Pecewse Lnear Cost Functons Common practce many utltes prefer to represent ther generator cost functons as sngle-

More information

Inductance Calculation for Conductors of Arbitrary Shape

Inductance Calculation for Conductors of Arbitrary Shape CRYO/02/028 Aprl 5, 2002 Inductance Calculaton for Conductors of Arbtrary Shape L. Bottura Dstrbuton: Internal Summary In ths note we descrbe a method for the numercal calculaton of nductances among conductors

More information

A METHOD FOR DETECTING OUTLIERS IN FUZZY REGRESSION

A METHOD FOR DETECTING OUTLIERS IN FUZZY REGRESSION OPERATIONS RESEARCH AND DECISIONS No. 2 21 Barbara GŁADYSZ* A METHOD FOR DETECTING OUTLIERS IN FUZZY REGRESSION In ths artcle we propose a method for dentfyng outlers n fuzzy regresson. Outlers n a sample

More information

Calculation of time complexity (3%)

Calculation of time complexity (3%) Problem 1. (30%) Calculaton of tme complexty (3%) Gven n ctes, usng exhaust search to see every result takes O(n!). Calculaton of tme needed to solve the problem (2%) 40 ctes:40! dfferent tours 40 add

More information

Lectures - Week 4 Matrix norms, Conditioning, Vector Spaces, Linear Independence, Spanning sets and Basis, Null space and Range of a Matrix

Lectures - Week 4 Matrix norms, Conditioning, Vector Spaces, Linear Independence, Spanning sets and Basis, Null space and Range of a Matrix Lectures - Week 4 Matrx norms, Condtonng, Vector Spaces, Lnear Independence, Spannng sets and Bass, Null space and Range of a Matrx Matrx Norms Now we turn to assocatng a number to each matrx. We could

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

Module 14: THE INTEGRAL Exploring Calculus

Module 14: THE INTEGRAL Exploring Calculus Module 14: THE INTEGRAL Explorng Calculus Part I Approxmatons and the Defnte Integral It was known n the 1600s before the calculus was developed that the area of an rregularly shaped regon could be approxmated

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

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

Numerical Solution of Ordinary Differential Equations

Numerical Solution of Ordinary Differential Equations Numercal Methods (CENG 00) CHAPTER-VI Numercal Soluton of Ordnar Dfferental Equatons 6 Introducton Dfferental equatons are equatons composed of an unknown functon and ts dervatves The followng are examples

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

APPENDIX A Some Linear Algebra

APPENDIX A Some Linear Algebra APPENDIX A Some Lnear Algebra The collecton of m, n matrces A.1 Matrces a 1,1,..., a 1,n A = a m,1,..., a m,n wth real elements a,j s denoted by R m,n. If n = 1 then A s called a column vector. Smlarly,

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

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

Workshop: Approximating energies and wave functions Quantum aspects of physical chemistry

Workshop: Approximating energies and wave functions Quantum aspects of physical chemistry Workshop: Approxmatng energes and wave functons Quantum aspects of physcal chemstry http://quantum.bu.edu/pltl/6/6.pdf Last updated Thursday, November 7, 25 7:9:5-5: Copyrght 25 Dan Dll (dan@bu.edu) Department

More information

More metrics on cartesian products

More metrics on cartesian products More metrcs on cartesan products If (X, d ) are metrc spaces for 1 n, then n Secton II4 of the lecture notes we defned three metrcs on X whose underlyng topologes are the product topology The purpose of

More information

The lower and upper bounds on Perron root of nonnegative irreducible matrices

The lower and upper bounds on Perron root of nonnegative irreducible matrices Journal of Computatonal Appled Mathematcs 217 (2008) 259 267 wwwelsevercom/locate/cam The lower upper bounds on Perron root of nonnegatve rreducble matrces Guang-Xn Huang a,, Feng Yn b,keguo a a College

More information

CHAPTER-5 INFORMATION MEASURE OF FUZZY MATRIX AND FUZZY BINARY RELATION

CHAPTER-5 INFORMATION MEASURE OF FUZZY MATRIX AND FUZZY BINARY RELATION CAPTER- INFORMATION MEASURE OF FUZZY MATRI AN FUZZY BINARY RELATION Introducton The basc concept of the fuzz matr theor s ver smple and can be appled to socal and natural stuatons A branch of fuzz matr

More information

CS 331 DESIGN AND ANALYSIS OF ALGORITHMS DYNAMIC PROGRAMMING. Dr. Daisy Tang

CS 331 DESIGN AND ANALYSIS OF ALGORITHMS DYNAMIC PROGRAMMING. Dr. Daisy Tang CS DESIGN ND NLYSIS OF LGORITHMS DYNMIC PROGRMMING Dr. Dasy Tang Dynamc Programmng Idea: Problems can be dvded nto stages Soluton s a sequence o decsons and the decson at the current stage s based on the

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

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

Heuristic Algorithm for Finding Sensitivity Analysis in Interval Solid Transportation Problems

Heuristic Algorithm for Finding Sensitivity Analysis in Interval Solid Transportation Problems Internatonal Journal of Innovatve Research n Advanced Engneerng (IJIRAE) ISSN: 349-63 Volume Issue 6 (July 04) http://rae.com Heurstc Algorm for Fndng Senstvty Analyss n Interval Sold Transportaton Problems

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

DUE: WEDS FEB 21ST 2018

DUE: WEDS FEB 21ST 2018 HOMEWORK # 1: FINITE DIFFERENCES IN ONE DIMENSION DUE: WEDS FEB 21ST 2018 1. Theory Beam bendng s a classcal engneerng analyss. The tradtonal soluton technque makes smplfyng assumptons such as a constant

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

Min Cut, Fast Cut, Polynomial Identities

Min Cut, Fast Cut, Polynomial Identities Randomzed Algorthms, Summer 016 Mn Cut, Fast Cut, Polynomal Identtes Instructor: Thomas Kesselhem and Kurt Mehlhorn 1 Mn Cuts n Graphs Lecture (5 pages) Throughout ths secton, G = (V, E) s a mult-graph.

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

Global Sensitivity. Tuesday 20 th February, 2018

Global Sensitivity. Tuesday 20 th February, 2018 Global Senstvty Tuesday 2 th February, 28 ) Local Senstvty Most senstvty analyses [] are based on local estmates of senstvty, typcally by expandng the response n a Taylor seres about some specfc values

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

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

A New Refinement of Jacobi Method for Solution of Linear System Equations AX=b

A New Refinement of Jacobi Method for Solution of Linear System Equations AX=b Int J Contemp Math Scences, Vol 3, 28, no 17, 819-827 A New Refnement of Jacob Method for Soluton of Lnear System Equatons AX=b F Naem Dafchah Department of Mathematcs, Faculty of Scences Unversty of Gulan,

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

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

Outline and Reading. Dynamic Programming. Dynamic Programming revealed. Computing Fibonacci. The General Dynamic Programming Technique

Outline and Reading. Dynamic Programming. Dynamic Programming revealed. Computing Fibonacci. The General Dynamic Programming Technique Outlne and Readng Dynamc Programmng The General Technque ( 5.3.2) -1 Knapsac Problem ( 5.3.3) Matrx Chan-Product ( 5.3.1) Dynamc Programmng verson 1.4 1 Dynamc Programmng verson 1.4 2 Dynamc Programmng

More information

Applied Mathematics Letters

Applied Mathematics Letters Appled Matheatcs Letters 2 (2) 46 5 Contents lsts avalable at ScenceDrect Appled Matheatcs Letters journal hoepage: wwwelseverco/locate/al Calculaton of coeffcents of a cardnal B-splne Gradr V Mlovanovć

More information

COEFFICIENT DIAGRAM: A NOVEL TOOL IN POLYNOMIAL CONTROLLER DESIGN

COEFFICIENT DIAGRAM: A NOVEL TOOL IN POLYNOMIAL CONTROLLER DESIGN Int. J. Chem. Sc.: (4), 04, 645654 ISSN 097768X www.sadgurupublcatons.com COEFFICIENT DIAGRAM: A NOVEL TOOL IN POLYNOMIAL CONTROLLER DESIGN R. GOVINDARASU a, R. PARTHIBAN a and P. K. BHABA b* a Department

More information

Chapter 3 Differentiation and Integration

Chapter 3 Differentiation and Integration MEE07 Computer Modelng Technques n Engneerng Chapter Derentaton and Integraton Reerence: An Introducton to Numercal Computatons, nd edton, S. yakowtz and F. zdarovsky, Mawell/Macmllan, 990. Derentaton

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

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

Chapter 6. Supplemental Text Material

Chapter 6. Supplemental Text Material Chapter 6. Supplemental Text Materal S6-. actor Effect Estmates are Least Squares Estmates We have gven heurstc or ntutve explanatons of how the estmates of the factor effects are obtaned n the textboo.

More information

Lecture 5 Decoding Binary BCH Codes

Lecture 5 Decoding Binary BCH Codes Lecture 5 Decodng Bnary BCH Codes In ths class, we wll ntroduce dfferent methods for decodng BCH codes 51 Decodng the [15, 7, 5] 2 -BCH Code Consder the [15, 7, 5] 2 -code C we ntroduced n the last lecture

More information

Fuzzy Boundaries of Sample Selection Model

Fuzzy Boundaries of Sample Selection Model Proceedngs of the 9th WSES Internatonal Conference on ppled Mathematcs, Istanbul, Turkey, May 7-9, 006 (pp309-34) Fuzzy Boundares of Sample Selecton Model L. MUHMD SFIIH, NTON BDULBSH KMIL, M. T. BU OSMN

More information

Inner Product. Euclidean Space. Orthonormal Basis. Orthogonal

Inner Product. Euclidean Space. Orthonormal Basis. Orthogonal Inner Product Defnton 1 () A Eucldean space s a fnte-dmensonal vector space over the reals R, wth an nner product,. Defnton 2 (Inner Product) An nner product, on a real vector space X s a symmetrc, blnear,

More information

Supplementary Notes for Chapter 9 Mixture Thermodynamics

Supplementary Notes for Chapter 9 Mixture Thermodynamics Supplementary Notes for Chapter 9 Mxture Thermodynamcs Key ponts Nne major topcs of Chapter 9 are revewed below: 1. Notaton and operatonal equatons for mxtures 2. PVTN EOSs for mxtures 3. General effects

More information

Lecture 21: Numerical methods for pricing American type derivatives

Lecture 21: Numerical methods for pricing American type derivatives Lecture 21: Numercal methods for prcng Amercan type dervatves Xaoguang Wang STAT 598W Aprl 10th, 2014 (STAT 598W) Lecture 21 1 / 26 Outlne 1 Fnte Dfference Method Explct Method Penalty Method (STAT 598W)

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

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

Using T.O.M to Estimate Parameter of distributions that have not Single Exponential Family

Using T.O.M to Estimate Parameter of distributions that have not Single Exponential Family IOSR Journal of Mathematcs IOSR-JM) ISSN: 2278-5728. Volume 3, Issue 3 Sep-Oct. 202), PP 44-48 www.osrjournals.org Usng T.O.M to Estmate Parameter of dstrbutons that have not Sngle Exponental Famly Jubran

More information

The Expectation-Maximization Algorithm

The Expectation-Maximization Algorithm The Expectaton-Maxmaton Algorthm Charles Elan elan@cs.ucsd.edu November 16, 2007 Ths chapter explans the EM algorthm at multple levels of generalty. Secton 1 gves the standard hgh-level verson of the algorthm.

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

Linear Regression Analysis: Terminology and Notation

Linear Regression Analysis: Terminology and Notation ECON 35* -- Secton : Basc Concepts of Regresson Analyss (Page ) Lnear Regresson Analyss: Termnology and Notaton Consder the generc verson of the smple (two-varable) lnear regresson model. It s represented

More information

LECTURE 9 CANONICAL CORRELATION ANALYSIS

LECTURE 9 CANONICAL CORRELATION ANALYSIS LECURE 9 CANONICAL CORRELAION ANALYSIS Introducton he concept of canoncal correlaton arses when we want to quantfy the assocatons between two sets of varables. For example, suppose that the frst set of

More information

Finding Dense Subgraphs in G(n, 1/2)

Finding Dense Subgraphs in G(n, 1/2) Fndng Dense Subgraphs n Gn, 1/ Atsh Das Sarma 1, Amt Deshpande, and Rav Kannan 1 Georga Insttute of Technology,atsh@cc.gatech.edu Mcrosoft Research-Bangalore,amtdesh,annan@mcrosoft.com Abstract. Fndng

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

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

U.C. Berkeley CS294: Spectral Methods and Expanders Handout 8 Luca Trevisan February 17, 2016

U.C. Berkeley CS294: Spectral Methods and Expanders Handout 8 Luca Trevisan February 17, 2016 U.C. Berkeley CS94: Spectral Methods and Expanders Handout 8 Luca Trevsan February 7, 06 Lecture 8: Spectral Algorthms Wrap-up In whch we talk about even more generalzatons of Cheeger s nequaltes, and

More information

FUZZY FINITE ELEMENT METHOD

FUZZY FINITE ELEMENT METHOD FUZZY FINITE ELEMENT METHOD RELIABILITY TRUCTURE ANALYI UING PROBABILITY 3.. Maxmum Normal tress Internal force s the shear force, V has a magntude equal to the load P and bendng moment, M. Bendng moments

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

Temperature. Chapter Heat Engine

Temperature. Chapter Heat Engine Chapter 3 Temperature In prevous chapters of these notes we ntroduced the Prncple of Maxmum ntropy as a technque for estmatng probablty dstrbutons consstent wth constrants. In Chapter 9 we dscussed the

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

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

Graph Reconstruction by Permutations

Graph Reconstruction by Permutations Graph Reconstructon by Permutatons Perre Ille and Wllam Kocay* Insttut de Mathémathques de Lumny CNRS UMR 6206 163 avenue de Lumny, Case 907 13288 Marselle Cedex 9, France e-mal: lle@ml.unv-mrs.fr Computer

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

Yong Joon Ryang. 1. Introduction Consider the multicommodity transportation problem with convex quadratic cost function. 1 2 (x x0 ) T Q(x x 0 )

Yong Joon Ryang. 1. Introduction Consider the multicommodity transportation problem with convex quadratic cost function. 1 2 (x x0 ) T Q(x x 0 ) Kangweon-Kyungk Math. Jour. 4 1996), No. 1, pp. 7 16 AN ITERATIVE ROW-ACTION METHOD FOR MULTICOMMODITY TRANSPORTATION PROBLEMS Yong Joon Ryang Abstract. The optmzaton problems wth quadratc constrants often

More information

Application of B-Spline to Numerical Solution of a System of Singularly Perturbed Problems

Application of B-Spline to Numerical Solution of a System of Singularly Perturbed Problems Mathematca Aeterna, Vol. 1, 011, no. 06, 405 415 Applcaton of B-Splne to Numercal Soluton of a System of Sngularly Perturbed Problems Yogesh Gupta Department of Mathematcs Unted College of Engneerng &

More information

The Geometry of Logit and Probit

The Geometry of Logit and Probit The Geometry of Logt and Probt Ths short note s meant as a supplement to Chapters and 3 of Spatal Models of Parlamentary Votng and the notaton and reference to fgures n the text below s to those two chapters.

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

On the Repeating Group Finding Problem

On the Repeating Group Finding Problem The 9th Workshop on Combnatoral Mathematcs and Computaton Theory On the Repeatng Group Fndng Problem Bo-Ren Kung, Wen-Hsen Chen, R.C.T Lee Graduate Insttute of Informaton Technology and Management Takmng

More information

THE SUMMATION NOTATION Ʃ

THE SUMMATION NOTATION Ʃ Sngle Subscrpt otaton THE SUMMATIO OTATIO Ʃ Most of the calculatons we perform n statstcs are repettve operatons on lsts of numbers. For example, we compute the sum of a set of numbers, or the sum of the

More information

x = , so that calculated

x = , so that calculated Stat 4, secton Sngle Factor ANOVA notes by Tm Plachowsk n chapter 8 we conducted hypothess tests n whch we compared a sngle sample s mean or proporton to some hypotheszed value Chapter 9 expanded ths to

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

Vapnik-Chervonenkis theory

Vapnik-Chervonenkis theory Vapnk-Chervonenks theory Rs Kondor June 13, 2008 For the purposes of ths lecture, we restrct ourselves to the bnary supervsed batch learnng settng. We assume that we have an nput space X, and an unknown

More information