家庭生存臨界動力學：微觀-宏觀耦合經濟系統的相變理論（完整版）

**家庭生存臨界動力學：微觀-**宏觀耦合經濟系統的相變理論（完整版）

Household Survival Criticality Dynamics: Phase Transition Theory of Micro-Macro Coupled Economic Systems

作者: Neo.K & Theia
機構: EveMissLab
日期: 2026年4月
版本: v2.0-Comprehensive

§0 核心命題

定理0（生存臨界原理）： $\boxed{ \text{系統崩潰} \Leftrightarrow \rho\left(\frac{\partial Q}{\partial W} \leq 0\right) > \rho_c \land \text{SR}_{\text{med}} < \text{SR}c \land \frac{\text{CFI}}{W{1\%}} > \tau_{\text{逃離}} }$

三重判據：

努力無效密度超過臨界值
購物車生存比率跌破閾值
資本外逃率超過恐慌線

第一部分：理論架構（§1-6）

§1 微觀層：單一家庭動力學

1.1 狀態向量

定義1.1： $\boxed{ \mathbf{H}_i(t) = \begin{pmatrix} I_i(t) \ E_i(t) \ S_i(t) \ W_i(t) \ Q_i(t) \end{pmatrix}, \quad Q_i = \frac{P_i \cdot T_i}{\Psi_i} \cdot \tanh\left(\lambda \frac{S_i}{E_i}\right) }$

其中：

<![if !msEquation]><![if !vml]><![endif]><![endif]>：實際購買力
<![if !msEquation]><![if !vml]><![endif]><![endif]>：自由時間（h/月）
<![if !msEquation]><![if !vml]><![endif]><![endif]>：壓力指數
<![if !msEquation]><![if !vml]><![endif]><![endif]>

1.2 努力效率

定義1.2： $\boxed{ \eta_i(t) = \frac{\partial Q_i}{\partial W_i} = \frac{\partial}{\partial W_i}\left[\frac{I_i(W_i) \cdot (720-W_i)}{\Psi_i(W_i) \cdot \text{CPI}} \cdot \phi(S_i/E_i)\right] }$

臨界家庭： $\boxed{ \mathcal{F}_c(t) = \left\{i ,|, \eta_i(t) \leq 0 \land \Delta t > \tau_{\text{心理}} \sim 6\text{-}12\text{月}\right\} }$

1.3 動力學方程

方程1.1： $\boxed{ \begin{align} \frac{dI_i}{dt} &= \alpha_i W_i - \beta_i U(t) - \gamma_i \cdot \text{Age}_i \\ \frac{dS_i}{dt} &= I_i - E_i - \delta_i |S_i|_- \\ \frac{dW_i}{dt} &= \kappa_i \left(\frac{E_i - I_i}{I_i}\right)_+ - \mu_i(W_i - W_0) \\ \frac{d\Psi_i}{dt} &= \nu_i \left[\left(\frac{E_i}{I_i}-1\right)_+ + \frac{|S_i|_-}{E_i} + \sigma_{\text{市場}}\right] \end{align} }$

§2 中觀層：網絡耦合

2.1 網絡結構

定義2.1： $\boxed{ \mathcal{G} = (\mathcal{V}, \mathcal{E}, \mathbf{A}), \quad A_{ij} = \begin{cases} w_{ij}^{\text{就業}}, & i \text{ 僱 } j \\ w_{ij}^{\text{借貸}}, & i \text{ 借 } j \\ w_{ij}^{\text{供需}}, & i \text{ 消費 } j \end{cases} }$

2.2 級聯動力學

定義2.2（節點狀態）： $\boxed{ \sigma_i(t+\Delta t) = \mathbb{1}\left[S_i(t) + \sum_{j \in \mathcal{N}_i} A_{ij} \sigma_j(t) \Delta S_{ij} > -\theta_{\text{破產}}\right] }$

雪崩分佈： $\boxed{ P(s) \sim s^{-\tau} e^{-s/s_c}, \quad \tau \approx 1.5 }$

§3 宏觀層：聚合動力學

3.1 序參量定義

定義3.1（核心序參量）： $\boxed{ \Phi(t) = 1 - \frac{1}{N}\sum_i \mathbb{1}_{\eta_i \leq 0} = \frac{|{i ,|, \eta_i > 0}|}{N} }$

統計量： $\boxed{ \begin{align} I_{\text{med}}(t) &= \text{Median}{I_i} \\ Q_{\text{med}}(t) &= \text{Median}{Q_i} \\ S_{\text{bias}}(t) &= \frac{\text{Mean}{Q_i}}{\text{Median}{Q_i}} \cdot \frac{\text{Mean}{I_i}}{\text{Median}{I_i}} \end{align} }$

3.2 與貨幣流動性統一

定義3.2： $\boxed{ \begin{align} M_A(t) &= \sum_i M_{A,i} = \sum_i \max(I_i - E_i^{\text{必需}}, 0) \\ V_A(t) &= \frac{\sum_i E_i^{\text{必需}}}{M_A} \\ R_A(t) &= \frac{M_A}{M_{\text{總}}} \\ \text{DMR}(t) &= \frac{D_{\text{總}}}{M_{\text{總}}} \end{align} }$

§4 臨界條件與相變

4.1 朗道自由能

方程4.1： $\boxed{ \mathcal{L}[\Phi, T] = a(T - T_c)\Phi^2 + b\Phi^4 + c\Phi^6 - h\Phi }$

平衡條件： $\boxed{ \frac{\delta \mathcal{L}}{\delta \Phi} = 0 \Rightarrow 2a(T-T_c)\Phi + 4b\Phi^3 + 6c\Phi^5 = h }$

溫度定義： $\boxed{ T(t) = w_1 \cdot \text{DMR}(t) + w_2 \cdot (0.4 - R_A(t))_+ + w_3 \cdot U(t) }$

4.2 臨界指數

定理4.1（標度律）： $\boxed{ \begin{align} \Phi &\sim (T_c - T)^\beta, \quad \beta \approx 0.5 \\ \chi &= \frac{\partial \Phi}{\partial h} \sim |T - T_c|^{-\gamma}, \quad \gamma \approx 1.0 \\ \xi &\sim |T - T_c|^{-\nu}, \quad \nu \approx 0.6 \end{align} }$

4.3 早期預警

定義4.2（臨界減速）： $\boxed{ \text{CSD}(t) = \frac{\text{Var}[\Phi(t-\tau:t)]}{\text{Mean}[\Phi(t-\tau:t)]}, \quad \tau \sim 6\text{月} }$

預警閾值：

<![if !msEquation]><![if !vml]><![endif]><![endif]><![if !supportLineBreakNewLine]>
<![endif]>

§5 統一動力學方程

方程5.1（主方程）： $\boxed{ \frac{d\Phi}{dt} = -2a(T - T_c)\Phi - 4b\Phi^3 - 6c\Phi^5 + h + \sqrt{2D},\xi(t) }$

耦合系統： $\boxed{ \begin{align} \frac{d\Phi}{dt} &= f_\Phi(\Phi, R_A, V_A, \text{DMR}) + \eta_1(t) \\ \frac{dR_A}{dt} &= f_{R_A}(R_A, \Phi, h_{\text{政策}}) + \eta_2(t) \\ \frac{dV_A}{dt} &= f_{V_A}(V_A, \Phi, R_A) + \eta_3(t) \\ \frac{d(\text{DMR})}{dt} &= f_{\text{DMR}}(\text{DMR}, R_A, \Phi) + \eta_4(t) \end{align} }$

§6 六層完備性檢驗

層級映射： $\boxed{ \begin{align} E[\text{HSCT}] &= {\mathbf{H}_i}_{i=1}^N \cup {\mathbf{A}, \rho, T, h} \in \mathbb{R}^{5N+N^2+3} \\ C[\text{HSCT}] &= (\Phi, R_A, V_A, \text{DMR}, \text{CSD}) \in \mathbb{R}^5 \\ N[\text{HSCT}] &: \Phi \sim (T_c - T)^\beta \\ P[\text{HSCT}] &= {\Phi(t_0) \to \Phi(t_1) \to \cdots} \\ M[\text{HSCT}] &\approx 0.88 \quad \text{（與CDMS、NCAT高耦合）} \\ S[\text{HSCT}] &\approx 0.75 \quad \text{（Level 2湧現）} \end{align} }$

第二部分：雙重溫度計（§7）

§7 購物車指數與資本外逃

7.1 購物車生存指數

定義7.1（標準購物車）： $\boxed{ \text{BCI}(t) = \sum_{i \in \mathcal{B}} p_i(t) q_i^{\text{標準}} }$

**標準籃（週/**家庭）： $\boxed{ \mathcal{B} = \begin{cases} \text{主食} & 5\text{kg} \\ \text{蛋白質} & 2\text{kg肉} + 10\text{蛋} + 2\text{L奶} \\ \text{蔬果} & 3\text{kg蔬} + 2\text{kg果} \\ \text{油鹽} & 0.5\text{L油} + \text{調味} \\ \text{日用} & \text{衛生紙、洗衣粉} \end{cases} }$

生存比率： $\boxed{ \text{SR}_i(t) = \frac{I_i^{\text{週}}(t)}{\text{BCI}(t)}, \quad \text{SR}_c = 3.5 }$

與序參量映射： $\boxed{ \Phi(t) \approx \frac{1}{N}\sum_i \mathbb{1}_{\text{SR}_i > 3.5} }$

7.2 避稅天堂資本流入

定義7.2： $\boxed{ \text{CFI}_s(t) = \frac{d}{dt}\left[\sum_{h \in \mathcal{H}} D_{h \leftarrow s}(t)\right] }$

避稅天堂集： $\boxed{ \mathcal{H} = {\text{瑞士, 開曼, BVI, 盧森堡, 新加坡, 巴拿馬, \ldots}} }$

逃離率： $\boxed{ r_{\text{逃離}}(t) = \frac{\text{CFI}_s(t)}{W_{1\%,s}(t)} }$

臨界閾值： $\boxed{ \begin{align} r < 5\%/\text{年} &: \text{正常} \\ 5\% < r < 15\% &: \text{警告} \\ r > 15\% &: \text{恐慌} \end{align} }$

7.3 系統分裂溫度

定義7.3： $\boxed{ T_{\text{split}}(t) = \frac{r_{\text{逃離}}(t)}{|\Delta \text{SR}_{\text{med}}(t)|/\text{SR}_{\text{med}}(t)} }$

解讀： $\boxed{ \begin{align} T_{\text{split}} \gg 1 &: \text{富人先逃（領先指標）} \\ T_{\text{split}} \approx 1 &: \text{同步崩潰} \\ T_{\text{split}} \ll 1 &: \text{底層先爆} \end{align} }$

7.4 綜合判據

定理7.1（三重臨界條件）： $\boxed{ \begin{align} &\text{條件A}: \text{SR}_{\text{med}} < 2.5 \\ &\text{\_\_條件B}: r_{\text{逃離}} > 15\% \\ &\text{條件C}: T_{\text{split}} > 5 \\ &\Rightarrow P(\text{崩潰}|12\text{月}) > 0.75 \end{align} }$

第三部分：國際校準（§8）

§8 歷史危機參數反推

8.1 校準數據集

定義8.1： $\boxed{ \mathcal{D} = {(c, \mathbf{X}_c(t), t_{\text{崩}}^{(c)}, \Delta_c) ,|, c \in \mathcal{C}_{\text{危機}}} }$

樣本（<![if !msEquation]><![if !vml]><![endif]><![endif]>）：

高收入：希臘2010, 冰島2008, 西班牙2012

中等收入：阿根廷2001, 土耳其2018, 巴西2015

災難級：委內瑞拉2016, 津巴布韋2008

8.2 宏觀代理變量

映射8.1： $\boxed{ \begin{align} \Phi_{\text{宏}} &\approx 1 - \frac{U + P_{\text{貧困}}}{2} \\ R_A &\approx \frac{M1 - M0}{M2} \\ \text{DMR} &= \frac{D_{\text{政府}} + D_{\text{企業}} + D_{\text{家庭}}}{M2} \\ S_{\text{bias}} &\approx \frac{\text{GDP}_{\text{\_\_人均}}}{\text{\_\_收入}_{\text{中位數}}} \end{align} }$

8.3 參數優化

問題8.1： $\boxed{ \min_{\boldsymbol{\theta}} \sum_{c \in \mathcal{C}} \left[\left(t_{\text{預測}}^{(c)} - t_{\text{實際}}^{(c)}\right)^2 + \lambda|\Delta_{\text{預測}}^{(c)} - \Delta_{\text{實際}}^{(c)}|^2\right] }$

結果（<![if !msEquation]><![if !vml]><![endif]><![endif]>）：

第四部分：跨國風險地圖（§9）

§9 匿名風險空間

9.1 國家類型學

定義9.1： $\boxed{ \text{Type} = f(\text{收入}, \text{政體}, \text{貨幣地位}, \text{資源}) }$

分類：

Type-I：高收入 + 民主 + 儲備貨幣

特徵：高DMR可承受

代表：美國、英國

Type-II：高收入 + 民主 + 非儲備

特徵：DMR臨界、無特權

代表：南歐部分國家

Type-III：中高收入 + 威權 + 資本管制

特徵：高DMR、低R_A、高S_bias

風險：結構性高危

Type-IV：中等收入 + 大宗商品依賴

特徵：DMR波動、外生衝擊敏感

代表：拉美部分國家

9.2 相空間投影

投影9.1： $\boxed{ \begin{align} x &= \text{DMR} - \text{DMR}_c \\ y &= R_A - R_{A,c} \end{align} }$

相圖分區： $\boxed{ \begin{array}{c|c} \text{象限} & \text{狀態} \\ \hline (x>0, y>0) & \text{危險區} \\ (x<0, y>0) & \text{安全區} \\ (x<0, y<0) & \text{脆弱區} \\ (x>0, y<0) & \text{特權區} \end{array} }$

9.3 匿名標註（2024 Q2）

不點名，用參數：

點A (x=1.2, y=-0.15)：

DMR = 3.68, R_A = 0.23

SR_med = 3.2, CFI/W_1% = 1.8%

CSD = 2.1, T_split = 4.5

→ 風險：高（Pr_崩潰,24月 = 0.78）

點B (x=-0.5, y=0.25)：

DMR = 1.98, R_A = 0.63

SR_med = 6.1, CFI/W_1% = 0.3%

→ 風險：低

點C (x=0.8, y=0.05)：

DMR = 3.28, R_A = 0.43

SR_med = 4.2, CFI/W_1% = 1.2%

CSD = 1.5

→ 風險：中-高（Pr = 0.52）

9.4 軌跡演化

情景分析（點A）： $\boxed{ \begin{align} &\text{基線}(h=0): \quad \Phi(24\text{月}) = 0.25, , P_{\text{崩}} = 78\% \\ &\text{中等}(h=0.3): \quad \Phi(24\text{月}) = 0.35, , P_{\text{崩}} = 65\% \\ &\text{激進}(h=0.7): \quad \Phi(24\text{月}) = 0.65, , P_{\text{崩}} = 28\% \end{align} }$

第五部分：政策空間（§10）

§10 干預策略與最優控制

10.1 政策向量

定義10.1： $\boxed{ \mathbf{h} = (h_1, h_2, h_3, h_4, h_5)^T = (\text{財富稅}, \text{UBI}, \text{債務重組}, \text{資本管制}, \text{QE}) }$

效果矩陣（實證）： $\boxed{ \frac{\partial \mathbf{S}}{\partial \mathbf{h}} = \begin{pmatrix} 0.3 & 0.5 & 0.2 & 0.1 & -0.1 \ 0.2 & 0.4 & 0.1 & 0.3 & -0.3 \ -0.1 & 0.05 & -0.3 & 0 & 0.4 \ \vdots & \vdots & \vdots & \vdots & \vdots \end{pmatrix} }$

列：<![if !msEquation]><![if !vml]><![endif]><![endif]>
行：<![if !msEquation]><![if !vml]><![endif]><![endif]>

10.2 最優控制

問題10.1： $\boxed{ \min_{\mathbf{h}(t)} \int_0^T \left[(1-\Phi)^2 + \lambda|\mathbf{h}|^2\right] dt }$

約束： $\boxed{ \begin{align} \frac{d\mathbf{S}}{dt} &= \mathbf{F}(\mathbf{S}, \mathbf{h}) \\ \mathbf{h}_{\min} \leq \mathbf{h} &\leq \mathbf{h}_{\max} \\ \sum_i c_i h_i &\leq B_{\text{財政}} \end{align} }$

數值解（點A****類型）： $\boxed{ \mathbf{h}^* = \begin{pmatrix} 0.6 \ 0.8 \ 0.4 \ 0.2 \ 0 \end{pmatrix}, \quad \text{成本} \approx 9\% , \text{GDP/年} }$

效果：

<![if !msEquation]><![if !vml]><![endif]><![endif]><![if !supportLineBreakNewLine]>
<![endif]>

第六部分：Type-III**詳細分析（§11）**

§11 匿名高風險案例

11.1 參數特徵（2024 Q2）

不點名數據： $\boxed{ \begin{align} \text{DMR} &= 3.2\text{-}3.8 , (\gg \text{DMR}_c = 2.48) \\ R_A &= 0.20\text{-}0.28 , (< R_{A,c} = 0.38) \\ S_{\text{bias}} &= 2.8\text{-}3.5 , (\text{極高}) \\ \text{SR}_{\text{med}} &= 3.2 , (< \text{SR}_c = 3.5) \\ \text{CSD} &= 1.8\text{-}2.3 , (\text{\_\_上升}) \\ \text{CFI}/W_{1\%} &= 1.5\text{-}2.0\% , (\text{加速中}) \end{align} }$

11.2 購物車詳細分析

標準購物車成本演化： $\boxed{ \begin{align} \text{BCI}(2020) &= 220 , \text{貨幣單位} \\ \text{BCI}(2024) &= 350 , (+59\%) \\ I_{\text{med}}(2020) &= 1100/\text{週} \\ I_{\text{med}}(2024) &= 1200/\text{週} , (+9\%) \\ \text{SR}(2020) &= 5.0 \\ \text{SR}(2024) &= 3.4 , (\text{跌破臨界}) \end{align} }$

分佈分析： $\boxed{ \begin{align} P(\text{SR} < 2.5) &\approx 35\% , \text{（2024）} \\ P(\text{SR} < 3.5) &\approx 58\% \\ P(\text{SR} > 5) &\approx 18\% , (\text{僅頂層}) \end{align} }$

11.3 資本外逃加速

時間序列： $\boxed{ \begin{align} \text{CFI}(2022) &= 650\text{億USD} \\ \text{CFI}(2023) &= 850\text{億USD} , (+31\%) \\ \text{CFI}(2024, \text{年化}) &= 1400\text{億USD} , (+65\%) \end{align} }$

逃離率演化：

<![if !msEquation]><![if !vml]><![endif]><![endif]><![if !supportLineBreakNewLine]>
<![endif]>

11.4 系統診斷

綜合指標： $\boxed{ \begin{align} \Phi_{\text{估}} &\approx 0.38\text{-}0.42 , (\text{臨界邊緣}) \\ T &= 1.5 \cdot \text{DMR} + 3.0 \cdot (0.4-R_A) + 0.5 \cdot U \\ &\approx 5.1 + 0.36 + 0.06 = 5.52 \\ T_c &\approx 4.8 \\ T - T_c &\approx 0.72 , (\text{超臨界}) \end{align} }$

預警信號矩陣： $\boxed{ \begin{array}{l|c|c} \text{指標} & \text{當前值} & \text{狀態} \\ \hline \Phi & 0.40 & \color{red}{\text{✗ 臨界}} \\ \text{SR}_{\text{med}} & 3.2 & \color{red}{\text{\_\_✗\_ \_低於閾值}} \\ \text{DMR} & 3.5 & \color{red}{\text{\_\_✗\_ \_嚴重超標}} \\ R_A & 0.24 & \color{red}{\text{\_\_✗\_ \_流動性危機}} \\ \text{CSD} & 2.1 & \color{red}{\text{\_\_✗\_ \_臨界減速}} \\ r_{\text{逃離}} & 1.75\% & \color{orange}{\text{⚠ 警告}} \\ T_{\text{split}} & 4.5 & \color{red}{\text{✗ 高度分裂}} \end{array} }$

11.5 軌跡預測

無干預情景： $\boxed{ \begin{align} t = 6\text{月} &: \Phi = 0.33, , \text{CSD} = 2.5 \\ t = 12\text{月} &: \Phi = 0.28, , \text{SR}_{\text{med}} = 2.8 \\ t = 18\text{\_\_月} &: \Phi = 0.22, , r_{\text{逃離}} > 3\% \\ t = 24\text{月} &: P(\text{進入崩潰相}) = 78\% \end{align} }$

政策窗口： $\boxed{ \tau_{\text{有效}} \leq 12\text{月}, \quad h_{\min} \geq 0.6 }$

超過此窗口，所需政策強度指數增長：

<![if !msEquation]><![if !vml]><![endif]><![endif]><![if !supportLineBreakNewLine]>
<![endif]>

第七部分：理論總結（§12）

§12 形式化命題鏈

12.1 核心定理

定理12.1（相變必然性）： $\boxed{ \begin{align} &\text{給定} , T > T_c \land h < h_c, \\ &\exists , t^* < \infty: \Phi(t^_) < \Phi_c \land \frac{d\Phi}{dt}\bigg|_{t^_} < 0 \\ &\Rightarrow \lim_{t \to \infty} \Phi(t) = 0 \end{align} }$

定理12.2（購物車判據）： $\boxed{ \text{SR}_{\text{med}} < 2.5 \Rightarrow \Phi < 0.4 , \text{（概率95\%）} }$

定理12.3（資本外逃領先性）： $\boxed{ \frac{dr_{\text{逃離}}}{dt} > 0.5\%/\text{年}^2 \Rightarrow P(\text{崩潰}|18\text{月}) > 0.6 }$

12.2 統一泛函

作用量： $\boxed{ \mathcal{S} = \int dt , d\mathbf{H} \left[\sum_i \mathcal{L}_i[\mathbf{H}i] + \sum{i<j} \mathcal{L}_{ij}[A_{ij}] + \mathcal{L}_{\text{宏}}[\Phi, R_A, \text{DMR}]\right] }$

極值原理：

<![if !msEquation]><![if !vml]><![endif]><![endif]><![if !supportLineBreakNewLine]>
<![endif]>

12.3 時空尺度

時間尺度分離： $\boxed{ \begin{align} \tau_{\text{微觀}} &\sim 1\text{月（家庭決策）} \\ \tau_{\text{中觀}} &\sim 6\text{月（網絡傳播）} \\ \tau_{\text{宏觀}} &\sim 18\text{月（相變完成）} \\ \tau_{\text{預警}} &\sim 24\text{月（CSD到崩潰）} \end{align} }$

空間尺度： $\boxed{ \begin{align} N_{\text{微觀}} &\sim 10^7\text{家庭} \\ \langle k \rangle_{\text{網絡}} &\sim 10\text{-}50\text{（平均度）} \\ \xi_{\text{相關}} &\sim 10^3\text{-}10^4\text{（臨界時）} \end{align} }$

附錄A：API****規範

A.1 核心接口

python

class HSCT_Analyzer:

"""家庭生存臨界動力學分析器"""

def init(self, country_type: str):

self.type = country_type

self.params = self.load_calibrated_params()

def compute_state(self, data: dict) -> dict:

"""計算系統狀態向量

Parameters:

data : dict

{'dmr': float, 'r_a': float,

'unemployment': float, 'gini': float,

'bci': float, 'income_med': float,

'cfi': float, 'w_1pct': float}

Returns:

state : dict

{'Phi': float, 'SR_med': float,

'CSD': float, 'T': float, 'T_split': float}

"""

pass

def predict_crisis(self, horizon: int = 24) -> dict:

"""預測危機概率

Returns:

{'prob_collapse': float,

'time_to_critical': int (months),

'trajectory': array}

"""

pass

def policy_simulation(self, policies: dict) -> dict:

"""政策情景模擬"""

pass

A.2 數據源

yaml

required_data_sources:

macroeconomic:

World Bank API: GDP, CPI, GINI
OECD: unemployment, household_debt
IMF: money_supply, credit_data
BIS: cross_border_flows

micro_proxies:

Numbeo: cost_of_living (BCI代理)
National Statistics: income_distribution

capital_flight:

BIS: banking_statistics
Tax Havens: deposit_flows (部分公開)

附錄B：參數表

校準結果（基於15個歷史危機）：

參數

物理意義

估計值

95% CI