无限线性方程组闭式解析解求解及相关参考咨询

阿华AIGC实验室

2026-4-23

看起来你遇到的是带二阶跳跃的连续时间马尔可夫链平稳分布方程组，这类问题在概率建模（比如人口动力学、排队论、化学反应模型）里很常见，我来帮你梳理下思路和相关方向：

问题定位与模型匹配

你的方程组本质是连续时间马尔可夫链的平稳平衡条件$\bar{x}Q=0$（$Q$是转移率矩阵），结合归一化条件$\sum_{i=0}^\infty x_i=1$。其中：

状态$i$的出转移率包含三类：$R$（跳到$i+1$）、$w_i=\gamma i$（跳到$i-1$）、$v_i=\beta i(i-1)$（跳到$i-2$）
入转移率则来自$i-1$（速率$R$）、$i+1$（速率$w_{i+1}$）、$i+2$（速率$v_{i+2}$）

这种带二阶跳跃的结构，属于广义生灭过程的范畴，区别于常规的只在相邻状态转移的生灭过程。

闭式解的求解思路

要得到非递归的解析解，常用的方法有两种：

1. 生成函数法（最常用）

定义生成函数$G(s) = \sum_{i=0}^\infty x_i s^{i$，将每个平衡方程乘以$s}i$后求和，转化为关于$G(s)$的微分方程：

对$i=0,1$的方程直接代入生成函数的低阶项
对$i\geq2$的方程，利用$\sum_{i=0}^\infty i x_i s^i = sG'(s)$、$\sum_{i=0}^\infty i(i-1)x_i s^i = s^2G''(s)$这类导数关系，把求和项转化为$G(s)$的一阶、二阶导数

整理后你会得到一个二阶线性常微分方程，如果这个方程有解析解（通常会涉及超几何函数、合流超几何函数这类特殊函数），那么$x_i$就是$G(s)$在$s=0$处的泰勒展开系数，这就是闭式解（用特殊函数表示的形式）。

2. 递推式求解

从平衡方程中整理出关于$x_{i+2}$的递推关系：
$$
x_{i+2} = \frac{1}{\beta(i+2)(i+1)} \left[ (R + \beta i(i-1) + \gamma i)x_i - R x_{i-1} - \gamma(i+1)x_{i+1} \right]
$$
结合$i=0,1$的初始方程，你可以尝试假设解的形式为$x_i = C \cdot \rho^i \cdot \frac{(\alpha)_i}{i!}$（其中$(\alpha)_i$是上升阶乘），代入递推式求解参数$\rho$和$\alpha$，再用归一化条件确定常数$C$。

平稳分布存在性的前提

需要注意：只有当$\sum_{i=0}^\infty x_i < \infty$时，平稳分布才存在。你可以通过递推式分析$x_i$的衰减速度，比如当$i$很大时，$x_{i+2} \approx x_i$（由$\beta i(i-1)x_i$主导），这时候需要验证级数是否收敛——通常只要参数满足一定约束（比如$\gamma$或$\beta$足够大，使得状态不会无限扩散），级数就会收敛。