电力系统牛拉法潮流计算
电力系统中的潮流计算与优化方法
电力系统中的潮流计算与优化方法潮流计算是电力系统运行和规划中的重要环节，它用于计算电力系统中各节点的电压、相角、有功、无功功率以及线路、变压器等的潮流分布情况。
对电力系统进行潮流计算可以帮助电力系统运行人员了解系统的稳定性、可靠性以及容载能力，也可以为电力系统规划提供数据支持。
本文将介绍电力系统潮流计算的基本方法与优化技术。
一、潮流计算的基本方法1.1 普通潮流计算方法潮流计算的基本方法是牛顿-拉夫逊迭代法（Newton-Raphson Iteration Method）和高尔顿法（Gauss-Seidel Method）。
牛顿-拉夫逊迭代法主要是通过不断迭代求解雅可比矩阵的逆，直到迭代误差小于给定阀值时停止迭代；高尔顿法则是逐一更新所有节点的电压与相角，直至所有节点的迭代误差都小于给定阀值。
1.2 快速潮流计算方法在大型电力系统中，普通的潮流计算方法计算速度较慢。
因此，研究人员提出了一些针对快速潮流计算的方法，如快速牛顿-拉夫逊法（Fast Newton-Raphson Method）和DC潮流计算方法。
快速牛顿-拉夫逊法通过简化牛顿-拉夫逊法的迭代公式，减少计算量，提高计算速度；DC潮流计算方法则是将潮流计算问题转化为一个线性方程组的求解问题，进一步提升计算效率。
二、潮流计算的优化技术2.1 改进的潮流计算算法为了提高潮流计算的准确性和收敛速度，研究人员提出了一些改进的潮流计算算法。
其中，改进的牛顿-拉夫逊法（Improved Newton-Raphson Method）是一种结合牛顿-拉夫逊法和割线法的算法，通过混合使用这两种方法，实现在减小迭代误差的同时加快计算速度。
此外，基于粒子群优化算法（Particle Swarm Optimization）和遗传算法（Genetic Algorithm）的潮流计算算法也得到了广泛研究和应用。
2.2 潮流优化潮流计算不仅可以用于分析电力系统的工作状态，还可以作为优化问题的约束条件。
电力系统潮流计算机算法
电力系统潮流计算机算法电力系统潮流计算是电力系统分析中最基本的一项计算，其目的是确定电力系统中各母线电压的幅值和相角、各元件中的功率以及整个系统的功率损耗等。
随着计算机技术的发展，电力系统潮流计算算法也在不断更新和完善。
以下是电力系统潮流计算的一些常用算法：1. 牛顿-拉夫逊法（Newton-Raphson Method）：这是一种求解非线性方程组的方法，应用于电力系统潮流计算中。
该方法在多数情况下没有发散的危险，且收敛性较强，可以大大节约计算时间，因此得到了广泛的应用。
2. 快速迪科法（Fast Decoupled Method）：这是一种高效的电力系统潮流计算方法，将电力系统分为几个子系统进行计算，从而提高了计算速度。
3. 最小二乘法（Least Squares Method）：这是一种用于求解线性方程组的方法，通过最小化误差平方和来获得最优解。
在电力系统潮流计算中，可用于优化电压幅值和相角。
4. 遗传算法（Genetic Algorithm）：这是一种全局优化搜索算法，应用于电力系统潮流计算中，可以解决一些复杂和非线性问题。
5. 粒子群优化算法（Particle Swarm Optimization）：这是一种启发式优化算法，通过模拟鸟群觅食行为来寻找最优解。
在电力系统潮流计算中，可用于优化网络参数和运行条件。
6. 模拟退火算法（Simulated Annealing）：这是一种全局优化搜索算法，应用于电力系统潮流计算中，可以在较大范围内寻找最优解。
7. 人工神经网络（Artificial Neural Network）：这是一种模拟人脑神经网络的计算模型，可用于电力系统潮流计算。
通过训练神经网络，可以实现对电力系统中复杂非线性关系的建模和预测。
以上所述算法在电力系统潮流计算中起着重要作用，为电力系统运行、设计和优化提供了有力支持。
同时，随着计算机技术的不断发展，未来还将出现更多高效、精确的电力系统潮流计算算法。
牛顿拉夫逊法计算潮流步骤
牛顿拉夫逊法计算潮流步骤牛顿拉夫逊法（Newton-Raphson method）是一种用于求解非线性方程组的迭代方法，它可以用来计算电力系统潮流的解。
潮流计算是电力系统规划和运行中的重要任务，它的目标是求解电力系统中各节点的电压幅值和相角，以及线路的功率流向等参数，用于分析电力系统的稳定性和安全性，以及进行电力系统规划和调度。
下面是使用牛顿拉夫逊法计算潮流的一般步骤：步骤1：初始化首先，需要对电力系统的各个节点（包括发电机节点和负荷节点）的电压幅值和相角进行初始化，一般可以使用其中一种估计值或者历史数据作为初始值。
步骤2：建立潮流方程根据电力系统的潮流计算模型，可以建立节点电压幅值和相角的平衡方程，一般采用节点注入功率和节点电压的关系来表示。
潮流方程一般是一个非线性方程组，包含了各个节点之间的复杂关系。
步骤3：线性化方程组将潮流方程组进行线性化处理，一般采用泰勒展开的方法，将非线性方程组变为线性方程组。
线性化的过程需要计算雅可比矩阵，即方程组中的系数矩阵。
步骤4：求解线性方程组利用线性方程组的求解方法，比如高斯消元法或LU分解法等，求解线性方程组，得到电压幅值和相角的修正量。
步骤5：更新节点电压根据线性方程组的解，更新各个节点的电压幅值和相角，得到新的节点电压。
步骤6：检查收敛性判断节点电压的修正量是否小于设定的收敛阈值，如果满足收敛条件，则停止迭代；否则，返回步骤3，循环进行线性化方程组和线性方程组的求解。
步骤7：输出结果当潮流计算收敛时，输出最终的节点电压幅值和相角，以及线路的功率流向等参数。
牛顿拉夫逊法是一种高效、快速且收敛性良好的方法，广泛应用于电力系统潮流计算。
在实际应用中，可能会遇到多次迭代或者收敛性不好的情况，此时可以采用退火技术或其他优化算法进行改进。
此外，牛顿拉夫逊法的计算也可以并行化，利用多核处理器或者分布式计算集群来加速计算过程。
总之，牛顿拉夫逊法是一种重要的潮流计算方法，通过迭代计算逼近非线性方程组的解，可以得到电力系统中各节点的电压幅值和相角，用于分析电力系统的稳定性和安全性。
关于直角坐标牛拉法系统潮流分布计算的答辩
关于直角坐标牛拉法系统潮流分布计算的答辩关于直角坐标牛拉法系统潮流分布计算的答辩一、引言直角坐标牛拉法是一种常用的电力系统潮流计算方法，用于计算电力系统中各节点的电压和功率分布。
本文将从算法原理、计算步骤、应用场景等方面进行阐述和答辩。
二、算法原理直角坐标牛拉法基于功率平衡方程和节点电压方程，通过迭代求解的方式，逐步逼近系统的潮流分布。
其核心思想是将电压和功率分别表示为实部和虚部，通过复数运算来求解未知量。
具体而言，直角坐标牛拉法将电流和导纳分别表示为复数形式，利用复数的乘法和除法运算，将节点电流和导纳联系起来，从而得到节点电压和功率的计算结果。
三、计算步骤直角坐标牛拉法的计算步骤包括以下几个部分：1. 初始化：给定电网拓扑结构、节点导纳和负荷信息，初始化节点电压和功率。
2. 潮流计算：根据功率平衡方程和节点电压方程，通过迭代计算节点电压和功率。
具体而言，每次迭代中，首先根据节点电压和导纳计算节点电流；然后，根据节点电流和导纳计算节点电压；再根据节点电压和导纳计算节点功率。
通过多次迭代，直到收敛为止。
3. 收敛判断：判断节点电压和功率的迭代计算是否收敛。
一般来说，可以通过判断节点电压和功率的变化量是否小于设定的收敛阈值来进行判断。
若满足收敛条件，则停止迭代；否则，继续迭代。
4. 输出结果：输出最终的节点电压和功率分布结果。
根据需要，还可以输出其他相关信息，如潮流方向、线路功率损耗等。
四、应用场景直角坐标牛拉法广泛应用于电力系统潮流计算和分析。
具体而言，它可以用于以下几个方面：1. 网络规划：通过潮流计算，可以评估电力系统的稳定性和可靠性，为电网规划提供依据。
例如，可以通过潮流计算来确定新建变电站的容量和位置，优化电网结构。
2. 运行调度：在电力系统的日常运行中，潮流计算可以用于实时监测和调度。
通过潮流计算，可以了解各节点的电压和功率情况，及时发现问题并采取措施，确保电力系统的安全稳定运行。
3. 短路分析：在电力系统发生短路故障时，潮流计算可以用于分析故障电流的分布情况，确定故障点和故障线路，为故障处理和保护调整提供参考。
电力系统专业课程设计牛顿拉夫逊法潮流计算
课程设计阐明书电力系统分析系（部）专业（班级）姓名指引教师起止日期电力系统分析课程设计任务书系（部）：专业：指引教师：系（部）主 管领导意 见年 月曰教研室 意见目录一、潮流计算基本原理1.1潮流方程基本模型1. 2潮流方程讨论和肖点类型划分1.3、潮流计算意义二、牛顿一拉夫逊法2. 1牛顿-拉夫逊法基本原理2.2节点功率方程2.3修正方程2. 4牛顿法潮流计算重要流程三、收敛性分析四、算例分析总结参照文献电力系统分析潮流计算一、潮流计算基本原理1.1潮流方程基本模型电力系统是由发电机、变压器、输电线路及负荷等构成，苴中发电机及负荷是非线性元件，但在进行潮流计算时，普通可以用接在相应节点上一种电流注入量来代表。
因而潮流计算所用电力网络系由变压器、输电线路、电容器、电抗器等静止线性元件所构成，并用集中参数表达串联或并联等值支路来模仿。
结合电力系统特点，对这样线性网络进行分析，普通采用是节点法，巧点电压与节点电流之间关系j = YV(1-1) 其展开式为Yn’j将式(1-6).式(1-7)代入以导纳矩阵为基本式(1 一4),并将实部与虚某些开，可以得到如下 两种形式潮流方程。
潮流方程直角坐标形式为P,=勺工9“勺 _◎/；)+£工(G/ + 竝勺)(Z = 1,2,3,…屮)(1-8)Q = f,工(G 泸j - 3』)f 为(G/ + 坊勺)a = 1,2,3,…(1-9)潮流方程极坐标形式为A 岭匕(i = 1,2,3,…，”) ;-1(1-2)在工程实际中，已经节点注入量往往不是节点电流而是节点功率，为此必要应用联系 肖点电流和节点功率关系式I, = P,\，Q,(，= 1,2,3,…，”)V/(1-3)将式(1-3)代入式(1-2)得到竺理=乞型(山123,…V/ >='(1-4)交流电力系统中复数电压变量可以用两种极坐标来表达(1-5)或而复数导纳为y i =e i +jf i(1 —6)(1-7)R =匕》匕(G» cos© + B. sin 0 ) (/ = 1,2,3, • •(1 一10)Qi =匕》V； (G, sin 0.- B3 cosq) (/ = 1,2,3,…,”)(1-11) 以上各式中，j门表达工号后标号门丫点必要直接和肖点r相联，并涉及)=,状况。
牛顿拉夫逊法潮流计算
牛顿拉夫逊法潮流计算牛顿-拉夫逊法（Newton-Raphson method）是一种用于求解非线性方程的数值方法。
它通过迭代逼近根的方式，将非线性方程转化为一系列的线性方程来求解。
在电力系统中，潮流计算用于确定电力网中节点的电压幅值和相角。
潮流计算是电力系统分析的重要基础，可以用于计算电力系统的潮流分布、功率损耗、节点电压稳定度等参数，为电力系统的规划、运行和控制提供参考依据。
牛顿-拉夫逊法是一种常用的潮流计算方法，它的基本思想是通过不断迭代来逼近电网的潮流分布，直到满足一定的收敛条件。
下面将对牛顿-拉夫逊法的具体步骤进行详细介绍。
首先，我们需要建立电力网络的节点潮流方程，即功率方程。
对于每一个节点i，其节点功率方程可以表示为：Pi - Vi * (sum(Gij * cos(θi - θj)) - sum(Bij * sin(θi -θj))) = 0Qi - Vi * (sum(Gij * sin(θi - θj)) + sum(Bij * cos(θi -θj))) = 0其中，Pi和Qi分别为节点i的有功功率和无功功率，Vi和θi分别为节点i的电压幅值和相角，Gij和Bij分别为节点i和节点j之间的导纳和电纳。
接下来，我们需要对每个节点的电压幅值和相角进行初始化。
一般情况下，可以将电压幅值设置为1，相角设置为0。
然后，我们可以开始进行迭代计算。
在每一轮迭代中，我们需要计算每个节点的雅可比矩阵和功率残差，然后更新电压幅值和相角。
雅可比矩阵可以通过对节点功率方程进行求导得到，具体如下：dPi/dVi = -sum(Vj * (Gij * sin(θi - θj) + Bij * cos(θi - θj)))dPi/dθi = sum(Vj * (Gij * Vi * cos(θi - θj) - Bij * Vi * sin(θi - θj)))dQi/dVi = sum(Vj * (Gij * cos(θi - θj) - Bij * sin(θi - θj)))dQi/dθi = sum(Vj * (Gij * Vi * sin(θi - θj) + Bij * Vi * cos(θi - θj)))功率残差可以通过将节点功率方程代入得到，如下：RPi = Pi - Vi * (sum(Gij * cos(θi - θj)) - sum(Bij *sin(θi - θj)))RQi = Qi - Vi * (sum(Gij * sin(θi - θj)) + sum(Bij *cos(θi - θj)))最后，我们可以使用牛顿-拉夫逊法的迭代公式来更新电压幅值和相角，具体如下：Vi(new) = Vi(old) + ΔViθi(new) = θi(old) + Δθi其中，ΔVi和Δθi分别为通过求解线性方程组得到的电压幅值和相角的增量。
(整理)电力系统分析课程设计——牛顿拉夫逊潮流计算
目录电力系统编程潮流计算 (2)1 设计任务及初步分析 (2)1.1 设计任务 (2)1.2 初步分析 (2)2 牛顿-拉夫逊法简介 (3)2.1概述 (3)2.2 一般概念 (3)2.3 潮流计算的修正方程 (4)2.4 直角坐标表示的修正方程 (4)3 程序设计 (8)3.1 程序流程图 (8)3.2 源程序 (8)3.3 结果及分析 (16)4 小结 (17)参考文献 (18)电力系统编程潮流计算1 设计任务及初步分析1.1 设计任务条件：节点数：3 支路数：3 计算精度：0.00010支路1： 0.0300+j0.09001┠—————□—————┨2支路2： 0.0200+j0.09002┠—————□—————┨3支路3： 0.0300+j0.09003┠—————□—————┨1节点1：PQ 节点，S （1）=-0.5000-j0.2000 节点2：PQ 节点，S （2）=-0.6000-j0.2500 节点3：平衡节点，U （3）=1.0000∠0.0000 要求：编写程序计算潮流1.2 初步分析潮流计算在数学上可归结为求解非线性方程组，其数学模型简写如下：⎪⎪⎩⎪⎪⎨⎧=⋅⋅⋅⋅⋅⋅=⋅⋅⋅=⋅⋅⋅0)(0)(0)(21212211n n n n x x x f x x x f x x x f ，，，，，，，，，2 牛顿-拉夫逊法简介2.1概述牛顿－拉夫逊法是目前求解非线性方程最好的一种方法。
这种方法的特点就是把对非线性方程的求解过程变成反复对相应的线性方程求解的过程，通常称为逐次线性化过程，就是牛顿－拉夫逊法的核心。
牛顿-拉夫逊法的基本原理是在解的某一邻域内的某一初始点出发，沿着该点的一阶偏导数——雅可比矩阵，朝减小方程的残差的方向前进一步，在新的点上再计算残差和雅可矩阵继续前进，重复这一过程直到残差达到收敛标准，即得到了非线性方程组的解。
因为越靠近解，偏导数的方向越准，收敛速度也越快，所以牛顿法具有二阶收敛特性。
第四章复杂电力系统潮流计算-牛顿-拉夫逊潮流计算
ΔX
(k )
?
Yes
收敛结束
极坐标形式的潮流方程
* * I U Y U Pi jQi Ui i i ij j * j 1 n
U i U i i，U j U j j Yij Gij jBij , ij i j
电压相量用 极坐标表示
极坐标下有功功率和无功功率方程
n Pi U iU j (Gij cos ij Bij sin ij ) j 1 n Q U U (G sin B cos ) ij ij ij i i j ij j 1
i 1, 2, , n
泰勒级数展开忽略步时的修正方程组为修正量修正方程的矩阵形式其中函数fx的jocabi雅可比矩阵收敛结束yesijijij极坐标下有功功率和无功功率方程电压相量用极坐标表示次迭代时pq节点
§3-4 牛顿-拉夫逊法潮流计算 （Newton-Raphson迭代法）
牛顿-拉夫逊法
单变量非线性方程
真解
f ( x ) 0,
j 1 n
n
U i U j [(Gij cos ij Bij sin ij ) j (Gij sin ij Bij cos ij )
j 1 n
U i U j (Gij cos ij Bij sin ij ) jU i U j (Gij sin ij Bij cos ij )
x( 0 )
1
x( 0 )
1
x1 x( 0 ) n x2 f n xn (0) x n x n f 1 x n
迭代至第 k 步时的修正方程组为
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
10
直角坐标下潮流方程
直角坐标下待求变量
直角坐标下功率方程
P 1 Pn Q 1 f ( x) Qn r V 2 n r 1 V 2 n
11
e1 en x f1 fn
G cos G sin B cos B sin
28
定Jacobian算法
考虑到正常情况下，ij 很小 节点自导纳要远大于节点注入功率
B G J J0 G B
'
则定Jacobian矩阵的潮流计算修正方程为
BH GM GN V P / V BL V Q / V
P T e Q f J T T x e 2 V eT
V 2 f T P f T Q f T
22
极坐标下牛顿-拉夫逊方法
P SP P(V , ) P(V , ) f ( x) SP Q ( V , ) Q Q ( V , )
电力系统潮流计算（1）
华北电力大学电气与电子工程学院 孙英云 Email: sunyy@ 办公室：教五 C204
1
问题
什么是潮流计算？
什么是潮流？ 什么是计算？ 电力系统状态不可直接测量 潮流和电力系统运行状态的关系 电力系统分析、计算的需要
为什么要进行潮流计算？
忽略高阶项，则有
f ( x ) f ( x )x
0
18
牛顿-拉夫逊法潮流计算
牛顿法的几何意义
19
牛顿-拉夫逊法潮流计算
牛顿法计算流程 1 初始化，形成节点导纳阵，给出初值 2 令k=0 进入迭代循环
(k )
x (0)
2.1 计算函数值 f ( x(k ) ) ，判断是否收敛 f ( x ) 2.2 计算Jacobian矩阵 f ( x( k ) ) 2.3 计算修正量 x( k ) (f ( x( k ) ))1 f ( x( k ) ) 2.4 对变量进行修正 x( k 1) x( k ) x(k ) ，k=k+1返回 2.1
H J M R N L S
极坐标
Pi Vi V j (Gij cos ij Bij sin ij ) ji Qi Vi V j (Gij sin ij Bij cos ij ) ji
H J M
N L
3 输出计算结果
20
牛顿-拉夫逊法潮流计算
牛顿法可写成如下简单迭代格式
x( k 1) x(k ) ( J ( x(k ) ))1 f ( x(k ) ) ( x(k ) )
( x) J 1 J 1 1 ( x) ( x) I T f ( x) J T f ( x) T T x x x x
7
极坐标功率平衡方程
V 如果将节点电压用直角坐标表示，即令 V i i i 则有：
Pi jQi Vi i (Gij jBij )V j j
ji
=Vi (Gij jBij )(cos ij j sin ij )
ji
i 1, 2, N
极坐标潮流方程的已知量和待求量？
13
潮流方程的解法
潮流方程是一组高维非线性方程组 所有能用于求解非线性方程组的方法都可以用 于求解潮流方程
Gauss法（简单迭代法） Newton法（包括其变形算法） 割线法 拟牛顿法 ……
14
以Gauss法为基础的潮流方 程解法
待求方程 高斯迭代法
P T J Q T
P V T V Q V T V
P T Q T
P V T P V Q V Q V V V T
26
潮流计算速度
目前的主流潮流计算算法都是迭代算法
计算时间=迭代次数×每次迭代所需计算时间
提高计算速度的两条思路
减少迭代次数 高阶收敛性算法 减少每次迭代所需时间 定Jacobian方法
27
定Jacobian算法
2 V 将极坐标Jacobian矩阵中的 移出矩阵
' V H P ' VQ M
Pi Vi V j (Gij cos ij Bij sin Bij ) ji Qi Vi V j (Gij sin ij Bij cos Bij ) ji
i 1, 2, N i 1, 2, N
8
潮流方程的讨论和节点类型 的划分
对于电力系统来讲，每个节点有四个运行变量 （电压×2，功率×2），两个功率平衡方程 （有功、无功） 负荷节点
N ' VP P V ' Q VQ L V
B sin Q P G sin P Q
H' J ' M
N ' B cos ' L G cos
17
牛顿-拉夫逊法潮流计算
牛顿法的历史 牛顿法基本原理
对于非线性方程 f ( x) 0 给定初值 x (0) 用Talor级数展开，有：
(0) x f ( x (0) x (0) ) f ( x (0) ) f ' ( x (0) )x (0) f '' ( x (0) ) 2! 0 (0) ' (0) (0)
令 则有
I YV Yn V n n s s
Yn L + D + U
= D-1 (I -YV - LV - UV ) V n n s s n n
ˆ i 1 n 1 S ( k ) ( k ) Y V YV i YisV s ij j ij j ˆ Yii V (k ) j 1 j i 1 i i 1, 2,, n
( x )的谱半径趋近于0，因此 随着迭代的进行， 越接近收敛点，牛顿法收敛越快，具备局部二 阶收敛性
21
直角坐标下牛顿-拉夫逊方 法
P(e, f ) P SP P(e, f ) SP f ( x) Q(e, f ) Q Q(e, f ) V 2 (e, f ) (V SP ) 2 V 2 (e, f )
直角坐标下潮流方程Pi P Nhomakorabea SP (ei ai f i bi ) 0 Qi QiSP ( fi ai ei bi ) 0 Vi 2 (Vi SP ) 2 (ei2 fi 2 ) 0
直角坐标潮流方程的已知量和待求量？
12
极坐标潮流方程
Pi Vi V j (Gij cos ij Bij sin Bij ) ji Qi Vi V j (Gij sin ij Bij cos Bij ) ji
f ( x) 0
x ( x)
x(0) x0
x( k 1) ( x( k ) )
当矩阵的谱半径小于1时收敛，谱半径越小， 收敛性越好
( x* ) ( x) xT x x*
15
基于节点导纳矩阵的高斯迭 代法
I Yn Ys V n n YT Y ss Vs s Is
16
( k 1) V i
高斯法的讨论
高斯法可分为基于节点导纳阵的高斯法和基于 阻抗阵的高斯法两种 高斯法的改进 高斯-赛德尔法 高斯法的PV节点处理较为困难
具体可参见 Kusic G L. Computer-aided power systems analysis. Prentice Hall, 1986
P T J Q T
P V T Q V T
23
极坐标下牛顿-拉夫逊法
为了使Jacobian矩阵中对电压的偏导项恢复为 关于V的二次函数，在对V的偏导项处乘以一 个V，在V的修正项中除以一个V，则有
x V V
3
简单电力系统等值电路（实例）
输电线路
发电机 升压变压器
配电线路
降压变压器
G
T1
L1
T2
L2
降压变压器
T3
负荷
K1ZT1
G
Z110 Z120
ZL1
K2ZT2
ZL2
K3ZT3
PD+jQD
YL1/2
YL1/2
Z210
Z220
YL2/2
YL2/2
Z310
Z320
4
电力系统稳态模型
发电机
出力可调，机端电压可控：PV或平衡节点
ji
6
直角坐标功率平衡方程
e jf 如果将节点电压用直角坐标表示，即令 V i i i 则有：
Pi jQi (ei jfi ) (Gij jBij )(e j jf j )
ji
(ei jfi )(ai jbi )
i 1, 2, N
ai (Gij e j Bij f j ) ji bi (Gij f j Bij e j ) ji i 1, 2, N Pi ei ai fi bi i 1, 2, N Qi fi ai ei bi