clear all
clc
ydata=zeros(1,41);
for i=1:41;
% | number | From | To | X |
% | of line | Bus | Bus | pu |
A = [ 1 1 2 0.0575
2 1 3 0.1852
3 2 4 0.1737
4 3 4 0.0379
5 2 5 0.1983
6 2 6 0.1763
7 4 6 0.0414
8 5 7 0.1160
9 6 7 0.0820
10 6 8 0.0420
11 6 9 0.2080
12 6 10 0.5560
13 9 11 0.2080
14 9 10 0.1100
15 4 12 0.2560
16 12 13 0.1400
17 12 14 0.2559
18 12 15 0.1304
19 12 16 0.1987
20 14 15 0.1997
21 16 17 0.1932
22 15 18 0.2185
23 18 19 0.1292
24 19 20 0.0680
25 10 20 0.2090
26 10 17 0.0845
27 10 21 0.0749
28 10 22 0.1499
29 21 22 0.0236
30 15 23 0.2020
31 22 24 0.1790
32 23 24 0.2700
33 24 25 0.3292
34 25 26 0.3800
35 25 27 0.2087
36 27 28 0.3960
37 27 29 0.4153
38 27 30 0.6027
39 29 30 0.4533
40 8 28 0.2000
41 6 28 0.0599 ];
k=max(ydata(:,1)+1);
k=A(:,1);
a=A(:,2);
b=A(:,3);
d=A(:,4);
ydata(k,1)=k;
ydata(k,2)=a;
ydata(k,3)=b;
ydata(k,4)=d;
%ydata
% i=i+1;
j=sqrt(-1); i = sqrt(-1);
nl = ydata(:,2);
nr = ydata(:,3);
X = ydata(:,4);
nbr=length(ydata(:,2));
nbus = max(max(nl), max(nr));
Z = j*X;
Y= ones(nbr,1)./Z;
% formulation of incedence matrix
A=zeros(nbr,nbus);
for j=1:nbus
for i=1:nbr
if (ydata(i,2)==j)
A(i,j)=1;
end
if (ydata(i,3)==j+1)
A(i,j+1)=-1;
end
end
end
A
A';
y=zeros(nbr,nbr);
for m=1:nbr
y(m,m)=Y(m);
end
break;
end
% end
Ybus=A'*y*A;
Ybus=-Ybus*sqrt(-1);
T=-sqrt(-1)*y;
S=1;
%S= Enter slack bus number
Yred=[];
Yred=Ybus;
Yred(:,S) = [];
Yred(S,:) = [];
Pbus=zeros(1,1);
for q=1:29
%w = input('1.Enter bus power one by one in sequence 2.finished ');
% | Pin |
w= [ .3586
-.024
-.076
-.938
0
-.228
.05
0
-.058
.1793
-.112
.1691
-.062
-.082
-.035
-.09
-.032
-.095
-.022
-.175
0
-.032
-.087
0
-.035
0
0
-.024
-.106 ];
%if w == 1
%e=input('bus power(pu) ');
Pbus=w;
q=q+1;
%end
%if w == 2
% break;
%end
end
Pbus=w;
Ared=A;
Ared(:,S)=[];
% Ared(S,:)=[];
Ared
theta=[];
theta=-inv(Yred)*Pbus;
theta
P=[];
K=[];
K=T*Ared;
P=K*theta;
disp('Power flow ');
P
dcloadflow.zip_DC FLOW_flow
版权申诉
193 浏览量
2022-07-14
09:14:53
上传
评论
收藏 1KB ZIP 举报
周楷雯
- 粉丝: 80
- 资源: 1万+
最新资源
- 2-08 http协议和nginx入门
- 基于Java的ADmobile ADSuyi-Flutter插件广告聚合SDK设计源码
- mkl-2021.1.1-py2.py3-none-manylinux1-x86-64.whl
- Drivingnotmessage.apk
- mkl-2021.4.0-py2.py3-none-manylinux1-x86-64.whl
- wiliwili-yoga.zip
- 项目计划书.doc
- mkl-2021.2.0-py2.py3-none-manylinux1-x86-64.whl
- 情侣恋爱相册记录-游玩生活
- 毕业设计项目源码-使用Matlab基于蚁群算法的三维路径规划算法实现.zip
资源上传下载、课程学习等过程中有任何疑问或建议,欢迎提出宝贵意见哦~我们会及时处理!
点击此处反馈