1、1,第2章 電力系統(tǒng)潮流計算,概述 2.1 潮流計算的數(shù)學(xué)模型 2.2 潮流計算的導(dǎo)納矩陣和阻抗矩陣迭代法 2.3 潮流計算的N-R法 2.4 潮流計算的P-Q分解法 2.5 保留非線性的潮流算法 2.6 最優(yōu)乘子法（病態(tài)潮流算法） 2.7 潮流迭代算法的收斂性問題,2,概 述 作為研究電力系統(tǒng)穩(wěn)態(tài)運(yùn)行情況的一種基本電氣計算，電力系統(tǒng)常規(guī)潮流計算的任務(wù)是根據(jù)給定的網(wǎng)絡(luò)結(jié)構(gòu)及運(yùn)行條件（網(wǎng)絡(luò)結(jié)構(gòu)包括線路、變電站、電源點(diǎn)的位置等；運(yùn)行條件是指負(fù)荷的大小及電源出力等），求出整個網(wǎng)絡(luò)的運(yùn)行狀態(tài)，其中包括各母線的電壓、網(wǎng)絡(luò)中的功率分布以及功率損耗等等。,3,潮流計算的應(yīng)用分為離線和在線應(yīng)用兩類。 離線應(yīng)用

2、主要有： 規(guī)劃及運(yùn)行規(guī)劃研究 靜態(tài)及暫態(tài)穩(wěn)定計算 故障分析及優(yōu)化計算 在線應(yīng)用主要有： 隨著現(xiàn)代化的調(diào)度控制中心的建立，為了對電力系統(tǒng)進(jìn)行實(shí)時安全監(jiān)控，需要根據(jù)實(shí)時數(shù)據(jù)庫所提供的信息，隨時判斷系統(tǒng)當(dāng)前的運(yùn)行狀態(tài)并對預(yù)想事故進(jìn)行安全分析，這就需要進(jìn)行廣泛的潮流計算，并且對計算速度等還提出了更高的要求，從而產(chǎn)生了潮流的在線計算。,4,由上可見，潮流計算是電力系統(tǒng)中應(yīng)用最為廣泛、最基本和最重要的一種電氣計算 潮流計算問題在數(shù)學(xué)上一般是屬于多元非線性代數(shù)方程組的求解問題，必須采用迭代計算方法。 自從50年代中期開始利用電子計算機(jī)進(jìn)行潮流計算以來，潮流計算是電力系統(tǒng)各種問題中投入研究力量最多的領(lǐng)域之一，

3、出現(xiàn)了大量的研究成果。這些成果： 開拓了各種特殊性質(zhì)的潮流計算問題 更多的是屬于為了提高計算性能而陸續(xù)提出的各種具體算法。,5,對于一個潮流算法，其基本要求可歸納成以下四個方面 (1)計算速度； (2)計算機(jī)內(nèi)存占用量； (3)算法的收斂可靠性； (4)程序設(shè)計的方便性以及算法擴(kuò)充移植 等的通用靈活性。 這四點(diǎn)要求也成為本章后面評價各種潮流算法性能時所依據(jù)的主要標(biāo)準(zhǔn)。,6,本章在對潮流計算問題的數(shù)學(xué)模型進(jìn)行簡單的回顧以后，將首先轉(zhuǎn)入三種最基本的潮流算法： 高斯一塞德爾法 牛頓法 快速解耦法的討論 這三種算法的基本原理在大學(xué)本科的電力系統(tǒng)分析教材中已作過介紹，但鑒于這些方法的重要性，將在大學(xué)本科

4、電力系統(tǒng)分析教材的基礎(chǔ)上作進(jìn)一步的討論。,7,牛頓法的特點(diǎn)是將非線性方程線性化。70年代后期，有人提出采用更精確的模型，即將泰勒級數(shù)的高階項也包括進(jìn)來，希望以此提高算法的性能，這便產(chǎn)生了保留非線性的潮流算法。 為了解決病態(tài)潮流計算，出現(xiàn)了將潮流計算表示為一個無約束非線性規(guī)劃問題的模型，并稱之為最小化潮流計算法。,8,一些實(shí)際用于生產(chǎn)的潮流程序往往在上述基本潮流的框架內(nèi)再加入模擬實(shí)際系統(tǒng)運(yùn)行控制特點(diǎn)的自動調(diào)整計算功能，如潮流控制，分接頭調(diào)整等，這部分內(nèi)容將在本章簡要介紹。,9,60年代中期，結(jié)合電力系統(tǒng)經(jīng)濟(jì)調(diào)度工作的開展，針對經(jīng)典的經(jīng)濟(jì)調(diào)度方法的不足，開辟了一個新的研究領(lǐng)域，稱之為最優(yōu)潮流問題。

5、這種以非線性規(guī)劃作為計算模型的潮流問題能夠統(tǒng)籌兼顧電力系統(tǒng)的經(jīng)濟(jì)性、安全性和電能質(zhì)量，因而受到很大的重視，發(fā)展很快，其應(yīng)用領(lǐng)域正在不斷擴(kuò)大。我們將在本章第3章中簡要討論。,10,2.1 潮流計算問題的數(shù)學(xué)模型,計算網(wǎng)絡(luò)的結(jié)構(gòu)： 網(wǎng)絡(luò)只包含L、T等線性元件，用 Y 或 Z 描述； Load、G 等非線性元件引出網(wǎng)絡(luò)之外，用注入網(wǎng)絡(luò)的功率 或 電流 描述； 聯(lián)絡(luò)節(jié)點(diǎn)視為帶有零注入功率的負(fù)荷。,11,2.1.1 節(jié)點(diǎn)方程和節(jié)點(diǎn)分類,2.1 潮流計算問題的數(shù)學(xué)模型,(1) Y 描述的節(jié)點(diǎn)電壓方程 I=YV,(2) Z 描述的節(jié)點(diǎn)電壓方程 V=ZI,12,2.1.1 節(jié)點(diǎn)方程和節(jié)點(diǎn)分類,2.1 潮流計算

6、問題的數(shù)學(xué)模型,(3) 節(jié)點(diǎn)分類,13,2.1.2 節(jié)點(diǎn)功率方程,2.1 潮流計算問題的數(shù)學(xué)模型,一般形式,14,2.1.2 節(jié)點(diǎn)功率方程,2.1 潮流計算問題的數(shù)學(xué)模型,(1) 極坐標(biāo)形式的節(jié)點(diǎn)功率方程,15,2.1.2 節(jié)點(diǎn)功率方程,2.1 潮流計算問題的數(shù)學(xué)模型,(2) 直角坐標(biāo)形式的節(jié)點(diǎn)功率方程,16,2.1.3 潮流計算中的幾個概念問題,2.1 潮流計算問題的數(shù)學(xué)模型,(1) 潮流計算需解決的基本問題和途徑,需解決的基本問題： 已知系統(tǒng)條件網(wǎng)絡(luò)結(jié)構(gòu)、負(fù)荷與出力水平及其分布， 確定系統(tǒng)運(yùn)行狀態(tài)節(jié)點(diǎn)電壓、支路潮流、系統(tǒng)網(wǎng)損等 解決的問題的途徑：求解潮流方程，即 節(jié)點(diǎn)功率平衡方程：,高階、

7、非線性代數(shù)方程組！ f(x)=0 迭代求解,(2) 潮流迭代的收斂條件,17,2.1.3 潮流計算中的幾個概念問題,2.1 潮流計算問題的數(shù)學(xué)模型,(3) 支路潮流分布的計算方法,18,2.1.3 潮流計算中的幾個概念問題,2.1 潮流計算問題的數(shù)學(xué)模型,(4) 網(wǎng)損的計算方法,線路阻抗支路功率損耗,變壓器阻抗支路功率損耗,全網(wǎng)功率損耗,(5) 給定節(jié)點(diǎn)負(fù)荷功率的處理2種基本情況 a) 負(fù)荷為恒定功率特性： Pis=Pis0=const ; Qis=Qis0=const b) 計及負(fù)荷的靜態(tài)電壓特性：,19,2.2.1 導(dǎo)納矩陣迭代法,2.2 導(dǎo)納矩陣迭代法和阻抗矩陣迭代法,(1) Gauss

8、-Seidel 導(dǎo)納矩陣迭代法的迭代格式,YV=I F(V)=YV-I=0,D-1F(V)=D-1(YV-I)=0,V=(E-D-1Y)V+D-1I,概述： Y、Z矩陣迭代法是求解非線性代數(shù)方程組最簡便的迭代算法；其中簡單迭代法計算流程簡單，但收斂性差，故常用的是其改進(jìn)迭代算法Gauss-Seidel 迭代算法。,20,2.2.1 導(dǎo)納矩陣迭代法,(2) 不同類型節(jié)點(diǎn)的處理,(i) PQ節(jié)點(diǎn)直接計算：,(ii) PV節(jié)點(diǎn) Q I V V：,迭代開始時：,21,2.2.1 導(dǎo)納矩陣迭代法,(2) 不同類型節(jié)點(diǎn)的處理,(iii) 平衡節(jié)點(diǎn) 只要在迭代收斂后計算節(jié)點(diǎn)注入功率,(3) 迭代步驟,22,

9、2.2.2 阻抗矩陣迭代法,(1) 簡單迭代法 基本公式：,主要問題：收斂慢,23,2.2.2 阻抗矩陣迭代法,(2) 簡單迭代法的改進(jìn) 塞德爾迭代法 目的：加速收斂,24,2.2.2 阻抗矩陣迭代法,(3) 不同類型節(jié)點(diǎn)的處理, PQ節(jié)點(diǎn)已知：Pis、Qis迭代過程中直接計算 計算待求狀態(tài)量 V, PV 節(jié)點(diǎn)已知：Pis、Vis 節(jié)點(diǎn)電壓(相位)、Q為待求變量 a) 求得 V 之后，必須依已知的 Vis 進(jìn)行修正 b) 由修正后的 V 計算電流變化量，進(jìn)而計算節(jié)點(diǎn)Q c) 為改善收斂性，I 變化量中，可計及平衡節(jié)點(diǎn) (R) 電流變化的影響 d) 迭代開始時，PV 節(jié)點(diǎn)之 Q 值可假定，通常取

10、 (1/2)Pis , V 節(jié)點(diǎn)(R)已知： 每次迭代必須計算節(jié)點(diǎn)電流IR 正是其不斷變化來平衡系統(tǒng)功率：,25,2.2.2 阻抗矩陣迭代法,(4) 完整計算流程Gauss-Seidel 阻抗矩陣迭代法,26,2.2.3 算法性能分析,(1) Gauss-Seidel 迭代法的一般描述,簡單迭代法的一般描述 非線性代數(shù)方程組： F(X)=0，求解 X 令 g(X)=X-F(X) 則 X=g(X) X(t+1)=g(X(t)迭代格式簡單迭代法 展開，即有,特點(diǎn)：簡單，但收斂慢改進(jìn)以加速收斂Gauss Seidel 迭代法,27,2.2.3 算法性能分析,(1) Gauss-Seidel 迭代法的

11、一般描述,(b) Gauss-Seidel 迭代法的一般描述,28,2.2.3 算法性能分析,(2) 收斂性分析,29,2.2.3 算法性能分析,(3) 阻抗矩陣迭代與導(dǎo)納矩陣迭代的性能比較,導(dǎo)納矩陣迭代法：簡單、每次迭代計算量小，內(nèi)存占用少； 收斂性差，迭代次數(shù)多,阻抗矩陣迭代法：收斂性比Y矩陣迭代好； 每次迭代計算量較大，內(nèi)存占用多；,30,(1) 單變量非線性代數(shù)方程的牛頓迭代算法 修正方程 設(shè) f (x)=0 x(0) 為近似解 與真解的誤差：x(0) 真解： x=x(0) -x(0) 則 f (x(0) - x(0)=0 f (x) 在 x(0) 附近泰勒展開：,2.3 潮流計算的牛

12、頓-拉夫遜法,2.3.1 牛頓法的基本概念,忽略高次項,31,(1) 單變量非線性代數(shù)方程的牛頓迭代算法 迭代格式： k=0 x(0)：迭代初值； x(0) ：初始修正量 修正方程： 修正量： 新的近似解：,2.3.1 牛頓法的基本概念,32,(2) n 階非線性方程組的牛純法： 修正方程：F(x) =f1(x) , f2(x) , , fn(x)T =0 ； x =x1 , x2 , , xnT =0,2.3.1 牛頓法的基本概念,33,(2) n 階非線性方程組的牛純法： 修正方程：矩陣形式,34,(2) n 階非線性方程組的牛純法： 迭代格式： 修正方程： F(x(k)-J (k) x(

13、k)=0 Or F(x(k) J (k) x(k) 修正量： x(k) J(k)-1 F(x(k) 新的近似解： x(k+1)x(k) - x(k) 其中 F(x(k)f1(x(k) , f2(x(k) , , fn(x(k)T x(k)x1(k), x2(k) , , xn(k)T x(k)x1(k), x2(k) , , xn(k)T x(k+1)x1(k+1), x2(k+1) , , xn(k+1)T 雅可比矩陣：,35,(2) n 階非線性方程組的牛純法：,2.3 潮流計算的牛頓-拉夫遜法,2.3.1 牛頓法的基本概念,36,(1) 節(jié)點(diǎn)功率平衡方程,2.3 潮流計算的牛頓-拉夫遜法

14、,2.3.2 直角坐標(biāo)形式的 N-R 算法,PQ節(jié)點(diǎn)： i=1,2, , m,PV節(jié)點(diǎn)： i=m+1, , n-1,(2) 修正方程矩陣形式, W = J V,設(shè)： 節(jié)點(diǎn)總數(shù) n、PV 節(jié)點(diǎn) 數(shù) r 、平衡節(jié)點(diǎn)數(shù) 1PQ節(jié)點(diǎn)數(shù) m=n-r-1 節(jié)點(diǎn)編號：PQ(1,2,m)、PV(m+1,n-1)、平衡節(jié)點(diǎn)(n),37,2.3.2 直角坐標(biāo)形式的 N-R 算法,(2) 修正方程各向量(矩陣)結(jié)構(gòu),PQ節(jié)點(diǎn)的Jij (i=1,2,m; j=1,2, n-1),PV節(jié)點(diǎn)的Jij (i=m+1, ,n-1; j=1,2, n-1),38,2.3.2 直角坐標(biāo)形式的 N-R 算法,(2) 修正方程各向量

15、(矩陣)結(jié)構(gòu),39,ij 非對角元素： i=j 對角元素：,2.3.2 直角坐標(biāo)形式的 N-R 算法,(3) 雅可比矩陣元素的計算,40,2.3.2 直角坐標(biāo)形式的 N-R 算法,(3) 雅可比矩陣特點(diǎn),41,(4) 潮流計算步驟：,2.3.2 直角坐標(biāo)形式的 N-R 算法,42,(5) 提高計算精度的措施,2.3.2 直角坐標(biāo)形式的 N-R 算法,原修正方程中，J 矩陣對應(yīng) 節(jié)點(diǎn) i 的對角子塊：,改變 V 的排序,改變 W 的排序,43,2.3.3 極坐標(biāo)形式的 N-R 潮流算法,(1) 節(jié)點(diǎn)功率平衡方程：,設(shè)： 節(jié)點(diǎn)總數(shù) n、PV 節(jié)點(diǎn) 數(shù) r 、平衡節(jié)點(diǎn)數(shù) 1PQ節(jié)點(diǎn)數(shù) m=n-r-1

16、 節(jié)點(diǎn)編號：PQ(1,2,m)、PV(m+1,n-1)、平衡節(jié)點(diǎn)(n), 有功平衡方程 n-1個，不包含平衡節(jié)點(diǎn), 無功平衡方程m=n-r-1 個，僅對應(yīng)PQ節(jié)點(diǎn),44,2.3.3 極坐標(biāo)形式的 N-R 潮流算法,(2) 修正方程：,修正方程的簡單形式：,45,(3) 雅可比矩陣結(jié)構(gòu)及其元素計算：,2.3.3 極坐標(biāo)形式的 N-R 潮流算法,46, 階數(shù)2PQ節(jié)點(diǎn)數(shù) + PV節(jié)點(diǎn)數(shù)n-1 + m 比較：“直角坐標(biāo)” J的階數(shù)2(PQ節(jié)點(diǎn)數(shù) + PV節(jié)點(diǎn)數(shù))=2(n-1) 極坐標(biāo)形式 比 直角坐標(biāo)形式 階數(shù)減小n-m-1 ) J元素是節(jié)點(diǎn)電壓的函數(shù)，數(shù)值在迭代過程中變化 J 與Yij有關(guān)， if

17、 Yij0 then 對應(yīng)元素0 J的每一分快陣(H、N、K、L) 與Y有相同的稀疏結(jié)構(gòu) JJT,2.3.3 極坐標(biāo)形式的 N-R 潮流算法,(3) 雅可比矩陣結(jié)構(gòu)及其元素J矩陣的特點(diǎn),47,2.3.3 極坐標(biāo)形式的 N-R 潮流算法,48,2.4 潮流計算的P-Q 分解法,2.4.1 P-Q 分解法的基本原理,49,(3) H、L的常數(shù)化表示：,2.4.1 P-Q 分解法的基本原理,50,2.4.2 P-Q 分解法的修正方程,B,B,51,2.4.3 P-Q 分解法 的計算過程,52,2.4.3 P-Q 分解法的計算過程,53,2.4.4 P-Q 分解法的應(yīng)用注意,B、B 的計算 為加速收斂

18、，在 B 中忽略對有功、電壓相角影響較小的因素； 在 B 中忽略對無功、電壓幅值影響較小的因素。 兩種基本的常用方法 BX 法、XB 法,BX 法：,XB 法：,54,2.4.4 P-Q 分解法的應(yīng)用注意,(2) 算法的收斂性 P-Q 分解法按幾何級數(shù)收斂， N-R 法為平方收斂 值得注意： 對元件參數(shù) r/x 較高的系統(tǒng) (如配電網(wǎng)絡(luò))，P-Q分解法可能不收斂，需要對算法改進(jìn)。,55,2.5 保留非線性的潮流算法,2.5.1 基本思路,直角坐標(biāo)潮流方程是關(guān)于節(jié)點(diǎn)電壓的二次方程泰勒展開并保留非線性項，即為無截斷誤差的精確表達(dá)式,56,2.5.2 保留非線性潮流算法的模型描述,(1) 二次齊次方

19、程,57,2.5.2 保留非線性潮流算法的模型描述,設(shè) 函數(shù)向量 y(x) 及其 給定的平衡點(diǎn)xs 的 函數(shù)向量值 ys ,設(shè) f(x) 為任意解 x 對應(yīng)的 不平衡量： f(x)= y(x)-ys 對于平衡點(diǎn)x=xs，應(yīng)當(dāng)滿足: y(x) = ys 即，在解點(diǎn)應(yīng)有： f(x)=y(x) - ys =0,rs AiAi 的行展開,(2) 二次齊次方程組及其平衡方程,58,2.5.2 保留非線性潮流算法的模型描述,將 yi(x) 在 x(0) 展開：,(3) 二次齊次方程組的泰勒展開,59,2.5.2 保留非線性潮流算法的模型描述,(4) 迭代公式,60,2.5.2 保留非線性潮流算法的模型描述

20、,(5) 注意的若干問題, x 是 y(x)的函數(shù)只能迭代求得滿足y(x)-ys=0的解； J對應(yīng)于x(0)，在迭代過程中保持不變； y(x)由y(x)表達(dá)式直接計算，且通常令x (0)=0以計算y(x(0); 收斂準(zhǔn)則：,61,2.5.2 保留非線性潮流算法的模型描述,(5) 注意的若干問題, 與N-R的比較： 比N-R法的計算速度快幾倍至10倍！,62,2.5.2 保留非線性潮流算法的模型描述,(5) 注意的若干問題, 與潮流迭代的對應(yīng)關(guān)系關(guān)鍵問題！,63,2.6 病態(tài)條件下的潮流計算最優(yōu)乘子法,2.6.1 問題的提出,(1) 正常條件與病態(tài)條件 正常條件 (Well-condition)

21、 潮流迭代收斂 病態(tài)條件(ill-condition) 潮流迭代發(fā)散 常規(guī)潮流算法不能解決問題 (2) 非線性規(guī)劃潮流算法的基本思路與特點(diǎn) 將潮流計算問題當(dāng)作非線性規(guī)劃問題，用最優(yōu)化技術(shù)求解 特點(diǎn)：a) 迭代過程不會發(fā)散 b) 可以明確潮流問題是否有解“有解可求、無解可判” if 潮流問題 有解目標(biāo)函數(shù)很快趨近 零 if 潮流問題無解目標(biāo)函數(shù)穩(wěn)定在一非零 值,64,2.6.2 潮流問題的非線性規(guī)劃模型,(1) 模型描述 潮流問題即為求解非線性方程組（節(jié)點(diǎn)功率（電壓）偏移量方程）：,Let,Then,f(x)=0，潮流問題有解,f(x) 0，潮流問題無解,65,2.6.2 潮流問題的非線性規(guī)劃模

22、型,(2) 描述潮流問題的非線性規(guī)劃模型,給定初值 x(0)，并置迭代次數(shù) k=0； 確定搜索方向(尋優(yōu)方向) x(k) ； 沿搜索方向確定最優(yōu)步長*(k)：,(3) 最優(yōu)乘子法的基本求解步驟,(d) 得到新的近似解： (e) 收斂判斷,66,2.6.3 最優(yōu)乘子潮流算法,(1) 搜索方向,() 的表達(dá)：,(2) 求解最優(yōu)步長因子*(k),一維搜索問題：,67,2.6.3 最優(yōu)乘子潮流算法,* 的確定：,(2) 求解最優(yōu)步長因子*(k),2.6 病態(tài)條件下的潮流計算最優(yōu)乘子法,68,(1) 線性方程式的迭代收斂性,2.7.1 迭代算法收斂性的基本概念,2.7 潮流算法的收斂性比較,ax-b=0

23、 x=x-ax+b=(1-a)x+b x(k+1)=x(k)+b x(k+1)-x(k)= k(x(1)-x(0) 結(jié)論：if |1 then k時， x(k+1)-x(k) 0，迭代收斂 |越小，收斂越快，收斂性越好！,(2) 線性方程組的迭代收斂性,Ax=b F(x)=Ax-b=0 x(k+1)=x(k)-F(x(k)=x(k)-(Ax(k)-b)=(E-A)x(k)+b x(k+1)=x(k)+b ; =E-A 結(jié)論：if |=q1 then k時， x(k+1)-x(k) 0，迭代收斂 |越小，收斂越快，收斂性越好！,69,2.7.1 迭代算法收斂性的基本概念,2.7 潮流算法的收斂性

24、比較,(2) 線性方程組的迭代收斂性關(guān)于|的定義：,|ij| 愈小，收斂愈快，收斂性愈好； 對 Ax=b , A 愈接近 E，收斂性愈好。,70,2.7.1 迭代算法收斂性的基本概念,2.7 潮流算法的收斂性比較,(3) 非線性方程組的迭代收斂性 對 非線性方程組 F(x)=0，泰勒展開：,J(0) 愈接近 E，收斂性愈好 !,71,2.7.1 迭代算法收斂性的基本概念,2.7 潮流算法的收斂性比較,(4) 非線性方程組的迭代格式 參照 Ax=b 的迭代格式：,F(x)展開式：,F(x)=Axb=0,構(gòu)造非線性方程組的迭代格式：,72,2.7.2 潮流算法迭代收斂性的比較,2.7 潮流算法的收

25、斂性比較,(1) 導(dǎo)納矩陣迭代法的收斂性 YV=I F(V)=YV-I=0V=(E-D-1Y)V+D-1I,迭代格式： V(k+1)=V(k)-D-1YV(k)-I(V(k) i.e. V(k+1)=V(k)-D-1F(V(k) ),Let I=const Then Jv-Y=D-1F(V(k) vD-1Y,Y 對角優(yōu)勢越強(qiáng)， JV 越接近 E 算法收斂性越好,73,2.7.2 潮流算法迭代收斂性的比較,2.7 潮流算法的收斂性比較,(2) 阻抗矩陣迭代法的收斂性 在Y迭代法中，F(xiàn)(V)=YV-I=0J=F(V)V YJ-1 Z 構(gòu)造 V(k+1)=V(k) - J-1 F(V(k) =V(k) Z F(V(k)=V(k)-ZYV(k)-I(V(k) Then, 阻抗矩陣迭代法 的收斂性 優(yōu)于 導(dǎo)納矩陣迭代法！,比較Y迭代法： Jv-Y=D-1F(V(k)vD-1Y,阻抗矩陣迭代法的迭代格式,V(k+1)=Z I(V(k) =V(k) Z F(V(k),JV-Z=Z F(V(k)v ZY=E,74,2.7.2 潮流算法迭代收斂性的比較,(3) N-R法 的收斂性 迭代格式：V(k+1)=V(k) J(k)-1 F(V(k), N-R迭代法的收斂性 優(yōu)于 Z 迭代法！,JV-NR= J(k)-1 F(V(k)v = J(k)-1 J(k) =E,比較Z迭