電力系統(tǒng)分析課程設(shè)計

1、電力系統(tǒng)分析課程設(shè)計學(xué) 院 電氣工程學(xué)院 班 級 學(xué) 號 姓 名 指導(dǎo)教師 時 間 前言 電力系統(tǒng)潮流計算是研究電力系統(tǒng)穩(wěn)態(tài)運行情況的一種計算，它根據(jù)給定的運行條件及系統(tǒng)接線情況確定整個電力系統(tǒng)各部分的運行狀態(tài)：各母線的電壓，各元件中流過的功率，系統(tǒng)的功率損耗等等。在電力系統(tǒng)規(guī)劃的設(shè)計和現(xiàn)有電力系統(tǒng)運行方式的研究中，都需要利用潮流計算來定量地分析比較供電方案或運行方式的合理性?？煽啃院徒?jīng)濟性。此外，電力系統(tǒng)潮流計算也是計算系統(tǒng)動態(tài)穩(wěn)定和靜態(tài)穩(wěn)定的基礎(chǔ)。所以潮流計算是研究電力系統(tǒng)的一種很重要和很基礎(chǔ)的計算。 隨著科學(xué)技術(shù)的發(fā)展，電力系統(tǒng)變得越來越復(fù)雜，電氣工程師掌握一種好的能對電力系統(tǒng)進行仿真的

2、軟件是學(xué)習(xí)和研究的需要。與眾多專門的電力系統(tǒng)仿真軟件相比，MATLAB軟件具有易學(xué)、功能強大和開放性好，是電力系統(tǒng)仿真研究的有力工具。目錄1 設(shè)計題目31.1系統(tǒng)圖的確定31.2各節(jié)點的初值及阻抗參數(shù)42 潮流計算52.1潮流計算概述與發(fā)展52.2復(fù)雜電力系統(tǒng)潮流計算52.3 MATLAB概述62.4牛頓-拉夫遜法原理72.5牛頓-拉夫遜法解決潮流計算問題82.6計算機潮流計算的步驟92.7計算機潮流計算流程圖113 手算潮流計算123.1確定節(jié)點類型123.2求節(jié)點導(dǎo)納矩陣Yb123.3計算各節(jié)點功率的修正方程的初始值（不平衡量）134 計算機算法潮流計算154.1計算機法潮流計算過程154

3、.2計算機法潮流計算結(jié)果234.3系統(tǒng)功率分布圖25個人心得29附錄：源程序30參考文獻401 設(shè)計題目1.1系統(tǒng)圖的確定選擇六節(jié)點、環(huán)網(wǎng)、兩電源和多引出的電力系統(tǒng)，簡化電力系統(tǒng)圖如圖1-1所示，等值導(dǎo)納圖如圖1-2所示。運用以直角坐標表示的牛頓-拉夫遜計算如圖1-1所示系統(tǒng)中的潮流分布。計算精度要求各節(jié)點電壓的誤差或修正量不大于。 圖1-1 電力系統(tǒng)圖圖1-2 電力系統(tǒng)等值導(dǎo)納圖1.2各節(jié)點的初值及阻抗參數(shù)該系統(tǒng)中，節(jié)點為平衡節(jié)點，保持=1.05+j0為定值，節(jié)點為PV節(jié)點，其他四個節(jié)點都是PQ節(jié)點。給定的注入電壓標幺值如表1-1、線路阻抗標幺值如表1-2、輸出功率標幺值如表1-3。表1-1

4、 各節(jié)點電壓標幺值參數(shù)UUUUUU1.051.001.001.001.001.05表1-2 線路、變壓器阻抗標幺值 線路L2L3L4L5T1T2Y/2阻抗0.06+j0.250.04+j0.250.08+j0.300.1+j0.35j0.03j0.015j0.25表1-3 節(jié)點輸出功率節(jié)點功率2+j11.8+j0.41.6+j0.83.7+j1.35注：各PQ節(jié)點的電壓取1是為了方便計算和最后驗證程序的正確性。2 潮流計算2.1潮流計算概述與發(fā)展 電力系統(tǒng)潮流計算也分為離線計算和在線計算兩種，前者主要用于系統(tǒng)規(guī)劃設(shè)計和安排系統(tǒng)的運行方式，后者則用于正在運行系統(tǒng)的經(jīng)常監(jiān)視及實時控制。 利用電子數(shù)

5、字計算機進行電力系統(tǒng)潮流計算從50年代中期就已經(jīng)開始。在這20年內(nèi)，潮流計算曾采用了各種不同的方法，這些方法的發(fā)展主要圍繞著對潮流計算的一些基本要求進行的。 牛頓-拉夫遜法作為一種實用的，有競爭力的電力系統(tǒng)潮流計算方法，是在應(yīng)用了稀疏矩陣技巧和高斯消去法求修正方程后。牛頓-拉夫遜法是求解非線性代數(shù)方程有效的迭代計算。2.2復(fù)雜電力系統(tǒng)潮流計算 電力系統(tǒng)潮流計算是對復(fù)雜電力系統(tǒng)正常和故障條件下穩(wěn)態(tài)運行狀態(tài)的計算。潮流計算的目標是求取電力系統(tǒng)在給定運行方式下的節(jié)點電壓和功率分布，用以檢查系統(tǒng)各元件是否過負荷、各點電壓是否滿足要求、功率的分布和分配是否合理以及功率損耗等。對現(xiàn)有電力系統(tǒng)的運行和擴建，

6、對新的電力系統(tǒng)進行規(guī)劃設(shè)計都是以潮流計算為基礎(chǔ)。 潮流計算結(jié)果的用途，例如用于電力系統(tǒng)穩(wěn)定研究、安全估計或最優(yōu)潮流等也對潮流計算的模型和方法有直接影響。 節(jié)點類型：（1）PV節(jié)點：柱入有功功率P為給定值，電壓也保持在給定數(shù)值。（2）PQ節(jié)點：諸如有功功率和無功功率是給定的。（3）平衡節(jié)點：用來平衡全電網(wǎng)的功率。選一容量足夠大的發(fā)電機擔(dān)任平衡全電網(wǎng)功率的職責(zé)。平衡節(jié)點的電壓大小與相位是給定的，通常以它的相角為參考量，即取其電壓相角為0。一個獨立的電力網(wǎng)中只設(shè)一個平衡點。 基本步驟：（1）形成節(jié)點導(dǎo)納矩陣；（2）將各節(jié)點電壓設(shè)初值U；（3）將節(jié)點初值代入相關(guān)求式，求出修正方程式的常數(shù)項向量；（4）

7、將節(jié)點電壓初值代入求式，求出雅可比矩陣元素；（5）求解修正方程，求修正向量；（6）求取節(jié)點電壓的新值；（7）檢查是否收斂，如不收斂，則以各節(jié)點電壓的新值作為初值自第3步重新開始進行狹義次迭代，否則轉(zhuǎn)入下一步；（8）計算支路功率分布，PV節(jié)點無功功率和平衡節(jié)點柱入功率。2.3 MATLAB概述 目前電子計算機已廣泛應(yīng)用于電力系統(tǒng)的分析計算，潮流計算是其基本應(yīng)用之一?，F(xiàn)有很多潮流計算方法。對潮流計算方法有五方面的要求：（1）計算速度快；（2）內(nèi)存需要少；（3）計算結(jié)果有良好的可靠性和可信性；（4）適應(yīng)性好，亦即能處理變壓器變比調(diào)整、系統(tǒng)元件的不同描述和與其它程序配合的能力強；（5）簡單。 MATL

8、AB是一種交互式、面向?qū)ο蟮某绦蛟O(shè)計語言，廣泛應(yīng)用于工業(yè)界與學(xué)術(shù)界，主要用于矩陣運算，同時在數(shù)值分析、自動控制模擬、數(shù)字信號處理、動態(tài)分析、繪圖等方面也具有強大的功能。 MATLAB程序設(shè)計語言結(jié)構(gòu)完整，且具有優(yōu)良的移植性，它的基本數(shù)據(jù)元素是不需要定義的數(shù)組。它可以高效率地解決工業(yè)計算問題，特別是關(guān)于矩陣和矢量的計算。MATLAB與C語言和FORTRAN語言相比更容易被掌握。通過M語言，可以用類似數(shù)學(xué)公式的方式來編寫算法，大大降低了程序所需的難度并節(jié)省了時間，從而可把主要的精力集中在算法的構(gòu)思而不是編程上。2.4牛頓-拉夫遜法原理假設(shè)有n個聯(lián)立的非線性代數(shù)方程： 假設(shè)以給出各變量的初值，令其分

9、別為個變量的修正量，使?jié)M足以上方程，所以： 將上式中的n個多元函數(shù)在初始值附近分別展開成泰勒級數(shù)，并略去含有，的二次及以上階次的各項，便得：方程可寫成： 以上方程是對于修正量，的線性方程組，稱為牛頓法的修正方程，可解出，。對初始近似解進行修正： （i=1,2,，n）反復(fù)迭代，在進行k+1次迭代時，從求解修正方程式：得到修正量,對各量進行修正 （i=1,2,，n）迭代過程一直進行到滿足收斂判據(jù)2.5牛頓-拉夫遜法解決潮流計算問題節(jié)點總數(shù)為n；PQ節(jié)點有m，；PV節(jié)點有n-m-1，平衡節(jié)點有1個，節(jié)點編號按照先PQ節(jié)點，再PV節(jié)點，最后平衡節(jié)點的順序進行編號，即:1，2，m為PQ節(jié)點；m+1，m+

10、2，n-1為PV節(jié)點；n為平衡節(jié)點。可形成結(jié)點導(dǎo)納矩陣。導(dǎo)納矩陣元素可表示為，本文中節(jié)點電壓以直角坐標形式表示，即。由此下列公式可求出Pi，Qi假設(shè)系統(tǒng)中的第1，2，m號節(jié)點為PQ節(jié)點，第i個節(jié)點的給定功率為和，對該節(jié)點可列方程： 假設(shè)系統(tǒng)中的第m+1，m+2，n-1號節(jié)點為PV節(jié)點，則對其中每一個節(jié)點可列方程: 第n號節(jié)點為平衡節(jié)點，其電壓為是給定的，故不參加迭代。 修正方程可寫成分塊矩陣的形式： 通過反復(fù)求解修正方程，解出各節(jié)點的未知量，再通過收斂判據(jù)判定是否已為真值。從而求得PQ節(jié)點的電壓V及相角的真值，PV節(jié)點的Q、真值，平衡節(jié)點的P、Q真值，以上即為牛頓-拉夫遜迭代法的潮流計算過程，

11、其優(yōu)點為計算精確，運行速度快。其中的各個環(huán)節(jié)都可通過MATLAB程序來實現(xiàn)。2.6計算機潮流計算的步驟（1）對電力網(wǎng)絡(luò)的所有參數(shù)設(shè)初值,包括電壓、相角、有功、無功等。（2）處理非標準變比支路,使其變成標準變比為1的變壓器支路。（3）形成節(jié)點導(dǎo)納矩陣Y。（4）計算有功功率的不平衡量Pi,從而求出。（5）根據(jù)節(jié)點的類型形成J。（6）解修正方程式,求各節(jié)點的電壓的變化量ei(i=1,2,3.n,is)（7）求各節(jié)點相角的新值ei=ei+ei (i=1,2,3. n,is)（8）計算無功功率的不平衡量Qi,從而求出 (i=1,2,3.n,is)（9）解修正方程式,求各節(jié)點的電壓大小的變化量 (i=1,

12、2,3.,n,is)。（10）求各節(jié)點的電壓大小的新值 (i=1,2,3.,n,is)。（11）運用個節(jié)點的電壓的新值自第四步開始下一次迭代。 計算平衡節(jié)點的功率和線路功率。其中平衡節(jié)點的功率的計算公式為線路上的功率為: 從而線路上的損耗的功率為:2.7計算機潮流計算流程圖選用牛頓-拉夫遜法進行計算機潮流計算。牛頓-拉夫遜法潮流計算流程圖如圖2-1 。圖2-1 牛頓-拉夫遜法潮流計算流程圖3 手算潮流計算3.1確定節(jié)點類型節(jié)點為平衡節(jié)點節(jié)點為PV節(jié)點，其他四個節(jié)點都是PQ節(jié)點。3.2求節(jié)點導(dǎo)納矩陣Yb 由圖1可知，該系統(tǒng)以串聯(lián)支路的導(dǎo)納標幺值和對地并聯(lián)導(dǎo)納標幺值得等值電路如圖2所示?？傻脠D2的

13、等值節(jié)點導(dǎo)納矩陣。對角線上的元素為：非對角線上的元素為：所以節(jié)點導(dǎo)納矩陣為：Y=3.3計算各節(jié)點功率的修正方程的初始值（不平衡量）設(shè)各節(jié)點電壓初始值為：根據(jù)上述N-R法的求解過程，將各節(jié)點電壓的初始值代入（1）式和（2）式，進行潮流計算，得節(jié)點功率和節(jié)點電壓：雅可比矩陣：J=修正各節(jié)點電壓：4 計算機算法潮流計算4.1計算機法潮流計算過程導(dǎo)納矩陣 Y=0 -30.2343i 0 +31.7460i 0 0 0 0 0 +31.7460i 14.8252 -42.6506i -14.2012 + 5.9172i 0 -0.6240 + 3.9002i 0 0 -14.2012+5.9172i 1

14、5.0311-8.5292i -0.8299+3.1120i 0 0 0 0 -0.8299 + 3.1120i 1.5846 - 5.5035i -0.7547 + 2.6415i 0 0 -0.6240 + 3.9002i 0 -0.7547 + 2.6415i 1.3787 -72.9583i 0 +63.4921i0 0 0 0 0 +63.4921i 0 -60.4686i功率方程第(1)次差值： Columns 1 through 12 0 0 4.2619 -2.1000 0.1000 -1.8000 -0.5500 -1.6000 8.4738 -3.7000 0 5.0000

15、形成的第(1)次Jacobi矩陣： Columns 1 through 11 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -37.3887 14.8252 5.9172 -14.2012 0 0 3.9002 -0.6240 0 0 0 -14.8252 -47.9126 14.2012 5.9172 0 0 0.6240 3.9002 0 0 0 5.9172 -14.2012 -8.0292 15.0311 3.1120 -0.8299 0 0 0 0 0 14.2012 5.9172 -15.0311 -9.0292 0.8299 3.

16、1120 0 0 0 0 0 0 0 3.1120 -0.8299 -5.2535 1.5846 2.6415 -0.7547 0 0 0 0 0 0.8299 3.1120 -1.5846 -5.7535 0.7547 2.6415 0 0 0 3.9002 -0.6240 0 0 2.6415 -0.7547 -63.1845 1.3787 63.4921 0 0 0.6240 3.9002 0 0 0.7547 2.6415 -1.3787 -82.7321 0 0 0 0 0 0 0 0 0 0 0 -2.4000 0 0 0 0 0 0 0 0 0 76.1905 0 Columns

17、 12 through 13 0 0 0 0 0 4.2619 0 -2.1000 0 0.1000 0 -1.8000 0 -0.5500 0 -1.6000 0 8.4738 63.4921 -3.7000 0 0 -63.4921 5.0000Jacobi矩陣第(1)次回代運算 Columns 1 through 12 0 0 -0.1490 0.1103 -0.0024 0.1657 0.1015 0.3342 -0.1435 0.0297 0 -0.0431各個節(jié)點電壓模 0 1.1543 1.0160 0.9586 1.1439 1.2008功率方程第(2)次差值： Columns

18、 1 through 12 0 0 -0.9196 -0.3298 -0.2847 -0.0127 -0.5652 0.0384 -1.5185 0.1960 -0.0019 -0.3928形成的第(2)次Jacobi矩陣： Columns 1 through 11 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -47.1567 23.2568 5.2333 -16.9700 0 0 4.4126 -1.1470 0 0 0 -20.2170 -47.5883 16.9700 5.2333 0 0 1.1470 4.4126 0 0 0 3.5

19、786 -15.2150 -5.6605 18.1967 2.9819 -1.3474 0 0 0 0 0 15.2150 3.5786 -14.7623 -6.4582 1.3474 2.9819 0 0 0 0 0 0 0 2.5188 -1.7857 -3.5899 4.7795 2.1211 -1.5609 0 0 0 0 0 1.7857 2.5188 -1.7466 -5.2407 1.5609 2.1211 0 0 0 4.4413 -0.8294 0 0 2.9982 -0.9415 -83.4895 7.1529 72.6032 0 0 0.8294 4.4413 0 0 0

20、.9415 2.9982 -0.3333 -83.2844 1.8854 0 0 0 0 0 0 0 0 0 0 -2.4000 0 0 0 0 0 0 0 0 -2.7375 76.1905 -1.8854 Columns 12 through 13 0 0 0 0 0 -0.9196 0 -0.3298 0 -0.2847 0 -0.0127 0 -0.5652 0 0.0384 -1.8854 -1.5185 72.6032 0.1960 -0.0862 -0.0019 -72.6032 -0.3928Jacobi矩陣第(2)次回代運算 Columns 1 through 12 0 0

21、0.0441 0.0155 0.0700 -0.0210 0.1790 -0.0275 0.0288 0.0422 -0.001 0.0487各個節(jié)點電壓模 0 1.1120 0.9435 0.7821 1.1170 1.2010功率方程第(3)次差值： Columns 1 through 12 0 0 -0.0417 0.0017 -0.0185 -0.0126 -0.1340 -0.0217 -0.0367 0.0646 -0.0024 -0.0916形成的第(3)次Jacobi矩陣： Columns 1 through 11 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0

22、 0 0 0 0 0 0 0 0 -44.1917 23.5230 4.7526 -16.4350 0 0 4.2309 -1.1798 0 0 0 -19.9621 -46.3315 16.4350 4.7526 0 0 1.1798 4.2309 0 0 0 3.4621 -14.0963 -5.0872 17.0579 2.7814 -1.2240 0 0 0 0 0 14.0963 3.4621 -13.4379 -6.4674 1.2240 2.7814 0 0 0 0 0 0 0 1.9845 -1.5515 -1.8991 4.3504 1.6691 -1.3531 0 0 0

23、0 0 1.5515 1.9845 -1.3054 -5.0484 1.3531 1.6691 0 0 0 4.3026 -0.9762 0 0 2.8902 -1.0313 -79.8824 10.0760 70.7751 0 0 0.9762 4.3026 0 0 1.0313 2.8902 -3.4952 -82.5737 4.5677 0 0 0 0 0 0 0 0 0 0 -2.4019 0 0 0 0 0 0 0 0 0.3535 76.2524 -4.5677 Columns 12 through 13 0 0 0 0 0 -0.0417 0 0.0017 0 -0.0185 0

24、 -0.0126 0 -0.1340 0 -0.0217 -4.5677 -0.0367 70.7751 0.0646 0.0111 -0.0024 -70.7751 -0.0916Jacobi矩陣第(3)次回代運算 Columns 1 through 12 0 0 0.0078 0.0043 0.0191 -0.0063 0.0768 -0.0055 0.0056 0.0156 0.0011 0.0181各個節(jié)點電壓模 0 1.1048 0.9237 0.7098 1.1126 1.2001功率方程第(4)次差值： Columns 1 through 12 0 0 -0.0004 0.000

25、2 -0.0002 -0.0012 -0.0276 -0.0058 0.0006 0.0016 -0.0003 -0.0053形成的第(4)次Jacobi矩陣： Columns 1 through 11 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -43.7449 23.5927 4.6455 -16.3504 0 0 4.1980 -1.1918 0 0 0 -20.0302 -45.9892 16.3504 4.6455 0 0 1.1918 4.1980 0 0 0 3.4389 -13.7881 -4.9899 16.7680 2.7

26、273 -1.1885 0 0 0 0 0 13.7881 3.4389 -13.0467 -6.4292 1.1885 2.7273 0 0 0 0 0 0 0 1.7502 -1.4706 -1.1217 4.2480 1.4705 -1.2806 0 0 0 0 0 1.4706 1.7502 -1.1039 -4.9984 1.2806 1.4705 0 0 0 4.2712 -1.0337 0 0 2.8637 -1.0684 -79.3737 11.1438 70.4222 0 0 1.0337 4.2712 0 0 1.0684 2.8637 -4.6945 -82.2282 5

27、.5608 0 0 0 0 0 0 0 0 0 0 -2.3998 0 0 0 0 0 0 0 0 1.5031 76.1844 -5.5608 Columns 12 through 13 0 0 0 0 0 -0.0004 0 0.0002 0 -0.0002 0 -0.0012 0 -0.0276 0 -0.0058 -5.5608 0.0006 70.4222 0.0016 0.0473 -0.0003 -70.4222 -0.0053Jacobi矩陣第(4)次回代運算 Columns 1 through 12 0 0 0.0021 0.0014 0.0056 -0.0020 0.025

28、9 -0.0018 0.0016 0.0048 0.0002 0.0053各個節(jié)點電壓模 0 1.1029 0.9178 0.6857 1.1114 1.2000功率方程第(5)次差值： Columns 1 through 12 0 0 -0.0000 0.0000 0.0000 -0.0001 -0.0032 -0.0007 -0.0000 0.0001 -0.0000 -0.0005形成的第(5)次Jacobi矩陣： Columns 1 through 11 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -43.6300 23.6210 4

29、.6137 -16.3286 0 0 4.1889 -1.1958 0 0 0 -20.0561 -45.8841 16.3286 4.6137 0 0 1.1958 4.1889 0 0 0 3.4337 -13.6963 -4.9688 16.6803 2.7114 -1.1776 0 0 0 0 0 13.6963 3.4337 -12.9312 -6.4137 1.1776 2.7114 0 0 0 0 0 0 0 1.6711 -1.4437 -0.8564 4.2162 1.4034 -1.2564 0 0 0 0 0 1.4437 1.6711 -1.0344 -4.9843 1

30、.2564 1.4034 0 0 0 4.2621 -1.0517 0 0 2.8559 -1.0801 -79.2374 11.4911 70.3226 0 0 1.0517 4.2621 0 0 1.0801 2.8559 -5.0505 -82.1223 5.8687 0 0 0 0 0 0 0 0 0 0 -2.3993 0 0 0 0 0 0 0 0 1.8419 76.1690 -5.8687Columns 12 through 13 0 0 0 0 0 -0.0000 0 0.0000 0 0.0000 0 -0.0001 0 -0.0032 0 -0.0007 -5.8687

31、-0.0000 70.3226 0.0001 0.0580 -0.0000 -70.3226 -0.0005Jacobi矩陣第(5)次回代運算 Columns 1 through 12 0 0 0.0003 0.0002 0.0008 -0.0003 0.0039 -0.0003 0.0002 0.0007 0.0000 0.0008各個節(jié)點電壓模 0 1.1026 0.9169 0.6820 1.1113 1.2000功率方程第(6)次差值： 1.0e-004 * Columns 1 through 12 0 0 -0.0026 0.0007 0.0032 -0.0182 -0.7419 -

32、0.1647 -0.0016 0.0184 -0.0063 -0.1036形成的第(6)次Jacobi矩陣： Columns 1 through 11 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -43.6128 23.6253 4.6089 -16.3254 0 0 4.1876 -1.1964 0 0 0 -20.0599 -45.8683 16.3254 4.6089 0 0 1.1964 4.1876 0 0 0 3.4330 -13.6825 -4.9658 16.6670 2.7090 -1.1760 0 0 0 0 0 13.6

33、825 3.4330 -12.9139 -6.4113 1.1760 2.7090 0 0 0 0 0 0 0 1.6591 -1.4396 -0.8161 4.2114 1.3932 -1.2528 0 0 0 0 0 1.4396 1.6591 -1.0238 -4.9822 1.2528 1.3932 0 0 0 4.2607 -1.0543 0 0 2.8548 -1.0818 -79.2170 11.5431 70.3078 0 0 1.0543 4.2607 0 0 1.0818 2.8548 -5.1036 -82.1067 5.9148 0 0 0 0 0 0 0 0 0 0

34、-2.3993 0 0 0 0 0 0 0 0 1.8924 76.1670 -5.9148 Columns 12 through 13 0 0 0 0 0 -0.0000 0 0.0000 0 0.0000 0 -0.0000 0 -0.0001 0 -0.0000 -5.9148 -0.0000 70.3078 0.0000 0.0596 -0.0000 -70.3078 -0.0000Jacobi矩陣第(6)次回代運算 1.0e-004 * Columns 1 through 12 0 0 0.0764 0.0496 0.2048 -0.0728 0.9503 -0.0641 0.0564 0.1753 0.0074 0.1923各