A scalar mean-square stability and variance test for a discrete feedback loop with random gain and delay

技术摘要：具有随机增益与延迟的离散反馈回路标量均方稳定性与方差测试

问题陈述
本文研究了一个离散时间单输入反馈回路在两种同时存在的随机不确定性源（随机整数延迟 $D_k$ 和随机回路增益 $G_k$）下的均方稳定性与性能分析。该系统由以下方程建模：
$$e_{k+1} = e_k - \kappa G_k e_{k-D_k} + w_k$$
其中 $w_k$ 是加性白噪声干扰。虽然具有随机系数的离散系统的均方稳定性在经典上可以通过二阶矩算子的谱半径（$\rho(E[A_k \otimes A_k]) < 1$）来表征，但这种方法需要构建一个阶数为 $n^2$ 的矩阵，其中 $n$ 是增广状态的维度（由最大延迟决定）。对于具有显著延迟的系统，该矩阵维数会迅速增长，使得谱半径的计算变得昂贵，并且掩盖了增益方差与延迟离散化对稳定性裕度的具体贡献。

方法论
作者针对由单输入产生的秩一不确定性通道，提出了一种针对既定二阶矩判据的“无损约减”（lossless reduction）。其方法如下：

1. 增广状态表述： 将系统转化为增广状态空间形式 $x_{k+1} = A_k x_k + b w_k$，其中不确定性通过一个秩一项进入。
2. 秩一约减： 利用不确定性的结构，论文推导出了一个等价于谱半径条件的标量不等式。这依赖于引理 1，该引理确立了：对于受秩一项扰动的正算子，其谱半径小于 1 当且仅当涉及 Lyapunov 方程解的特定迹条件成立。
3. 脉冲响应表征： 该 Lyapunov 解被识别为“平均回路”（即用期望值代替随机变量后的回路）脉冲响应（$h$）的自相关托普利茨矩阵（Toeplitz matrix）。
4. “预算”的推导： 稳定性条件被简化为一个涉及“均方预算”（$B$）的标量不等式。预算定义为：
   $$B = \kappa^2 \left( E[G_k^2] \|h\|^2 - \|v * h\|^2 \right)$$
   其中 $\|h\|^2$ 是平均回路脉冲响应的能量，$v$ 代表增益-延迟联合分布的一阶矩，$*$ 表示离散卷积。
5. 性能分析： 论文将该标量指标扩展到加性白噪声下的性能分析，推导出了受控偏差稳态协方差的闭式表达式。

核心贡献与结果

• 精确标量判据： 论文证明了对于具有联合随机增益和延迟的单输入回路，均方稳定的充分必要条件是 $B < 1$（假设平均回路是 Schur 稳定的）。这用标量递归和两次内积取代了 $O(n^2)$ 的特征值计算。
• 不确定性成本的分解： 在增益与延迟相互独立的假设下，预算 $B$ 精确地拆分为两个非负项：
  • $B_G = \kappa^2 \sigma_G^2 \|h\|^2$：由增益方差引起的成本。
  • $B_D = \gamma^2 (\|h\|^2 - \|p * h\|^2)$：由延迟变异性引起的成本。
    论文证明了 $B_D \geq 0$，且当延迟为确定性时，$B_D$ 等于零。这量化了“延迟离散化带来的惩罚”，即脉冲响应被延迟分布平滑时所损失的能量。
• 方差放大因子： 预算 $B$ 被证明是一个精确的方差放大因子。受控偏差的稳态方差由 $\sigma_w^2 \|h\|^2 / (1 - B)$ 给出。这使二元稳定性测试转变为定量设计度量，揭示了增加控制强度虽能降低响应能量（$\|h\|^2$）但会增加预算（$B$），从而存在一个使稳态方差最小化的最优控制强度。
• 计算效率： 数值验证表明，对于中等延迟（例如 $n=65$ 时），标量法比直接矩阵计算快几个数量级（约 660 倍）。
• 验证： 结果通过直接二阶矩判据、迭代协方差递归以及数百种随机配置下的蒙特卡洛轨迹模拟进行了验证，显示出与机器精度一致的结果。

意义与主张
论文强调其贡献并非一种新的通用稳定性理论，而是针对特定结构化系统（单输入、秩一不确定性）的一种“无损约减”。其意义体现在三个方面：

1. 计算可行性： 它消除了 $n^2$ 的计算负担，使得具有较大延迟的回路进行稳定性分析成为可能。
2. 诊断清晰度： 它明确地将稳定性裕度的消耗归因于增益方差或延迟离散化，而全二阶矩矩阵的谱半径无法提供这种区分。
3. 设计实用性： 通过将稳定性直接与通过预算因子联系起来的稳态方差相连，它提供了一个直接的选择控制强度的规则。论文指出，仅基于确定性稳定性边界选择强度，会导致稳态方差高于通过标量预算得出的最优选择（例如在提供的示例中高出约三分之一）。

作者明确指出了局限性：结果假设增益和延迟在时间上是独立同分布（i.i.d.）的对（尽管允许步内相关），且采用单输入，并假设平均矩阵是 Schur 稳定的。若不进行进一步修改，该约减方法不适用于时间相关的系数、马尔可夫切换模型或多输入通道。