电力网络计算中的稀疏技术.pptxVIP

电力网络计算中的稀疏技术;电力网络的稀疏性 稀疏存储技术 稀疏矩阵的因子表分解 稀疏线性方程组的求解 稀疏矢量技术 节点优化编号技术;稀疏矩阵：只有少量非零元素 稀疏度：对m*n阶矩阵A，其中有τ个非零元素，则可定义稀疏度 线性方程组的求解主要涉及大量的四则运算，而0的四则运算相当特殊 稀疏计算亦可称作“排零”计算;以求解节点电流-电压线性方程为例： 非线性的潮流方程本质相同，且也需在迭代过程中求解线性方程 系数矩阵为节点导纳矩阵 对角元：与相应节点相连的所有支路导纳之和，称自导纳 非对角元：与相应行列对应的节点间所有支路导纳之和的相反数，称互导纳 节点导纳矩阵为对称矩阵 只有电力网络中存在支路，相应非对角元才不为0;设有1000条母线的电力系统，母线出线度平均为10，其稀疏度为 实际电网仅有非常少量的枢纽变电站存在出线度为10左右的母线 大量母线出线度仅为1~2 发电机机端母线 终端负荷母线 联络母线 …… 算法是否采用排零操作可影响计算速度几十上百倍;核心：不存储零元素?仅保留非零元素在原矩阵中的数值及位置信息?应在必要时轻易恢复成满阵存储格式 上世纪六七十年代提出 当前计算机硬件速度和容量已发生了翻天覆地的变化，还要考虑稀疏存储吗？ 需要分析的电力系统规模也显著扩大 要求计算的速度也更快（如在线分析） 节省计算机内存占有量 尽量减少检索矩阵元素所耗时间;a11;a11;;a12;矩阵化为上三角矩阵的初等变换过程 等价的矩阵计算 图论解释;;;;;;;;;;节点号从小到大依次操作 对每个节点分为两步： 规格化运算 消去运算;用对角元去除各非零非对角元 图上的操作： 修正由p发出的互边 新边权为原边权除以p自边边权;app正下方若有非零元素，则相应行需做消去运算 消去某行时app右方非零元素对应列对本行有影响;;;;;;;;;;;;;;;;;;;;节点1的点位为0，对相关节点电位无影响;;;;;;;;;;前述内容研究的是A稀疏时的排零技术 若b中也只有少量非零元素，如何提高计算速度？ 若只需计算少数的x，如何提高计算速度？ 核心内容：分析出前代回代过程中那些计算步是必不可少的，哪些是可以忽??的;定义：有向因子图上，每个节点发出边中对侧节点编号最小的边构成的子图 为何一定是树？ 连通图 节点数比支路数多1;定义：有向因子图上，每个节点发出边中对侧节点编号最小的边构成的子图 为何一定是树？ 连通图 节点数比支路数多1;定义：有向因子图上，每个节点发出边中对侧节点编号最小的边构成的子图 为何一定是树？ 连通图 节点数比支路数多1;定义：道路树上点沿道路树到树根的路径 右边分别为节点2和节点4的道路树 相关性质：原有向因子图上任一点发出的边的终点必在该点的道路上 树支 连支;点集中所有点的道路的并集 右边为点集{2,4}的路集 相关性质 有向因子图上任一边两个端点的路集即为小号节点的道路 若任一点集中节点对之间都有边，则点集的路集即为最小号节点的道路 ……;前代只在非零元素对应点集的路集上进行（前代只直接影响该边的终点） 设b 中只有2号、4号两个元素非零;路集上任一点的前代必须在路集上比该点编号小且其道路经过该点的点的前代完成后才能进行 以10号节点为例;回代只在待解元素对应点集的路集上进行（回代只直接影响该边的起点） 设x 中只有2号、4号两个元素需求解;路集上任一点的回代必须在路集上比该点编号大的点的回代完成后才能进行 以10号节点为例;改善稀疏矩阵计算的性能 改善稀疏矢量计算的性能;目标：因子表分解过程中尽可能少出现新的非零元素 Tinney-1 Tinney-2 Tinney-3;;;;;;;;;;;;;;;;不同的节点顺序，高斯消去时新增的非零元情况不同 如何确定节点顺序，使得高斯消去时新增非零元尽可能少？ 静态优化：按原网络各节点所连接的支路数由少到多 半动态优化：按每次消去后网络各节点所连的支路数由少到多 全动态优化（不实用，不讲）;按有向A图中节点出度由小到大编号 称静态节点优化编号 缺点：某点做完消去后，其发出的边对后续消去无影响，不应计入;先选出度最小节点进行编号 得到消去该节点的新网络结构，该点及发出的边均去掉，在剩下的图中重复刚才过程 称半动态节点优化编号 可大大减少新出现的非零元素个数 缺点：出线度小并不一定意味着新增非零元素个数少（若原来已有非零元素，则只改其值，不会出现新非零元素）;与Tinney-2基本相同 编号原则不是出线度最小，而是新增边最少 效果最好 缺点：最耗机时 最广泛应用的是Tinney-2;目标：使生成的道路树尽量矮 较复杂，不讲 实际应用时应与改善稀疏矩阵计算性能的方法综合考虑