JOURNAL OF ECONOMIC LITERATURE · COMPUTATIONAL ECONOMICS

Deep Learning for Solving Economic Models

深度学习求解经济学模型:理论范式、算法架构与启示评述

Jesús Fernández‑Villaverde | Journal of Economic Literature 64(3), 2026, 829–875

VIDEO EXPLAINER · 中文视频讲解

Deep Learning 中文有声讲解

7 分 2 秒 · 854×480 压缩视频

AUDIO EXPLAINER · 中文音频讲解

深度学习破解经济学维数灾难

18 分 27 秒 · 压缩音频

论文英文标题 Deep Learning for Solving Economic Models
论文作者 Jesús Fernández-Villaverde (University of Pennsylvania, NBER, and CEPR)
发表期刊与卷期 Journal of Economic Literature 2026, 64(3), pp. 829–875
文献DOI https://doi.org/10.1257/jel.20261794
JEL 分类号 C45, C61, D83, O41

一、论文 JEL 分类号及其学术内涵解析

本文在《Journal of Economic Literature》(JEL)正式发表,论文标注了四个标准的JEL分类代码。这些分类代码精准地界定了论文横跨现代计算经济学、数理规划、信息经济学与宏观经济增长的理论交叉属性:

C45: Econometric and Statistical Methods: Artificial Intelligence / Neural Networks / Machine Learning(计量经济与统计方法:人工智能、神经网络与机器学习)

【学科对应】:本类别代表人工智能、神经网络与现代机器学习在经济学中的方法论研究。在本文中,对应将多层前向神经网络、反向传播算法、随机梯度下降(SGD)、Adam优化器以及隐式正则化理论系统化应用于经济学函数逼近的核心议题。

C61: Mathematical Methods; Programming Models; Dynamic Programming(数学方法;规划模型;动态规划)

【学科对应】:本类别涵盖动态规划、贝尔曼方程(Bellman Equation)、连续时间哈密顿-雅可比-贝尔曼(HJB)方程、最优控制理论及KKT优化条件的数值解法。本文的核心正是解决高维状态空间下动态规划的“维度灾难”,建立离散与连续时间下的非线性动态规划深度求解范式。

D83: Search, Learning, Information and Knowledge; Communication; Belief; Unawareness(搜寻、学习、信息与知识;沟通;信念;未觉知)

【学科对应】:本类别关注经济主体的信息结构、学习过程以及算法对未知函数的“学习(Learning)”。在本文中,一方面体现为神经网络对未知策略函数和价值函数的自适应表示学习(Representation Learning);另一方面体现为经济主体在不确定性环境下形成理性预期与主观信念的数值逼近。

O41: One, Two, and Multisector Growth Models(一部门、两部门及多部门宏观经济增长模型)

【学科对应】:本类别属于宏观经济学中的新古典增长理论(Neoclassical Growth Model)。本文通篇以确定性与随机新古典增长模型、连续时间新古典增长模型为基准示范平台,展示深度学习如何求解资本积累、跨期消费欧拉方程及稳态遍历分布。

二、论文的核心理论与方法论贡献

Jesús Fernández-Villaverde 教授的这篇长篇权威综述,不仅是对近年来深度学习在经济学计算领域前沿进展的系统梳理,更是对定量宏观经济学与微观结构模型计算范式的一场深刻重塑。其核心学术贡献体现在以下五个维度:

1. 突破传统数值方法的“维度灾难”,彻底打开高维经济模型的求解空间

自20世纪70年代理性预期革命以来,动态一般均衡模型大多缺乏闭式解析解。传统全局求解方法(如基于网格的值函数迭代VFI、基于切比雪夫多项式或有限元的正交投影法)在状态变量超过4至5维时即遭遇指数级算力爆炸(若每维取100个网格点,d维需100^d个点);而局部扰动法(Perturbation)虽计算迅速,但在处理偶发绑定约束(如零下限ZLB、流动性约束)与极端尾部灾难时严重失真。本文系统性论证并确立了深度学习作为破解高维状态空间(数十至数百维甚至更高)非线性动态模型全局解的基石地位,使经济学家不再因算力限制而被迫在模型中施加不符合现实的过度简化假设。

2. 阐明“自适应学习基函数(Learned Bases)”超越“先验固定基函数(Fixed Bases)”的本质机理

传统投影法要求研究者在见到具体问题前先验设定正交基函数(如Chebyshev多项式张量积),在高维空间中面临基函数组合爆炸和严重的多重共线性。深度学习通过多层复合架构,利用反向传播将复杂非线性全部“封装”在简单的单变量激活函数(如ReLU、Tanh)中,使得输入特征被自动投影、扭曲和折叠,自发生成一个信息紧凑的最优隐空间基底(Hidden Units)。在经过训练的隐空间中,高度复杂的未知经济学目标函数转变为隐单元的近乎线性组合,从根本上兼顾了强大的非线性表达力与极佳的外推泛化能力。

3. 首次清晰形式化并提出破解经济学独特难题——“均衡循环(Equilibrium Loop)”的自洽算法

计算机科学中的常规深度学习依赖外生固定的数据集与标签进行离线监督学习。然而,动态经济模型的求解存在本质不同的内在循环:训练网络所需的评估点应当来自经济体自身的遍历分布(Ergodic Distribution),但该遍历分布本身又由经济主体尚未求解出的最优策略函数内生决定!这一“先有鸡还是先有蛋”的悖论被作者高度概括为“均衡循环”。论文提出了基于轮次(Epoch)的策略模拟与自洽重采样迭代机制,使得算法在训练过程中动态自适应聚焦于系统实际漫游的经济有效状态区域,彻底摆脱了高维超立方体无效空间的网格算力浪费。

4. 构建涵盖离散/连续时间、多控制变量与非线性互补约束的大一统计算框架

论文超越了简单的无约束规划,构建了一个具有极高通用性的动态规划大一统计算框架:无论是离散时间的随机贝尔曼方程,还是连续时间的Hamilton-Jacobi-Bellman(HJB)方程;无论是等式约束、状态转移律,还是带有偶发绑定约束(如借贷受限、名义利率零下限)的不等式约束,均可被统一编码为多输出神经网络系统(同时拟合价值函数、决策规则及KKT拉格朗日乘子)。论文通过构建包含主可行性残差、对偶可行性残差和互补松弛残差的加权复合损失函数,给出了完备的优化表述。

5. 规范现代深度学习工程生态的落地标准与稳健实践指南

结合NVIDIA GPU高度并行化硬件生态与PyTorch/JAX自动微分软件生态,论文通过连续时间新古典增长模型的实际案例,给出了诸如:将消费参数化为Sigmoid产出比例以天然保证预算可行与非负性、利用Softplus保证边际效用严格正向、针对HJB残差进行绝对值动态归一化、以及运用梯度裁剪(Gradient Clipping)与余弦退火学习率调度等一系列高价值工程技巧。在普通笔记本电脑上仅需1分钟即可实现欧拉方程误差低于0.1%的超高精度全局解,确立了经济学模型深度学习求解的可重复性与透明度标杆。

三、核心概念、数学定义与模型符号详析

为了深入理解深度学习如何切入动态经济学求解,本节对论文中涉及的全部数学对象、系统符号、结构定义与算法概念进行模块化分解与详细解释。

1. 基础特征表示与目标函数映射

在经济学模型的函数逼近设定中,目标是在特征空间上寻找一个未知的连续或离散映射。定义输入特征向量(即模型的状态变量):

X=(1,x1,x2,...,xN)∈X⊆ℝN+1X = (1,\ x_{1},\ x_{2},\ ...,\ x_{N}) \in X \subseteq {\mathbb{R}}^{N + 1}

其中常数项 11 被形式化纳入状态向量首位,用于将仿射变换中的常数偏置项(Bias)直接吸收进矩阵乘法中。未知目标函数记为 y=f(X)y = f(X),其中 y∈ℝy \in {\mathbb{R}} 为标量经济学变量(如价值函数值、最优消费水平、资产价格)。若目标为向量值(如包含多个资产的投资组合策略向量),则通过堆叠多个标量输出构成多头网络(Multi-head Network)。

在随机动态博弈或含观测噪声的场景中,实际观测到的目标值满足:

y=f(X)+ϵ,𝔼[ϵ|X]=0y = f(X) + \epsilon,\ \ \ \ {\mathbb{E}}\lbrack\epsilon|X\rbrack = 0

此时 ϵ\epsilon 反映了对手混合策略的随机实现或不可观测的状态扰动,但期望收益函数 f(X)f(X) 本身依然是确定性的映射。

2. 经典多项式基函数与切比雪夫正交投影

传统数值经济学通过先验设定的多项式基底构建参数化近似族:

(1)线性近似模型:仅保留一阶项,无法刻画动态规划中由于边际效用递减所必然导致的价值函数严格凹性:

f(X)≈gθ(X)=θ0+∑n=1Nθnxnf(X) \approx g_{\theta}(X) = \theta_{0} + \sum_{n = 1}^{N}\theta_{n}x_{n}

(2)二次单项式展开:引入二次平方项与双变量交叉乘积项,但高阶单项式在有界闭区间上高度共线,导致法方程矩阵严重病态(Ill-conditioned):

gθQA(X)=θ0+∑n=1Nθn1xn+∑n=1Nθn2xn2+∑n=1N−1∑m=n+1Nθn,mxnxmg_{\theta}^{QA}(X) = \theta_{0} + \sum_{n = 1}^{N}\theta_{n}^{1}x_{n} + \sum_{n = 1}^{N}\theta_{n}^{2}x_{n}^{2} + \sum_{n = 1}^{N - 1}\sum_{m = n + 1}^{N}\theta_{n,m}x_{n}x_{m}

(3)第一类切比雪夫多项式(Chebyshev Polynomials of the First Kind)投影:利用在权重 (1−x2)−1/2\left( 1 - x^{2} \right)^{- 1/2} 下的正交性改善条件数:

gθCP(X)=θ0+∑m∈ℳθm∏j=1NTmj(xj)g_{\theta}^{CP}(X) = \theta_{0} + \sum_{m \in \mathcal{M}}\theta_{m}\prod_{j = 1}^{N}T_{m_{j}}\left( x_{j} \right)

其中 m=(m1,...,mN)m = (m_{1},\ ...,\ m_{N}) 为多重指标(Multi-index)。若采用张量积集 ℳ={0,...,D}N\mathcal{M} = \{ 0,\ ...,\ D\}^{N},项数达到 (D+1)N(D + 1)^{N};若采用完全度集(Total-degree Set)ℳ={m∈ℕ0N:||m||1≤D}\mathcal{M} = \{ m \in {\mathbb{N}}_{0}^{N}:\ ||m||_{1} \leq D\},项数仍随 NN 呈组合级剧烈膨胀。

3. 神经网络结构:从单隐层到深度网络架构

神经网络采取完全不同的函数构造逻辑。单隐层神经网络(Single-hidden-layer Neural Network)通过三步实现输入到输出的映射:

【步骤1:预激活计算(Preactivations)】输入特征通过权重矩阵 W∈ℝM×(N+1)W \in {\mathbb{R}}^{M \times (N + 1)} 进行仿射线性投影,得到 MM 维预激活向量 aa:

a=WX,am=wm,0+∑n=1Nwm,nxna = WX,\ \ \ \ a_{m} = w_{m,0} + \sum_{n = 1}^{N}w_{m,n}x_{n}

其中 WW 的首列 wm,0w_{m,0} 为偏置(Bias),MM 为网络宽度(Width)。该操作本质是简单的矩阵乘法,高度适配现代GPU的并行计算流。

【步骤2:隐单元与非线性封装(Hidden Units / Neurons)】对预激活向量逐元素应用非线性激活函数 ϕ(⋅)\phi( \cdot ),生成隐单元向量 z∈ℝMz \in {\mathbb{R}}^{M}:

z=ϕ(a)=ϕ(WX)z = \phi(a) = \phi(WX)

zz 所处的空间被称为【表示空间(Representation Space)】。此处体现了深度学习的核心思想——【封装非线性(Encapsulated Nonlinearity)】:整个模型中所有的非线性弯曲完全被锁入单变量标量激活函数 ϕ(⋅)\phi( \cdot ) 之中,矩阵 WW 负责旋转、平移和缩放坐标轴,自发为目标函数构造最贴合的自适应基底。

【步骤3:仿射输出层(Linear Output Layer)】对隐空间基底进行最终的线性组合:

gθNN(X)=v0+v⊤z=(r∘ϕ∘W)(X)g_{\theta}^{NN}(X) = v_{0} + v^{\top}z = (r \circ \phi \circ W)(X)

其中 r(z)=v0+v⊤zr(z) = v_{0} + v^{\top}z 为仿射输出映射,所有网络权重收集在元组 θ=(v0,v,W)\theta = (v_{0},\ v,\ W) 中。

【深度神经网络(Deep Neural Networks, DNN)的推广】:若包含 JJ 个隐层,前向传播(Forward Propagation)形成多层复合嵌套映射:

z1=ϕ(W1X),zj+1=ϕ(Wj+1zj),gθDL(X)=(r∘ϕJ∘WJ∘...∘ϕ1∘W1)(X)z^{1} = \phi(W^{1}X),\ \ \ \ z^{j + 1} = \phi(W^{j + 1}z^{j}),\ \ \ \ g_{\theta}^{DL}(X) = (r \circ \phi^{J} \circ W^{J} \circ ... \circ \phi^{1} \circ W^{1})(X)

深层网络允许较低层抽取原始状态的局部相关性,较高层将低层特征递归复合为复杂的全局几何表征,实现特征的层次化抽象与极致压缩。

4. 激活函数谱系(Activation Functions)

激活函数必须具备计算廉价性与高效导数求取特性。论文中系统比较了以下六类激活函数,其数学定义与经济学特征如下:

• 线性恒等函数:ϕ(z)=z\phi(z) = z。退化为经典多元线性回归,验证了神经网络是对传统计量模型的直接推广。

• Sigmoid 逻辑斯蒂函数:ϕ(z)=11+e−z\phi(z) = \frac{1}{1 + e^{- z}}。值域位于 (0,1)(0,\ 1),常用于拟合离散选择概率、政策份额或将经济学变量约束在有界正区间内。

• 双曲正切函数:tanh(z)=e2z−1e2z+1tanh(z) = \frac{e^{2z} - 1}{e^{2z} + 1}。值域为 (−1,1)( - 1,\ 1),零中心化且平滑无限阶可微。适合已知目标函数高度平滑的连续时间动态规划模型(如HJB方程涉及价值函数一阶及高阶导数计算)。

• 线性整流单元:ReLU(z)=max(0,z)ReLU(z) = max(0,\ z)。分段线性函数,导数仅为 0 或 1,计算极其迅速,彻底规避深层反向传播中的梯度消失问题。具备【局部自适应性(Local Adaptability)】,在非平滑、存在断点或偶发绑定约束的模型中表现优异。

• 平滑增广函数:Softplus(z)=log(1+ez)Softplus(z) = log(1 + e^{z})。ReLU 的处处光滑可微近似,不仅保持梯度在正轴上的顺畅流动,且能严格保证输出恒大于零,常用于强加价值函数单调递增约束 V′(k)>0V\prime(k) > 0。

• Softmax 多项概率映射:πk(X)=esk(X)∑j=1Kesj(X)\pi_{k}(X) = \frac{e^{s_{k}(X)}}{\sum_{j = 1}^{K}e^{s_{j}(X)}}。将实数向量 logits 映射到离散概率单纯形上,用于拟合多主体博弈中的混合策略选择概率。

5. 经济学残差系统与损失函数设计

深度学习求解经济模型,本质是将泛函方程算子 𝒯(d)=0\mathcal{T}(d) = 0 转化为评估点上的残差平方和最小化问题:

(1)贝尔曼方程残差(Bellman Equation Residual):直接基于动态规划的递推自洽条件:

ℛ[VθV](X)=VθV(X)−max{u(c)+β𝔼VθV(X′)}\mathcal{R}\lbrack V_{\theta}^{V}\rbrack(X) = V_{\theta}^{V}(X)\ - \ max\{ u(c)\ + \ \beta\ {\mathbb{E}}\ V_{\theta}^{V}(X\prime)\}

(2)欧拉方程决策残差(Euler Equation Residual):当市场均衡非帕累托最优,或直接求解一阶条件比求解值函数更便捷时,拟合决策规则 ct=d(Xt)c_{t} = d(X_{t}):

ℛ[dθ](Xt)=u′(dθ(Xt))−β𝔼t{u′(dθ(Xt+1))(1+αezt+1kt+1α−1−δ)}\mathcal{R}\lbrack d_{\theta}\rbrack(X_{t}) = u\prime(d_{\theta}(X_{t}))\ - \ \beta\ {\mathbb{E}}_{t}\{ u\prime(d_{\theta}(X_{t + 1}))(1\ + \ \alpha\ e^{z_{t + 1}}k_{t + 1}^{\alpha - 1}\ - \ \delta)\}

(3)参数化期望残差(PEA Residual):当政策函数在临界点存在硬折角(如利率零下限ZLB),直接拟合平滑的条件期望项 d(Xt)=𝔼t[⋅]d(X_{t}) = {\mathbb{E}}_{t}\lbrack \cdot \rbrack,避免非光滑政策函数的近似困难。进而通过 ct=(u′)−1(βd(Xt))c_{t} = (u\prime)^{- 1}(\beta d(X_{t})) 倒求消费。

(4)KKT 系统复合约束残差:在存在不等式约束 G(X,α)≤0G(X,\ \alpha) \leq 0 与等式约束 H(X,α)=0H(X,\ \alpha) = 0 时,构建包含四大误差分量的复合损失,并均严格使用原生数学公式表征:

• 主可行性残差(Primal Feasibility Residuals):

errPF1,θ(X)=max{0,G(X,α̂θα(X))},errPF2,θ(X)=H(X,α̂θα(X)){err}_{PF1,\theta}(X)\ = \ max\{ 0,\ G(X,\ {\alpha\hat{}}_{\theta}^{\alpha}(X))\},\ \ \ \ \ \ \ \ {err}_{PF2,\theta}(X)\ = \ H(X,\ {\alpha\hat{}}_{\theta}^{\alpha}(X))

• 对偶可行性残差(Dual Feasibility Residual):保证不等式约束拉格朗日乘子非负:

errDF,θ(X)=max{0,−μ̂θμ(X)}{err}_{DF,\theta}(X)\ = \ max\{ 0,\ - {\mu\hat{}}_{\theta}^{\mu}(X)\}

• 互补松弛残差(Complementary Slackness Residual):确保乘子与约束松弛度互补:

errCS,θ(X)=(μ̂θμ(X))⊤G(X,α̂θα(X)){err}_{CS,\theta}(X)\ = \ ({\mu\hat{}}_{\theta}^{\mu}{(X))}^{\top}G(X,\ {\alpha\hat{}}_{\theta}^{\alpha}(X))

• 单评估点加权复合损失函数(Per-evaluation-point Loss):

ll(θ)=∑i=16ωi||erri,θ(Xl)||22l_{l}(\theta)\ = \ \sum_{i = 1}^{6}\omega_{i}||{err}_{i,\theta}(X_{l}{)||}_{2}^{2}

通过稳态基准值对各方程残差实施无量纲归一化,平衡不同经济学量纲在多目标优化中的梯度贡献。

四、核心理论命题与推论的直觉解释与生动案例

本节梳理支撑深度学习求解经济学模型的五大核心数学命题与理论推论。严格按照要求,不展开复杂的数理证明,而是深入浅出地剖析其经济与数学内涵,将命题中的数学符号与算例均转换为Word原生公式呈现,并配备生动具体的经济学直觉案例。

命题一:通用近似定理(Universal Approximation Theorem)

【文献溯源】:Cybenko (1989), Hornik, Stinchcombe and White (1989), Kidger and Lyons (2020)

【命题核心内涵】:对于定义在有限维实数紧集上的任意 Borel 可测函数 f(X)f(X)(即便该函数不处处可微,或在局部存在离散跳跃断点),只要隐含层包含足够多的神经元,单隐层前向神经网络即可在任意预先指定的逼近精度 ϵ>0\epsilon > 0 下一致逼近该目标函数。

【直觉解释与生动案例】:

• 直觉比喻:用标准模块化乐高积木拼装任意复杂几何雕塑。论文第4.3节给出了一个极其生动的算例:目标函数为 f(x)=x3+x2−x−1f(x) = x^{3} + x^{2} - x - 1,在区间 [−2,1.5]\lbrack - 2,\ 1.5\rbrack 上包含局部极大、局部极小及拐点。如果只用 1 个 ReLU 神经元,网络拟合函数为 gθ1(x)=−max(0,−7.7−5x)g_{\theta}^{1}(x) = - max(0,\ - 7.7 - 5x),只能拟合一段平截折线;当神经元逐步增加到 6 个时,每个 ReLU 单元凭借其自身的阈值偏移与斜率变换,像 6 段互相咬合的铰链活动臂,在各自激活的区间内独立调整,而在钝化区域输出严格为 0,互不干扰,最终几乎无缝重合了原目标曲线(见论文图5)。

• 经济学应用场景:考虑存在流动性约束或借贷上限的家庭生命周期消费问题。在家庭资产触及借贷边界(at=0a_{t} = 0)的一瞬间,消费函数会出现明显的尖锐拐折(Kink Point),边际消费倾向出现结构性跃升。传统高阶平滑切比雪夫多项式在遇到这种尖锐拐折时,会出现灾难性的“龙格现象(Runge Phenomenon)”——为了拟合局部折角,多项式会在整个平滑区域引起剧烈且虚假的波浪状振荡。相反,神经网络的ReLU单元可以恰好将一个神经元的折点锚定在借贷边界上:边界左侧输出为零,边界右侧线性增长,不仅在折点处实现了绝对平稳的拟合,而且丝毫不会污染其他健康资产区间的平滑消费决策。

命题二:巴伦近似误差界定理与破除维度灾难(Barron's Theorem on Approximation Bounds)

【文献溯源】:Barron (1993), Jacot, Choi, and Wen (2024)

【命题核心内涵】:对于属于巴伦函数类(Barron Class,即其傅里叶变换具备良好的一阶谱衰减性质,这是绝大多数经济学连续决策函数的自然属性)的目标函数,单隐层神经网络的均方积分逼近误差收敛速率为 𝒪(1/M)\mathcal{O}(1/M)(MM 为隐单元神经元数量)。最关键的是:该收敛速度完全独立于输入状态变量的维度 NN!相比之下,传统 NN 维正交级数/切比雪夫多项式逼近的误差收敛速率为 𝒪(D−2s/N)\mathcal{O}(D^{- 2s/N})(DD 为基函数项数,ss 为目标函数的光滑阶数)。

【直觉解释与生动案例】:

• 直觉比喻:“网格盲目地毯式排查”对比“智能雷达自适应锁定关键方位”。传统多项式方法为了保证逼近精度,必须在 NN 维超立方体空间的每个维度上均匀布设节点;由于维度 NN 赫然出现在收敛阶数的指数分母上(−2s/N- 2s/N),维度每增加一维,为了维持相同精度所需的节点数就要呈几何倍数暴增,这就是经典的“维度灾难”。

• 经济学应用场景:假设某主权基金需要求解包含 100 种金融资产(N=100N = 100)的跨期动态投资组合模型。如果采用切比雪夫多项式网格法,即使每个资产维度仅粗糙地取 3 个离散网格点,总网格数量也高达 3100≈5.15×10473^{100} \approx 5.15 \times 10^{47} 个,这远超地球上所有海滩沙粒的总和,没有任何超级计算机能够存储和运算。然而,现实金融世界中 100 种资产的价格联动态势,绝大多数由几个共同宏观因子(如基准利率、经济增长率、全要素生产率)所支配。巴伦定理证明:神经网络根本不需要在 100 维空间的无边荒原上地毯式铺点,而是通过学习隐层权重矩阵 WW,自动把 100 维高维状态压缩投影到几组最关键的线性组合方向上。在实践中,只需将隐单元数量从 50 个增加到 500 个,均方误差就能呈 1/M1/M 的确定性速率平稳衰减,成功粉碎了高维状态空间对收敛速度的指数诅咒!

命题三:深度网络的指数级表达能力定理(Expressive Power of Deep vs. Shallow Networks)

【文献溯源】:Goujon, Etemadi, and Unser (2024), Poggio et al. (2017)

【命题核心内涵】:采用分段线性激活函数(如ReLU)的深度神经网络,其能够对输入空间划分的独立线性超平面区域(Linear Regions)数量随网络层数(深度 JJ)呈【指数级增长】,而网络的总参数量和计算复杂度仅随层数呈【线性增长】。与浅层单隐层网络相比,深度架构在逼近具备组合结构或层级对称性的函数时展现出无可比拟的高效性。

【直觉解释与生动案例】:

• 直觉比喻:神奇的折纸实验。拿一张平整的白纸代表输入状态空间。浅层网络相当于拿着一把剪刀,每增加一个神经元就在纸上剪出一道笔直的切口,区域划分能力随神经元数量线性增加;而深度网络则是先将白纸对折(第1层),再对折(第2层),连续对折 JJ 次。此时只需用剪刀轻轻剪一刀,展开后就会在纸面上瞬间呈现出 2J2^{J} 个错综复杂的多边形划分区域!

• 经济学应用场景:考虑一个包含复杂产业链分工的宏观生产网络模型。最终消费品依赖汽车、电子、机械等若干部门,而每个部门又依赖上游的钢铁、塑料、稀土,再向上追溯至采矿与电力能源,构成一个典型的跨行业深层嵌套生产函数。浅层网络如果想拟合这一巨型生产网络,必须强行把所有层级的状态一次性拍平,需要动用海量的神经元才能勉强捕捉不同行业间的非线性替代弹性;而深度神经网络结构天生契合这种现实经济生产的层级抽象:第1层提取能源与原材料对初级零部件的非线性约束,第2层将零部件融合成功能模块,第3层拟合整车装配,层层递进、特征复用,以极少的参数规模实现了对巨型复杂生产网络政策函数的高保真刻画。

命题四:超参数化与梯度下降的隐式正则化命题(Implicit Regularization and Minimum-Norm Interpolation)

【文献溯源】:Wilson (2025), Soudry et al. (2018), Du et al. (2019)

【命题核心内涵】:当神经网络的参数量远大于训练评估点数量(dim(θ)≫Ldim(\theta) \gg L),系统跨过“插值阈值(Interpolation Threshold)”进入高度过参数化区间时,尽管存在无穷多个能够实现零训练损失 ℒL(θ)=0\mathcal{L}_{L}(\theta) = 0 的参数解,高维梯度下降类优化器(如SGD、Adam)会自动偏向收敛至泛函半范数 ||gθNN||ψ||g_{\theta}^{NN}||_{\psi}(如Sobolev半范数,衡量函数各阶导数平方积分的曲率平滑度)最小的那个特解。即梯度算法自发偏好“最平滑、最自然”的曲面。

【直觉解释与生动案例】:

• 直觉比喻:穿过若干指定钉子的橡皮筋。传统计量经济学教科书时刻告诫学者“自由度陷阱”——参数一旦多于样本,回归就会严重过拟合,在数据点之间疯狂震荡。但深度学习打破了这一旧教条。在论文第6.3节图7中,作者用12个数据点拟合一条复杂曲线:34个参数时严重欠拟合;91个参数时刚刚实现数据点拟合;而当参数激增到301个(严重过参数化)时,测试集上的均方根误差(RMSE)反而从0.108大幅断崖式下降至0.048!曲线在未见过的测试区间内展现出极其优雅顺畅的拟合轮廓。

• 经济学应用场景:在求解随机动态一般均衡模型的资本积累路径时,研究者通常只能生成有限长度的蒙特卡洛模拟轨迹作为训练集。如果数值方法出现伪震荡,会导致资本在模拟点之间剧烈起伏,甚至触发资本爆炸式积累或发散,直接违背经济学中的横截性条件(Transversality Condition, TVC: limT→∞βTcT−1kT+1=0\lim_{T \rightarrow \infty}\beta^{T}c_{T}^{- 1}k_{T + 1} = 0)。得益于隐式正则化机制,梯度下降算法在超参数化空间中寻找解时,天然倾向于惩罚剧烈的导数变化,在数据点之间架起曲率最小的平滑曲面,使得求得的宏观政策函数天然收敛于有界的稳定鞍点路径,甚至无需人工硬性编写复杂的横截性边界惩罚项!

命题五:均衡循环与遍历测度重采样收敛机制(The Equilibrium Loop and Iterative Self-Consistent Sampling)

【文献溯源】:Fernández-Villaverde et al. (2020), Maliar, Maliar, and Winant (2021)

【命题核心内涵】:在动态随机经济模型求解中,理论最优的损失函数应当是基于经济系统内生未知遍历分布 ξ*(X)\xi^{*}(X) 加权的期望残差;通过在每个训练时期(Epoch)利用当前网络蕴含的最优决策规则向前模拟生成新的状态样本,并交替执行“模拟采样—网络重训”,该迭代更新过程能够使得训练点集动态逼近真实的遍历测度支撑集,实现计算精度的自洽收敛。

【直觉解释与生动案例】:

• 直觉比喻:“以车修路,车路协同”的动态闭环。在荒无人烟的未开辟区域,传统网格法就像不分青红皂白在所有悬崖、森林、沼泽上每隔1米铺设一根水泥桩,结果99.9%的水泥桩汽车一辈子也不会路过;但如果不修路,又不知道汽车真正开起来会在哪里形成车流。均衡循环算法的做法是:先基于经验常识规划一条简化的泥土小路(比如用一阶扰动法求出的线性化政策规则);让车在泥土路上跑一圈(模拟采样),记录汽车轮胎压出来的车辙印(生成状态样本);针对这些被频繁碾压的路线集中加铺高标准柏油路面(训练神经网络减少欧拉方程残差);路修得越平顺,汽车跑得越稳健,又引导汽车进一步探索更广泛的安全驾驶边界。如此轮番迭代,车流最终完全稳定在平滑的高速路网上。

• 经济学应用场景:央行在评估包含偶发金融危机和流动性冻结的宏观模型时,整个经济状态空间极其辽阔(涵盖所有极端的失业率、恶性通胀和天文数字杠杆)。但经济体在95%的正常年份都平稳运行在稳态附近的狭小遍历流形上,只有极少数年份会在外生大冲击下短暂滑入危机走廊。如果采用静态固定网格,算力将被浩瀚且永远无法到达的无意义状态淹没;而均衡循环机制使神经网络从第一轮开始,就动态追随模拟经济体真实的繁荣—衰退路径,把宝贵的算力100%倾注在稳态常态区以及由冲击内生触发的金融危机轨迹上,实现了计算资源的最优化倾斜配置。

五、连续时间新古典增长模型求解范式与数值表现

为了给经济学者提供一个直观、透明且高度可复现的完整操作样板,论文第9节以连续时间确定性新古典增长模型(Continuous-Time Deterministic Neoclassical Growth Model)为例,全流程展示了深度学习如何实现高精度全局求解。

1. 经济学模型与 HJB 方程系统构建

社会规划者在无限跨期视野内最大化代表性家庭的折现贴现效用流:

max∫0∞e−ρtu(c)dt,u(c)=c1−γ1−γmax\ \int_{0}^{\infty}e^{- \rho t}u(c)dt,\ \ \ \ u(c) = \frac{c^{1 - \gamma}}{1 - \gamma}

资本积累状态转移微分方程约束为:k̇=Akα−δk−ck\dot{}\ = \ Ak^{\alpha}\ - \ \delta k\ - \ c,且资本满足非负下界 k>0k > 0。对应的连续时间 Hamilton-Jacobi-Bellman (HJB) 泛函方程为:

ρV(k)=max{u(c)+V′(k)[Akα−δk−c]}\rho V(k)\ = \ max\ \{ u(c)\ + \ V\prime(k)\lbrack A\ k^{\alpha}\ - \ \delta k\ - \ c\rbrack\}

联立求解出消费一阶条件:u′(c)=V′(k)⇒c=(u′)−1(V′(k))u\prime(c) = V\prime(k) \Rightarrow c = (u\prime)^{- 1}\left( V\prime(k) \right)。作者在此设定两个独立的深度神经网络分别逼近价值函数 V̂θV(k){V\hat{}}_{\theta}^{V}(k) 与消费政策函数 ĉθc(k){c\hat{}}_{\theta}^{c}(k)。训练目标由两个核心残差构成,均严格以原生公式呈现:

• HJB 泛函残差:

errHJB,θ(k)=ρV̂θV(k)−u(ĉθc(k))−∂V̂θV(k)∂k[Akα−δk−ĉθc(k)]{err}_{HJB,\theta}(k)\ = \ \rho{V\hat{}}_{\theta}^{V}(k)\ - \ u({c\hat{}}_{\theta}^{c}(k))\ - \ \frac{\partial\ {V\hat{}}_{\theta}^{V}(k)}{\partial k}\left\lbrack Ak^{\alpha}\ - \ \delta k\ - \ {c\hat{}}_{\theta}^{c}(k) \right\rbrack

• 一阶决策残差:

errc,θ(k)=(u′)−1(∂V̂θV(k)∂k)−ĉθc(k){err}_{c,\theta}(k)\ = \ (u\prime)^{- 1}\left( \frac{\partial\ {V\hat{}}_{\theta}^{V}(k)}{\partial k} \right)\ - \ {c\hat{}}_{\theta}^{c}(k)

2. 架构设计细节与数值调优工程实践

论文详细披露了14条兼具理论美感与工程实效的调优规范,为学术界提供了极佳的范式参考:

(结构参数标准校准):设相对风险厌恶系数 γ=2\gamma = 2,时间贴现率 ρ=0.04\rho = 0.04,生产率水平 A=0.5A = 0.5,资本产出弹性 α=0.36\alpha = 0.36,资本折旧率 δ=0.05\delta = 0.05。解析求得模型内生稳态资本为 kss=2.95k_{ss} = 2.95。

(网络拓扑与隐藏单元):选用 3 个隐含层(Depth = 3),每层包含 64 个神经元(Width = 64)。实践表明,“层数稍多、每层神经元适中”的深而窄架构,在计算效率与泛化性上普遍优于“层数极少但极宽”的浅层架构。

(激活函数严格匹配光滑性):全连接隐层统一采用双曲正切激活函数 tanh(x)tanh(x)。由于HJB残差必须通过自动微分精确求取一阶导数 V′(k)V\prime(k),而分段线性的ReLU的二阶导数处处为0,光滑的 tanh\tanh 能确保高阶梯度顺畅流动。

(状态定义域的经济学裁剪):将资本评估区间严格限定在 [k_,k‾]=[0.1,10]\lbrack k_{\_},\ k\bar{}\rbrack = \lbrack 0.1,\ 10\rbrack,分别对应稳态资本的 3.4% 与 339%,既确保了超宽的宏观全局探索范围,又排除了没有经济学意义的虚假数值区域。

(物理先验约束参数化设计):将消费决策强制参数化为产出比例形态:c(k)=σ(net(k))⋅Y(k)c(k) = \sigma(net(k)) \cdot Y(k),其中 σ(⋅)\sigma( \cdot ) 为 Sigmoid 函数,Y(k)=AkαY(k) = Ak^{\alpha} 为总产出。该结构在数学上恒等保证消费 c∈(0,Y(k))c \in (0,\ Y(k)),杜绝了数值迭代中可能出现的负消费或超额透支,并在低资本极端区域彻底避免了梯度消失。

(严格施加边际效用正向公理):对网络输出的原始导数实施 Softplus 变换:V′(k)=log(1+ex)>0V\prime(k) = log(1 + e^{x}) > 0,在架构层面原生内嵌了“效用随资本单调递增”的经济学公理。

(HJB 残差动态自适应归一化):在计算损失时,将HJB残差除以动态因子 1+|ρV(k)|1 + |\rho V(k)|。由于低资本区域效用绝对值 |u(c)||u(c)| 与 |V(k)||V(k)| 极大,若不归一化,优化器会片面专注于拟合极低资本点,导致常态区域精度严重恶化。

(创新混合采样分布 Mixture Sampling):每个轮次抽取 128 个训练点,由三部分独立分布混合构成:(1)以稳态 kssk_{ss} 为中心、标准差为 1 的正态分布,强化稳态邻域精度;(2)在 [0.1,10]\lbrack 0.1,\ 10\rbrack 上的均匀分布,确保全局覆盖;(3)对数均匀分布(Log-uniform),使 30% 的样本点落在高曲率的低资本区间 [0.1,1]\lbrack 0.1,\ 1\rbrack,集中攻克非线性最强区域。

(自适应优化与退火调度):采用 Adam 优化算法,初始学习率设为 10−410^{- 4},配合余弦退火(Cosine Annealing)学习率调度机制,在训练末期将步长优雅收缩至趋近于 0;同时施加阈值为 1.0 的梯度范数裁剪(Gradient Clipping),严防偶发梯度爆炸。

3. 数值结果与收敛精度评估

该算法在普通个人笔记本电脑上仅耗时约 1 分钟即可完成全部 100,000 个 Epoch 的训练。数值评估揭示了以下三个核心结论:

(1)与高精度有限差分基准高度吻合:论文图10将神经网络求得的价值函数 V(k)V(k) 与消费政策函数 c(k)c(k)(虚线)同基于超密网格的传统有限差分法(Finite Difference Method, 实线)进行了对比,二者在整个定义域 [0.1,6.0]\lbrack 0.1,\ 6.0\rbrack 内完全重合,几乎无法分辨任何肉眼可见的偏差。

(2)残差呈阶梯式平稳收敛:论文图11记录了优化轨迹,随着训练推进,HJB残差与消费FOC残差从最初的 10−110^{- 1} 迅速发生雪崩式下降,最终双双稳定收敛在 10−710^{- 7} 的极低微观量级。

(3)欧拉方程误差优于千分之一基准:论文图12绘制了全定义域上的无量纲对数欧拉方程误差 log10|1−c̃/c|\log_{10}|1\ - \ c\tilde{}/c|(其中 c̃=[V′(k)]−1/γc\tilde{} = \lbrack V\prime(k)\rbrack^{- 1/\gamma} 为根据边际价值倒推的消费)。结果显示,在绝大多数资本区间内,误差对数值稳定保持在 -4 至 -7 之间,远低于 -3 的学术公认高精度红线(即相对误差小于消费总量的 0.1%),仅在边界处出现极其微弱的网格自由法常见轻微边缘扰动。

六、全文结论的经济学与管理学启示总结

Jesús Fernández-Villaverde 教授在本文中不仅展示了一项前沿计算工具的成熟,更宣告了定量经济学在建模理念上的重大跃迁。突破“维度灾难”这一技术瓶颈,对于宏观经济治理、微观商业管理与公共政策制定产生了深远的影响:

1. 宏观经济调控与央行货币政策的新范式:从“代表性主体假象”到“全景异质性评估”

长期以来,中央银行与财政部门高度依赖代表性主体(Representative Agent)假设或低阶线性化的DSGE模型。然而,现实经济充斥着财富不平等、收入风险与借贷流动性约束。近年来兴起的异质性主体新凯恩斯模型(HANK)虽然在理论上能刻画贫富差距对货币政策传导的影响,但其状态变量包含全体居民的财富财富连续分布,维度无穷。深度学习使得经济学者能够直接以神经网络表征宏观财富分布的动态演化律,并在模型中无缝嵌入名义利率零下限(ZLB)的偶发绑定特征。这使宏观调控当局首次能够精准模拟:一项加息政策究竟如何非对称地冲击高负债贫困家庭,进而如何通过边际消费倾向差异反馈至总体总需求,彻底改变了宏观调控政策制定的颗粒度。

2. 极端尾部风险与金融系统性危机的全局防范

传统宏观金融模型严重依赖在确定性稳态附近的局部泰勒展开(Perturbation),这导致其天然预设经济系统具有强回归性,从而系统性低估了金融危机的爆发概率与传导烈度。深度学习作为一种真正的全局非线性解法,能够高精度刻画资产负债表螺旋破裂、非线性去杠杆以及“罕见灾难(Rare Disasters)”对跨期投资预期的永久性冲击。对于金融监管机构而言,这为更严谨地设计宏观审慎资本充足率要求、跨市场流动性压力测试以及系统重要性金融机构的风险传染链阻断提供了全新的量化评估底座。

3. 气候经济学与跨期绿色转型的长程最优规划

气候综合评估模型(Integrated Assessment Models, IAMs)要求将复杂的全球大气碳循环物理模型与多区域非平稳经济增长模型联立。这种系统状态维度极高,且伴随着不可逆的气候临界点(Tipping Points)与极端灾难非平稳动态。传统数值方法几乎完全瘫痪。深度学习的引入使经济学者能够在兼顾地理空间高维异质性的前提下,精准求解兼顾跨代公平、减排成本与能源技术内生演化的全球最优动态碳税路径,为实现“双碳”目标提供了可信的计算理论支撑。

4. 企业高维动态定价、供应链协调与竞合博弈策略

在产业组织与现代企业管理领域,拥有多品类资产和跨国生产线的大型平台企业面临着海量竞争对手交互的高维动态定价与库存控制问题。传统的启发式经验法则往往只能捕捉局部市场变化。深度学习方法能够将数以千计的互补品/替代品需求、竞对动态反应以及多级供应链在途库存状态转化为统一的动态规划网络进行全局求解。企业管理层可以据此构建抵御地缘政治冲击、关键零部件断供与原材料价格暴涨的高韧性动态协同决策系统。

5. 不完全信息下的机制设计与高维最优契约制定

在激励理论、拍卖设计与平台经济学中,当代理人的私人信息属于高维类型空间(High-dimensional Types,例如消费者的多维度偏好向量或供应商的多维度成本结构)时,传统的解析契约设计往往遭遇瓶颈。借助深度学习对最优拍卖机制和非线性契约的直接数值求解,平台运营者与政府规制机构能够设计出激发动机自洽、抗共谋且社会总福利最大化的前沿交易制度。

6. 组织数字化转型与决策范式的底层重塑

这篇综述给现代企业数字化战略最根本的启示在于:AI在企业中的最高阶应用,绝非仅仅停留在客服聊天、文本撰写或简单的销售预测等表面感知层,而在于**“微观感知—动态规划—全局最优闭环决策”**的认知升维。通过将深度学习融入运筹学优化与跨期动态资源分配,企业能够从“被动应对突发事件”转向“面向未来不确定性的全局动态策略规划”,实现管理科学的真正数智化跃升。

七、精选核心参考文献

以下为论文中引用的奠基性理论文献、前沿计算经济学成果及机器学习经典文献(按作者姓氏英文字母升序排列):

Abreu, Dilip, David Pearce, and Ennio Stacchetti. 1990. 'Toward a Theory of Discounted Repeated Games with Imperfect Monitoring.' Econometrica 58 (5): 1041–63.

Aggarwal, Charu C. 2023. Neural Networks and Deep Learning: A Textbook. 2nd ed. Springer.

Athey, Susan, and Guido W. Imbens. 2019. 'Machine Learning Methods That Economists Should Know About.' Annual Review of Economics 11: 685–725.

Auclert, Adrien, Bence Bardóczy, Matthew Rognlie, and Ludwig Straub. 2021. 'Using the Sequence-Space Jacobian to Solve and Estimate Heterogeneous-Agent Models.' Econometrica 89 (5): 2375–408.

Azinovic, Marlon, Luca Gaegauf, and Simon Scheidegger. 2022. 'Deep Equilibrium Nets.' International Economic Review 63 (4): 1471–525.

Barron, Andrew R. 1993. 'Universal Approximation Bounds for Superpositions of a Sigmoidal Function.' IEEE Transactions on Information Theory 39 (3): 930–45.

Bellman, Richard. 1957. Dynamic Programming. Princeton University Press.

Brumm, Johannes, and Simon Scheidegger. 2017. 'Using Adaptive Sparse Grids to Solve High-Dimensional Dynamic Models.' Econometrica 85 (5): 1575–612.

Cybenko, George. 1989. 'Approximation by Superpositions of a Sigmoidal Function.' Mathematics of Control, Signals and Systems 2 (4): 303–14.

Dell, Melissa. 2025. 'Deep Learning for Economists.' Journal of Economic Literature 63 (1): 5–58.

den Haan, Wouter J., and Albert Marcet. 1990. 'Solving the Stochastic Growth Model by Parameterizing Expectations.' Journal of Business and Economic Statistics 8 (1): 31–34.

Du, Simon, Jason Lee, Haochuan Li, Liwei Wang, and Xiyu Zhai. 2019. 'Gradient Descent Finds Global Minima of Deep Neural Networks.' In Proceedings of the 36th ICML, Vol. 97, 1675–85.

Fernández-Villaverde, Jesús, Samuel Hurtado, and Galo Nuño. 2023. 'Financial Frictions and the Wealth Distribution.' Econometrica 91 (3): 869–901.

Fernández-Villaverde, Jesús, Joël Marbet, Galo Nuño, and Omar Rachedi. 2024. 'Inequality and the Zero Lower Bound.' Journal of Econometrics 249 (C): 105819.

Fernández-Villaverde, Jesús, Galo Nuño, George Sorg-Langhans, and Maximilian Vogler. 2020. 'Solving High-Dimensional Dynamic Programming Problems Using Deep Learning.' Unpublished.

Fernández-Villaverde, Jesús, Juan Francisco Rubio-Ramírez, and Frank Schorfheide. 2016. 'Solution and Estimation Methods for DSGE Models.' In Handbook of Macroeconomics, Vol. 2, 527–724. North-Holland.

Goujon, Alexis, Arian Etemadi, and Michael Unser. 2024. 'On the Number of Regions of Piecewise Linear Neural Networks.' Journal of Computational and Applied Mathematics 441: 115667.

Hornik, Kurt, Maxwell Stinchcombe, and Halbert White. 1989. 'Multilayer Feedforward Networks Are Universal Approximators.' Neural Networks 2 (5): 359–66.

Jacot, Arthur, Seok Hoan Choi, and Yuxiao Wen. 2024. 'How DNNs Break the Curse of Dimensionality: Compositionality and Symmetry Learning.' Preprint, arXiv:2407.05664.

Judd, Kenneth L. 1998. Numerical Methods in Economics. MIT Press.

Kahou, Mahdi Ebrahimi, Jesús Fernández-Villaverde, Jesse Perla, and Arnav Sood. 2021. 'Exploiting Symmetry in High-Dimensional Dynamic Programming.' NBER Working Paper 28981.

Kingma, Diederik P., and Jimmy Ba. 2017. 'Adam: A Method for Stochastic Optimization.' Preprint, arXiv:1412.6980.

Maliar, Lilia, Serguei Maliar, and Pablo Winant. 2021. 'Deep Learning for Solving Dynamic Economic Models.' Journal of Monetary Economics 122: 76–101.

Poggio, Tomaso, Hrushikesh Mhaskar, Lorenzo Rosasco, Brando Miranda, and Qianli Liao. 2017. 'Why and When Can Deep-But Not Shallow-Networks Avoid the Curse of Dimensionality: A Review.' International Journal of Automation and Computing 14 (5): 503–19.

Prince, Simon J. D. 2023. Understanding Deep Learning. MIT Press.

Rumelhart, David E., Geoffrey E. Hinton, and Ronald J. Williams. 1986. 'Learning Representations by Back-Propagating Errors.' Nature 323 (6088): 533–36.

Soudry, Daniel, Elad Hoffer, Mor Shpigel Nacson, Suriya Gunasekar, and Nathan Srebro. 2018. 'The Implicit Bias of Gradient Descent on Separable Data.' Journal of Machine Learning Research 19: 1–57.

Stokey, Nancy L., Robert E. Lucas Jr., and Edward C. Prescott. 1989. Recursive Methods in Economic Dynamics. Harvard University Press.

Wilson, Andrew Gordon. 2025. 'Deep Learning Is Not So Mysterious or Different.' Preprint, arXiv:2503.02113.