基于前推回代法的连续潮流计算IEEE33节点算例matlab程序

基于前推回代法的连续潮流计算IEEE33节点算例matlab程序
1 理论介绍
在图11-1(a)所示的网络中,供电点A通过馈电干线向负荷节点b、c和d供电,各负荷节点功率已知。如果节点d的电压也给定,就可以从节点d开始,利用同一点的电压和功率计算第三段线路的电压降落和功率损耗,得到节点c的电压,并算出第二段线路末端的功率,然后依次计算第二段线路和第一段线路的电压降落和功率损耗,一次性地求得解答。但是实际的情况并不这么简单,多数的情况是已知电源点电压和负荷节点的功率,要求确定各负荷点电压和网络中的功率分布。在这种情况下,可以采取近似的方法通过迭代计算求得满足一定精度的解答。
基于前推回代法的连续潮流计算IEEE33节点算例matlab程序_第1张图片
在进行电压和功率分布计算以前,先要对网络的等值电路(见图11-1(b))作些简化处理。具体的做法是,将输电线等值电路中的电纳支路都分别用额定电压U、下的充电功率代替,这样,对每段线路的首端和末端的节点都分别加上该段线路充电功率的一半。
在这里插入图片描述
基于前推回代法的连续潮流计算IEEE33节点算例matlab程序_第2张图片
习惯上称S、S。和S:为电力网的运算负荷。这样,我们就把原网络简化为由三个集中的阻抗元件相串联,而在四个节点(包括供电点)接有集中负荷的等值网络(见图11-1©)。针对图11-1©所示的等值网络将按以下两个步骤进行电压和功率分布的计算。第一步,从离电源点最远的节点d开始,利用线路额定电压,逆着功率传送的方向依次算出各段线路阻抗中的功率损耗和功率分布。对于第三段线路
基于前推回代法的连续潮流计算IEEE33节点算例matlab程序_第3张图片
基于前推回代法的连续潮流计算IEEE33节点算例matlab程序_第4张图片
2 算例
1)IEEE33节点算例结构
基于前推回代法的连续潮流计算IEEE33节点算例matlab程序_第5张图片
2)负荷时序值
基于前推回代法的连续潮流计算IEEE33节点算例matlab程序_第6张图片
3 程序运行结果
基于前推回代法的连续潮流计算IEEE33节点算例matlab程序_第7张图片
4 matlab程序

clc
clear
close all

%% 日前负荷时序值
load=[0.450000000000000,0.337500000000000,0.300000000000000,0.287500000000000,0.275000000000000,0.337500000000000,0.362500000000000,0.375000000000000,0.412500000000000,0.475000000000000,0.500000000000000,0.525000000000000,0.587500000000000,0.637500000000000,0.562500000000000,0.537500000000000,0.575000000000000,0.650000000000000,0.837500000000000,0.962500000000000,1,0.900000000000000,0.812500000000000,0.675000000000000];

figure
plot(load)
xlabel('时间/h')
ylabel('负荷/pu')
title('负荷时序值')


%% 前推回代法算潮流
v=zeros(33,24);
for t=1:24
%% 系统负荷参数 节点  有功  无功
Bus=[1,0,0   ;
2,100, 60;
3,90,40;
4,120,80;
5,60,30;
6,60,20;
7, 200, 100 ;
8,200,  100 ;
9, 60, 20;
10,60, 20;
11,45,30;
12,60, 35;
13,60, 35;
14,120,80;
15,60, 10;
16,60, 20;
17,60, 20;
18,90, 40;
19,90, 40;
20,90, 40;
21,90, 40;
22,90, 40; 
23,90, 40; 
24,420, 200;
25,420, 200 ;
26,60, 25;
27,60 ,25;
28,60, 20;
29,120, 70;
30,200,600 ;
31,150, 70;
32,210, 100 ;
33,60,40;];
%% 负荷更新
  Bus(:,2)=load(t)*Bus(:,2);
  Bus(:,3)=load(t)*Bus(:,3);
  %% 计算潮流
[voltage]=back_function(Bus);
v(:,t)=voltage;
end
%% 标幺化
v=v/12.66;
%% 输出
figure
ylabel('时间/h')
xlabel('节点号')
zlabel('电压幅值/pu')
title('电压分布三维图')
。。。。。。。。略

你可能感兴趣的:(电气工程,电力系统,潮流计算,matlab,算法)