【网络科学1】无时延固定拓扑网络同步性例题证明

最新推荐文章于 2024-07-12 17:06:20 发布

ssss1s2s1s3s8

最新推荐文章于 2024-07-12 17:06:20 发布

阅读量185

点赞数

分类专栏：网络科学导论文章标签：算法

本文链接：https://blog.csdn.net/qq_43710355/article/details/127420904

版权

网络科学导论专栏收录该内容

2 篇文章 0 订阅

订阅专栏

0绪论

本文主要参考书有：网络科学导论、现代控制理论、线性代数等。
主要解决的是《网络科学导论》一书中第11章网络同步与控制，课后习题11-1。如有证明任何问题欢迎指出并讨论，由于后面有矩阵运算，因此使用katex编辑会出现问题，因此部分公式使用截图。

1问题的提出

节点 $i$ 的初始状态演化方程为：
${\dot x_i}(t) = \sum\limits_{i = 1}^N {{a_{ij}}({x_j}(t) - {x_i}(t))}$
其中 $i = 1, 2, 3, 4..., N$

节点数为 $N$ 个，初始状态为任意给定，网络拓扑固定。证明：
$\mathop {\lim }\limits_{t \to \infty } {x_i}(t) = {1 \over N}\sum\limits_{j = 1}^N {{x_j}(0)}$

2问题的解决（两步走）

步骤一：将 ${x_i}(t)$ 求转化为求系数% ${e^{ - Lt}}$ 的极限：

由《网络科学导论》的366页11-4的结论可得
${l_{ij}} = - {a_{ij}}(i \ne j)$
邻接矩阵行和即每个节点的度为：
${l_{ii}} = \sum\limits_{j = 1}^N {aij} = {k_i}$
将节点度改写为对角线矩阵：
$D = diag({k_1},{k_2},...,{k_N})$
其中矩阵A为邻接矩阵
$A = ({a_{ij}})$

其中 ${l_{ij}}$ 是拉普拉斯矩阵 $L$ 的元素，即可得出 $L = D - A$ ，所以题目中所给演化方程可以改写为关于拉普拉斯矩阵 $L$ 和的函数 ${x_i}(t)$ （结合实例理解）：
${\dot X_i}(t) = - L{X_i}(t)$
而后使用一阶微分方程求解的方法：
两边积分
$\int {{{d{X_i}(t)} \over {{X_i}(t)}}} = \int { - Ldt}$
得到积分后的等式
$ln ({X_i}(t)) = - Lt + C$
两边取对数
${X_i}(t) = {e^{ - Lt}}*{X_i}(0)$
而后问题就转换为求时间趋于无穷的系数极限即可：
$\mathop {\lim }\limits_{t \to \infty } {e^{ - Lt}} = ?$
步骤二：求解系数 ${e^{ - Lt}}$ 的极限：
为解决这一问题首先在这陈述一下相关应用的定理：凯莱哈密顿定理

介绍凯莱哈密顿定理：
矩阵满足自身的特征方程，若阶矩阵的特征多项式为：
$f(\lambda ) = [\lambda I - B] = {\lambda ^n} + {b_{n - 1}}{\lambda ^{n - 1}} + ... + {b_1}\lambda + {b_0} （公式1.1）$
则有：
$f(B) = {B^n} + {b_{n - 1}}{B^{n - 1}} + ... + {b_1}B + {b_0}I = 0（公式1.2）$
矩阵指数 ${e^{Bt}}$ 可表示为 $B$ 的 $(n - 1)$ 阶多项式：
${e^{Bt}} = \sum\limits_{m = 0}^{n - 1} {{b_m}(t){B^m}} （公式1.3）$
凯莱哈密顿定理所得矩阵 $B$ 满足自己的特征方程，对于上式用 $B$ 的特征值 ${\lambda _i},i = 1,2,...,n$ 代替 $B$ 后仍然满足：
${e^{{\lambda _i}t}} = \sum\limits_{j = 0}^{n - 1} {{b_j}(t){\lambda _i}^j}（公式1.4）$
利用该式可求得待定系数 ${a_j}(t)$
所得的 $n$ 元一次方程组为：公式1.5
在这里插入图片描述
这里插入一个定理：有关实对称阵的性质：
任何一个实对称矩阵都可以正交对角化，即存在一个正交矩阵 $S$ （ ${S^T} = {S^{ - 1}}$ ）使得 ${S^{ - 1}}BS = {S^T}BS = \Lambda$ 且正交矩阵 $S$ 的列向量不相关、可单位化。

正式求解：
由上述现代控制理论中求状态转移矩阵使用的凯莱哈密顿定理公式（1.5）可得：
在这里插入图片描述
由于系数矩阵行列式必不为0因此可逆，所以两边左乘系数矩阵的逆矩阵得：

再由公式（1.3）可得下述公式：
在这里插入图片描述
由于 $B = - L$ 矩阵的特征值分别为 ${\lambda _1} = 0, - {\lambda _2}, - {\lambda _3},..., - {\lambda _n} < 0$ 因此极限：

即

因此设矩阵 $- L = B$ 转换为求 $\mathop {\lim }\limits_{t \to \infty } {e^{Bt}} = ?$ 由矩阵B为实对称阵可将B矩阵转化为：
$S\Lambda {S^T}$
矩阵S为正交矩阵
$S{S^T} = I$
将原极限转化为下式：
$\mathop {\lim }\limits_{t \to \infty } {e^{Bt}} = I - B - ... - {B^{{\rm{n}} - 1}} = S(I - \Lambda - {\Lambda ^2} - .. - {\Lambda ^{n - 1}}){S^T}$
矩阵S为矩阵B的特征向量，因此第一列为
在这里插入图片描述
题中所求转换为：

其中：

其中：矩阵S为

将S带入上述矩阵乘法可得：