4.第四章潮流计算中的特殊问题
[理学]4 潮流计算
![[理学]4 潮流计算](https://img.taocdn.com/s3/m/5f0f7c30a8114431b90dd82b.png)
方程数： n 1 m n m 1 2(n 1)
i
未知量： ei , f i , i ( PQ PV ) , 2(n 1)
第四讲 电力系统潮流计算
（2）极坐标下的数学方程
将 Vi Vi e
j i
和 Yij Gij jBij 代入 Pi jQi Vi Yij V j
n个节点，m个PQ节点，n m 1个PV节点
编号
1 m
PQ
m 1 n 1
PV
n
平衡点
I YV
I i Yij V j
j 1 n
Vi I i Vi Y ij V j
j 1
n
Pi jQi Vi Yij V j
x(k )
x
讨论收敛区域和收敛条件。又称切线法。
第四讲 电力系统潮流计算
2、设初始点 xo , f ( xo ) 0
f ( xo x ) 0 df f ( xo ) dx f ( xo ) x df dx 1 d2 f x 2 dx 2 x 0
x0
x 2
( 0) x1 0 ( 0) x 2
(1) ( 0) ( 0) x1 x1 x1 (1) ( 0) ( 0) x2 x2 x2
第四讲 电力系统潮流计算
3、多维非线性方程组的迭代公式
基于同样的思想，我们可以得到 n维非线性方程 —牛顿 拉夫逊迭代公式
j 1
j i j i ( G jB ) V e ij ij i j 1 n
n
*
第四章 电力系统潮流计算
第四章 电力系统潮流分析与计算电力系统潮流计算是电力系统稳态运行分析与控制的基础，同时也是安全性分析、稳定性分析电磁暂态分析的基础（稳定性分析和电磁暂态分析需要首先计算初始状态，而初始状态需要进行潮流计算）。
其根本任务是根据给定的运行参数，例如节点的注入功率，计算电网各个节点的电压、相角以及各个支路的有功功率和无功功率的分布及损耗。
潮流计算的本质是求解节点功率方程，系统的节点功率方程是节点电压方程乘以节点电压构成的。
要想计算各个支路的功率潮流，首先根据节点的注入功率计算节点电压，即求解节点功率方程。
节点功率方程是一组高维的非线性代数方程，需要借助数字迭代的计算方法来完成。
简单辐射型网络和环形网络的潮流估算是以单支路的潮流计算为基础的.本章主要介绍电力系统的节点功率方程的形成，潮流计算的数值计算方法,包括高斯迭代法、牛顿拉夫逊法以及PQ 解藕法等。
介绍单电源辐射型网络和双端电源环形网络的潮流估算方法.4—1 潮流计算方程——节点功率方程1。
支路潮流所谓潮流计算就是计算电力系统的功率在各个支路的分布、各个支路的功率损耗以及各个节点的电压和各个支路的电压损耗.由于电力系统可以用等值电路来模拟，从本质上说,电力系统的潮流计算首先是根据各个节点的注入功率求解电力系统各个节点的电压，当各个节点的电压相量已知时，就很容易计算出各个支路的功率损耗和功率分布.假设支路的两个节点分别为k 和l ，支路导纳为kl y ,两个节点的电压已知，分别为kV 和l V ，如图4—1所示。
图4-1 支路功率及其分布那么从节点k 流向节点l 的复功率为（变量上面的“－”表示复共扼）:)]([lk kl k kl k kl V V y V I V S -== （4—1） 从节点l 流向节点k 的复功率为:)]([kl kl l lk l lk V V y V I V S -== （4—2） 功率损耗为：2)()(klkl l k kl l k lk kl kl V y V V y V V S S S ∆=--=+=∆ （4—3)因此，潮流计算的第一步是求解节点的电压和相位,根据电路理论，可以采用节点导纳方程求解各个节点的电压。
电力网络分析学后总结-推荐下载
电力网络分析是电力系统分析的关键环节。
随着国民经济的不断提高，社会对电能质量的需求也越来越高。
电力系统分析的作用至关重要。
高等电力网络分析是通过归纳、总结、提升，抽象出电网分析中的共性问题，从更基础的层面来描述和解决电网分析问题。
此书把电力网络分为两部分来研究。
第一部分为基础篇，介绍电力网络分析的基本原理。
第二部分为应用篇，介绍潮流计算和故障分析。
第一部分 电力网络分析基本原理一、电力网络分析的一般方法1.1 网络分析概述电力网络包含两个要素：电气元件及其联接方式。
电力网络的运行特性的约束和元件之间联接关系的约束（拓扑约束）共同决定。
元件的特性约束由欧姆定律来描述：, , .Ri =u dL dt=u ∫1C idt =u 网络的拓扑约束由基尔霍夫定律来描述：基尔霍夫电流定律：. 基尔霍夫电压定律：.∑I =0∑V =0有关电力系统分析计算问题包括状态估计、潮流计算、经济调度、故障分析、稳定计算等，这些问题既相互关联，又各有侧重点。
如状态估计可以为潮流计算提供良好的初值，而潮流计算则是经济调度、故障分析、稳定计算与系统控制的出发点。
网络分析是解决这些所有问题的共同基础。
研究一个特定的电力系统运行问题应当包括四个基本步骤：1、建立电力网络元件的数学模型；2、建立电力网络的数学模型；3、选择合理的数值计算方法；4、电力网络问题的计算机求解。
网络分析中常用的关联矩阵有：节-支关联矩阵、回-支关联矩阵、割-支关联矩阵。
1.2 电力网络支路特性的约束 一般支路如图：图1：一般支路元件的约束特性可用以下支路方程来表示：=或=V k+E k z k(I k+I sk)I k+I sk y k(V k+E k)把网络内所有支路方程集合在一起，引入电动势矢量和电流源矢量.E S,I S可以得到网络的支路方程=或=V b+E s z b(I b+I s)I b+I s y b(V b+E s)为原始导纳矩阵和原始阻抗矩阵，若网络内所有的支路间不存在互感，z b，y b是对角阵，对角线元素既是相应的支路阻抗和支路导纳；若存在互感则z b，y b在相应于互感支路相关的位置上存在非零非对角线元素。
第四章复习题新
第四章复习题一、选择题1、计算机求解潮流方程组的方法是( B ).A．解析法B．数值方法C．手算法D．对数法2、在电力系统潮流计算中，PV节点的待求量是（）A．Q、δB．P、Q C．V、δD．P、V3、计算机解潮流方程时，经常采用的方法是（）A．递推法B．迭代法C．回归法D．替代法4、潮流计算中的P-Q分解法是在哪一类方法的基础上简化来的（）A．阻抗法B．直角坐标形式的牛顿——拉夫逊法C．高斯——赛德尔法D．极坐标形式的牛顿——拉夫逊法5、电力系统潮流计算时某物理量的单位为Mvar，则该量是( )A.有功功率B.无功功率C.视在功率D.有功电量6、在电力系统中平衡节点的数量( )A.必有一个B.是大量的C.少量或没有D.数量最多7、一般潮流分析中将节点分为几类( )A.四类B.三类C.二类D.一类8、对于含有n-1个PQ节点的n节点网络，雅可比矩阵的阶数为（）（A）n阶（B）n-1阶（C）2n阶（D）2（n-1）阶9、用牛顿—拉夫逊法进行潮流计算时，线性修正方程求解的是（）A．线路的功率B．节点的注入功率C．节点的电压值D．节点电压的修正量10、节点导纳矩阵为方阵，其阶数等于（）A．网络中所有节点数B．网络中除参考节点以外的节点数C．网络中所有节点数加1 D．网络中所有节点数加211、电力系统潮流计算的方程属于( )A．线性代数方程组B．微分方程C．非线性代数方程组D．微分方程组二、判断题1、用牛顿—拉夫逊法进行潮流计算时，线性修正方程求解的是节点的电压值。
（）2、同样的迭代次数，牛顿-拉夫逊法比PQ分解法精度高。
（）3、在收敛判据相同的情况下，牛顿-拉普逊法和PQ分解法潮流计算的精度相同。
（）4、PQ分解法潮流计算过程中，由于简化了很多因素，所以比牛顿-拉普逊法的收敛速度快。
（）5、PQ分解法潮流计算过程中，由于简化了很多因素，所以比牛顿-拉普逊法的计算速度快。
（）6、在潮流的计算机算法中，由于节点导纳矩阵是对称矩阵，所以对应的雅克比矩阵也是对称矩阵，可以只计算上三角的元素。
五、最小化潮流计算及潮流计算中的自动调整
limit
Qi Qilim it Qi
(5－5) (5－6)
Qi 和Ui之间的灵敏度关系由下式给出：
Ui Rii Qi
Rii是增广的B”的逆矩阵中和节点i相对应的对角线元素。根据需调整的Qi ， 用（5－6）式算出 Ui，最后将节点i 的给定电压调整到新值
Uinew UiSP Ui
©版权所有
8
电力系统稳态分析
计算步长公式的推导（续）
目标函数改写为 F ( x) fi ( x) ai bi 2ci ( ) (5－3)
2 2 i 1 i 1
n
n
将F(x)对μ求导，并令其等于零，求得μ*
2 d ( ) d n 2 ai bi ci d d i 1 2 2 a b ci bi 2 ci 0 i i i 1 n
ys－ y(x(k＋1)) ＝ ys－ y(x(k)＋μ(k)x(k)) ＝ys－[y(x(k))+μ(k)J(x(k))∆x(k)+(μ(k))2y(∆x(k))] =a(k)+ μ(k)b(k)+ (μ(k))2c(k)
©版权所有
11
电力系统稳态分析
具体应用，三种情况
从一定的数值出发，原来的潮流问题有解。
©版权所有
13
电力系统稳态分析
自动调整的两类方法 s 按照所要保持的系统状态量y 和当前的计算值y的大 小，不断地在一次次迭代中间改变某一个控制参数 x的大小。x大小的调整按照偏差反馈的原理进行。 ∆x＝a(ys－y) 改变原来潮流方程的构成。
©版权所有
潮流计算问题
潮流计算问题潮流计算的定义（课后题）各种潮流计算模型和算法的特点、适用范围以及相互之间的区别和联系（课后题）影响潮流收敛性的因素，以及如何改善潮流计算的收敛性（课后题）通过功率方程说明为什么潮流计算的数学模型是非线性的应该采用什么样的数学方法求解（03A、05A）电力系统的潮流计算有哪些常规算法有哪些扩展算法（05B）潮流计算的目的是什么其数学模型是什么有何特点（06B）简要说明潮流计算的概念、模型及计算方法。
（07B）高斯赛德尔迭代法和牛顿拉夫逊迭代法是常规的潮流计算方法，请介绍一下最优潮流（OPF）算法的原理及其应用。
（04电科院）潮流计算的目的：常规潮流计算的目的是在已知电力网络参数和各节点的注入量的条件下，求解各节点电压。
目的1：1.在电网规划阶段，通过潮流计算，合理规划电源容量和接入点，合理规划网架，选择无功补偿方案，满足规划水平年的大小方式下潮流交换控制、调峰、调相、调压的要求。
2.在编制年运行方式，在预计复合增长及新设备投运基础上，选择典型方式进行潮流计算，发现电网中的薄弱环节，供调度人员异常调度控制参考，并对规划、基建部门提出改进网架结构，加快基建进度的建议。
3.正常检修及特殊运行方式下的潮流计算，用于日常运行方式的编制，指导发电厂开机方式，有功、无功调整方案及负荷调整方案，满足线路、变压器热稳定要求及电压质量要求。
4.预想事故、设备退出运行对静态安全分析的影响及做出预想的运行方式调整方案。
目的2：A.检查电力系统各元件是否过负荷；B.检查电力系统各节点的电压是否满足电压质量的要求；C.根据对各种运行方式的潮流分布计算，可以正确的选择系统接线方式，合理调整负荷，以保证电力系统安全、可靠地运行，向用户供给高质量的电能；D. 根据功率分布，可以选择电力系统的电气设备和导线截面积，可以为电力系统继电保护整定计算提供必要的数据等；E. 为电力系统扩建和规划提供依据；F. 为调压计算、经济运行计算、短路计算等提供必要的数据。
潮流计算中的特殊问题
Vti fd
sin
Q
Vt It
sin
xad xs
Vti fd
cos
Vt2 xs
消去角度，有
P2
Q
Vt 2 xs
xad xs
Vti fd
2
i fd i fd max
Q
0,
Vt 2 xs
xad xs
Vti fd max
P
发电机无功输出的转子导体末端过热限制
欠励时电枢电流 产生的磁通和励 磁电流产生的磁 通同向叠加，造 成转子导体末端 过热（Kundur）
与kT相关的雅可比矩阵元素表达式 （直角坐标）
Pp
K T
KT
ep
eqG pq
f q B pq
f p eq B pq f qG pq
2
e
2 p
f
2 p
G pq
KT
Q p
KT
KT
f p eqG pq f q B pq
e p eq B pq f qG pq
2
e
2 p
f
2 p
B pq
Pq
KT
KT
eq e pG pq f p B pq
f q e p B pq f pG pq
KT
KT
f q e pG pq f p B pq
eq e p B pq f pG pq
K
T
直流潮流——应用场合
只关心有功潮ቤተ መጻሕፍቲ ባይዱ分布 不关心节点电压情况 对计算速度要求高
kT），节点q除已知P、Q外，还已知V，故称PQV 节点 潮流计算电压变量少一个 潮流计算多了一个变压器变比的变量 方程和变量个数仍相等，可求解
国网考试之电力系统分析：第四章复习题---3页
第四章复习题一、选择题1、n节点电力系统中，PQ节点的数目为m，平衡节点的数目应为（）。
A．n-m-1 B．n-m-2 C．1 D．可以为02、在电力系统潮流计算中，PV节点的待求量是（）A．Q、δB．P、Q C．V、δD．P、V3、计算机解潮流方程时，经常采用的方法是（）A．递推法B．迭代法C．回归法D．替代法4、潮流计算中的P-Q分解法是在哪一类方法的基础上简化来的（）A．阻抗法B．直角坐标形式的牛顿——拉夫逊法C．高斯——赛德尔法D．极坐标形式的牛顿——拉夫逊法5、电力系统潮流计算时某物理量的单位为Mvar，则该量是( )A.有功功率B.无功功率C.视在功率D.有功电量6、在电力系统中平衡节点的数量( )A.必有一个B.是大量的C.少量或没有D.数量最多7、一般潮流分析中将节点分为几类( )A.四类B.三类C.二类D.一类8、潮流计算时，下列节点中，哪一类节点一般只有一个，且必需有一个( )A.PQ节点B.PU节点C.平衡节点D.QU节点9、用牛顿—拉夫逊法进行潮流计算时，线性修正方程求解的是（）A．线路的功率B．节点的注入功率C．节点的电压值D．节点电压的修正量10、节点导纳矩阵为方阵，其阶数等于（）A．网络中所有节点数B．网络中除参考节点以外的节点数C．网络中所有节点数加1 D．网络中所有节点数加211、潮流方程是( )A．代数方程B．微分方程C．代数方程组D．微分方程组12、计算潮流时牛顿—拉夫逊法与高斯—塞德尔法相比的主要优点是（）A．对初值要求低B．占用内存少C．收敛性好，计算速度快D．复杂13、解潮流方程的方法是( )A．解析法B．数值方法C．手算法D．对数法二、判断题1、用牛顿—拉夫逊法进行潮流计算时，线性修正方程求解的是节点的电压值。
（）2、同样的迭代次数，牛顿-拉夫逊法比PQ分解法精度高。
（）三、填空题1、在计算机算法中，若PV节点的无功功率越限，则该节点应________________________。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
第四章 潮流计算中的特殊问题第一节 负荷的静态特性负荷的功率是系统频率和电压的函数。
在潮流计算中可以认为频率变化不大。
但由于发电机或输电设备的开断会引起电压较大的变化，在潮流计算中计及负荷的静态电压特性是合理的。
负荷的电压静态特性就是负荷的有功和无功功率与电压大小的关系，一般表达如下：⎪⎪⎪⎩⎪⎪⎪⎨⎧⎥⎥⎦⎤⎢⎢⎣⎡+⎪⎪⎭⎫ ⎝⎛+⎪⎪⎭⎫ ⎝⎛=⎥⎥⎦⎤⎢⎢⎣⎡+⎪⎪⎭⎫ ⎝⎛+⎪⎪⎭⎫ ⎝⎛=Qi is i Qi is i Qi Di Di Pi is i Pi is i Pi Di Di c V V b V V a Q Q c V V b V V a P P 2)0(2)0( （4-1） 式中系数满足11=++=++Qi Qi Qi Pi Pi Pi c b a c b a)0(Di P 、)0(Di Q 是在设定电压is V 下的负荷值。
组成负荷的三部分被分别看做恒定阻抗部分、恒定电流部分和恒定功率部分，所以（4-1）称为负荷的ZIP 模型。
当0=Pi a 、0=Qi a 时，忽略电压的二次项。
潮流计算中计及负荷的静态电压特性的方法：1、节点功率的不平衡量计算:⎪⎪⎪⎩⎪⎪⎪⎨⎧-⎥⎥⎦⎤⎢⎢⎣⎡+⎪⎪⎭⎫ ⎝⎛+⎪⎪⎭⎫ ⎝⎛-=--=∆-⎥⎥⎦⎤⎢⎢⎣⎡+⎪⎪⎭⎫ ⎝⎛+⎪⎪⎭⎫ ⎝⎛-=--=∆),(),(),(),(2)0(2)0(θθθθV Q c V V b V V a Q Q V Q Q Q Q V P c V V b V V a P P V P P P P i Qi is i Qi is i Qi Di Gi i Di Gi i i Pi is i Pi is i Pi Di Gi i Di Gi i (4-2)2、牛顿法雅可比矩阵子矩阵N 和L 的对角线元素要增加i i V P ∂∆∂和ii V Q ∂∆∂ 3、P-Q 分解法，Q-V 迭代的系数矩阵B ''的对角线元素也应增加i i V Q ∂∆∂，这样B ''不再是常数了。
为了节省计算量，ii V Q ∂∆∂也可取为常数，如忽略二次项取0=Qi a ，或不改变B ''，但功率不平衡量要按（4-2）计算。
负荷电压静态特性模型的指数形式⎪⎪⎩⎪⎪⎨⎧⎪⎪⎭⎫ ⎝⎛=⎪⎪⎭⎫ ⎝⎛=βαis i Di Di is i Di Di V V Q Q V V P P )0()0( （4-3） 8.1~5.0=α、6~5.1=β在动态潮流计算中，不能不考虑频率的变化。
考虑频率变化时式（4-1）、（4-3）变为。
⎪⎪⎪⎩⎪⎪⎪⎨⎧⎪⎪⎭⎫ ⎝⎛-+⎥⎥⎦⎤⎢⎢⎣⎡+⎪⎪⎭⎫ ⎝⎛+⎪⎪⎭⎫ ⎝⎛=⎪⎪⎭⎫ ⎝⎛-+⎥⎥⎦⎤⎢⎢⎣⎡+⎪⎪⎭⎫ ⎝⎛+⎪⎪⎭⎫ ⎝⎛=002)0(002)0(11f f f k c V V b V V a Q Q f f f k c V V b V V a P P Qi Qi is i Qi is i Qi Di Di Pi Pi is i Pi is i Pi Di Di ⎪⎪⎩⎪⎪⎨⎧⎪⎪⎭⎫ ⎝⎛-+⎪⎪⎭⎫ ⎝⎛=⎪⎪⎭⎫ ⎝⎛-+⎪⎪⎭⎫ ⎝⎛=00)0(00)0(11f f f k V V Q Q f f f k V V P P Qi is i Di Di Pi is i Di Di βα 当考虑频率变化时，频率也是待求的未知量，应出现在潮流方程中。
模型中系数的选取属于负荷建模的问题，仍未得到很好的解决。
第二节 节点类型的相互转换一、PV 节点转换为PQ 节点当在迭代过程中出现PV 节点无功功率越限时，可以再迭代几次，如果无功仍越限，说明PV 节点电压设置不合理，应进行调整：如果无功功率越下限，检查是否电压设置过低如是可适当提高电压设定值，或转换为PQ 节点，无功定值置下限值。
如果无功功率越上限，说明节点无功功率不能支持设定的电压，可适当调低电压设定值，或转换为PQ 节点，无功定值取上限值。
PV 节点转换为PQ 节点的处理方法：1、直角坐标方式的节点不平衡量由2i V ∆变为i Q ∆；2、牛顿法极坐标方式的修正方程加1个Q ∆方程；3、P-Q 分解法，θ-P 迭代不变，V Q -迭代的系数矩阵有两种处理方法：（1）B ''增加一行一列，如增加到最后：⎥⎦⎤⎢⎣⎡''=''ii T i B B B B B i ~ （4-4） 新的矩阵的因子表可由右下角加边的因子表修正法求出。
（2）B ''的对角元加大数在形成B ''时包含PV 节点对应的导纳，但PV 节点的对角元加一个很大的数。
这样在正常Q-V 迭代时，PV 节点的电压修正零接近于0，不会影响其他节点的电压修正量。
当PV 转换为PQ 节点时，将加的大数去掉。
ΔB B B -''=''~（4-5）采用因子表秩1修正法得到新的因子表二、PQ 节点转换为PV 节点当在迭代过程中出现PQ 节点电压越限时，可以再迭代几次，如果电压仍越限，说明PQ 节点无功设置不合理，应进行调整：如果电压越下限，说明无功设置较低，可适当提高无功设定值，或转换为PV 节点，电压定值取下限值。
如果电压越上限，说明节点无功设定偏高，可适当调低无功设定值，或转换为PV 节点，电压定值取上限值。
PQ 节点转换为PV 节点的处理方法：1、直角坐标方式的节点不平衡量由i Q ∆变为2i V ∆；2、牛顿法极坐标方式的修正方程减1个Q ∆方程；3、P-Q 分解法，θ-P 迭代不变，V Q -迭代的系数矩阵有两种处理方法：（1）在B ''中划去将要转换为PV 节点的节点所在的行和列，重新形成因子表。
（2）在B ''中将要转换为PV 节点的节点对应的对角元加一个很大的数，用因子表秩1修正法得到新的因子表三、因子表修正方法1、因子表秩1修正法设系数矩阵A 已因子化为如下的形式LDU A = （4-6）由于某种原因，A 变化为：A A N M A A ∆+=+=T a ~ （4-7）其中M 和N 为1⨯n 的列矢量，a 为标量。
新矩阵A ~的因子表为：U D L A ~~~~= （4-8）将（4-8）、（4-6）代入（4-7）有：T M aN LDU U D L +=~~~ （4-9）为了求出U D L ~~~中的各元素，将U D L ~~~和LDU 各矩阵的第一行和第一列单独列出，并写成分块矩阵的形式：⎥⎦⎤⎢⎣⎡=11L l L 1 ⎥⎦⎤⎢⎣⎡=1D D 1d ⎥⎦⎤⎢⎣⎡=1U u U 11 (4-10) 和⎥⎦⎤⎢⎣⎡=11L l L ~~1~ ⎥⎦⎤⎢⎣⎡=1D D ~~~1d ⎥⎦⎤⎢⎣⎡=1U u U ~~1~1 (4-11) 及⎥⎦⎤⎢⎣⎡=11M M m ⎥⎦⎤⎢⎣⎡=11N N n （4-12） 将（4-10）代入（4-6），A 矩阵可写为：⎥⎦⎤⎢⎣⎡+=11111111111U D L u l l u A d d d d （4-13） 将（4-11）代入（4-8），A ~矩阵可写为：⎥⎦⎤⎢⎣⎡+=11111111111~~~~~~~~~~~~U D L u l l u A d d d d （4-14） 将（4-12）代入T MaN A =∆⎥⎦⎤⎢⎣⎡=∆T T a an a m an m 11111111N M M N A （4-15） 将（4-13）、（4-14）、（4-15）代入（4-7）有⎥⎦⎤⎢⎣⎡+⎥⎦⎤⎢⎣⎡+=⎥⎦⎤⎢⎣⎡+T T a an a m an m d d d d d d d d 111111111111111111111111111111~~~~~~~~~~~N M M N U D L u l l u U D L u l l u （4-16） 根据等号两端矩阵对应元素相等，可得：（1） 1111~an m d d += （4-17）（2） T a m d d 111111~~N u u +=将（4-17）变为1111~an m d d -=代入，有 T a m d 111111~~~N u u -+= （4-18）其中1111~u N N n T T -= （4-19） （3） 111111~~an d d M l l +=将（4-17）变为1111~an m d d -=代入，有111111~~~-+=d an M l l （4-20）其中1111~m l M M -= （4-21）由上（1）、（2）、（3）可计算出新矩阵因子表上三角矩阵第一行元素、下三角矩阵第一列元素和对角线矩阵第一个元素。
（4） T a d d 11111111111111~~~~~~N M U D L u l U D L u l ++=+ 重写为111111111111111111~~~~~~A U D L N M U D L u l u l U D L ∆+=++-=T a d d（4-22） 其中T a d d 111111111~~~N M u l u l A +-=∆将（4-18）、（4-20）代入得 T T T T T T T T T T T T T T T T T T T T T T T T a a m d an a a m d an n a m a m d an an m a m a m d an a an a m an m a m d an a an m an an m a m an m a m d an a an m n a m an m a a m d an an a m an m a a m d d d an d d an a m d d d an m d a a m d d d an an m d a d d 111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111~~~~)~(~~~~)()(~~~)()(~~~~~~~~~)()(~~~~~~~~~~~~~~~~~~)~~(~)~~()~(~~~N M N M N M u N l M N M u l M N l M N M N M u M N l u l N M N M u l u M u l N l u l N M N M u l M u N l u l N M N M u M N l u l N M N M u M N l u l u l u l N M N u M l u l N M u l u l A =-=---=----=-+--=-++-+--=-+-----=+----=+-----=+++--=+-=∆------------- （4-23）其中)~(~1111a m d an a a --= （4-24）因此，（4-22）可写为如下的形式T a 11111111~~~~~~N M U D L U D L += （4-25） （4-25）与（4-9）有同样的形式，可用（1）、（2）、（3）的方法分别求出111~~~U D L 矩阵的第一行第一列元素。