@@ -109,13 +109,13 @@ class Neumann:
109109 assert (self.A >= 0).all() and (self.B >= 0).all(), 'The input and \
110110 output matrices must have only non-negative entries!'
111111
112- # (1) Check whether Assumption 1 is satisfied:
112+ # (1) Check whether Assumption I is satisfied:
113113 if (np.sum(B, 0) <= 0).any():
114114 self.AI = False
115115 else:
116116 self.AI = True
117117
118- # (2) Check whether Assumption 2 is satisfied:
118+ # (2) Check whether Assumption II is satisfied:
119119 if (np.sum(A, 1) <= 0).any():
120120 self.AII = False
121121 else:
@@ -328,13 +328,13 @@ We use the following notation.
328328$\mathbf{0}$ denotes
329329a vector of zeros.
330330
331- We call an $n$-vector positive and write
331+ We call an $n$-vector ** positive** and write
332332$x\gg \mathbf{0}$ if $x_i>0$ for all $i=1,2,\dots,n$.
333333
334- We call a vector non-negative and write $x\geq \mathbf{0}$ if $x_i\geq 0$ for
334+ We call a vector ** non-negative** and write $x\geq \mathbf{0}$ if $x_i\geq 0$ for
335335all $i=1,2,\dots,n$.
336336
337- We call a vector semi-positive and written $x > \mathbf{0}$ if
337+ We call a vector ** semi-positive** and written $x > \mathbf{0}$ if
338338$x\geq \mathbf{0}$ and $x\neq \mathbf{0}$.
339339
340340For two conformable vectors $x$ and $y$, $x\gg y$,
@@ -389,18 +389,18 @@ Two key assumptions restrict economy $(A,B)$:
389389```
390390````
391391
392- A semi-positive * intensity* $m$-vector $x$ denotes levels at which
392+ A semi-positive * intensity* $m$-vector $x$ denotes levels at which
393393activities are operated.
394394
395395Therefore,
396396
397- - vector $x^TA $ gives the total amount of * goods used in
397+ - vector $x^T A $ gives the total amount of * goods used in
398398 production*
399- - vector $x^TB $ gives * total outputs*
399+ - vector $x^T B $ gives * total outputs*
400400
401401An economy $(A,B)$ is said to be * productive* , if there exists a
402402non-negative intensity vector $x \geq 0$ such
403- that $x^T B > x^TA $.
403+ that $x^T B > x^T A $.
404404
405405The semi-positive $n$-vector $p$ contains prices assigned to
406406the $n$ goods.
@@ -412,21 +412,24 @@ The $p$ vector implies *cost* and *revenue* vectors
412412
413413Satisfaction or a property of an input-output pair $(A,B)$ called * irreducibility*
414414(or indecomposability) determines whether an economy can be decomposed
415- into multiple "sub-economies".
415+ into multiple "sub-economies".
416416
417417``` {prf:definition}
418418For an economy $(A,B)$, the set of goods
419- $S\subset \{1,2,\dots,n\}$ is called an *independent subset* if
419+ $S\subset \{1,2,\dots,n\}$ is called an ** independent subset* * if
420420it is possible to produce every good in $S$ without consuming
421- goods from outside $S$. Formally, the set $S$ is independent if
421+ goods from outside $S$.
422+
423+ Formally, the set $S$ is independent if
422424$\exists T\subset \{1,2,\dots,m\}$ (a subset of activities) such
423- that $a_{i,j}=0$ $\forall i\in T$ and $j\in S^c$ and
424- for all $j\in S$, $\exists i\in T$ for which $b_{i,j}>0$.
425+ that $a_{i,j}=0$, $\forall i\in T$ and $j\in S^c$ and
426+ for all $j\in S$, $\exists i\in T$ for which $b_{i,j}>0$.
427+
425428The economy is **irreducible** if there are no proper independent
426429subsets.
427430```
428431
429- We study two examples, both in Chapter 9.6 of Gale {cite}` gale1989theory `
432+ We study two examples, both in Chapter 9.6 of Gale {cite}` gale1989theory `
430433
431434``` {code-cell} ipython3
432435# (1) Irreducible (A, B) example: α_0 = β_0
@@ -511,7 +514,7 @@ We follow John von Neumann in studying “balanced growth”.
511514Let $./$ denote an elementwise division of one vector by another and let
512515$\alpha >0$ be a scalar.
513516
514- Then * balanced growth* is a situation in which
517+ Then ** balanced growth* * is a situation in which
515518
516519$$
517520x_{t+1}./x_t = \alpha , \quad \forall t \geq 0
@@ -553,17 +556,17 @@ relationship between technological and valuation characteristics of
553556the economy:
554557
555558``` {prf:definition}
556- The *technological expansion problem* (TEP) for the economy
559+ The ** technological expansion problem* * (TEP) for the economy
557560$(A,B)$ is to find a semi-positive $m$-vector $x>0$
558561and a number $\alpha\in\mathbb{R}$ that satisfy
559- ```
560562
561563$$
562564\begin{aligned}
563565 &\max_{\alpha} \hspace{2mm} \alpha\\
564566 &\text{s.t. }\hspace{2mm}x^T B \geq \alpha x^T A
565567 \end{aligned}
566568$$
569+ ```
567570
568571Theorem 9.3 of David Gale’s book {cite}` gale1989theory ` asserts that if {prf: ref }` assumption1 ` and {prf: ref }` assumption2 ` are
569572both satisfied, then a maximum value of $\alpha$ exists and that it is
@@ -596,7 +599,7 @@ and the *economical expansion problem* are both linearly homogeneous,
596599the optimality of $x_0$ and $p_0$ are defined only up to a
597600positive scale factor.
598601
599- For convenience (and to emphasize a close connection to zero-sum games), we normalize both vectors
602+ For convenience (and to emphasize a close connection to zero-sum games), we normalize both vectors
600603$x_0$ and $p_0$ to have unit length.
601604
602605A standard duality argument (see Lemma 9.4. in (Gale, 1960) {cite}` gale1989theory ` ) implies
628631```
629632
630633``` {prf:proof} (Sketch)
634+
631635{prf:ref}`assumption1` and {prf:ref}`assumption2` imply that there exist $(\alpha_0,
632- x_0)$ and $(\beta_0, p_0)$ that solve the TEP and EEP, respectively. If
633- $\gamma^*>\alpha_0$, then by definition of $\alpha_0$, there cannot
636+ x_0)$ and $(\beta_0, p_0)$ that solve the TEP and EEP, respectively.
637+
638+ If $\gamma^*>\alpha_0$, then by definition of $\alpha_0$, there cannot
634639exist a semi-positive $x$ that satisfies $x^T B \geq \gamma^{* }
635- x^T A$. Similarly, if $\gamma^*<\beta_0$, there is no semi-positive
640+ x^T A$.
641+
642+ Similarly, if $\gamma^*<\beta_0$, there is no semi-positive
636643$p$ for which $Bp \leq \gamma^{* } Ap$. Let $\gamma^{*
637644}\in[\beta_0, \alpha_0]$, then $x_0^T B \geq \alpha_0 x_0^T A \geq
638- \gamma^{* } x_0^T A$. Moreover, $Bp_0\leq \beta_0 A p_0\leq \gamma^* A
645+ \gamma^{* } x_0^T A$.
646+
647+ Moreover, $Bp_0\leq \beta_0 A p_0\leq \gamma^* A
639648p_0$. These two inequalities imply $x_0\left(B - \gamma^{* } A\right)p_0
640649= 0$.
641650```
@@ -655,7 +664,7 @@ be unused.
655664
656665Therefore, the conditions stated in {prf: ref }` theorem1 ` ex encode all equilibrium conditions.
657666
658- So {prf: ref }` theorem1 ` essentially states that under {prf: ref }` assumption1 ` and {prf: ref }` assumption2 ` there
667+ So {prf: ref }` theorem1 ` essentially states that under {prf: ref }` assumption1 ` and {prf: ref }` assumption2 ` there
659668always exists an equilibrium $\left(\gamma^{* }, x_0, p_0\right)$
660669with balanced growth.
661670
@@ -682,14 +691,15 @@ To compute the equilibrium $(\gamma^{*}, x_0, p_0)$, we follow the
682691algorithm proposed by Hamburger, Thompson and Weil (1967), building on
683692the key insight that an equilibrium (with balanced growth) can be
684693solves a particular two-player zero-sum game.
694+
685695First, we introduce some notation.
686696
687697Consider the $m\times n$ matrix $C$ as a payoff matrix,
688698with the entries representing payoffs from the ** minimizing** column
689699player to the ** maximizing** row player and assume that the players can
690700use mixed strategies. Thus,
691701
692- * the row player chooses the $m$-vector $x > \mathbf{0}$ subject to $\iota_m^T x = 1$
702+ * the row player chooses the $m$-vector $x > \mathbf{0}$ subject to $\iota_m^T x = 1$
693703* the column player chooses the $n$-vector $p > \mathbf{0}$ subject to $\iota_n^T p = 1$.
694704
695705``` {prf:definition}
703713\end{aligned}
704714$$
705715
706- The number $V(C)$ is called the *value* of the game.
716+ The number $V(C)$ is called the ** value* * of the game.
707717```
708718
709719From the above definition, it is clear that the value $V(C)$ has
729739
730740### Connection with Linear Programming (LP)
731741
732- Nash equilibria of a finite two-player zero-sum game solve a linear programming problem.
742+ Nash equilibria of a finite two-player zero-sum game solve a linear programming problem.
733743
734- To see this, we introduce
735- the following notation
744+ To see this, we introduce the following notation
736745
737746* For a fixed $x$, let $v$ be the value of the minimization problem: $v \equiv \min_p x^T C p = \min_j x^T C e^j$
738747* For a fixed $p$, let $u$ be the value of the maximization problem: $u \equiv \max_x x^T C p = \max_i (e^i)^T C p$
@@ -776,6 +785,7 @@ $x_0^T B \gg \mathbf{0}$, where $x_0$ is a maximizing
776785vector. Since $B$ is non-negative, this requires that each
777786column of $B$ has at least one positive entry, which is
778787{prf:ref}`assumption1`.
788+
779789* $\Leftarrow$ From {prf:ref}`assumption1` and the fact
780790that $p>\mathbf{0}$, it follows that $Bp > \mathbf{0}$.
781791This implies that the maximizing player can always choose $x$
@@ -797,9 +807,11 @@ calculating the solution of the game implies
797807 $\exists j\in\{ 1, \dots, n\} $, s.t.
798808 $[ x^T M(\gamma)] _ j < 0$ implying
799809 that $V(M(\gamma)) < 0$.
810+
800811- If $\gamma < \beta_0$, then for all $p>0$, there
801812 $\exists i\in\{ 1, \dots, m\} $, s.t.
802813 $[ M(\gamma)p] _ i > 0$ implying that $V(M(\gamma)) > 0$.
814+
803815- If $\gamma \in \{ \beta_0, \alpha_0\} $, then (by {prf: ref }` theorem1 ` ) the
804816 optimal intensity and price vectors $x_0$ and $p_0$
805817 satisfy
@@ -897,6 +909,7 @@ Compute $\alpha_0$ and $\beta_0$
897909 1 . If $V(M(\gamma)) \geq 0$, then set $LB = \gamma$,
898910 otherwise let $UB = \gamma$.
899911 1 . Iterate on 1. and 2. until $|UB - LB| < \epsilon$.
912+
900913- Finding $\beta_0$
901914 1 . Fix $\gamma = \frac{UB + LB}{2}$ and compute the solution
902915 of the two-player zero-sum game associated.
@@ -905,6 +918,7 @@ Compute $\alpha_0$ and $\beta_0$
905918 1 . If $V(M(\gamma)) > 0$, then set $LB = \gamma$,
906919 otherwise let $UB = \gamma$.
907920 1 . Iterate on 1. and 2. until $|UB - LB| < \epsilon$.
921+
908922- * Existence* : Since $V(M(LB))>0$ and $V(M(UB))<0$ and
909923 $V(M(\cdot))$ is a continuous, nonincreasing function, there is
910924 at least one $\gamma\in[ LB, UB] $, s.t. $V(M(\gamma))=0$.
@@ -1049,7 +1063,7 @@ The latter shows that $1/\alpha_0$ is a positive eigenvalue of
10491063$A$ and $x_0$ is the corresponding non-negative left
10501064eigenvector.
10511065
1052- The classic result of ** Perron and Frobenius* * implies
1066+ The classic result of * Perron and Frobenius* implies
10531067that a non-negative matrix has a non-negative
10541068eigenvalue-eigenvector pair.
10551069
0 commit comments