1、Matlab形成節(jié)點(diǎn)導(dǎo)納矩陣學(xué)號(hào)： 姓名：侯成濱引言：電力網(wǎng)的運(yùn)行狀態(tài)可用節(jié)點(diǎn)方程或回路方程來描述。節(jié)點(diǎn)導(dǎo)納矩陣是以系統(tǒng)元件的等值導(dǎo)納為基礎(chǔ)所建立的、描述電力網(wǎng)絡(luò)各節(jié)點(diǎn)電壓和注入電流之間關(guān)系的線性方程。導(dǎo)納矩陣計(jì)算是電力系統(tǒng)分析最基本的計(jì)算。除它自身的重要作用之外，還是網(wǎng)損計(jì)算、靜態(tài)安全分析、暫態(tài)穩(wěn)定計(jì)算、小干擾靜態(tài)穩(wěn)定計(jì)算、短路計(jì)算、靜態(tài)和動(dòng)態(tài)等值計(jì)算的基礎(chǔ)。本次任務(wù)是用計(jì)MATLAB語言編寫程序求出潮流計(jì)算中要用到的導(dǎo)納矩陣。為了確定結(jié)果是否正確，與一個(gè)手工計(jì)算比較運(yùn)算結(jié)果， 驗(yàn)證程序是否正確。一、分析網(wǎng)絡(luò)等效電路此電力系統(tǒng)是一個(gè)6節(jié)點(diǎn)，7支路的電力網(wǎng)絡(luò)?？梢园严到y(tǒng)等值網(wǎng)絡(luò)畫出來，如圖1-

2、1。圖1-1 某電力系統(tǒng)的等值網(wǎng)絡(luò)在計(jì)算電力系統(tǒng)網(wǎng)絡(luò)的潮流分布時(shí)，我們需要把變壓器轉(zhuǎn)化成變壓器的型等值電路來進(jìn)行計(jì)算器等效導(dǎo)納，根據(jù)等效電路的等效原則，可以把上圖等效成如圖1-2導(dǎo)納等值網(wǎng)絡(luò)圖。對(duì)導(dǎo)納等值網(wǎng)絡(luò)圖簡(jiǎn)化電路圖，可以得到圖1-3簡(jiǎn)化導(dǎo)納等值電路圖，方便潮流計(jì)算中導(dǎo)納矩陣的計(jì)算。圖1-2電力系統(tǒng)網(wǎng)絡(luò)的導(dǎo)納等值電路圖1-3電力系統(tǒng)簡(jiǎn)化等值電路圖二、MATLAB程序形成導(dǎo)納矩陣導(dǎo)納矩陣的計(jì)算總結(jié)如下：1）導(dǎo)納矩陣的階數(shù)等于電力系統(tǒng)網(wǎng)絡(luò)的節(jié)點(diǎn)數(shù)；2）導(dǎo)納矩陣各行非對(duì)角元素中非零元素的個(gè)數(shù)等于對(duì)應(yīng)節(jié)點(diǎn)所連的不接地支路數(shù)；3）導(dǎo)納矩陣的對(duì)角元素，即各節(jié)點(diǎn)的自導(dǎo)納等于相應(yīng)節(jié)點(diǎn)所連之路的導(dǎo)納之和：Y

3、ij=jiyij其中，yij為節(jié)點(diǎn)i與節(jié)點(diǎn)支路阻抗Zij的倒數(shù)，符號(hào)ji表示j屬于i或與i相連的j，即內(nèi)只包括與節(jié)點(diǎn)i直接相連的節(jié)點(diǎn)j。當(dāng)節(jié)點(diǎn)i有接地支路時(shí)，還應(yīng)包括j=0的情況。4）導(dǎo)納矩陣非對(duì)角元素等于節(jié)點(diǎn)i與節(jié)點(diǎn)j之間的導(dǎo)納的負(fù)數(shù)。2.1 MATLAB程序及其運(yùn)行節(jié)點(diǎn)導(dǎo)納程序如下：N=input(請(qǐng)輸入節(jié)點(diǎn)數(shù): N=);L=input(請(qǐng)輸入支路數(shù): L=);B=input(請(qǐng)輸入支路信息: B=);X=input(請(qǐng)輸入由節(jié)點(diǎn)號(hào)及其對(duì)地阻抗形成的矩陣：X=);Y=zeros(N);for n=1:N;if X(n,2)=0; p=X(n,1); Y(p,p)=1./X(n,2);endendfor n=1:L if B(n,6)=0 p=B(n,1);q=B(n,2); else p=B(n,2);q=B(n,1); end Y(p,q)=Y(p,q)-1./(B(n,3)*B(n,5); Y(q,p)=Y(p,q); Y(q,q)=Y(q,q)+1./(B(n,3)*B(n,5)2)+B(n,4)./2; Y(p,p)=Y