39IEEE隐式梯形积分程序

是用来计算39节点的iEEE的隐式梯形积分法，用来算故障的。

%动态仿真作业（隐式梯形积分）
%说明: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