function?Y?=?Admittance_matrix（nbAB）
%形成導納矩陣（節(jié)點數(shù)，支路數(shù)，支路參數(shù)矩陣A，電容對地導納矩陣B）
%函數(shù)說明
%{
????1.支路信息不可重復輸入，即已經(jīng)輸入i-j支路的信息，就不可輸入j-i支路的信息
????2.含變壓器的支路第一個節(jié)點編號默認為變壓器一次側(cè)（即變壓器等值電路中的阻抗歸算側(cè)）
????3.變比為1表示線路；變比不為1表示變壓器
???*4.考慮到課本例題4-1，節(jié)點可能接有電容，輸入矩陣B作為節(jié)點電容對地導納參數(shù)。如果沒有可以置為：
???????B=[1???0.00;
??????????2???0.00;???????????????
??????????3???0.00;
??????????...?...?];
%}

????Y=zeros（n）;%生成導納矩陣
????y1=zeros（n）;%串聯(lián)支路導納矩陣
????y2=zeros（n）;%節(jié)點處并聯(lián)支路導納矩陣

????%考慮到節(jié)點可能接有電容（課本例題4-1），計算節(jié)點電容對地導納參數(shù)
????for?p=1:n
????????if?B（p2）~=0;
????????????t=B（p1）;
????????????y2（tt）=B（p2）*i;
????????end
????end
????%計算串聯(lián)支路導納
????for?p=1:n
????????if?A（p6）~=1????%變比不為1時
????????????y1（A（p1）A（p2））=（1/（A（p4）*i））/A（p6）;
????????????y1（A（p2）A（p1））=y1（A（p