39. 输入-输出模型#

39.1. 概述#

在我们学习输入-输出模型之前，本讲座需要以下导入和安装。

!pip install quantecon_book_networks
!pip install quantecon
!pip install pandas-datareader

Show output

Hide code cell output

Collecting quantecon_book_networks

  Downloading quantecon_book_networks-1.5-py2.py3-none-any.whl.metadata (1.7 kB)
Requirement already satisfied: numpy in /home/runner/miniconda3/envs/quantecon/lib/python3.12/site-packages (from quantecon_book_networks) (1.26.4)
Requirement already satisfied: scipy in /home/runner/miniconda3/envs/quantecon/lib/python3.12/site-packages (from quantecon_book_networks) (1.13.1)
Requirement already satisfied: pandas in /home/runner/miniconda3/envs/quantecon/lib/python3.12/site-packages (from quantecon_book_networks) (2.2.2)
Requirement already satisfied: matplotlib in /home/runner/miniconda3/envs/quantecon/lib/python3.12/site-packages (from quantecon_book_networks) (3.9.2)
Requirement already satisfied: networkx in /home/runner/miniconda3/envs/quantecon/lib/python3.12/site-packages (from quantecon_book_networks) (3.3)
Requirement already satisfied: wbgapi in /home/runner/miniconda3/envs/quantecon/lib/python3.12/site-packages (from quantecon_book_networks) (1.0.12)
Requirement already satisfied: contourpy>=1.0.1 in /home/runner/miniconda3/envs/quantecon/lib/python3.12/site-packages (from matplotlib->quantecon_book_networks) (1.2.0)
Requirement already satisfied: cycler>=0.10 in /home/runner/miniconda3/envs/quantecon/lib/python3.12/site-packages (from matplotlib->quantecon_book_networks) (0.11.0)
Requirement already satisfied: fonttools>=4.22.0 in /home/runner/miniconda3/envs/quantecon/lib/python3.12/site-packages (from matplotlib->quantecon_book_networks) (4.51.0)
Requirement already satisfied: kiwisolver>=1.3.1 in /home/runner/miniconda3/envs/quantecon/lib/python3.12/site-packages (from matplotlib->quantecon_book_networks) (1.4.4)
Requirement already satisfied: packaging>=20.0 in /home/runner/miniconda3/envs/quantecon/lib/python3.12/site-packages (from matplotlib->quantecon_book_networks) (24.1)
Requirement already satisfied: pillow>=8 in /home/runner/miniconda3/envs/quantecon/lib/python3.12/site-packages (from matplotlib->quantecon_book_networks) (10.4.0)
Requirement already satisfied: pyparsing>=2.3.1 in /home/runner/miniconda3/envs/quantecon/lib/python3.12/site-packages (from matplotlib->quantecon_book_networks) (3.1.2)
Requirement already satisfied: python-dateutil>=2.7 in /home/runner/miniconda3/envs/quantecon/lib/python3.12/site-packages (from matplotlib->quantecon_book_networks) (2.9.0.post0)

Requirement already satisfied: pytz>=2020.1 in /home/runner/miniconda3/envs/quantecon/lib/python3.12/site-packages (from pandas->quantecon_book_networks) (2024.1)
Requirement already satisfied: tzdata>=2022.7 in /home/runner/miniconda3/envs/quantecon/lib/python3.12/site-packages (from pandas->quantecon_book_networks) (2023.3)
Requirement already satisfied: requests in /home/runner/miniconda3/envs/quantecon/lib/python3.12/site-packages (from wbgapi->quantecon_book_networks) (2.32.3)
Requirement already satisfied: PyYAML in /home/runner/miniconda3/envs/quantecon/lib/python3.12/site-packages (from wbgapi->quantecon_book_networks) (6.0.1)
Requirement already satisfied: tabulate in /home/runner/miniconda3/envs/quantecon/lib/python3.12/site-packages (from wbgapi->quantecon_book_networks) (0.9.0)
Requirement already satisfied: six>=1.5 in /home/runner/miniconda3/envs/quantecon/lib/python3.12/site-packages (from python-dateutil>=2.7->matplotlib->quantecon_book_networks) (1.16.0)
Requirement already satisfied: charset-normalizer<4,>=2 in /home/runner/miniconda3/envs/quantecon/lib/python3.12/site-packages (from requests->wbgapi->quantecon_book_networks) (3.3.2)
Requirement already satisfied: idna<4,>=2.5 in /home/runner/miniconda3/envs/quantecon/lib/python3.12/site-packages (from requests->wbgapi->quantecon_book_networks) (3.7)
Requirement already satisfied: urllib3<3,>=1.21.1 in /home/runner/miniconda3/envs/quantecon/lib/python3.12/site-packages (from requests->wbgapi->quantecon_book_networks) (2.2.3)
Requirement already satisfied: certifi>=2017.4.17 in /home/runner/miniconda3/envs/quantecon/lib/python3.12/site-packages (from requests->wbgapi->quantecon_book_networks) (2024.8.30)

Downloading quantecon_book_networks-1.5-py2.py3-none-any.whl (514 kB)

Installing collected packages: quantecon_book_networks

Successfully installed quantecon_book_networks-1.5

Requirement already satisfied: quantecon in /home/runner/miniconda3/envs/quantecon/lib/python3.12/site-packages (0.10.1)
Requirement already satisfied: numba>=0.49.0 in /home/runner/miniconda3/envs/quantecon/lib/python3.12/site-packages (from quantecon) (0.60.0)
Requirement already satisfied: numpy>=1.17.0 in /home/runner/miniconda3/envs/quantecon/lib/python3.12/site-packages (from quantecon) (1.26.4)
Requirement already satisfied: requests in /home/runner/miniconda3/envs/quantecon/lib/python3.12/site-packages (from quantecon) (2.32.3)
Requirement already satisfied: scipy>=1.5.0 in /home/runner/miniconda3/envs/quantecon/lib/python3.12/site-packages (from quantecon) (1.13.1)
Requirement already satisfied: sympy in /home/runner/miniconda3/envs/quantecon/lib/python3.12/site-packages (from quantecon) (1.13.2)
Requirement already satisfied: llvmlite<0.44,>=0.43.0dev0 in /home/runner/miniconda3/envs/quantecon/lib/python3.12/site-packages (from numba>=0.49.0->quantecon) (0.43.0)
Requirement already satisfied: charset-normalizer<4,>=2 in /home/runner/miniconda3/envs/quantecon/lib/python3.12/site-packages (from requests->quantecon) (3.3.2)
Requirement already satisfied: idna<4,>=2.5 in /home/runner/miniconda3/envs/quantecon/lib/python3.12/site-packages (from requests->quantecon) (3.7)
Requirement already satisfied: urllib3<3,>=1.21.1 in /home/runner/miniconda3/envs/quantecon/lib/python3.12/site-packages (from requests->quantecon) (2.2.3)
Requirement already satisfied: certifi>=2017.4.17 in /home/runner/miniconda3/envs/quantecon/lib/python3.12/site-packages (from requests->quantecon) (2024.8.30)
Requirement already satisfied: mpmath<1.4,>=1.1.0 in /home/runner/miniconda3/envs/quantecon/lib/python3.12/site-packages (from sympy->quantecon) (1.3.0)

Requirement already satisfied: pandas-datareader in /home/runner/miniconda3/envs/quantecon/lib/python3.12/site-packages (0.10.0)
Requirement already satisfied: lxml in /home/runner/miniconda3/envs/quantecon/lib/python3.12/site-packages (from pandas-datareader) (5.2.1)
Requirement already satisfied: pandas>=0.23 in /home/runner/miniconda3/envs/quantecon/lib/python3.12/site-packages (from pandas-datareader) (2.2.2)
Requirement already satisfied: requests>=2.19.0 in /home/runner/miniconda3/envs/quantecon/lib/python3.12/site-packages (from pandas-datareader) (2.32.3)
Requirement already satisfied: numpy>=1.26.0 in /home/runner/miniconda3/envs/quantecon/lib/python3.12/site-packages (from pandas>=0.23->pandas-datareader) (1.26.4)
Requirement already satisfied: python-dateutil>=2.8.2 in /home/runner/miniconda3/envs/quantecon/lib/python3.12/site-packages (from pandas>=0.23->pandas-datareader) (2.9.0.post0)
Requirement already satisfied: pytz>=2020.1 in /home/runner/miniconda3/envs/quantecon/lib/python3.12/site-packages (from pandas>=0.23->pandas-datareader) (2024.1)
Requirement already satisfied: tzdata>=2022.7 in /home/runner/miniconda3/envs/quantecon/lib/python3.12/site-packages (from pandas>=0.23->pandas-datareader) (2023.3)
Requirement already satisfied: charset-normalizer<4,>=2 in /home/runner/miniconda3/envs/quantecon/lib/python3.12/site-packages (from requests>=2.19.0->pandas-datareader) (3.3.2)
Requirement already satisfied: idna<4,>=2.5 in /home/runner/miniconda3/envs/quantecon/lib/python3.12/site-packages (from requests>=2.19.0->pandas-datareader) (3.7)
Requirement already satisfied: urllib3<3,>=1.21.1 in /home/runner/miniconda3/envs/quantecon/lib/python3.12/site-packages (from requests>=2.19.0->pandas-datareader) (2.2.3)
Requirement already satisfied: certifi>=2017.4.17 in /home/runner/miniconda3/envs/quantecon/lib/python3.12/site-packages (from requests>=2.19.0->pandas-datareader) (2024.8.30)
Requirement already satisfied: six>=1.5 in /home/runner/miniconda3/envs/quantecon/lib/python3.12/site-packages (from python-dateutil>=2.8.2->pandas>=0.23->pandas-datareader) (1.16.0)

import numpy as np
import networkx as nx
import matplotlib.pyplot as plt
import quantecon_book_networks
import quantecon_book_networks.input_output as qbn_io
import quantecon_book_networks.plotting as qbn_plt
import quantecon_book_networks.data as qbn_data
import matplotlib as mpl
from matplotlib.patches import Polygon

quantecon_book_networks.config("matplotlib")
mpl.rcParams.update(mpl.rcParamsDefault)

FONTPATH = "fonts/SourceHanSerifSC-SemiBold.otf"
mpl.font_manager.fontManager.addfont(FONTPATH)
plt.rcParams['font.family'] = ['Source Han Serif SC']

下图展示了从美国经济分析局2021年输入-输出账户数据中的15个生产部门之间的联系网络。

_images/a09cc809a4350b89f92871d5d0fe9e24fc415d8db079c0b0e9e1481709251a91.png — Fig. 39.1 美国15个生产部门之间的联系网络#

标签	行业	标签	行业	标签	行业
ag	农业	wh	批发	pr	专业服务
mi	采矿业	re	零售	ed	教育与健康
ut	公用事业	tr	运输	ar	艺术与娱乐
co	建筑	in	信息	ot	其他服务（不包括政府）
ma	制造业	fi	金融	go	政府

从 $i$ 到 $j$ 的箭头表示 $i$ 行业的一些产出作为 $j$ 行业生产的投入。

经济的特征是存在许多这样的联系。

分析这些联系的基本框架是列昂惕夫的投入产出模型。

在介绍投入产出模型之后，我们将描述它与线性规划讲座的一些联系。

39.2. 投入产出分析#

设

$x_{0}$ 为单一外生生产要素的数量，例如劳动力
$x_{j}, j = 1, \dots n$ 为最终产品 $j$ 的总产出
$d_{j}, j = 1, \dots n$ 为可用于最终消费的最终产品 $j$ 的净产出
$z_{i j}$ 为用作生产产品 $j$ 所投入的产品 $i$ 的数量， $i = 1, \dots n$ ， $j = 1, \dots n$
$z_{0 j}$ 为用于生产产品 $j$ 的劳动力数量
$a_{i j}$ 为生产一单位产品 $j$ 所需的产品 $i$ 的单位数， $i = 0, \dots, n, j = 1, \dots n$
$w > 0$ 为外生劳动力工资，以每单位劳动力的美元计
$p$ 为 $n \times 1$ 的生产品 $i = 1, \dots, n$ 价格向量

生产产品 $j \in {1, \dots, n}$ 的技术由以下列昂惕夫函数描述

x_{j} = min_{i \in {0, \dots, n}} (\frac{z_{i j}}{a_{i j}})

39.2.1. 两种产品#

为说明起见，我们首先设 $n = 2$ 并制定以下网络。

_images/c9c282cdc521360a6f44f1e3ec33143bdb1a166a775e621dbd013479195b960f.png

可行分配必须满足

\begin{aligned} (1 - a_{11}) x_{1} - a_{12} x_{2} & \geq d_{1} \\ - a_{21} x_{1} + (1 - a_{22}) x_{2} & \geq d_{2} \\ a_{01} x_{1} + a_{02} x_{2} & \leq x_{0} \end{aligned}

这可以用图形表示如下。

_images/2dc6782b8763b78563d6ecf7745c8733c4aecfff5466cf1798cb4e29736c172b.png

生产的约束条件是

(39.1)#

\begin{aligned} (I - A) x & \geq d \\ a_{0}^{⊤} x & \leq x_{0} \end{aligned}

其中 $A$ 是一个 $n \times n$ 的矩阵，其每个值为 $a_{i j}$ ，而 $a_{0}^{⊤} = [\begin{matrix} a_{01} & \dots & a_{0 n} \end{matrix}]$ 。

如果我们解 (39.1) 的第一组方程，得到总产出 $x$ ，我们得到

(39.2)#

x = (I - A)^{- 1} d \equiv L d

其中矩阵 $L = (I - A)^{- 1}$ 有时被称为 列昂惕夫逆矩阵。

为确保 (39.2) 的解 $X$ 是一个正向量，以下 霍金斯-西蒙条件 就足够了：

\begin{array}{r} \begin{matrix} det (I - A) > 0 且 \\ (I - A)_{i j} > 0 对所有 i = j \end{matrix} \end{array}

Example 39.1

例如，一个包含两种商品的经济体

(39.3)#

\begin{array}{r} A = [\begin{array}{c} 0.1 & 40 \\ 0.01 & 0 \end{array}] 和 d = [\begin{array}{c} 50 \\ 2 \end{array}] \end{array}

A = np.array([[0.1, 40],
             [0.01, 0]])
d = np.array([50, 2]).reshape((2, 1))

I = np.identity(2)
B = I - A
B

array([[ 9.e-01, -4.e+01],
       [-1.e-02,  1.e+00]])

让我们检查一下霍金斯-西蒙条件

np.linalg.det(B) > 0 # 检查霍金斯-西蒙条件

True

现在我们计算列昂惕夫逆矩阵

L = np.linalg.inv(B) # 得到列昂惕夫逆矩阵
L

array([[2.0e+00, 8.0e+01],
       [2.0e-02, 1.8e+00]])

x = L @ d   # 求解毛产出
x

array([[260. ],
       [  4.6]])

39.3. 生产可能性边界#

(39.1)的第二个等式可以写成

a_{0}^{⊤} x = x_{0}

或

(39.4)#

A_{0}^{⊤} d = x_{0}

其中

A_{0}^{⊤} = a_{0}^{⊤} (I - A)^{- 1}

对于 $i \in {1, \dots, n}$ ， $A_{0}$ 的第 $i$ 个分量是生产一单位第 $i$ 种商品的最终产出所需的劳动量。

等式 (39.4) 描绘了一个生产可能性边界，表示在给定外生劳动投入 $x_{0}$ 的情况下，可以生产的最终消费束 $d$ 。

Example 39.2

考虑(39.3)中的例子。

假设我们现在给出

a_{0}^{⊤} = [\begin{matrix} 4 & 100 \end{matrix}]

然后我们可以通过以下方式找到 $A_{0}^{⊤}$ ：

a0 = np.array([4, 100])
A0 = a0 @ L
A0

array([ 10., 500.])

因此，这个经济体的生产可能性边界是：

10 d_{1} + 500 d_{2} = x_{0}

39.4. 价格#

[Dorfman et al., 1958] 认为， $n$ 种商品的相对价格必须满足

\begin{array}{r} \begin{matrix} p_{1} = a_{11} p_{1} + a_{21} p_{2} + a_{01} w \\ p_{2} = a_{12} p_{1} + a_{22} p_{2} + a_{02} w \end{matrix} \end{array}

更一般地，

p = A^{⊤} p + a_{0} w

这表明每种最终商品的价格等于生产的总成本，包括中间投入品的成本 $A^{⊤} p$ 和劳动力成本 $a_{0} w$ 。

这个方程可以写成

(39.5)#

(I - A^{⊤}) p = a_{0} w

这意味着

p = (I - A^{⊤})^{- 1} a_{0} w

注意 (39.5) 与 (39.1) 通过相互转置的算子形成了一个共轭对。

这种联系在经典线性规划及其对偶问题中再次出现。

39.5. 线性规划#

原始问题是

min_{x} w a_{0}^{⊤} x

约束条件为

(I - A) x \geq d

39.6. 列昂惕夫逆矩阵#

我们已经讨论过，总产出 $x$ 由公式 (39.2) 给出，其中 $L$ 被称为列昂惕夫逆矩阵。

回顾诺伊曼级数引理，它指出如果谱半径 $r (A) < 1$ ，则 $L$ 存在。事实上

L = \sum_{i = 0}^{\infty} A^{i}

39.6.1. 需求冲击#

现在我们考虑需求冲击 $Δ d$ 的影响，它将需求从 $d_{0}$ 转变为 $d_{1} = d_{0} + Δ d$ 。总产出从 $x_{0} = L d_{0}$ 转变为 $x_{1} = L d_{1}$ 。

如果 $r (A) < 1$ ，则存在解，且

Δ x = L Δ d = Δ d + A (Δ d) + A^{2} (Δ d) + \dots

这说明 $L$ 的一个元素 $l_{i j}$ 显示了对商品 $j$ 的需求单位变化对部门 $i$ 的总影响。

39.7. 图论的应用#

我们可以通过应用图论来进一步研究投入产出网络。

投入产出网络可以通过邻接矩阵 $A$ 导出的加权有向图来表示。

节点集 $V = [n]$ 是部门列表，边集由以下给出：

E = {(i, j) \in V \times V : a_{i j} > 0}

在 Fig. 39.1 中，权重由箭头的宽度表示，与相应的投入产出系数成正比。

现在我们可以使用中心性度量来对部门进行排序，并讨论它们相对于其他部门的重要性。

39.7.1. 特征向量中心性#

节点 $i$ 的特征向量中心性由以下公式衡量：

\begin{array}{r} e_{i} = \frac{1}{r (A)} \sum_{1 \leq j \leq n} a_{i j} e_{j} \end{array}

我们为 Fig. 39.1 中表示的部门绘制了基于枢纽的特征向量中心性条形图。

_images/c1866a9f90e53156bf20259ce984cf6a49a80688532625a21404b971388a6f0e.png

较高的指标表示作为供应商的重要性更高。

因此，大多数行业的需求冲击将显著影响具有高特征向量中心性的行业的活动。

上图表明制造业是美国经济中最主导的行业。

39.7.2. 产出乘数#

在投入产出网络中对行业进行排名的另一种方法是通过产出乘数。

行业 $j$ 的产出乘数，记为 $μ_{j}$ ，通常定义为行业 $j$ 需求单位变化所产生的整个行业范围内的总影响。

早些时候在讨论需求冲击时，我们得出结论：对于 $L = (l_{i j})$ ，元素 $l_{i j}$ 表示行业 $j$ 需求单位变化对行业 $i$ 的影响。

因此，

μ_{j} = \sum_{j = 1}^{n} l_{i j}

这可以写成 $μ^{⊤} = 1^{⊤} L$ 或

μ^{⊤} = 1^{⊤} (I - A)^{- 1}

请注意，这里我们用 $1$ 表示一个全为1的向量。

在这个指标中排名较高的行业是中间品的重要购买者。

这些行业的需求冲击将对整个生产网络造成巨大影响。

下图显示了 Fig. 39.1 中表示的各行业的产出乘数。

_images/45a743574eb03ec287795b1400e12436839c80eb15983e262bf4ea566c47e484.png

我们观察到制造业和农业是排名最高的行业。

39.8. 练习#

Exercise 39.1

[Dorfman et al., 1958] 第9章讨论了一个具有以下参数设置的例子：

\begin{array}{r} A = [\begin{array}{c} 0.1 & 1.46 \\ 0.16 & 0.17 \end{array}] 和 a_{0} = [\begin{array}{c} .04 & .33 \end{array}] \end{array}

\begin{array}{r} x = [\begin{array}{c} 250 \\ 120 \end{array}] 和 x_{0} = 50 \end{array}

\begin{array}{r} d = [\begin{array}{c} 50 \\ 60 \end{array}] \end{array}

描述他们如何从以下关于农业和制造业的假设”数据”中推断出 $A$ 和 $a_{0}$ 中的投入-产出系数：

\begin{array}{r} z = [\begin{array}{c} 25 & 175 \\ 40 & 20 \end{array}] 和 z_{0} = [\begin{array}{c} 10 & 40 \end{array}] \end{array}

其中 $z_{0}$ 是每个行业使用的劳动服务的向量。

Solution to Exercise 39.1

For each $i = 0, 1, 2$ and $j = 1, 2$

a_{i j} = \frac{z_{i j}}{x_{j}}

Exercise 39.2

推导上一练习中所描述经济的生产可能性边界。

Solution to Exercise 39.2

A = np.array([[0.1, 1.46],
              [0.16, 0.17]])
a_0 = np.array([0.04, 0.33])

I = np.identity(2)
B = I - A
L = np.linalg.inv(B)

A_0 = a_0 @ L
A_0

array([0.16751071, 0.69224776])

因此，生产可能性边界由以下给出

0.17 d_{1} + 0.69 d_{2} = 50