Logo

线性系统的应用

1. 经济学中的列昂惕夫模型

1.1 问题背景

在复杂的经济体系中,各生产部门通过相互交换资源实现协作。例如,煤炭部门需要电力支持生产,电力部门消耗煤炭发电,钢铁部门则依赖两者的资源。这种相互依赖关系要求收入与支出达到平衡,即每个部门的总收入(由其产出的价格决定)必须等于购买其他部门资源的总支出。列昂惕夫的“投入-产出”模型通过线性方程组将这种平衡关系数学化,从而求解出使所有部门收支相等的均衡价格。

1.2 线性方程组的作用

通过经济交换表,我们可以将各部门的产出分配比例转化为线性方程,构建齐次系统。通过求解该系统,可找到一组价格向量,使得每个部门的价格(收入)等于其购买其他部门产出的总成本(支出)。由于系统是齐次的(方程右侧全为零),解需通过自由变量参数化,体现价格比例关系而非绝对值。

1.3 示例

 1 ~1~:假设一个经济体由煤炭、电力、钢铁个部门组成,各部门的产出分配比例如下表所示:

产出分配来源煤炭电力钢铁
煤炭部门00.40.6
电力部门0.60.10.2
钢铁部门0.40.50.2

对于表格含义,以表格第二行为例说明:电力部门生成的电力 60% ~60\%~分配给了煤炭部门、10% 10\%~用于自身运营、20% 20\%~分配给了钢铁部门。目标:求均衡价格 pC, pE, pS ~p_C,~p_E,~p_S~,使得每个部门的总收入等于其总支出。

  • 建立线性方程组

    • 煤炭部门

      • 设收入: pC ~p_C~

      • 支出:购买电力部门 40% ~40\%~的产出( 0.4pE ~0.4p_E~)和钢铁部门 60% ~60\%~的产出( 0.6pS ~0.6p_S~)。

      • 方程:
        pC=0.4pE+0.6pSpC0.4pE0.6pS=0\begin{align*} & p_C = 0.4p_E + 0.6p_S\\[2ex] \Rightarrow \quad &p_C - 0.4p_E - 0.6p_S = 0\end{align*}
    • 电力部门

      • 设收入: pE ~p_E~

      • 支出:购买煤炭部门 60% ~60\%~的产出( 0.6pC ~0.6p_C~)、自身 10% ~10\%~的产出( 0.1pE ~0.1p_E~)、钢铁部门 20% ~20\%~的产出( 0.2pS ~0.2p_S~)。

      • 方程:
        pE=0.6pC+0.1pE+0.2pS0.6pC+0.9pE0.2pS=0\begin{align*} & p_E = 0.6p_C + 0.1p_E + 0.2p_S\\[2ex] \Rightarrow \quad &-0.6p_C + 0.9p_E - 0.2p_S = 0\end{align*}
    • 钢铁部门

      • 设收入: pS ~p_S~

      • 支出:购买煤炭部门 40% ~40\%~的产出( 0.4pC ~0.4p_C~)、电力部门 50% ~50\%~的产出( 0.5pE ~0.5p_E~)、自身 20% ~20\%~的产出( 0.2pS ~0.2p_S~)。

      • 方程:
        pS=0.4pC+0.5pE+0.2pS0.4pC0.5pE+0.8pS=0\begin{align*} & p_S = 0.4p_C + 0.5p_E + 0.2p_S\\[2ex] \Rightarrow \quad &-0.4p_C - 0.5p_E + 0.8p_S = 0\end{align*}
  • 行化简求解

    • 线性方程组:
      {pC0.4pE0.6pS=00.6pC+0.9pE0.2pS=00.4pC0.5pE+0.8pS=0\left\{ \begin{array}{l} p_C - 0.4p_E - 0.6p_S = 0 \\[2ex] -0.6p_C + 0.9p_E - 0.2p_S = 0 \\[2ex] -0.4p_C - 0.5p_E + 0.8p_S = 0 \end{array} \right.
    • 对应的增广矩阵:
      [10.40.600.60.90.200.40.50.80]\left[ \begin{array}{ccc|c} 1 & -0.4 & -0.6 & 0 \\ -0.6 & 0.9 & -0.2 & 0 \\ -0.4 & -0.5 & 0.8 & 0 \end{array} \right]
    • 化简得最终行阶梯形:
      [100.940010.8500000]\left[ \begin{array}{ccc|c} 1 & 0 & -0.94 & 0 \\ 0 & 1 & -0.85 & 0 \\ 0 & 0 & 0 & 0 \end{array} \right]
      解的参数形式:
      {pC=0.94pSpE=0.85pSpS为自由变量\left\{ \begin{array}{l} p_C = 0.94p_S \\[2ex] p_E = 0.85p_S \\[2ex] p_S \quad \footnotesize{为自由变量} \end{array} \right.
      参数化向量:
      P=[pCpEpS]=[0.94pS0.85pSpS]=[0.940.851]pS\mathbf{P} =\begin{bmatrix}p_C \\ p_E \\ p_S\end{bmatrix} = \begin{bmatrix} 0.94p_S \\ 0.85p_S \\ p_S \end{bmatrix} = \begin{bmatrix} 0.94 \\ 0.85 \\ 1 \end{bmatrix}p_S

参数化的均衡价格向量揭示了经济部门间的相对价值关系,钢铁部门作为基准价格,煤炭和电力部门的价格分别为其 94% ~94\%~ 85% ~85\%~,这种比例关系确保了各部门的收支平衡。

2. 化学方程式的平衡

 2 ~2~:平衡以下丙烷燃烧的化学反应式,要求系数为最小整数:
(x1)C3H8+(x2)O2(x3)CO2+(x4)H2O(x_1)C_3H_8 + (x_2)O_2 \rightarrow (x_3)CO_2 + (x_4)H_2O
求解过程如下:

  1. 定义原子守恒方程

    • 根据反应物和产物中碳、氢、氧的原子数守恒,建立向量方程:

      • 碳原子守恒:3x1=x33x_1 = x_3

      • 氢原子守恒:8x1=2x48x_1 = 2x_4

      • 氧原子守恒:2x2=2x3+x42x_2 = 2x_3 + x_4

      将方程整理为齐次线性方程组:
      {3x1x3=08x12x4=02x22x3x4=0\begin{cases} 3x_1 - x_3 = 0 \\ 8x_1 - 2x_4 = 0 \\ 2x_2 - 2x_3 - x_4 = 0 \end{cases}
  2. 构造增广矩阵并化简

    将方程组表示为增广矩阵:
    [301008002002210]\left[ \begin{array}{cccc|c} 3 & 0 & -1 & 0 & 0 \\ 8 & 0 & 0 & -2 & 0 \\ 0 & 2 & -2 & -1 & 0 \end{array} \right]
    行化简后得:
    [101/4000105/400013/40]\left[ \begin{array}{cccc|c} 1 & 0 & -1/4 & 0 & 0\\ 0 & 1 & 0 & -5/4 & 0\\ 0 & 0 & 1 & -3/4 & 0 \end{array} \right]
    解的参数形式:
    {x1=14x4x2=54x4x3=34x4x4为自由变量\left\{ \begin{array}{l} x_1 = \frac{1}{4}x_4 \\[2ex] x_2 = \frac{5}{4}x_4 \\[2ex] x_3 = \frac{3}{4}x_4 \\[2ex] x_4 \quad \footnotesize{为自由变量} \end{array} \right.
    参数化向量:
    x=[1/45/43/41]x4\mathbf{x} = \begin{bmatrix} 1/4 \\ 5/4 \\ 3/4 \\ 1 \end{bmatrix} x_4
     x4=4 ~x_4 = 4~得最小整数解:
    C3H8+5O23CO2+4H2OC_3H_8 + 5O_2 \rightarrow 3CO_2 + 4H_2O

3. 网络流分析

 3 ~3~:某城市的交通网络如图所示,包含四个交叉口(节点 A,B,C,D ~A,B,C,D~)和五条单向道路(分支)。已知以下流量数据:
  • 从外部进入网络的流量:A A~点流入 200 ~200~辆/小时,D D~点流入 300 ~300~辆/小时。

  • 从网络流出的流量:B B~点流出 150 ~150~辆/小时,C C~点流出 250 ~250~辆/小时。

  • 各分支的流量方向及未知流量标记为 x1,x2,x3,x4,x5 ~x_1,x_2,x_3,x_4,x_5~

求解:
  1. 建立流量守恒方程组
    交叉口流入 (Flow in)流出 (Flow out)
    A300+500300 + 500x1+x2x_1 + x_2
    Bx2+x4x_2 + x_4300+x3300 + x_3
    C100+400100 + 400x4+x5x_4 + x_5
    Dx1+x5x_1 + x_5600600
  2. 构造增广矩阵并通过行化简求解各分支流量

    • 流量守恒方程组为:
      x1+x2+0x3+0x4+0x5=8000x1+x2x3+x4+0x5=3000x1+0x2+0x3+x4+x5=500x1+0x2+0x3+0x4+x5=6000x1+0x2+x3+0x4+0x5=400\begin{alignat*}{6} x_1 & {}+{} & x_2 & \textcolor{dddddd}{{}+{}} & \textcolor{dddddd}{0x_3} & \textcolor{dddddd}{{}+{}} & \textcolor{dddddd}{0x_4} & \textcolor{dddddd}{{}+{}} & \textcolor{dddddd}{0x_5} & \hspace{1ex}{}={}\hspace{1ex} & 800 \\ \textcolor{dddddd}{0x_1} & \textcolor{dddddd}{{}+{}} & x_2 & {}-{} & x_3 & {}+{} & x_4 & \textcolor{dddddd}{{}+{}} & \textcolor{dddddd}{0x_5} & \hspace{1ex}{}={}\hspace{1ex} & 300 \\ \textcolor{dddddd}{0x_1} & \textcolor{dddddd}{{}+{}} & \textcolor{dddddd}{0x_2} & \textcolor{dddddd}{{}+{}} & \textcolor{dddddd}{0x_3} & \textcolor{dddddd}{{}+{}} & x_4 & {}+{} & x_5 & \hspace{1ex}{}={}\hspace{1ex} & 500 \\ x_1 & \textcolor{dddddd}{{}+{}} & \textcolor{dddddd}{0x_2} & \textcolor{dddddd}{{}+{}} & \textcolor{dddddd}{0x_3} & \textcolor{dddddd}{{}+{}} & \textcolor{dddddd}{0x_4} & \textcolor{dddddd}{{}+{}} & x_5 & \hspace{1ex}{}={}\hspace{1ex} & 600 \\ \textcolor{dddddd}{0x_1} & \textcolor{dddddd}{{}+{}} & \textcolor{dddddd}{0x_2} & \textcolor{dddddd}{{}+{}} & x_3 & \textcolor{dddddd}{{}+{}} & \textcolor{dddddd}{0x_4} & \textcolor{dddddd}{{}+{}} & \textcolor{dddddd}{0x_5} & \hspace{1ex}{}={}\hspace{1ex} & 400 \end{alignat*}
    • 增广矩阵为:
      [1100080001110300000115001000160000100400][10001600010012000010040000011500000000]\left[\begin{array}{ccccc|c} 1 & 1 & 0 & 0 & 0 & 800 \\ 0 & 1 & -1 & 1 & 0 & 300 \\ 0 & 0 & 0 & 1 & 1 & 500 \\ 1 & 0 & 0 & 0 & 1 & 600 \\ 0 & 0 & 1 & 0 & 0 & 400 \end{array}\right] \sim \left[\begin{array}{ccccc|c} 1 & 0 & 0 & 0 & 1 & 600 \\ 0 & 1 & 0 & 0 & -1 & 200 \\ 0 & 0 & 1 & 0 & 0 & 400 \\ 0 & 0 & 0 & 1 & 1 & 500 \\ 0 & 0 & 0 & 0 & 0 & 0 \end{array}\right]
    • 解的参数形式:
      {x1=600x5x2=200+x5x3=400x4=500x5x5为自由变量\left\{ \begin{array}{l} x_1 = 600 - x_5 \\[2ex] x_2 = 200+ x_5 \\[2ex] x_3 = 400 \\[2ex] x_4 = 500 - x_5 \\[2ex] x_5 \quad \footnotesize{为自由变量} \end{array} \right.
  3. 确定自由变量的取值范围,确保所有流量非负
    • 非负性约束(因街道为单行道,流量不可为负):
      {x10x5600x20x5200x40x5500(自动满足)\begin{cases} x_1 \geq 0 & \Rightarrow x_5 \leq 600 \\[2ex] x_2 \geq 0 & \Rightarrow x_5 \geq -200 \\[2ex] x_4 \geq 0 & \Rightarrow x_5 \leq 500 \end{cases} \quad (\text{自动满足})
      综合得: 0x5500 ~0 \leq x_5 \leq 500~,这对应实际交通流量的合理分配。
    • 例如,当自由变量 x5 ~x_5~接近其上限 500 ~500~时,得
      x1=100,x2=700,x4=0x_1 = 100, \quad x_2 = 700, \quad x_4 = 0
      从交叉口 C ~C~ D ~D~的分支流量 x4 ~x_4~趋近于零,而交叉口 A ~A~ B ~B~的分支流量 x2 ~x_2~显著增加至 700 ~700~辆/小时,可能导致路段 B ~B~拥堵,同时反映出路段 C ~C~ D ~D~未被充分利用,需调整交通管理策略以优化路网效率。