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