運(yùn)用MATLAB軟件進(jìn)行潮流計(jì)算論文

摘要本文運(yùn)用MATLAB軟件進(jìn)行潮流計(jì)算，對(duì)給定題目進(jìn)行分析計(jì)算，再應(yīng)用DDRTS軟件，構(gòu)建系統(tǒng)圖進(jìn)行仿真，最終得到合理的系統(tǒng)潮流。潮流計(jì)算是電力系統(tǒng)最基本最常用的計(jì)算。根據(jù)系接線及元件參數(shù)，通過(guò)潮流計(jì)算可以確定各母線的電壓幅值和相角，各元件流過(guò)的功率，整個(gè)系統(tǒng)的功率損耗。潮流計(jì)算是實(shí)現(xiàn)電力系統(tǒng)安全經(jīng)濟(jì)發(fā)供電的必要手段和重要工作環(huán)節(jié)。因此，潮流計(jì)算在電力系統(tǒng)的規(guī)劃計(jì)算，生產(chǎn)運(yùn)行，調(diào)度管理及科學(xué)計(jì)算中都有著廣泛的應(yīng)用。拉夫遜Newton-Raphson法以及MATLAB軟件進(jìn)行計(jì)算對(duì)給定系統(tǒng)圖進(jìn)行了四種不同負(fù)荷下的潮流計(jì)算，經(jīng)過(guò)調(diào)節(jié)均得到符合電壓限制及功率限制的潮流分布。經(jīng)過(guò)四到五次迭代后均能收斂。根據(jù)運(yùn)算結(jié)果，分析各支路損耗最后，應(yīng)用DDRTS軟件，構(gòu)建系統(tǒng)圖，對(duì)給定負(fù)荷重新進(jìn)行分析，潮流計(jì)算后的結(jié)果也能滿足相應(yīng)的參數(shù)要求。關(guān)鍵詞：牛頓-拉夫遜法MATLABDDRTS潮流計(jì)算目錄 2 …5 …6 6 7 8 10 …10 …………10 ……12 15 19 29一、題目原始資料：變電所1變電所2變電所3變電所4一次側(cè)電壓220kV母線3一次側(cè)電壓220kV和20MW；發(fā)電廠二總裝機(jī)容量為（200MW）。（二）變電所的負(fù)荷分別為：（四）變電所1和變電所3分別配有兩臺(tái)容量為75MVA的變壓器，短路損耗414KW，短路電壓（%）=16.7；變的變壓器，短路損耗為245KW，短路電壓（%）=10.5；二、課程設(shè)計(jì)容及要求(一)課設(shè)容：1.對(duì)給定的網(wǎng)絡(luò)查找潮流計(jì)算所需的各元件等值參數(shù)，畫出等值電路圖。況下的潮流計(jì)算，并對(duì)計(jì)算結(jié)果進(jìn)行分析。3跟.隨變電所負(fù)荷按一定比例發(fā)生變化，進(jìn)行潮流計(jì)算分析。3）1和4號(hào)變電所的負(fù)荷同時(shí)以2的%比例下降，而2和3號(hào)變電所的負(fù)荷同時(shí)以2的%比例上升；壓35KV要求調(diào)整圍在35-36之間）結(jié)果的比較。7.最終形成課程設(shè)計(jì)成品說(shuō)明書。(二)課設(shè)要求：3.對(duì)不同的負(fù)荷變化，分析潮流分布，寫出分析說(shuō)明。表表示調(diào)節(jié)控制的參數(shù)變化。5.打印利用DDRTS進(jìn)行潮流分析繪制的系統(tǒng)圖，以及潮流分布圖。三、題目分析根據(jù)系統(tǒng)圖可知此系統(tǒng)為兩端供電網(wǎng)路，將母線1，2設(shè)為節(jié)點(diǎn)1，2，將變電所1、2、3、4的高低壓側(cè)分別設(shè)為節(jié)點(diǎn)3、4、5、6、7、8、9、10。并且，將節(jié)點(diǎn)1設(shè)為平衡節(jié)點(diǎn)，將節(jié)點(diǎn)2設(shè)為PV節(jié)點(diǎn)，其余節(jié)點(diǎn)設(shè)為PQ節(jié)點(diǎn)。根據(jù)題目原始資料，計(jì)算發(fā)電廠、變壓器及線路的參數(shù)。（1）運(yùn)用下列公式計(jì)算變壓器參數(shù)：采用變壓器參數(shù)為折算至高壓側(cè)的數(shù)值，其變比K<1，其中，NN變壓器有5個(gè)抽頭，電壓調(diào)節(jié)圍為U±2*2.5%，U對(duì)應(yīng)的分接頭開始時(shí)NN設(shè)變壓器高壓側(cè)接主接頭。Z=R+jX=(r+jx)L（3）變電所負(fù)荷分別為：變電所1S=50+j30.98722變電所2LL變電所3L=40+j24.78977變電所4L=70+j43.38210將參數(shù)整理如下表：首端號(hào)11123357948.5+j20.14.25+j10.058.5+j20.1電納有名值5678.79利用牛頓拉夫遜法進(jìn)行求解，用MATLAB軟件編程，可以求解系統(tǒng)潮流分布根據(jù)題目的不同要求對(duì)參數(shù)進(jìn)行調(diào)整，通過(guò)調(diào)節(jié)變壓器變比和發(fā)電廠的電壓，求解出合理的潮流分布，最后用DDRTS進(jìn)行潮流分析，將兩者進(jìn)行比較。四、題目求解求PQ節(jié)點(diǎn)的ΔP(k),ΔQ(k)，求PV節(jié)點(diǎn)的ΔP(k)，ΔU(k)iiii是否ijijijijijij求解修正方程，得Δei(k)，Δfi(k)否是B1矩陣：2、末端號(hào)；3、支路阻抗；5、支路的變比；6、支路首端處于K側(cè)為1，1側(cè)為0；7、線路/變壓器標(biāo)識(shí)（0/1）變壓器參數(shù)當(dāng)支路首端處于K側(cè)標(biāo)識(shí)為1時(shí)歸B2矩陣：2、該節(jié)點(diǎn)負(fù)荷功率；79900000000000020002002000200由于各節(jié)點(diǎn)電壓不應(yīng)超過(guò)各自的允許圍，由題目要求知，變壓器低壓側(cè)電壓有名值變化圍如下表：電壓下限（KV）電壓上限(KV)46表明線路有功功率損耗越小，越滿足實(shí)際要求。由于程序較復(fù)雜，故將其放入附錄中，具體程序和結(jié)果見附錄一由附錄一的結(jié)果可以觀察到各個(gè)節(jié)點(diǎn)電壓有名值分別為:電壓162738495電壓由此觀察到節(jié)點(diǎn)4、6、8、10的電壓都較正常圍偏高，因此調(diào)節(jié)變壓器分接頭和發(fā)電廠電壓，此系統(tǒng)進(jìn)行多次調(diào)整，最終得到合理結(jié)果。由于調(diào)節(jié)的次數(shù)很多，過(guò)程較為復(fù)雜，只選取其中有代表性的兩組數(shù)據(jù)如下所述。調(diào)節(jié)方法：電廠2電壓電廠2電壓5第一次230第二次230調(diào)節(jié)結(jié)果：12512538.5/225.554979662第一次32807343第二次15765516總損耗0.48834-27.7819i0.263394-13.9481.76971-21.9077i6.529-13.411i0.12987+3.9287i0.10514+2.8384i0.081433+2.4635i2.58404-16.6331i0.21589+5.828112.1668-78.6239i第一次0.490064-27.7753i0.264126-13.9466i1.77418-21.8951i6.52693-13.413i0.13705+4.146i0.1108+2.9912i0.085847+2.597i2.58541-16.6258i0.2279+6.1524i12.2023-77.7692i第二次0.491865-27.7685i0.264885-13.9444i1.77882-21.8821i6.52482-13.415i0.14448+4.3707i0.11664+3.1488i0.090404+2.7349i2.58684-16.6181i0.24034+6.4883i12.2391-76.8854i由電壓和功率損耗的有名值可以觀察到，最后一次的電壓和有功損耗與題目要求較接近，但仍超出要求的圍。第二次調(diào)變比的具體的潮流分布如下：各條支路的末端功率SjS(3,1)=-50.1445-35.358i各條支路的末端功率SjS(3,1)=-50.1445-35.358iS(1,3)=50.6363+7.58949iS(1,5)=50.3815+20.1921iS(1,7)=-78.77876+64.83188i S(2,9)=200-39.35203i S(3,4)=50.1445+35.358i S(5,6)=50.1166+34.1365i S(7,8)=40.0904+27.5247i S(7,9)=-120.648+59.18933iS(9,10)=70.2403+49.8704iS(5,1)=-50.1166-34.1365iS(7,1)=80.55759-86.71398iS(9,2)=-193.4752+25.93705iS(4,3)=-50-30.9872iS(6,5)=-50-30.9877i S(8,7)=-40-24.7898iS(9,7)=123.2348-75.80746iS(10,9)=-70-43.3821i各節(jié)點(diǎn)的功率S為(節(jié)點(diǎn)號(hào)從小到大排列)：得到的圖像如下：.由于調(diào)節(jié)變壓器分接頭達(dá)到極限仍未達(dá)到題目要求的圍所，以接下來(lái)調(diào)節(jié)發(fā)節(jié)方法和結(jié)果見下表：調(diào)節(jié)方法：調(diào)節(jié)方法電廠1電壓廠2電壓第一次第二次第三次調(diào)節(jié)結(jié)果：第一次第二次第三次231814454305898866777785531916556382(1,3)支路(1,5)支路(1,7)支路(2,9)支路(3,4)支路(5,6)支路(7,8)支路(7,9)支路總損耗0.491865-27.7685i0.264885-13.94441.77882-21.88216.52482-13.4150.14448+4.3707i0.11664+3.1488i0.090404+2.7349i2.58684-16.6181i0.24034+6.4883i12.2391-76.8854i第一次0.497649-27.4970.267686-13.8091.8-21.5991i6.57847-13.0289i0.14593+4.4148i0.11775+3.1788i0.091303+2.7621i2.60918-16.355i0.24283+6.5555i12.3508-75.3792i第二次0.509516-26.9583i0.273416-13.5422i1.84334-21.0337i6.68799-12.2547i0.14891+4.5049i0.12003+3.2403i0.093144+2.8178i2.65474-15.8293i0.24794+6.6934i12.579-72.3618i第三次0.521798-26.4226i0.279324-13.2764i1.88803-20.4691i6.80057-11.4776i0.15199+4.5981i0.12238+3.3037i0.095044+2.8753i2.7015-15.3042i0.25322+6.836i12.8139-69.3369i由電壓和功率損耗的有名值可以觀察到，最后一次的電壓符合題目要求，可各條支路的首端功率SiS(1,3)=50.6738+9.16273iS(1,5)=50.4017+21.015i-78.26164+67.14246iS(2,9)=200-36.51027iS(3,4)=50.152+35.5853iS(5,6)=50.1224+34.2914iS(7,8)=40.095+27.665i-120.2447+59.94649iS(9,10)=70.2532+50.2181i各條支路的末端功率SjS(3,1)=-50.152-35.5853iS(5,1)=-50.1224-34.2914iS(7,1)=80.14967-87.61153iS(9,2)=-193.1994+25.03265iS(4,3)=-50-30.9872iS(6,5)=-50-30.9877iS(8,7)=-40-24.7898iS(9,7)=122.9462-75.25073iS(10,9)=-70-43.3821i各節(jié)點(diǎn)的功率S為(節(jié)點(diǎn)號(hào)從小到大排列)：.得到的圖像如下：與DDRTS比較：潮流計(jì)算方式w4個(gè)變電所的負(fù)荷同時(shí)以2%比例增大；為了使電壓符合題目要求，因此將K4變比由11/231調(diào)為11/225.5各節(jié)點(diǎn)的電壓模值大小V為(節(jié)點(diǎn)號(hào)從小到大排列)：各節(jié)點(diǎn)的電壓相角sida為(節(jié)點(diǎn)號(hào)從小到大排列)：各節(jié)點(diǎn)的功率S為(節(jié)點(diǎn)號(hào)從小到大排列)：1.0e+002*各節(jié)點(diǎn)電壓各節(jié)點(diǎn)相角各節(jié)點(diǎn)功率225.0002.2500225.009.9672220.92-0.9616節(jié)點(diǎn)435.18-4.612.2160+0.3894i2.2089-0.0371i222.74-0.45062.2274-0.0175i節(jié)點(diǎn)6-3.09221.342.2891節(jié)點(diǎn)835.60-0.582.2116+0.0884i220.145.10162.1927+0.1958i節(jié)點(diǎn)7.與DDRTS比較：.潮流計(jì)算方式節(jié)點(diǎn)4節(jié)點(diǎn)6節(jié)點(diǎn)8節(jié)點(diǎn)各225．0225.0221.135.29222.8221.535.69220.3節(jié)點(diǎn)電壓各0-0.93-4.42-0.43-2.962.371-0.38節(jié)點(diǎn)相角769各0.18662.000-0.49-0.49-0.39-0.68節(jié)+點(diǎn)0.9540功率i.潮流計(jì)算方式⑶變電所1、4負(fù)荷同時(shí)以2的%比例下降，變電所2、3負(fù)荷同時(shí)以2的%比例上升：各節(jié)點(diǎn)電壓各節(jié)點(diǎn)相角225.000225.001221.15-0.9302節(jié)點(diǎn)435.29-4.42222.74-4.4206節(jié)點(diǎn)6-3.09221.422.3451節(jié)點(diǎn)835.62-0.53220.315.2062節(jié)點(diǎn)5.2.0000-0.3723i0.0000+0.0000i-0.490.0000+2.0000-0.3723i0.0000+0.0000i-0.490.0000+0.0000i-0.40-0.0000.0000i-0.68點(diǎn)0.9755率潮流計(jì)算方式⑴斷開①③支路后，參數(shù)變化及其調(diào)節(jié)方法：為了使電壓符合題目要求，因此將K1變比由38.5/231調(diào)為38.5/225.5支路矩陣B1第一行發(fā)生變化：各節(jié)點(diǎn)的電壓模值大小V為(節(jié)點(diǎn)號(hào)從小到大排列)：各節(jié)點(diǎn)的電壓相角sida為(節(jié)點(diǎn)號(hào)從小到大排列)：各節(jié)點(diǎn)的功率S為(節(jié)點(diǎn)號(hào)從小到大排列)：1.0e+002*節(jié)點(diǎn)4節(jié)點(diǎn)6節(jié)點(diǎn)8節(jié)點(diǎn)各225.00225.00215.4935.18222.80221.4235.65220.24節(jié)點(diǎn)電壓各0-1.797-5.37-0.442-3.032.3309-0.485.1697.337點(diǎn)相角3372.0000-0.3651i0.0000+0.0000i-0.50-0.0002.0000-0.3651i0.0000+0.0000i-0.50-0.0000.0000i-0.500.0000+0.0000i-0.40-0.0000.0000i-0.700.433點(diǎn)1.1318率與DDRTS比較：.潮流計(jì)算方式⑵斷開①⑤支路后，參數(shù)變化及其調(diào)節(jié)方法：為了使電壓符合題目要求，因此將K2變比由11/231調(diào)為11/225.5支路矩陣B第二行發(fā)生變化：節(jié)點(diǎn)各225.00225.00221.0435.24220.2210.41221.4235.65220.24節(jié)點(diǎn)電壓各節(jié)點(diǎn)相角010.044-0.945-4.51-0.861-3.382.3309-0.485.16973各0.23142.0000-0.000-0.50-0.000-0.500.0000-0.40-0.000-0.70.點(diǎn)1.0506率-0.3651i0.00000.3090.00000.3090.0000潮流計(jì)算方式⑶斷開①⑦支路后，參數(shù)變化及其調(diào)節(jié)方法：為了使電壓符合題目要求，因此將K4變比由11/231調(diào)為11/225.5支路矩陣B第三行發(fā)生變化：225.00218.93218.87225.00221.04222.80節(jié)點(diǎn)435.24節(jié)點(diǎn)835.21節(jié)點(diǎn)6節(jié)點(diǎn)各.點(diǎn)電壓各0-0.945-4.51-0.442-3.034.57147.33973.709節(jié)點(diǎn)相角各0.241172.000090.0000-0.507-0.000-0.500.00004-0.40-0.0003-0.70節(jié)+---點(diǎn)0.98450.21600.00000.00000.00000.0000功率iiiiii潮流計(jì)算方式.（4）斷開⑦⑨支路后，參數(shù)變化及其調(diào)節(jié)方法：為了使電壓符合題目要求，因此將K4變比由11/231調(diào)為11/225.5支路矩陣B第八行發(fā)生變化：各節(jié)點(diǎn)電壓各節(jié)點(diǎn)相角各節(jié)點(diǎn)功率225.0000.2531+i225.0022.0000-0.2432i221.04-0.9459-0.0000.0000i節(jié)點(diǎn)435.24-4.51-0.50222.80-0.440.0000+0.0000i節(jié)點(diǎn)6-3.03-0.50220.972.3240節(jié)點(diǎn)835.57-0.50-0.0000.0000i-0.40219.127.9387-0.0000.0000i節(jié)點(diǎn)4.3171-0.700.433潮流計(jì)算方式（5）斷開②⑨支路后，參數(shù)變化及其調(diào)節(jié)方法：為了使電壓符合題目要求，因此將K4變比由11/231調(diào)為11/225.5支路矩陣B第四行發(fā)生變化：各節(jié)點(diǎn)電壓各節(jié)點(diǎn)相角各節(jié)點(diǎn)功率225.0000.3029+i225.0052.0000-0.2610i221.04-0.9459-0.0000.0000i節(jié)點(diǎn)435.24-4.51-0.50222.80-0.4427-0.0000.0000i節(jié)點(diǎn)6-3.03-0.50219.412.3909節(jié)點(diǎn)835.30-0.48-0.0000.0000i-0.40216.515.33710.0000+0.0000i節(jié)點(diǎn)8-0.700.433.與DDRTS比較：潮流計(jì)算方式進(jìn)行潮流計(jì)算，得到潮流圖如下：.將DDRTS仿真得到的潮流圖中各個(gè)點(diǎn)電壓的有名值與潮流S節(jié)點(diǎn)1225.0節(jié)點(diǎn)2225.0節(jié)點(diǎn)3221.0400節(jié)點(diǎn)4節(jié)點(diǎn)5222.8節(jié)點(diǎn)6節(jié)點(diǎn)7221.4節(jié)點(diǎn)8節(jié)點(diǎn)9220.2453節(jié)點(diǎn)由上述比較可觀察到:應(yīng)用仿真和MATLAB計(jì)算的到的結(jié)果大致相同，可以認(rèn)為此潮流計(jì)算是正確的。的基礎(chǔ)知識(shí)，也知道電力系統(tǒng)在日后的學(xué)習(xí)工作中占據(jù)很重要的地位。特別是其中的三大計(jì)算：潮流計(jì)算、斷路計(jì)算、功率計(jì)算，更為重要。所以我們花了兩周的時(shí)間進(jìn)行潮流計(jì)算的課程設(shè)計(jì)，課程設(shè)計(jì)的過(guò)程我也是受益匪淺。課程設(shè)計(jì)的過(guò)程比基礎(chǔ)知識(shí)的學(xué)習(xí)更為復(fù)雜，在短短兩周的時(shí)間要應(yīng)用好MATLAB和DDRTS技術(shù)并進(jìn)行潮流計(jì)算對(duì)我來(lái)說(shuō)也是一種挑戰(zhàn)。課設(shè)的過(guò)程在老力系統(tǒng)安全經(jīng)濟(jì)發(fā)展的必要手段和重要環(huán)節(jié)。是我們學(xué)電會(huì)使用了MATLAB和DDRTS的最基本的技術(shù)。在潮流計(jì)算的過(guò)程中也是坎坷不斷，做了三次潮流計(jì)算的更改才形成最終這篇論文，但也是頗值都采用的標(biāo)幺值，后來(lái)覺得不合理又重新做了一次；第壓再調(diào)的變比，但是考慮到電廠電壓是不可輕易調(diào)的，要考慮實(shí)后采用先調(diào)變比再調(diào)電廠電壓的方法，符合了要求。到更多關(guān)于電力系統(tǒng)的知識(shí)，更深入的了解潮流計(jì)算以及其他計(jì)算分析。而我需感在課設(shè)過(guò)程中幫助我的老師們和同學(xué)，讓我在課設(shè)中能夠有很多收獲!附錄一：%本程序的功能是用牛頓——拉夫遜法進(jìn)行%B1矩陣：1、支路首端號(hào)；2、末端號(hào)；3、支路阻抗；%5、支路的變比；6、支路首端處于K側(cè)為1，1側(cè)為0；%B2矩陣：1、該節(jié)點(diǎn)發(fā)電機(jī)功率；2、該節(jié)點(diǎn)負(fù)荷功率；pr=1e-5;%input('請(qǐng)輸入誤差精度：pr=');n=10;%input('請(qǐng)輸入節(jié)點(diǎn)數(shù)：n='nl=9;%input('請(qǐng)輸入支路數(shù)：nl='%---------------------------------------------070+43.38210i1002]%i%n=3;%input('請(qǐng)輸入節(jié)點(diǎn)數(shù)：n='%nl=3;%input('請(qǐng)輸入支路數(shù)：nl=');%isb=1;%input('請(qǐng)輸入平衡母線節(jié)點(diǎn)號(hào)：isb=');%pr=1e-5;%input('請(qǐng)輸入誤差精度：pr=');%010+6i113.803]%inY=zeros(n);e=zeros(1,n);f=zeros(1,n);V=zeros(1,n);sida=zeros(1,n);S1=zero%%%-----------求導(dǎo)納矩陣------------------------ifB1(i,7)==1%-----Y(p,q)=Y(p,q)-1./(B1(i,3)*B1(%左節(jié)點(diǎn)(首端)處于1側(cè)%左節(jié)點(diǎn)(首端)處于K側(cè)Y(q,q)=Y(q,q)+1./(B1(i,3)*B1(i,5)^2);%對(duì)角元K側(cè)Y(p,p)=Y(p,p)+1./B1(i,3)+B1(i,4);%對(duì)角元1側(cè)+勵(lì)磁導(dǎo)納else%------------否則為線路支路-----------------Y(p,q)=Y(p,q)-1./B1(i,3);%非Y(q,p)=Y(p,q);%非Y(q,q)=Y(q,q)+1./B1(i,3)+B1(i,4)./2.0;%對(duì)角元j側(cè)+線路電納的一半Y(p,p)=Y(p,p)+1./B1(i,3)+B1(i,4)./2.0;%對(duì)角元i側(cè)+線路電納的一半disp('導(dǎo)納矩陣Y=');disp(%-----------給定各節(jié)點(diǎn)初始電壓及給定各節(jié)點(diǎn)注入功率--------------------------%==================用牛頓-拉夫遜法迭代求解非線性代數(shù)方程（功率方程）=======================P=real(S);Q=imag(S);%分解出各節(jié)點(diǎn)注入的有功和無(wú)whileIT2~=0%----------------求取各個(gè)節(jié)點(diǎn)fori=1:n%n個(gè)節(jié)點(diǎn)2n行(每節(jié)點(diǎn)兩個(gè)方程P和Q或U)p=2*i-1;m=p+1;C(i)=0;DC(i)=C(i)+G(i,j1)*e(j1)-B(i,j1)*f(j1);%Σ(Gij*ej-Bij*fj)D(i)=D(i)+G(i,j1)*f(j1)+B(i,j1)*e(j1);%Σ(Gij*fj+BP1=C(i)*e(i)+f(i)*D(i);%節(jié)點(diǎn)功率P計(jì)%非平衡節(jié)點(diǎn)(PQ或PV節(jié)點(diǎn))J(p,N1)=V(i)^2-V2;%Pend%(ifi~=isb)非平衡節(jié)點(diǎn)(PQ或PV節(jié)end%(fori=1:n)n個(gè)節(jié)點(diǎn)2n行(每節(jié)點(diǎn)兩個(gè)方程P和Q或U)fork=3:N0除%去平衡節(jié)點(diǎn)1、2號(hào)以外的所有節(jié)點(diǎn)%--------------------以上為求取各個(gè)節(jié)點(diǎn)的功率及功率偏差及PV節(jié)點(diǎn)的電壓偏差%=================求取Jacobi矩陣形成修正方程===================fori=2:n%n個(gè)節(jié)點(diǎn)2n行(每節(jié)點(diǎn)兩個(gè)方程P和Q或U)ifB2(i,5)~=3%下面是針對(duì)PQ節(jié)點(diǎn)來(lái)求取Jacobi矩陣的元C(i)=C(i)+G(i,j1)*e(j1)-B(i,j1)*f(j1);%Σ(Gij*ej-Bij*fj)D(i)=D(i)+G(i,j1)*f(j1)+B(i,j1)*e(j1);%Σ(Gij*fj+Bforj1=2:n%第i行共n列(2n個(gè)Jacobi矩陣元素dP/de及dP/df或dQ/de及dQ/df)X1=-G(i,j1)*e(i)-B(i,j1)*f(i);%X1=dP/de=-dQ/dfX2=B(i,j1)*e(i)-G(i,j1)*f(i);%X2=dP/df=dQ/dX3=X2;%X2=dp/dfX3X4=-X1;%X1=dP/deX4J(p,q)=X4;J(m,q)=X2;%X1=-C(i)-G(i,i)*e(i)-B(i,i)*f(i);%X2=-D(i)+B(i,i)*e(i)-G(i,i)*f(i);%X3=D(i)+B(i,i)*e(i)-G(i,i)*f(i);%X4=-C(i)+G(i,i)*e(i)+B(i,i)*f(i);%p=2*i-1;q=2*j1-1;J(p,q)=X3;%擴(kuò)展列△QJ(p,N)=DQ;m=p+1;J(m,q)=X1;q=q+1;J(p,q)=X4;%擴(kuò)展列△Pelse%ifB2(i,5)~=3否則（即為PV節(jié)點(diǎn)）%===============下面是針對(duì)PV節(jié)點(diǎn)來(lái)求取Jacobi矩陣的元素===========X1=-G(i,j1)*e(i)-B(i,j1)*f(i);%X2=B(i,j1)*e(i)-G(i,j1)*f(i);%m=p+1;X1=-C(i)-G(i,i)*e(i)-B(i,i)*f(i);%X2=-D(i)+B(i,i)*e(i)-G(i,i)*f(i);%m=p+1;end%(ifB2(i,5)~=end%(fori=1:n)n個(gè)節(jié)點(diǎn)2n行(每節(jié)點(diǎn)兩個(gè)方程P和Q或U)JZ0=['形成的第(',num2str(%===============================以上為形成完整的Jacobi矩陣%====下面用高斯消去法對(duì)由Jacobi矩陣形成的修正方程進(jìn)行求解(按列消去、回代)==========fork1=k+1:N1%從k+1列的Jacobi元素到擴(kuò)展列的△P、△Q或△UJ(k,k1)=J(k,k1)./J(k,k);%用K行K列對(duì)角元素去除K行K列后的非對(duì)角元素進(jìn)行規(guī)格化J(k,k)=1;%對(duì)角元規(guī)格化K行K列對(duì)角元素賦1%==================按列消去運(yùn)算==================================fork2=k+1:N0%從k+1行到2*n最后一行fork3=k+1:N1%從k2+1列到擴(kuò)展列消去k+1行后各行下三角元素J(k2,k3)=J(k2,k3)-J(k2,kend%用當(dāng)前行K3列元素減去當(dāng)前行K列元素乘以第k行K3列元素%JZ=['Jacobi矩陣前代第(',num2str(a),')次'];JZ1=['Jacobi矩陣回代第(',num2str(a),')次%====================按列回代運(yùn)算=======================================%disp(JZ1);disp(J);disp(%----------------------------------修改節(jié)點(diǎn)電壓-------------------------------e(L)=e(L)-J(k,N1);%修改節(jié)%==============================結(jié)束一次迭代==============================%**********************下面為迭代計(jì)算結(jié)束后的有關(guān)輸出過(guò)程*****************disp('迭代次數(shù)：');disp('沒(méi)有達(dá)到精度要求的個(gè)數(shù)：');V(k)=sqrt(e(k)^2+f(k)^2);%計(jì)算各節(jié)點(diǎn)電壓sida(k)=atan(f(k)./e(k))*180./pi;E(k)=e(k)+f(k)*j;%將各節(jié)點(diǎn)電壓用復(fù)%===============計(jì)算各輸出量===========================disp(E);%顯示各節(jié)點(diǎn)的實(shí)際電壓值E用復(fù)數(shù)表示disp('--------------------------------------------------disp('各節(jié)點(diǎn)的電壓模值大小V為(節(jié)點(diǎn)號(hào)從小到大排列)：');disp(V);%顯示各節(jié)點(diǎn)的電壓大小V的模值disp('--------------------------------------------------disp(sida);%顯示各節(jié)點(diǎn)的C(p)=C(p)+conj(Y(p,q))*conj(E(q));%計(jì)算各節(jié)點(diǎn)注入電流的S(p)=E(p)*C(p);%計(jì)算各節(jié)點(diǎn)的功率S=電壓X注入電流的共軛值disp('各節(jié)點(diǎn)的功率S為(節(jié)點(diǎn)號(hào)從小到大排列)：');disp(S);%顯示各節(jié)點(diǎn)的disp('--------------------------------------------------disp