33节点潮流计算

用于IEEE33节点潮流计算。潮流计算，电力学名词，指在给定电力系统网络拓扑、元件参数和发电、负荷参量条件下，计算有功功率、无功功率及电压在电力网中的分布。潮流计算是根据给定的电网结构、参数和发......

%% IEEE-33节点系统潮流计算
%输入的支路信息是有名值
%输入Bus矩阵，第一列节点号，第二列首节点/有功，第三列尾节点/无功
%输入Branch矩阵第一列存支路号，第2列存支路首节点，
%%%%%%%%%%第3列存支路尾节点，第4存支路电阻，第5列存支路电抗）;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%format long;
clc
clear all
tic
h=waitbar（0‘please wait‘）;

for i=1:1000
    %computation here%
    str=[‘运行中...‘num2str（i/1000*100）‘%‘];
    waitbar（i/1000hstr）
end
delete（h）;
Bus=[100;210060;39040;412080;56030;66020;7200100;8200100;96020;106020;114530;126035;136035;1412080;156010;166020;176020;189040;199040;209040;219040;229040;239040;24420200;25420200;266025;276025;286020;2912070;30200600;3115070;32210100;336040;];
Branch=[1120.09220.0407;2230.49300.2511;3340.36600.1864;4450.38110.1941;5560.81900.7070;6670.18720.6188;7780.71440.2351;8891.03000.7400;99101.04400.7400;1010110.19660.065;1111120.37440.1238;1212131.46801.1550;1313140.54160.7129;1414150.59100.5260;1515160.74630.5450;1616171.28901.7210;1717180.73200.5740;182190.16400.1565;1919201.50421.3554;2020210.40950.4784;2121220.70890.9373;223230.45120.3083;2323240.89800.7091;2424250.89600.7011;256260.20300.1034;2626270.28420.1447;2727281.05900.9337;2828290.80420.7006;2929300.50750.2585;3030310.97440.9630;3131320.31050.3619;3232330.34100.5302;];
[busnumrow]=size（Bus）;
[branchnumrow]=size（Branch）;
soubus=Branch（:2）;
mobus=Branch（:3）;
Vbus=ones（busnum1）;
Vbus（:1）=12.66;
Vbus1=Vbus;
Ploss=zeros（busnum1）;
Qloss=zeros（busnum1）;
e=1;
k=0;
Branch1=Branch;
n=1;
%%%%%%%%%%%%%%%支路重新排序，各个分支线同时进行计算
while ~isempty（Branch1）%%%%T