%%本程序為求解節點導納矩陣

originalData36
t1=menu（‘選擇系統是否含有變壓器；是選1，否選2；‘‘1‘‘2‘）
?switch?t1
?????case?1
?????????disp（‘系統含有變壓器‘）
?????????linedata=oriLineData（）;
?????????transformdata=transformdata（）;
?????????I?=?linedata（:2）;??????%?起節點
?????????J?=?linedata（:3）;??????%?終節點
?????????l?=?linedata（:4）;??????%線路長度
?????????nbus?=?max（max（I）max（J））;????%?節點數
?????????nbranch?=?length（J）;??????????%?支路數
?????????Y=zeros（nbus*3）;
?????????nline=size（linedata1）;
?????????ntrans=size（transformdata1）;
????????if?nline>=1
?????????for?m=1:nbranch
????????????I?=?linedata（m2）;??????%?起節點
????????????J?=?linedata（m3）;??????%?終節點
????????????l?=?linedata（m4）;??????%線路長度為英尺feet
????????????L=l/5280;???????????????%線路長度轉換成mile
????????????Zl=oriZl*L*SB/UB^2/1000;
????????????Yl=inv（Zl）;
???????????Y（3*（I-1）+1:3*（I-1）+33*（J-1）+1:3*（J-1）+3）=Y（3*（I-1）+1:3*（I-1）+33*（J-1）+1:3*（J-1）+3）-Yl;
???????????Y（3*（J-1）+1:3*（J-1）+33*（I-1）+1:3*（I-1）+3）=Y（3*（I-1）+1:3*（I-1）+33*（J-1）+1:3*（J-1）+3）;
???????????Y（3*（I-1）+1:3*（I-1）+33*（I-1）+1:3*（I-1）+3）=Y（3*（I-1）+1:3*（I-1）+33*（I-1）+1:3*（I-1）+3）+Yl;
???????????Y（3*（J-1）+1:3*（J-1）+33*（J-1）+1:3*（J-1）+3）=Y（3*（J-1）+1:3*（J-1）+33*（J-1）+1:3*（J-1）+3）+Yl;??
?????????end
????????
????if?ntrans>=1
????????typeT=transformdata（:4）;
???????for?m=1:ntrans?
?????????I?=?transformdata（m2）;??????%?起節點
?????????J?=?transformdata（m3）;??????%?終節點
?????????U1=transformdata（m6）;
?????????U2=transformdata（m7）;
?????????k=（U1*UB2）/（U2*UB）;??
?????????SN=transformdata（m5）;
?????????R_100=transformdata（m8）;
?????????X_100=transformdata（m9）;
?????????ZT=U1^2/SN*1000;
?????????RT=R_100/100*ZT*SB/UB^2/1000;
?????????XT=X_100/100*ZT*SB/UB^2/1000;??
?????????yt=1.0/（RT+j*XT）;
??????????if?typeT（m）==1??%類型-1表示變壓器為YNyn接線方式
?????????????y11=[yt/k^2?0?0?;0?yt/k^2?0;0?0?yt/k^2];
?????????????y12=[yt/k?0?0;0?yt/k?0;0?0?yt/k];
?????????????y22=[yt?0?0;0?yt?0;0?0?yt];
??????????end
??????????if?typeT（m）==2?%類型-2表示變壓器為YNd11接線方式
??????????????y11=[yt/k^2?0?0?;0?yt/k^2?0;0?0?yt/k^2];
??????????????y12=[yt/k/sqrt（3）?0?0;0?yt/k/sqrt（3）?0;0?0?yt/k/sqrt（3）];
??????????????y22=[yt/3?0?0;0?yt/3?0;0?0?yt/3];
??????????end
??????????????Y（3*（I-1）+1:3*（I-1）+33*（J-1）+1:3*（J-1）+3）=Y（3*（I-1）+1:3*（I-1）+33*（J-1）+1:3*（J-1）+3）-y12;
??????????????Y（3*（J-1）+1:3*（J-1）+33*（I-1）+1:3*（I-1）+3）=Y（3*（I-1）+1:3*（I-1）+33*（J-1）+1:3*（J-1）+3）;
??????????????Y（3*（I-1）+1:3*（I-1）+33*（I-1）+1:3*（I-1）+3）=Y（3*（I-1）+1:3*（I-1）+33*（I-1）+1:3*（I-1）+3）+y11;
??????????????Y（3*（J-1）+1:3*（J-1）+33*（J-1）+1:3*（J-1）+3）=Y（3*（J-1）+1:3*（J-1）+33*（J-1）+1:3*（J-1）+3）+y22;?????????
???????end
????end
?end
?
?????case?2
?????????disp（‘系統沒有變壓器‘）
?????????linedata=oriLineData（）;
?????????I?=?linedata（:2）;??????%?起節點
?????????J?=?linedata（:3）;??????%?終節點
?????????l?=?linedata（:4）;??????%線路長度
?????????nbus?=?max（max（I）max（J））;?%?節點數
?????????nbranch?=?length（J）;??%?支路數
?????????Y=zeros（nbus*3）;
????for?m=1:nbranch
???????????I?=?linedata（

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????4310??2013-05-08?12:54??yinshigaosi_method\jiediandaona.asv

?????文件???????3792??2013-05-15?21:39??yinshigaosi_method\jiediandaona.m

?????文件???????2911??2013-05-08?14:41??yinshigaosi_method\originalData12.asv

?????文件???????2716??2013-05-08?14:44??yinshigaosi_method\originalData12.m

?????文件???????6802??2013-05-08?15:19??yinshigaosi_method\originalData36.m

?????文件???????1972??2013-05-08?14:31??yinshigaosi_method\originalData6.asv

?????文件???????1955??2013-05-08?14:36??yinshigaosi_method\originalData6.m

?????文件???????3809??2013-05-15?15:20??yinshigaosi_method\yinshigaosifa.asv

?????文件???????3809??2013-05-15?21:44??yinshigaosi_method\yinshigaosifa.m

?????目錄??????????0??2013-05-16?13:33??yinshigaosi_method

-----------?---------??----------?-----??----

????????????????32076????????????????????10