电力系统潮流计算的MATLAB辅助程序设计潮流计算程序Word文件下载.docx
《电力系统潮流计算的MATLAB辅助程序设计潮流计算程序Word文件下载.docx》由会员分享，可在线阅读，更多相关《电力系统潮流计算的MATLAB辅助程序设计潮流计算程序Word文件下载.docx（26页珍藏版）》请在冰豆网上搜索。
节点信息输入格式为单行输入，输入的数据形成一个矩阵，叫做矩阵。
第一列为节点号；
第二列为节点类型；
第三列和第四列分别为节点电压幅值（标幺值）和相角（单位为度）；
第五列和第六列分别为负荷的有功功率和无功功率；
第七列到十列分别为发电机的有功功率、无功功率、最小无功出力和最大无功出力；
最后一列为并联电容器注入无功功率。
第二列的编码用0、1、2来区分节点、平衡节点和节点：
0表示节点，输入正的有功功率（）和无功功率（），并且要设定节点电压初始估计值，一般幅值和相角分别设为1和0，若已经给定初始值，则用其给定值来代替1和0。
1表示平衡节点，且已知该节点的电压幅值和相角。
2表示节点，要设定该节点的节点电压幅值和发电机的有功功率（），并设定发电机的无功最小出力和最大出力（）。
2.线路数据文件线路数据用节点对的方法来确定，数据包含在称为的矩阵中。
第一列和第二列为节点号码，第三列到第五列为线路电阻、电抗及该线路电纳值的一半，以标幺值表示。
最后一列为变压器分接头设定值，对线路来说，需要输入1。
线路输入为无输入顺序，对变压器来说，左侧的节点号设为分接头端。
3是线路数据输入变量，包括四项，前两项是节点编号，后两项是线路电阻和电抗，均以标幺值表示，函数返回节点导纳矩阵。
三、潮流计算的程序清单
1.程序清单
%
0;
0;
=0;
=（（:
1））;
[][];
[];
[];
1
（k,1）;
（n）（k,2）;
（n）（k,3）;
（n）（k,4）;
（n）（k,5）;
（n）（k,6）;
（n）（k,7）;
（n）=（k,8）;
（n）（k,9）;
（n）（k,10）;
（n）（k,11）;
（n）<
=0（n）=1.0;
V（n）=1+j*0;
（n）=180*（n）;
V（n）=（n）*（（（n））+j*（（n）））;
P（n）=（（n）（n））;
Q（n）=（（n）（n）+（n））;
S（n）=P（n）+j*Q（n）;
（n）=0;
=0;
=1;
=（,1）*（,1）;
（'
'
）1
=1.3;
=0.001;
100;
=100;
（,1）;
（（（k）（m））&
（（k）（m）））;
（m）=2;
10;
>
=&
<
=
1;
n=1;
=0*0;
L=1;
（（L）n&
（L）1）,（L）;
=+（）*V（k）;
（L）1）,（L）;
=（V（n））*（（）*V（n）+）;
=（）;
（n）=P（n）-（）;
（n）=Q（n）-（）;
（n）1
S（n）;
P（n）=（）;
Q（n）=（）;
（n）=0;
（n）=0;
（n）=V（n）;
（n）2
（n）0
=Q（n）*+（n）-（n）;
（（n））<
=.005&
=10
=0.045
（n）,
（n）=（n）+0.005;
（n）=（n）+.005;
（n）=（n）-0.005;
（n）（n）+.005;
（n）=（（S（n））（V（n））-）/（）;
V（n）=V（n）+*（（n）（n））;
=（（n））;
=（（n）^2-^2）;
（n）=+j*;
V（n）=V（n）+*（（n）（n））;
（（（（）））,（（（））））;
&
\:
'
）
）,（'
.\n\n'
\n'
,
）;
n=1
（n）=（V（n））;
（n）=（V（n））*180;
S（n）（n）*Q（n）;
（n）=P（n）*+（n）;
（n）=Q（n）*+（n）-（n）;
1;
（k）（n）;
（n）=（（n）-j*（n）*（n））/（*（n）^2）;
（:
3）'
;
（:
4）'
L
2程序清单
（-1）;
i=（-1）;
=（:
1）;
2）;
R=（:
3）;
X=（:
4）;
=j*（:
5）;
a=（:
6）;
（（:
=（（）,（））;
Z=R+j*X;
（,1）;
%支路导纳
a（n）<
=0a（n）=1;
（）;
%将初始化为0
%非对角元素的数值
（（k）（k））（（k）（k））（k）（k）;
（（k）（k））（（k）（k））;
%对角元素的数值
（k）
（）=（）（k）/（a（k）^2）+（k）;
（）=（）（k）（k）;
3.程序清单
%a
%.
（）
（'
=\n'
）
.=\n\n'
=['
..'
];
%5g'
n）,（'
%7.3f'
（n））,
%8.3f'
（n））,（'
%9.3f'
%9.3f'
%8.3f\n'
（n））
）,（'
）,
%9.3f\n\n'
4程序清单
\n'
\n\n'
&
0
%6g'
P（n）*）
%9.3f'
Q（n）*）,（'
%9.3f\n'
（S（n）*））
（L）k=（L）;
=（V（n）-a（L）*V（k））*y（L）（L）^2+（L）（L）^2*V（n）;
=（V（k）-V（n）（L））*y（L）+（L）*V（k）;
=V（n）*（）*;
=V（k）*（）*;
=+;
=（V（n）-V（k）（L））*y（L）+（L）*V（n）;
=（V（k）-a（L）*V（n））*y（L）（L）^2+（L）（L）^2*V（k）;
（L）|（L）
%12g'
k）,
（））,（'
（））
（））,
（L）&
a（L）1
a（L））
（'
=2;
5程序清单
%
（k）1,=1;
（k）2=1;
（k）=;
t=（）;
2*2*;
%雅可比矩阵
AJ
A（）=0;
%初始化雅可比矩阵
（n）;
（n）（n）;
J11=0;
J22=0;
J33=0;
J44=0;
（）1
（）n|（）n
（）n,l=（）;
J1111+（n）*（l）*（）*（t（）-（n）+（l））;
J3333+（n）*（l）*（）*（t（）-（n）+（l））;
（n）1
J2222+（l）*（）*（t（）-（n）+（l））;
J4444+（l）*（）*（t（）-（n）+（l））;
（n）1&
（l）1
=（l）（l）;
=l（l）;
%J1的非对角元素
A（,）（n）*（l）*（）*（t（）-（n）+（l））;
（l）0%J2的非对角元素
A（,）（n）*（）*（t（）-（n）+（l））
（n）0%J3的非对角元素
A（,）（n）*（l）*（）*（t（）-（n）（l））;
（n）0&
（l）0%J4的非对角元素
=（n）^2*（）*（t（））33;
=（n）^2*（）*（t（））11;
（n）1P（n）;
Q（n）=;
%P
（n）2Q（n）;
=7
2
（n）=（n）+0.01;
（n）=（n）-0.01
A（）=J11;
%J1对角元素
（）=P（n）;
A（）=2*（n）*（）*（t（））22;
%J2对角元素
A（）=J33;
%J3对角元素
A（）2*（n）*（）*（t（））44;
%J4对角元素
（）=Q（n）;
\'
（n）;
（n）（n）;
（n）=（n）（）;
（n）（n）（）;
（（））;
V=.*（）*.*（）;
180*;
S（n）=P（n）*Q（n）;
6程序清单
（n）1,=1;
（n）=;
（）0|（）2
B1（）（（））;
（）0
B2（）（（））;
B1（B1）;
B2=（B2）;
%开始迭代
=%检验不平衡功率
（）n,l=（）;
=20%1&
6
=10%
（n）,%.（n）
%0.05
（n）,%
（n）=（n）-0.005%.
（）=（P（n））（n）;
（）=（Q（n））（n）;
1\'
2\'
00;
（（（））（（）））;
A;
;
J11J33B1B1B2B2
四、30节点电力系统计算实例
潮流计算时，必须将前面的六个程序保存在目录下格式为的文件，然后在的命令窗口输入如下命令：
=1.8;
%30节点电力系统
%母线母线电压相角负载发电机注入功率
%编号节点幅值角度有功无功有功无功无功最小值无功最大值无功
[111.060.00.00.00.00.0000
221.0430.021.7012.740.00.0-40500
301.00.02.41.20.00.0000
401.060.07.61.60.00.0000
521.010.094.219.00.00.0-40400
601.00.00.00.00.00.0000
701.00.022.810.90.00.0000
821.010.030.030.00.00.0-30400
901.00.00.00.00.00.0000
1001.00.05.82.00.00.0-62419
1121.0820.00.00.00.00.0000
1201.0011.27.500000
1321.071000.000-6240
140106.21.600000
150108.22.500000
160103.51.800000
170109.05.800000
180103.20.900000
190109.53.400000
200102.20.700000
2101017.511.200000
2201000.000000
230103.21.600000
240108.76.700004.3
2501000.000000
260103.52.300000
2701000.000000
2801000.000000
290102.40.900000
3001010.61.900000];
%线路数据
%RX1/2B1
[120.01920.05750.026401
130.04520.18520.020401
240.05700.17370.018401
340.01320.03790.004201
250.04720.19830.020901
260.05810.17630.018701
460.01190.04140.004501
570.04600.11600.010201
670.02670.08200.008501
680.01200.04200.004501
690.00.20800.00.978
6100
升级会员 