前推回代潮流计算

通过前推回代计算潮流，固定节点编号，形成节点关联矩阵，计算潮流。

clc;
clear;
data=xlsread（‘1.xls‘）;      %读取线路参数

[vSdelt SI]=qthd（data）;     %潮流计算


fprintf（‘初始潮流计算线路情况：\n‘）
v                                        %各节点电压输出

S=S*1000                                 %单位改为 kVA
deltS=1000*deltS                         %各线段上的损耗

sumdS=sum（deltS）                         %所有线路损耗之和

I=1000*abs（I）                            %单位改为 A

IN=[595595595595595595595595595595...
    360360215215215215215 180215215...
    215215215215215530530]‘;           %分段线路额定载流量


Ibai=100*I./IN                           %各线路电流的百分数


r=zeros（271）;
for k=1:27
    a=real（S（k+1））;
    b=imag（S（k+1））;
    r（k）=a/（sqrt（a*a+b*b））;
end
r                                      %功率因数

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件        2245  2018-05-10 17:40  qthd.m
     文件         831  2018-05-10 17:41  pfc.m