branch_dscb=load（‘branch_desc.txt‘）;%導入支路描述表
node_dscb=load（‘node_desc.txt‘）;%導入節點描述表
branch_num=size（branch_dscb1）;%支路數
node_num=size（node_dscb1）;%節點數
pq_num=0;
for?k=1:node_num????%統計pq節點數
????if?node_dscb（k6）==1
?????????pq_num=pq_num+1;
????end
end

Y=zeros（node_num）;%定義節點導納矩陣

for?k=1:branch_num
????head=branch_dscb（k1）;%首節點寫入
????tail=branch_dscb（k2）;%末節點寫入
????r=branch_dscb（k3）;%電阻寫入
????x=branch_dscb（k4）;%電抗寫入
????ban=branch_dscb（k5）;%對地電納寫入
????K=branch_dscb（k6）;%變比寫入
????
????%==節點導納寫入==
????Y（headhead）=Y（headhead）+1/（r+1j*x）+1j*ban/2;
????Y（headtail）=Y（headtail）-1/（K*（r+1j*x））;
????Y（tailhead）=Y（tailhead）-1/（K*（r+1j*x））;
????Y（tailtail）=Y（tailtail）+1/（K*K*（r+1j*x））+1j*ban/2;

end
G=real（Y）;%實部為電導
B=imag（Y）;%虛部為電納
precision=1;%精度初始值設

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

?????文件????????642??2019-04-15?16:12??chaoliujisuan\branch_desc.txt

?????文件????????396??2019-04-15?16:12??chaoliujisuan\jiedian.txt

?????文件???????4042??2019-04-15?16:11??chaoliujisuan\loadflowcal.m

?????目錄??????????0??2019-04-15?16:13??chaoliujisuan

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

?????????????????5080????????????????????4