資源描述：
《節(jié)點導(dǎo)納矩陣的建立》由會員上傳分享，免費在線閱讀，更多相關(guān)內(nèi)容在行業(yè)資料-天天文庫。
1、如上圖所示的簡單電力系統(tǒng)中，網(wǎng)絡(luò)各元件參數(shù)的標(biāo)幺值如下：z12=0.10+j0.40y120=y210=j0.01528z13=j0.3,k=1.1z14=0.12+j0.50y140=y410=j0.01920z24=0.08+j0.40y240=y420=j0.01413系統(tǒng)中節(jié)點1、2為PQ節(jié)點，節(jié)點3為PV節(jié)點，節(jié)點4為平衡節(jié)點。節(jié)點導(dǎo)納矩陣的運行程序如下：clcCleardisp('網(wǎng)絡(luò)各元件參數(shù)用標(biāo)幺值表示?。?！');N0=input('請輸入節(jié)點數(shù):N0=');n1=input('請輸入支路數(shù):n1=');l=input('請
2、輸入PQ節(jié)點的個數(shù)=');form=1:lc(m)=input(['請輸入第',num2str(m),'個PQ節(jié)點的節(jié)點號為：']);endt=input('請輸入PV節(jié)點的個數(shù)=');form=1:t???c(m)=input(['請輸入第',num2str(m),'個PV節(jié)點的節(jié)點號為：']);endb=input('請輸入平衡節(jié)點號:b=');%%由支路參數(shù)形成矩陣B1disp('各支路連接情況：')?i=1;form=1:n1?????symsYN???p=input(['第',num2str(m),'條支路的起始節(jié)點']);???
3、q=input(['第',num2str(m),'條支路的終止節(jié)點']);???mn=input(['第',num2str(m),'條支路是否有變壓器（請輸入‘Y’或‘N’）']);???y=0;k=1;???ifmn=='Y';???????k=input('請輸入變壓器變比(標(biāo)幺值):');???????z=input(['請輸入第',num2str(m),'條支路的線路阻抗']);???else????????z=input(['請輸入第',num2str(m),'條支路的線路阻抗：']);????????y=input(['請輸入第
4、',num2str(m),'條支路線路的對地阻抗：']);???end?????B1(i,1)=p;B1(i,2)=q;B1(i,3)=z;B1(i,4)=y;B1(i,5)=1/k;???i=i+1;enddisp('由支路參數(shù)形成的矩陣B1')B1%求節(jié)點導(dǎo)納矩陣Y=zeros(N0);e=zeros(1,N0);f=zeros(1,N0);fori=1:n1????????????????p=B1(i,1);q=B1(i,2);???Y(p,q)=Y(p,q)-1./(B1(i,3)*B1(i,5));???Y(q,p)=Y(p,q)
5、;????????????????????????Y(q,q)=Y(q,q)+1./(B1(i,3)*B1(i,5)^2)+B1(i,4)./2;????Y(p,p)=Y(p,p)+1./B1(i,3)+B1(i,4)./2;???Enddisp('導(dǎo)納矩陣Y=');disp(Y)程序運行結(jié)果：
