arxiv: v1 [math.oc] 10 Dec 2018 R. Jana a,1, A. K. Das b,2, Vishnu Narayan Mishra c,3

Transcription

1 Iterative Descet Method for Geeralized Leotief Model arxiv: v1 [mathoc] 10 Dec 2018 R Jaa a,1, A K Das b,2, Vishu Naraya Mishra c,3 a Jadavpur Uiversity, Kolkata b Idia Statistical Istitute, Kolkata c Idira Gadhi Natioal Tribal Uiversity, Amarkatak 1 rwitamjaaju@gmailcom 2 akdas@isicalaci 3 vishuarayamishra@gmailcom Abstract I this paper we cosider geeralized Leotief model We show that uder certai coditio the geeralized Leotief model is solvable by iterative descet method based o ifeasible iterior poit algorithm We prove the covergece of the method from strictly positive startig poit A umerical example is preseted to demostrate the performace of the algorithm Keywords: Leotief model, geeralized Leotief model, liear complemetarity problem, vertical liear complemetarity problem, ifeasible iterior poit algorithm 1 Itroductio The purpose of Leotief model [9] is to fid the iterrelatioship amog goods ad services for differet sectors of the ecoomy Leotief model cosiders productio of items withi some idustries where umber of idustries ad umber of products are equal I other words the model idicates a balace betwee demad ad supply The model is very useful to aalyze the atioal ecoomy of ay sector as each of the idustries uses iput from itself ad other idustries to produce a particular product Leotief model is classified as ope model ad closed model (see [2]) Ope model deals with fidig the productio level based o exteral demad whereas the closed model deals oly with iteral demad The iput-output model has wide applicatios i the area of regioal ecoomics [14], iteratioal trade [1], multi facility ivetory systems [13] The Leotief model describes a facilitated view of a ecoomical situatio The target of this model is to state the exact level of productio for each of various types of services or goods Suppose b j deotes the uits available 3 Correspodig author 1

2 or required at idustry j, ad a jk be the techical coefficiets represetig uits of output of sector j required per uit output of sector k The et output b j is ormally called the fial demad of the jth good Suppose the basic iput-output equatios are give as x j = a jk x k +b j (11) k=1 the we ca fid x j which is the output from idustry j There is a assumptio that each idustry or sector produces oly oe output I case if a idustry produces more tha oe output, the the aalysis is doe by aggregatio Several variats of the iput-output model are available i the literature [11] Leotief model cosiders sigle techology Ebiefug et al [6] itroduced geeralized Leotief model by cosiderig multiple techologies I this paper we cosider a ifeasible iterior poit method i lie with Kojima et al [7] ad show that geeralized Leotief model ca be solved usig this method by cotrollig step legths The paper is orgaized as follows Sectio 2 presets otatios ad some basic results related to the Leotief model ad geeralized Leotief model I sectio 3, we propose a ifeasible iterior poit method i lie with Kojima et al [7] to solve geeralized Leotief model We prove the covergece of the algorithm A umerical example is give i sectio 4 to illustrate the performace of the proposed algorithm 2 Prelimiaries We cosider matrices ad vectors with real etries Ay vector z R is a colum vector, z T deotes the traspose of z z k deotes the kth compoet of the vector z e deotes the vector of all 1 z deotes the euclidea orm of the vector z For ay matrix A R, A T deotes its traspose R deotes real -dimesioal space R + ad R ++ doote the oegative ad positive orthat i R respectively 21 Leotief Model Suppose there are idustries ad let b j be the demad for output j where j = 1,2,, b j 0 implies the requiremet at the ed of the period ad b j 0 represets the umber of uits available at the begiig of the period Geeralizig the Leotief model (11), the requiremet that each output at idustry j ca be writte as, x j b j + a jk x k (21) Agai Equatio 21 ca be rewritte as k=1 x 0 (22) b+(i A)x 0 (23) 2

3 x T [ b+(i A)x] = 0 (24) wherea = (a jk )adi istheidetitymatrixoforder HeceEquatio22, Equatio 23 ad Equatio 24 ca be writte as LCP(q,M) by takig q = b ad M = (I A) The complemetary coditio implies either x j > 0, the ( b+mx) j = 0 ie productio of output j at idustry j meets demad exactly or x j = 0 the from Equatio 21, we get b j > k=1 a jkx k ie o productio is ecessary at idustry j Now we cosider the followig result o Leotief model which will be used i the subsequet sectio Theorem 21 [2] I a ope Leotief model with iput matrix T ad let A = I T The the followig statemets are equivalet: 1 The model is feasible 2 The model is profitable 3 A is o sigular M matrix Proof If the model is feasible the by choosig demad vector d > 0, it follows that there exists a x 0 with Ax = d > 0 This coditio characterizes the property of o sigular M matrices Coversely, suppose A is o sigular M matrix, the A 1 0 Thus Ax = d has a oegative solutio x = A 1 d for each d 0 The equivalece of the model is profitable ad A is osigular M matrix follows from the fact that A is o sigular M matrix if ad oly if A T is o sigular M matrix There are several solutio procedure of Leotief model developed i the recet times Datzig [4] preseted a method to solve this model usig substitutio techique 22 Vertical Liear Complemetarity Problem Cosider a rectagular matrix N oforder m k with m k Suppose N ispartitioed row-wise ito k blocks i the form N 1 N 2 N = N k where each N j = ( j rs ) Rm j k with k j=1 m j = m The N is called vertical block matrix of type (m 1,m 2,,m k ) If m j = 1, j = 1,2,,k, the N is a square matrix The cocept of vertical block was itroduced by Cottle ad Datzig [3] i coectio with the geeralizatio of liear complemetarity problem Cottle- Datzig s geeralizatio ivolves a system w Nz = q, w 0, z 0, where N R m k, m k ad the variable w 1,w 2,,w m are partitioed ito k oempty sets S j,j = 1,2,,k Let T j = S j {z j }, j = 1,2,,k This problem is to fid a 3

4 solutio pair (w,z),w R m, z R k of the system such that atleast oe member of each set T j is o basic The formal statemet of this problem is as follows: Give a m k vertical block matrix A of type (m 1,m 2,,m k ) ad a give vector q R m where m = k j=1 m j, fid w R m ad z R k such that w Az = q, w 0, z 0 mj z j i=1 wj i = 0, for j = 1,2,,k This geeralizatio is kow as vertical liear complemetarity problem ad this problem is deoted as VLCP(q, A) I recet years, a umber of applicatios of the vertical liear complemetarity problem have bee reported i the literature Cottle-Datzig algorithm [3] is well kow for solvig vertical liear complemetarity problem For pivotal algorithms see [5] 221 Results i VLCP Theory We give some defiitios ad results which will be required i the ext sectio Defiitio 21 Let A be a vertical block matrix of type (m 1,m 2,,m k ) A submatrix of size k of A is called represetative submatrix if its jth row is draw from jth block A j of A Remark 21 A vertical block matrix A of type (m 1,m 2,,m k ) has atmost k j=1 m j distict represetative submatrices Defiitio 22 A vertical block matrix A of type (m 1,m 2,,m k ) is called vertical block P(P 0 ) matrix, if all its represetative submatrices are P(P 0 ) matrices Moha et al [10] cosider a vertical block matrix A of type (m 1,m 2,,m k ) where m j is the size of the jth block Costruct a matrix M by copyig A j ie jth colum of A, m j times for j = 1,2,,k This leads to a square matrix M of order m m from A M is said to be equivalet square matrix of A Now we state the followig theorem which will be required for our proposed algorithm Theorem 22 [10] Give the VLCP(q,A), let M be the equivalet square matrix of A The VLCP(q,A) has a solutio if ad oly if LCP(q,M) has a solutio 23 Geeralized Leotief Model Ebiefug ad Kostreva [6] exteded Leotief model cosiderig multiple techology Further they formulated the geeralized Leotief model as vertical liear complemetarity problem The formulatio is give below Cosider there are total m j 1 umber of differet techologies ad correspodig output with atleast oe m j > 1 for j = 1,2,, The geeralized Leotief model ca be writte as x j b j i + a j ik x k,i = 1,2, m j (25) k=1 4

5 For j = 1,2,,, b j i deotes the demad for output i at idustry j If bj i 0 the b j i represets quatity of goods i to be produced by idustry j ad if b j i 0 the b j i represets the umber of uits already available at the begiig of the period to satisfy demad Cosider a matrix E of order m k as follows: e e 2 0 E = 0 0 e k where e j be a colum vector of dimesio m j 1 with each compoet 1 Now by settig m = j=1 m j, we obtai E of dimesio m with m Cosider We write b = b j = (b j i )m j i=1, A j = (a j ik )m j, b 1 b 2 b i,k=1, A = A 1 A 2 A Note that A is a matrix of order m We write N = E A Here N is a vertical matrix of order m ad of type (m 1,m 2,,m ) The from the geeralized Leotief model 25, we write x 0 Nx b mj ad x j i=1 (Nj x b j ) i = 0, for j = 1,2,, Note that complemetary coditio states miimum cost requiremet Ebiefug [6] exteded the Chadrasekhara algorithm to gave a approach for solvig geeralized Leotief model 3 Results I this sectio we propose a iterative descet method based o ifeasible iterior poit algorithm to solve a geeralized Leotief model We defie the feasible regio of the LCP(q,M) as FEA(q,M) FEA(q,M) = {(z,w) R : z 0, w 0, w = q +Mz} ad iterior of the set FEA(q,M) as FEA + (q,m) = {(z,w) R : z > 0, w > 0} 5

6 Thealgorithmmovesfromthecurretiterate(z k,w k )tothesolutioofthelcp(q,m) by itroducig (z k+1,w k+1 ) defied as z k+1 = z k +α k d k z w k+1 = w k +α k d k w where α k is the suitable step legth of the algorithm The geerated sequece {(z k,w k )} is required to satisfy z k > 0, w k > 0 Now we defie cetral trajectory as the set of solutios (z,w) > 0 to the system of equatios w = q +Mz, ZWe = µe, for every µ > 0 Here the eighborhood N of the cetral trajectory be defied as N = {(z,w) > 0 : z i w i γ zt w (i = 1,2,,),zT w γ w Mz q or w Mz q ǫ} where ǫ > 0, γ > 0 ad γ (0,1) Startig from a strictly positive poit (z 0,w 0 ) the algorithm iteratively geerates a sequece {(z k,w k )} We defie a o-liear system with o-egative costraits [ ] Mz +w q F(z,w) = = 0 (31) ZWe µe where Z = diag(z i ), W = diag(w i ) ad e is the vector of all 1 s We obtai [ ] M I F (z,w) = W Z We apply Newto method to fid the search directio (d k z, dk w ) for the algorithm Now we solve [ M I W k Z k ][ d k z d k w where σ [0,1) ad µ k = σ zkt w k Hece Therefore we get ] = [ Mz k w k +q Z k W k e+µ k e Md k z +dk w = Mzk w k +q W k d k z +Zk d k w = Zk W k e+µ k e d k z = (Zk M +W k ) 1 (Z k q +Z k Mz k µ k e) (32) d k w = M(zk +d k z ) wk +q (33) Select the step legth α k suitably such that the algorithm geerates strictly positive poits (z k,w k > 0) i every step We cosider the merit fuctio as give i [12] φ(z,w) = w Mz q 2 + ZWe 2 We show that value of the merit fuctio at each step reduces ad algorithm stops whe φ(z,w) δ for some pre-determied δ > 0 Now we state our propose algorithm for solvig geeralized Leotief model give i (25) 6 ]

7 31 Algorithm Step I: Let M be a equivalet square matrix of the vertical matrix A usig the Theorem 22 of [10] Step II: Let k = 0 Let (z 0,w 0 ) > 0, δ > 0, β (0,1/2], γ (0,1), σ = 05 Compute the value of merit fuctio φ(z k,w k ) Step III: If φ(z k,w k ) δ, STOP ad (z k,w k ) is a approximate solutio to the LCP(q,M) Otherwise go to Step IV Step IV: Let µ k = σ zkt w k ad fid (d k z,d k w) from Equatio (32) ad (33) Step V: Compute α k so that (z k+1,w k+1 ) = (z k,w k )+α k (d k z,d k w) N, φ(z k+1,w k+1 ) φ(z k,w k ) α k β φ(z k,w k ) T (d k z,d k w) Step VI: Set k = k +1 ad go to the Step I To process geeralized Leotief model, we assume the matrix MZ + W is osigular ateach step of thealgorithmfor ay diagoal matrices Z ad W withstrictly positive elemets Remark 32 Note that we assume the o-sigularity of the matrix MZ + W to fid the solutio of geeralized Leotief model However i case of solvig Leotief model, similar assumptio is ot required Now we prove the followig theorem Theorem 31 I a Leotief model, det(mz + W) 0 for ay positive diagoal matrices Z ad W Proof It is kow that to for ay two strictly positive Z ad W the followig is true: [ ] M I (MZ +W) is o-sigular = W Z is o-sigular Now i LCP formulatio of Leotief model, [ M is a P-matrix ] by the M I Theorem21 AgaibyLemma(41)byKojimaetal [8], is o-sigular W Z if ad oly if M P 0 Hece det(mz +W) 0 7

8 32 Covergece Aalysis I this sectio we show that the algorithm preseted i the previous sectio coverges to the solutio of geeralized Leotief model We defie N = {(z,w) N : z T w ǫ, (z,w) ω } (34) where ǫ = mi{ǫ,γ ǫ} ad ω > 0 Firstly the ewto directio (d k z,dk w ) determied by the system of equatios is a cotiuous fuctio such that (z k,w k ) N As the matrix [ M I W k Z k ] is o-sigular for ay (z k,w k ) N Hece the Newto directio (d k z,d k w) is uiformly bouded for all (z k,w k ) over the set N Hece we ca fid positive costats η 1, η 2 such that (d k z,d k w) computed at every iteratios satisfies d k z i d k w i γd k zt d k w / η 1 (35) d k zt d k w η 2 (36) Firstly we show that the geerated sequece from the algorithm {z k,w k } satisfies the followig coditio Lemma 31 The geerated sequece {z k,w k } by the algorithm satisfies (i) Mz k+1 +w k+1 q = (1 α k )( Mz k +w k q) (ii) z kt d k w +w kt d k z = (1 σ)z kt w k Proof (i) From the Newto equatio we get [ ][ ] [ M I d k z Mz W k Z k = k w k +q Z k W k e+µ k e d k w where σ [0,1) ad µ k = σ zkt w k Hece, Md k z +dk w = q +Mzk w k Md k z Mz k +d k w +w k = q Mα k d k z Mα kz k +α k d k w +α kw k = α k q M(z k+1 z k ) Mα k z k +(w k+1 w k )+α k w k = α k q Mz k+1 +w k+1 q = Mz k +Mα k z k +w k α k w k q +α k q Mz k+1 +w k+1 q = ( Mz k +w k q) α k ( Mz k +w k q) Mz k+1 +w k+1 q = (1 α k )( Mz k +w k q) (ii) 2d part follows form the Newto equatio by restrictig µ k = σz kt w k We have W k d k z +Z k d k w = Z k W k e+µ k e z kt d k w +wkt d k z = zkt w k +σz kt w k z kt d k w +wkt d k z = (1 σ)zkt w k 8 ]

9 To specify the selectio of suitable step legth α k, we defie real valued fuctios g I (α) ad g II (α) as follows: gi I (α) = zk+1 i w k+1 i γ zk+1t w k+1 g II (α) = z k+1t w k+1 γ w k Mz k q where γ (0,1) The ext lemma will be required for the choice of α k The proof of the followig two lemmas are i lie with the idea give i [8] Lemma 32 If gi I i Proof (α) = zk+1 i w k+1 i γ zk+1t w k+1, the g I i (α) σ(1 γ)( ǫ/)α η 1α 2 gi I (α) = zk+1 i wi k+1 γ zk+1t w k+1 = zi kwk i +α(zk i dk w i +wi kdk z i )+α 2 d k z i d k w i γ zkt w k +α(z kt d k w+w kt d k z)+α 2 d k z = z k i wk i α(zk i wk i σzkt w k /)+α 2 d k z i d k w i γ zkt w k α(1 σ)z kt w k +α 2 d k z = (1 α)zi kwk i +ασzkt w k /+α 2 d k z i d k w i γ (1 α)zkt w k +ασz kt w k +α 2 d k z = (1 α)(zi kwk i γzkt w k /)+σ(1 γ)αz kt w k /+α 2 (d k z i d k w i γd k zt d k w /) σ(1 γ)( ǫ/)α η 1 α 2 (by(34), (35)) T d k w T d k w T d k w Lemma 33 If g II (α) = z k+1t w k+1 γ w k+1 Mz k+1 q, the g II (α) µ k α α 2 η 2 Proof g II (α) = z k+1t w k+1 γ w k+1 Mz k+1 q = (z k +αd k z )T (w k +αd k w ) γ(1 α) wk Mz k q = z kt w k +α[ (1 σ/)]z kt w k γ(1 α) w k Mz k q +α 2 d k zt d k w = (1 α)(z kt w k γ w k Mz k q )+µ k α+α 2 d k zt d k w µ k α α 2 d k zt d k w µ k α α 2 η 2 (by(36)) To fid step legth we have to fid α k to be the largest umber (0,1] for which g I i (α) 0 ad gii (α) 0 From Lemma (32) ad Lemma (33) we have g I i(α) 0, g II (α) 0 Hece we ca easily compute the lower boud of α k by solvig them for α Hece lettig α k = mi{1,σ ǫ(1 γ)/η 1,µ k /η 2 }, we obtai g I i (α) 0 g II (α) 0 9

10 hold for every α [0,α k ] ad i I [12], Simatiraki ad Shao state the idea to show the descet directios of the proposed algorithm We prove the followig theorem to show that the directios (d k z,dk w ) geerated by the algorithm is descet directio Theorem 32 The directios (d k z,d k w) geerated at the kth iteratio by the algorithm is a descet directio for the merit fuctio φ(z k,w k ) Proof Cosider p k = (z k,w k ), the where ē is the vector (0,e) T Hece Agai we have µ k = σ zkt w k Hece d k p = F (p k )[ F(p k )+µ k ē] φ(p k ) = 1 F(p k ) F (p k ) T F(p k ), φ(p k )d k p = 1 F(p k ) F(pk ) T F (p k )F (p k ) 1 [ F(p k )+µ k ē] = 1 F(p k ) ( F(pk ) 2 +µ k z kt w k ) = ( F(p k ) µ k z kt w k F(p k ) ) = (φ(p k ) µ k z kt w k φ(p k ) ) µ k z kt w k = σ (zkt w k ) 2 σ Z k W k e 2 σφ(p k ) 2 φ(p k ) T d k p φ(p k )(1 σ) 0 Hece (d k z,dk w ) is the decet directio for the algorithm Now we prove the covergece result of the proposed algorithm Theorem 33 Let the sequece p k = {z k,w k } be geerated by the algorithm The for ay { µ k z kt w k} (0,1) bouded away from zero, φ(pk ) coverges to zero Proof We have φ(p k+1 ) = α k β φ(p k ) T d k p φ(p k ) α k β(1 σ)φ(p k ) [1 α k β(1 σ)]φ(p k ) Hece, φ(p k+1 ) φ(p k ) 1 α k β(1 σ) Now α k s are bouded away from zero Hece correspodig subsequeces coverges to zero So lettig k, we get Hece {φ(p k )} coverges to zero φ(p k+1 ) φ(p k ) < 1 10

11 4 Numerical illustratio We cosider the example give by Ebiefug et al [6] where they cosider a ecoomy with three sectors The outputs of the sectors are pair of shoes, food ad light bulbs The techology matrix ad demads for the ecoomy are A = , h I = Here the colums ad rows are give i the order of shoes, food ad light bulbs i that order Cosider a ew techology for each of these sectors has bee itroduced i the market The ew techology matrix B ad demad vector h II are A = , h II = Now formulatig this problem uder as geeralized Leotief iput-output model, we ca get a vertical block matrix N of type (2,2,2) ad the demad vector b, A = , b =

12 Iteratio (k) z k w k d k z d k w µ φ(z k,w k )

13 Iteratio (k) z k w k d k z d k w µ φ(z k,w k ) Table 1: Summary of computatio for the proposed algorithm The algorithm rus o a HP PC with itel Core i5 processor 310 GHz 4 GB of RAM The proposed algorithm coverges to the solutio [415,0,54] T after 99 iteratios ad time take to reach the optimal solutio is s Ackowledgmet: The author R Jaa is thakful to the Departmet of Sciece ad Techology, Govt of Idia, INSPIRE Fellowship Scheme for fiacial support Refereces [1] Terece S Barker Foreig trade i multisectoral models Uiversity of Cambridge, Departmet of Applied Ecoomics, 1973 [2] Abraham Berma ad Robert J Plemmos Noegative matrices i the mathematical scieces, volume 9 Siam, 1994 [3] Richard W Cottle ad George B Datzig A geeralizatio of the liear complemetarity problem Joural of Combiatorial Theory, 8(1):79 90, 1970 [4] George B Datzig Optimal solutio of a dyamic leotief model with substitutio Ecoometrica: Joural of the Ecoometric Society, pages , 1955 [5] AK Das, R Jaa, et al Fiiteess of criss-cross method i complemetarity problem I Iteratioal Coferece o Mathematics ad Computig, pages Spriger, 2017 [6] Aieka A Ebiefug ad Michael M Kostreva The geeralized leotief iputoutput model ad its applicatio to the choice of ew techology Aals of Operatios Research, 44(2): , 1993 [7] Masakazu Kojima, Nimrod Megiddo, ad Shiji Mizuo A primaldual ifeasibleiterior-poit algorithm for liear programmig Mathematical programmig, 61(1-3): ,

14 [8] Masakazu Kojima, Nimrod Megiddo, Toshihito Noma, ad Akiko Yoshise A uified approach to iterior poit algorithms for liear complemetarity problems: A summary Operatios Research Letters, 10(5): , 1991 [9] Wassily Leotief Iput-output ecoomics Oxford Uiversity Press, 1986 [10] SR Moha, SK Neogy, ad R Sridhar The geeralized liear complemetarity problem revisited Mathematical Programmig, 74(2):197, 1996 [11] Ja Oosterhave O the plausibility of the supply-drive iput-output model Joural of Regioal Sciece, 28(2): , 1988 [12] Evagelia M Simatiraki ad David F Shao A ifeasible-iterior-poit method for liear complemetarity problems SIAM Joural o Optimizatio, 7(3): , 1997 [13] Arthur F Veiott Jr Miimum cocave-cost solutio of leotief substitutio models of multi-facility ivetory systems Operatios Research, 17(2): , 1969 [14] Thomas Wiedma A review of recet multi-regio iput output models used for cosumptio-based emissio ad resource accoutig Ecological Ecoomics, 69(2): ,