matlab潮流计算仿真,MATLAB潮流计算仿真

MATLAB潮流计算仿真

附录A MATLAB程序

%本程序的功能就是用牛顿——拉夫逊法进行潮流计算

% B1矩阵:1、支路首端号;2、末端号;3、支路阻抗;4、支路对地电纳

% 5、支路的变比;6、支路首端处于K侧为1,1侧为0

% B2矩阵:1、该节点发电机功率;2、该节点负荷功率;3、节点电压初始值

% 4、PV节点电压V的给定值;5、节点所接的无功补偿设备的容量

% 6、节点分类标号

clear;

n=10;%input('请输入节点数:n=');

nl=11;%input('请输入支路数:nl=');

isb=1;%input('请输入平衡母线节点号:isb=');

pr=0、00001;%input('请输入误差精度:pr=');

B1=[1 2 1、755e-2+4、155e-2i 0、26i 1 0;

1 4 3、159e-2+7、479e-2i 0、1215i 1 0;

1 6 3、159e-2+7、479e-2i 0、1215i 1 0;

2 3 3、68e-3+0、11135i 0 0、909 1;

4 5 3、68e-3+0、11135i 0 0、909 1;

4 6 2、808e-2+6、648e-2i 0、108i 1 0;

6 7 3、0865e-3+0、0833i 0 0、909 1;

6 8 3、159e-2+7、479e-2i 0、1215i 1 0;

6 10 2、457e-2+5、817e-2i 0、0945i 1 0;

8 9 3、0865e-3+0、08