高等电力系统分析 课件.ppt

逐行回代，得： 1、由于电力网络结构的特点，每个节点仅与3~5个节点相连，因此描述网络结构的矩阵是稀疏矩阵。 n*m的矩阵，非零元τ个，稀疏度等于τ/ n*m 如果系统有N=500个节点，平均每个节点与5条支路相连，则 稀疏度=5*500/（500*500）=1% 2、计算中，我们仅关心一部分的变量：稀疏矢量。 3、与稀疏矩阵和稀疏矢量相关的运算中，零元素不参与存储和计算——排零存储和排零计算 1.4 电力网络求解的稀疏技术：引入原因 特点：排零存储，即只存储其中的非零元和有关的检索信息。 要求： 节省内存 方便地检索和存取 考虑网络结构变化时能方便地对存储的信息加以修改 稀疏矢量：存储矢量中的非零元值和相应的下标 稀疏矩阵：考虑稀疏结构和所采用的算法 1.4 电力网络求解的稀疏技术：稀疏存储 VA a11 a12 a14 a21 a22 a23 a33 a42 a43 a44 IA 1 1 1 2 2 2 3 4 4 4 JA 1 2 4 1 2 3 3 2 3 4 散居格式 a13 1 3 常用存储方式 散居格式 按行（列）存储格式 三角检索存储格式 链表存储格式 例： VA a11 a12 a14 a21 a22 a23 a33 a42 a43 a44 JA 1 2 4 1 2 3 3 2 3 4 IA 1 4 7 8 按行存储格式 VA a11 a12 a13 a14 a21 a22 a23 a33 a42 a43 a44 JA 1 2 3 4 1 2 3 3 2 3 4 IA 1 5 8 9 修改后 三角检索存储格式 U a12 a14 a23 JU 2 4 3 IU 1 3 4 4 L a21 a42 a43 IL 2 4 4 JL 1 2 3 4 D a11 a22 a33 a44 U a12 a13 a14 a23 JU 2 3 4 3 IU 1 4 5 5 L a21 a42 a43 IL 2 4 4 JL 1 2 3 4 D a11 a22 a33 a44 修改后 链表存储格式 VA a11 a12 a14 a21 a22 a23 a33 a42 a43 a44 JA 1 2 4 1 2 3 3 2 3 4 LINK 2 3 0 5 6 0 0 9 10 0 IA 1 4 7 8 VA a11 a12 a14 a21 a22 a23 a33 a42 a43 a44 a13 JA 1 2 4 1 2 3 3 2 3 4 3 LINK 2 11 0 5 6 0 0 9 10 0 3 IA 1 4 7 8 修改后 小结 存储单元 检索 修改 散居格式 3τ 工作量大 易 按行存储格式 2τ+n 按行检索 难 三角检索存储格式 2τ+n 按行（列）检索 难 链表格式 3τ+n 按行检索 易 1、稀疏矩阵的因子分解 采用高斯消去法进行计算 1）按行规格化 2）消去运算 1.4 电力网络求解的稀疏技术：排零计算 采用三角检索存储格式时 例 U 7 -3 JU 2 4 IU 1 3 3 3 L 5 -2 IL 2 4 JL 1 2 3 3 D 2 4 5 6 U 3.5 -1.5 7.5 JU 2 4 4 IU 1 3 4 4 L 5 -2 IL 2 4 JL 1 2 3 3 D 2 -13.5 5 6 U 3.5 -1.5 -0.555 JU 2 4 4 IU 1 3 4 4 L 5 -2 IL 2 4 JL 1 2 3 3 D 2 -13.5 5 4.889 U 3.5 -1.5 -0.555 JU 2 4 4 IU 1 3 4 4 L 5 -2 IL 2 4 JL 1 2 3 3 D 2 -13.5 5 4.889 3 = p 2、利用稀疏矩阵因子表求解稀疏线性代数方程组 1）前代过程 计算流程 除法运算 回代运算 计算流程 采用三角检索存储格式时 采用按列存储格式 U u12 u24 u34 u15 u25 u45 IU 1 2 3 1 2 4 JU 1 2 2 4 7 例 U 3.5 -1.5 -0.555 IU 2 4 4 IU 1 3 4 4 L 5 -2 IL 2 4 JL 1 2 3 3 D 2 -13.5 5 4.889 A图 有向A图 赋权有向A图 A图：和矩阵A有相同拓扑结构的网络图 有向A图：对给定A图及节点编号，规定边的正方向由小号节点指向大号节点 赋权有向A图：在有向A图中，将A的非对角非零元的值赋给互边，将A的对角元素的值赋给自边 1、基本定义和术语 1.4 电力网络求解的稀疏技术：基于图论 因子图 有向因子图 赋权有向因子图 因子图：和因子表矩阵U有相同拓扑结构的网络图 有向因子图：在因子图上规定边的正方向由小号节点指向大号节点形成 赋权有向因子图：在有向