%動態仿真作業（隱式梯形積分）
%說明:1.潮流計算采用直角坐標形式保留二階項的快速潮流算法;
%?????2.暫態穩定分析采用改進歐拉法和直接解法;
%?????3.發電機模型為三階模型;
%?????4.勵磁是1階模型;

function?PowerFlow_TAnaly392
clear;
clc;
format?long;
i=sqrt（-1）;
%%==============================================讀取導納矩陣數據===================================================%
Y_S=xlsread（‘DAT39.xlsx‘‘Y_DAT‘‘A1:E1000‘）;????
[Y_S_n?Y_S_m]=size（Y_S）;????????????%確認共多少行
Ynum=39;?????????????????????????????%節點數
L_DAT=xlsread（‘DAT39.xlsx‘‘L_DAT‘‘A1:E40‘）;??
PL=L_DAT（2:Ynum+14）/100;
QL=L_DAT（2:Ynum+15）/100;
VL=L_DAT（2:Ynum+12）;
Y_Real=zeros（YnumYnum）;???????????????%?儲存導納矩陣的實部
Y_Imag=zeros（YnumYnum）;???????????????%?儲存導納矩陣的實部?
for?j=1:Y_S_n
????Yn=Y_S（j1）;
????Ym=Y_S（j2）;
????Y_Real（YnYm）=Y_S（j3）;
????Y_Imag（YnYm）=Y_S（j4）;
end
L_Real=PL./（

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

?????文件??????16509??2014-07-19?10:14??excel\DAT39.xlsx

?????文件??????35841??2014-07-23?09:25??excel\PowerFlow_TAnaly39_1step_excel.m

?????文件??????36478??2014-07-21?21:31??excel\PowerFlow_TAnaly39_3step_excel.m

?????目錄??????????0??2014-08-22?09:06??excel

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

????????????????88828????????????????????4