%本程序的功能是用牛頓——拉夫遜法進行潮流計算
%?B1矩陣：1、支路首端號；2、末端號；3、支路阻抗；4、支路對地電納
%?????????5、支路的變比；6、支路首端處于K側為1，1側為0
%?B2矩陣：1、該節點發電機功率；2、該節點負荷功率；3、節點電壓初始值
%?????????4、PV節點電壓V的給定值；5、節點所接的無功補償設備的容量
%?????????6、節點分類標號：1為平衡節點（應為1號節點）；2為PQ節點；
%????????????3為PV節點；
clear;
n=7;%input（‘請輸入節點數:n=‘）;
nl=7;%input（‘請輸入支路數:nl=‘）;
isb=1;%input（‘請輸入平衡母線節點號:isb=‘）;
pr=0.00001;%input（‘請輸入誤差精度:pr=‘）;
B1=[1?2?0.069+0.134i?0.027i?1?1;1?3?0.0281+0.068i?0.0135i?1?1;1?4?0.0380+0.115i?0.2316i?1?1;?1?5?0.052+0.103i?0.200i?1?1;1?6?0.098+0.2372i?0.02316i?1?1;6?7?0.00046+0.0011i?0.0135i?1?0;1?7?0.1123+0.271i?0.0136i?1?1;];%input（‘請輸入由各支路參數形成的矩陣:B1=‘）;
B2=[0?0?1?1?0?1;?0?0.25+0.1875i?1?0?0?2;0?0.32+0.24i?1?0?0?2;0?0.4+0.3i?1?0?0?2;0?0.24+0.18i?1?0?0?2;0?0.25+0.1875i?1?0?0?2;0?0.29+0.2175i?1?0?0?2;];%input（‘請輸入各節點參數形成的矩陣:B2=‘）

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

?????文件???????7855??2015-06-19?19:08??DG.m

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

?????????????????7855????????????????????1