>> %本程序的功能是用牛顿——拉夫逊法进行潮流计算
% B1矩阵：1、支路首端号；2、末端号；3、支路阻抗；4、支路对地电纳
% 5、支路的变比；6、支路首端处于K侧为1，1侧为0
% B2矩阵：1、该节点发电机功率；2、该节点负荷功率；3、节点电压初始值% 4、PV节点电压V的给定值；5、节点所接的无功补偿设备的容量% 6、节点分类标号：1为平衡节点（应为1号节点）；2为PQ节点；% 3为PV节点；
clear;
n=10;%input('请输入节点数：n=');
nl=10;%input('请输入支路数：nl=');
isb=1;%input('请输入平衡母线节点号：isb=');
pr=0.00001;%input('请输入误差精度：pr=');
B1=[1 2 0.03512+0.08306i 0.13455i 1 0;
2 3 0.0068+0.18375i 0 1.02381 1;
1 4 0.05620+0.13289i 0.05382i 1 0;
4 5 0.00811+0.24549i 0 1.02381 1;
1 6 0.05620+0.13289i 0.05382i 1 0;
4 6 0.04215+0.09967i 0.04037i 1 0;
6 7 0.0068+0.18375i 0 1.02381 1;
6 8 0.02810+0.06645i 0.10764i 1 0;
8 10 0.00811+0.24549i 0 1 1;
8 9 0.03512+0.08306i 0.13455i 1 0]
B2=[0 0 1.1 1.1 0 1;
0 0 1 0 0 2;
0 0.343+0.21256i 1 0 0 2;
0 0 1 0 0 2;
0 0.204+0.12638i 1 0 0 2;
0 0 1 0 0 2;
0 0.306+0.18962i 1 0 0 2;
0 0 1 0 0 2;
0.5 0 1.1 1.1 0 3;
0 0.343+0.21256i 1 0 0 2]
;%input('请输入各节点参数形成的矩阵：B2=');
Y=zeros(n);e=zeros(1,n);f=zeros(1,n);V=zeros(1,n);sida=zeros(1,n);S1=zeros(nl);
% % %---------------------------------------------------
for i=1:nl %支路数
if B1(i,6)==0 %左节点处于1侧
p=B1(i,1);q=B1(i,2);
else %左节点处于K侧
p=B1(i,2);q=B1(i,1);
end
Y(p,q)=Y(p,q)-1./(B1(i,3)*B1(i,5)); %非对角元
Y(q,p)=Y(p,q); %非对角元
Y(q,q)=Y(q,q)+1./(B1(i,3)*B1(i,5)^2)+B1(i,4)./2; %对角元K侧
Y(p,p)=Y(p,p)+1./B1(i,3)+B1(i,4)./2; %对角元1侧
end
%求导纳矩阵
disp('导纳矩阵Y=');
disp(Y)
%----------------------------------------------------------
G=real(Y);B=imag(Y); %分解出导纳阵的实部和虚部
for i=1:n %给定各节点初始电压的实部和虚部
e(i)=real(B2(i,3));
f(i)=imag(B2(i,3));
V(i)=B2(i,4); %PV节点电压给定模值
end
for i=1:n %给定各节点注入功率
S(i)=B2(i,1)-B2(i,2); %i节点注入功率SG-SL
B(i,i)=B(i,i)+B2(i,5); %i节点无功补偿量
end
%=================================================================== P=real(S);Q=imag(S); %分解出各节点注入的有功和无功功率
ICT1=0;IT2=1;N0=2*n;N=N0+1;a=0; %迭代次数ICT1、a；不满足收敛要求的节点数IT2 while IT2~=0 % N0=2*n 雅可比矩阵的阶数；N=N0+1扩展列IT2=0;a=a+1;
for i=1:n
if i~=isb %非平衡节点
C(i)=0;D(i)=0;
for j1=1:n
C(i)=C(i)+G(i,j1)*e(j1)-B(i,j1)*f(j1);%Σ(Gij*ej-Bij*fj)
D(i)=D(i)+G(i,j1)*f(j1)+B(i,j1)*e(j1);%Σ(Gij*fj+Bij*ej)
end
P1=C(i)*e(i)+f(i)*D(i);%节点功率P计算eiΣ(Gij*ej-Bij*fj)+fiΣ(Gij*fj+Bij*ej)
Q1=C(i)*f(i)-e(i)*D(i);%节点功率Q计算fiΣ(Gij*ej-Bij*fj)-eiΣ(Gij*fj+Bij*ej) %求i节点有功和无功功率P',Q'的计算值
V2=e(i)^2+f(i)^2; %电压模平方
%========= 以下针对非PV节点来求取功率差及Jacobi矩阵元素========= if B2(i,6)~=3 %非PV节点
DP=P(i)-P1; %节点有功功率差
DQ=Q(i)-Q1; %节点无功功率差
%=============== 以上为除平衡节点外其它节点的功率计算=================
%================= 求取Jacobi矩阵===================
for j1=1:n
if j1~=isb&j1~=i %非平衡节点&非对角元
X1=-G(i,j1)*e(i)-B(i,j1)*f(i); % dP/de=-dQ/df
X2=B(i,j1)*e(i)-G(i,j1)*f(i); % dP/df=dQ/de
X3=X2; % X2=dp/df X3=dQ/de
X4=-X1; % X1=dP/de X4=dQ/df
p=2*i-1;q=2*j1-1;
J(p,q)=X3;J(p,N)=DQ;m=p+1; % X3=dQ/de J(p,N)=DQ节点无功功率差
J(m,q)=X1;J(m,N)=DP;q=q+1; % X1=dP/de J(m,N)=DP节点有功功率差
J(p,q)=X4;J(m,q)=X2; % X4=dQ/df X2=dp/df
elseif j1==i&j1~=isb %非平衡节点&对角元
X1=-C(i)-G(i,i)*e(i)-B(i,i)*f(i);% dP/de
X2=-D(i)+B(i,i)*e(i)-G(i,i)*f(i);% dP/df
X3=D(i)+B(i,i)*e(i)-G(i,i)*f(i); % dQ/de
X4=-C(i)+G(i,i)*e(i)+B(i,i)*f(i);% dQ/df
p=2*i-1;q=2*j1-1;J(p,q)=X3;J(p,N)=DQ;%扩展列△Q
m=p+1;
J(m,q)=X1;q=q+1;J(p,q)=X4;J(m,N)=DP;%扩展列△P
J(m,q)=X2;
end
end
else
%=============== 下面是针对PV节点来求取Jacobi矩阵的元素=========== DP=P(i)-P1; % PV节点有功误差
DV=V(i)^2-V2; % PV节点电压误差
for j1=1:n
if j1~=isb&j1~=i %非平衡节点&非对角元
X1=-G(i,j1)*e(i)-B(i,j1)*f(i); % dP/de
X2=B(i,j1)*e(i)-G(i,j1)*f(i); % dP/df
X5=0;X6=0;
p=2*i-1;q=2*j1-1;J(p,q)=X5;J(p,N)=DV; % PV节点电压误差
m=p+1;
J(m,q)=X1;J(m,N)=DP;q=q+1;J(p,q)=X6; % PV节点有功误差
J(m,q)=X2;
elseif j1==i&j1~=isb %非平衡节点&对角元
X1=-C(i)-G(i,i)*e(i)-B(i,i)*f(i);% dP/de
X2=-D(i)+B(i,i)*e(i)-G(i,i)*f(i);% dP/df
X5=-2*e(i);
X6=-2*f(i);
p=2*i-1;q=2*j1-1;J(p,q)=X5;J(p,N)=DV; % PV节点电压误差
m=p+1;
J(m,q)=X1;J(m,N)=DP;q=q+1;J(p,q)=X6; % PV节点有功误差
J(m,q)=X2;
end
end
end
end
end
%========= 以上为求雅可比矩阵的各个元素及扩展列的功率差或电压差=====================
for k=3:N0 % N0=2*n （从第三行开始，第一、二行是平衡节点）
k1=k+1;N1=N; % N=N0+1 即N=2*n+1扩展列△P、△Q 或△U
for k2=k1:N1 % 从k+1列的Jacobi元素到扩展列的△P、△Q 或△U
J(k,k2)=J(k,k2)./J(k,k);% 用K行K列对角元素去除K行K列后的非对角元素进行规格化
end
J(k,k)=1; % 对角元规格化K行K列对角元素赋1
%==================== 回代运算=======================================
if k~=3 % 不是第三行k > 3
k4=k-1;
for k3=3:k4 % 用k3行从第三行开始到当前行的前一行k4行消去for k2=k1:N1 % k3行后各行上三角元素
J(k3,k2)=J(k3,k2)-J(k3,k)*J(k,k2);%消去运算（当前行k列元素消为0）
end %用当前行K2列元素减去当前行k列元素乘以第k行K2列元素
J(k3,k)=0; %当前行第k列元素已消为0
end
if k==N0 %若已到最后一行
break;
end
%================== 前代运算==================================
for k3=k1:N0 % 从k+1行到2*n最后一行
for k2=k1:N1 % 从k+1列到扩展列消去k+1行后各行下三角元素
J(k3,k2)=J(k3,k2)-J(k3,k)*J(k,k2);%消去运算
end %用当前行K2列元素减去当前行k列元素乘以第k行K2列元素
J(k3,k)=0; %当前行第k列元素已消为0
end
else %是第三行k=3
%====================== 第三行k=3的前代运算========================
for k3=k1:N0 %从第四行到2n行（最后一行）
for k2=k1:N1 %从第四列到2n+1列（即扩展列）
J(k3,k2)=J(k3,k2)-J(k3,k)*J(k,k2);%消去运算（当前行3列元素消为0）
end %用当前行K2列元素减去当前行3列元素乘以第三行K2列元
素
J(k3,k)=0; %当前行第3列元素已消为0
end
end
end
%====上面是用线性变换方式高斯消去法将Jacobi矩阵化成单位矩阵===== for k=3:2:N0-1
L=(k+1)./2;
e(L)=e(L)-J(k,N); %修改节点电压实部
k1=k+1;
f(L)=f(L)-J(k1,N); %修改节点电压虚部
end
%------修改节点电压-----------
for k=3:N0
DET=abs(J(k,N));
if DET>=pr %电压偏差量是否满足要求
IT2=IT2+1; %不满足要求的节点数加1
end
end
ICT2(a)=IT2; %不满足要求的节点数
ICT1=ICT1+1; %迭代次数
end
%用高斯消去法解"w=-J*V"
disp('迭代次数：');
disp(ICT1);
disp('没有达到精度要求的个数：');
disp(ICT2);
for k=1:n
V(k)=sqrt(e(k)^2+f(k)^2); %计算各节点电压的模值
sida(k)=atan(f(k)./e(k))*180./pi; %计算各节点电压的角度
E(k)=e(k)+f(k)*j; %将各节点电压用复数表示
end
%=============== 计算各输出量=========================== disp('各节点的实际电压标幺值E为(节点号从小到大排列)：');
disp(E); %显示各节点的实际电压标幺值E用复数表示
disp('-----------------------------------------------------');
disp('各节点的电压大小V为(节点号从小到大排列)：');
disp(V); %显示各节点的电压大小V的模值
disp('-----------------------------------------------------');
disp('各节点的电压相角sida为(节点号从小到大排列)：');
disp(sida); %显示各节点的电压相角
for p=1:n
C(p)=0;
for q=1:n
C(p)=C(p)+conj(Y(p,q))*conj(E(q)); %计算各节点的注入电流的共轭值end
S(p)=E(p)*C(p); %计算各节点的功率S = 电压X 注入电流的共轭值end
disp('各节点的功率S为(节点号从小到大排列)：');
disp(S); %显示各节点的注入功率
disp('-----------------------------------------------------');
disp('各条支路的首端功率Si为(顺序同您输入B1时一致)：');
for i=1:nl
p=B1(i,1);q=B1(i,2);
if B1(i,6)==0
Si(p,q)=E(p)*(conj(E(p))*conj(B1(i,4)./2)+(conj(E(p)*B1(i,5))...
-conj(E(q)))*conj(1./(B1(i,3)*B1(i,5))));
Siz(i)=Si(p,q);
else
Si(p,q)=E(p)*(conj(E(p))*conj(B1(i,4)./2)+(conj(E(p)./B1(i,5))...
-conj(E(q)))*conj(1./(B1(i,3)*B1(i,5))));
Siz(i)=Si(p,q);
end
disp(Si(p,q));
SSi(p,q)=Si(p,q);
ZF=['S(',num2str(p),',',num2str(q),')=',num2str(SSi(p,q))];
disp(ZF);
disp('-----------------------------------------------------');
end
disp('各条支路的末端功率Sj为(顺序同您输入B1时一致)：');
for i=1:nl
p=B1(i,1);q=B1(i,2);
if B1(i,6)==0
Sj(q,p)=E(q)*(conj(E(q))*conj(B1(i,4)./2)+(conj(E(q)./B1(i,5))...
-conj(E(p)))*conj(1./(B1(i,3)*B1(i,5))));
Sjy(i)=Sj(q,p);
else
Sj(q,p)=E(q)*(conj(E(q))*conj(B1(i,4)./2)+(conj(E(q)*B1(i,5))...
-conj(E(p)))*conj(1./(B1(i,3)*B1(i,5))));
Sjy(i)=Sj(q,p);
end
disp(Sj(q,p));
SSj(q,p)=Sj(q,p);
ZF=['S(',num2str(q),',',num2str(p),')=',num2str(SSj(q,p))];
disp(ZF);
disp('-----------------------------------------------------');
end
disp('各条支路的功率损耗DS为(顺序同您输入B1时一致)：');
for i=1:nl
p=B1(i,1);q=B1(i,2);
DS(i)=Si(p,q)+Sj(q,p);
disp(DS(i));
DDS(i)=DS(i);
ZF=['DS(',num2str(p),',',num2str(q),')=',num2str(DDS(i))];
disp(ZF);
disp('-----------------------------------------------------');
end
figure(1);
subplot(1,2,1);
plot(V);
xlabel('节点号');ylabel('电压标幺值');
grid on;
subplot(1,2,2);
plot(sida);
xlabel('节点号');ylabel('电压角度');
grid on;
figure(2);
subplot(2,2,1);
P=real(S);Q=imag(S);
bar(P);
xlabel('节点号');ylabel('节点注入有功');
grid on;
subplot(2,2,2);
bar(Q);
xlabel('节点号');ylabel('节点注入无功');
grid on;
subplot(2,2,3);
P1=real(Siz);Q1=imag(Siz);
bar(P1);
xlabel('支路号');ylabel('支路首端注入有功');
grid on;
subplot(2,2,4);
bar(Q1);
xlabel('支路号');ylabel('支路首端注入无功');
grid
on;
B1 =
1.0000
2.0000 0.0351 + 0.0831i 0 + 0.1346i
1.0000 0
2.0000
3.0000 0.0068 + 0.1838i 0 1.0238 1.0000
1.0000 4.0000 0.0562 + 0.1329i 0 + 0.0538i 1.0000 0
4.0000
5.0000 0.0081 + 0.2455i 0 1.0238 1.0000
1.0000 6.0000 0.0562 + 0.1329i 0 + 0.0538i 1.0000 0
4.0000 6.0000 0.0422 + 0.0997i 0 + 0.0404i 1.0000 0
6.0000
7.0000 0.0068 + 0.1838i 0 1.0238 1.0000
6.0000 8.0000 0.0281 + 0.0664i 0 + 0.1076i 1.0000 0
8.0000 10.0000 0.0081 + 0.2455i 0 1.0000 1.0000
8.0000 9.0000 0.0351 + 0.0831i 0 + 0.1346i 1.0000 0
B2 =
0 0 1.1000 1.1000 0 1.0000
0 0 1.0000 0 0 2.0000
0 0.3430 + 0.2126i 1.0000 0 0 2.0000
0 0 1.0000 0 0 2.0000
0 0.2040 + 0.1264i 1.0000 0 0 2.0000
0 0 1.0000 0 0 2.0000
0 0.3060 + 0.1896i 1.0000 0 0 2.0000
0 0 1.0000 0 0 2.0000
0.5000 0 1.1000 1.1000 0 3.0000
0 0.3430 + 0.2126i 1.0000 0 0 2.0000
导纳矩阵Y=
9.7177 -22.8591i -4.3185 +10.2135i 0 -2.6996 + 6.3834i 0 -2.6996 + 6.3834i 0 0 0 0
-4.3185 +10.2135i 4.5104 -15.3311i -0.1964 + 5.3083i 0 0 0 0 0 0 0
0 -0.1964 + 5.3083i 0.2011 - 5.4347i 0 0 0 0 0 0 0
-2.6996 + 6.3834i 0 0 6.4271 -18.7292i -0.1313 + 3.9744i -3.5993 + 8.5110i 0 0 0 0
0 0 0 -0.1313 + 3.9744i 0.1344 - 4.0690i 0 0 0 0 0
-2.6996 + 6.3834i 0 0 -3.5993 + 8.5110i 0 11.8891 -32.7444i -0.1964 + 5.3083i -5.3984 +12.7660i 0 0
0 0 0 0 0 -0.1964 + 5.3083i 0.2011 - 5.4347i 0 0 0
0 0 0 0 0 -5.3984 +12.7660i 0 9.8514 -26.9275i -4.3185 +10.2135i -0.1344 + 4.0690i
0 0 0 0 0 0 0 -4.3185 +10.2135i 4.3185 -10.1462i 0
0 0 0 0 0 0 0 -0.1344 + 4.0690i 0 0.1344 - 4.0690i
迭代次数：
4
没有达到精度要求的个数：
17 18 17 0
各节点的实际电压标幺值E为(节点号从小到大排列)：
1.1000 1.0757 - 0.0207i 1.0050 - 0.0780i 1.0772 - 0.0175i 1.0171 -
0.0631i 1.0762 - 0.0152i 1.0112 - 0.0666i 1.0778 - 0.0051i 1.0996 + 0.0304i
1.0177 - 0.0814i
-----------------------------------------------------
各节点的电压大小V为(节点号从小到大排列)：
1.1000 1.0759 1.0080 1.0773 1.0191 1.0763 1.0134 1.0778 1.1000 1.0209
-----------------------------------------------------
各节点的电压相角sida为(节点号从小到大排列)：
0 -1.1046 -4.4373 -0.9283 -3.5503 -0.8106 -3.7665 -0.2718
1.5822 -4.5707
各节点的功率S为(节点号从小到大排列)：
0.7165 + 0.2587i 0.0000 + 0.0000i -0.3430 - 0.2126i -0.0000 + 0.0000i -0.2040 - 0.1264i -0.0000 + 0.0000i -0.3060 - 0.1896i -0.0000 + 0.0000i 0.5000 + 0.0089i -0.3430 - 0.2126i
-----------------------------------------------------
各条支路的首端功率Si为(顺序同您输入B1时一致)：
0.3485 + 0.0932i
S(1,2)=0.3485+0.093157i
-----------------------------------------------------
0.3441 + 0.2420i
S(2,3)=0.34409+0.24201i
-----------------------------------------------------
0.1904 + 0.0760i
S(1,4)=0.19038+0.07599i
-----------------------------------------------------
0.2044 + 0.1400i
S(4,5)=0.20445+0.13999i
-----------------------------------------------------
0.1777 + 0.0895i
S(1,6)=0.17767+0.089525i
-----------------------------------------------------
-0.0163 - 0.0055i
S(4,6)=-0.016305-0.0054856i
-----------------------------------------------------
0.3069 + 0.2128i
S(6,7)=0.30686+0.21281i
-----------------------------------------------------
-0.1477 - 0.0234i
S(6,8)=-0.14767-0.02338i
-----------------------------------------------------
0.3443 + 0.2509i
S(8,10)=0.34427+0.25091i
-----------------------------------------------------
-0.4925 - 0.1508i
S(8,9)=-0.49251-0.15077i
-----------------------------------------------------
各条支路的末端功率Sj为(顺序同您输入B1时一致)：-0.3441 - 0.2420i
S(2,1)=-0.34409-0.24201i
-----------------------------------------------------
-0.3430 - 0.2126i
S(3,2)=-0.343-0.21256i
-----------------------------------------------------
-0.1881 - 0.1345i
S(4,1)=-0.18815-0.13451i
-----------------------------------------------------
-0.2040 - 0.1264i
S(5,4)=-0.204-0.12638i
-----------------------------------------------------
-0.1755 - 0.1482i
S(6,1)=-0.17551-0.14815i
-----------------------------------------------------
0.0163 - 0.0413i
S(6,4)=0.016326-0.041272i
-----------------------------------------------------
-0.3060 - 0.1896i
S(7,6)=-0.306-0.18962i
-----------------------------------------------------
0.1482 - 0.1001i
S(8,6)=0.14824-0.10014i
-----------------------------------------------------
-0.3430 - 0.2126i
S(10,8)=-0.343-0.21256i
-----------------------------------------------------
0.5000 + 0.0089i
S(9,8)=0.5+0.0089402i
-----------------------------------------------------
各条支路的功率损耗DS为(顺序同您输入B1时一致)：
0.0044 - 0.1488i
DS(1,2)=0.0044095-0.14885i
-----------------------------------------------------
0.0011 + 0.0294i
DS(2,3)=0.0010897+0.029445i
-----------------------------------------------------
0.0022 - 0.0585i
DS(1,4)=0.0022306-0.058518i
-----------------------------------------------------
0.0004 + 0.0136i
DS(4,5)=0.00044972+0.013613i
-----------------------------------------------------
0.0022 - 0.0586i
DS(1,6)=0.0021584-0.058629i
-----------------------------------------------------
0.0000 - 0.0468i
DS(4,6)=2.1344e-005-0.046758i
-----------------------------------------------------
0.0009 + 0.0232i
DS(6,7)=0.00085804+0.023186i
-----------------------------------------------------
0.0006 - 0.1235i
DS(6,8)=0.00056584-0.12352i
-----------------------------------------------------
0.0013 + 0.0384i
DS(8,10)=0.001267+0.038353i
-----------------------------------------------------
0.0075 - 0.1418i
DS(8,9)=0.0074931-0.14183i
-----------------------------------------------------
华中科技大学 信息工程学院课程设计报告书题目: 电力系统潮流计算 专业：电气工程及其自动化 班级： 学号： 学生姓名： 指导教师： 2015年 11 月 10 日
2015年11月12日
信息工程学院课程设计成绩评定表
摘要 电力系统稳态分析包括潮流计算和静态安全分析。本文主要运用的事潮流计算，潮流计算是电力网络设计与运行中最基本的运算，对电力网络的各种设计方案及各种运行方式进行潮流计算，可以得到各种电网各节点的电压，并求得网络的潮流及网络中的各元件的电力损耗，进而求得电能损耗。本位就是运用潮流计算具体分析，并有MATLAB仿真。 关键词：电力系统潮流计算 MATLAB仿真
Abstract Electric power system steady flow calculation and analysis of the static safety analysis. This paper, by means of the calculation, flow calculation is the trend of the power network design and operation of the most basic operations of electric power network, various design scheme and the operation ways to tide computation, can get all kinds of each node of the power grid voltage and seek the trend of the network and the network of the components of the power loss, and getting electric power. The standard is to use the power flow calculation and analysis, the specific have MATLAB simulation. Key words: Power system; Flow calculation; MATLAB simulation
电力系统计算程序设计（包含matlab源程序）
广西大学电气工程学院 2007年1月 第一章原始数据 电力系统原始数据是电力系统计算的基础。电力系统每个计算程序都要求输入一定的原始数据，这些数据可以反映电力网络结构、电力系统正常运行条件、电力系统各元件参数和特性曲线。不同的计算程序需要不用的原始数据。 第一节电力网络的描述 电力网络是由输电线路、电力变压器、电容器和电抗器等元件组成。这些元件一般用集中参数的电阻、电抗和电容表示。为了表示电力网络中各元件是怎样互相连接的，通常要对网络节点进行编号。电力网络的结构和参数由电力网络中各支路的特性来描述。 1.1.1 线路参数 在电力系统程序设计中，线路参数一般采用线路的Π型数学模型，即线路用节点间的阻抗和节点对地容性电纳来表示，由于线路的对地电导很小，一般可忽略不计。其等价回路如下： r+jx -jb/2 对于线路参数的数据文件格式一般可写为： 线路参数（序号，节点i，节点j，r，x，b/2） 1.1.2 变压器参数
在电力系统程序设计中，变压器参数一般采用Π型等值变压器模型，这是一种可等值地体现变压器电压变换功能的模型。在多电压级网络计算中采用这种变压器模型后，就可不必进行参数和变量的归算。双绕组变压器的等值回路如下： k Z T k:1 Z T (a)接入理想变压器后的等值电路(b) 等值电路以导纳表示 (c) 等值电路以导纳表示 三绕组变压器的等值回路如下： 综合所述，三绕组变压器的等值电路可以用两个双绕组变压器的等值电
路来表示。因此，对于变压器参数的数据文件格式一般可写为：变压器参数（序号，节点i，节点j，r，x，k0） 其中，k0表示变压器变比。 1.1.3对地支路参数 对地支路参数一般以导纳形式表示，其等价回路如下： i g-jb 对地支路参数的数据文件格式一般可写为： 接地支路参数（序号，节点i，g i，b i） 第二节电力系统运行条件数据 电力系统运行条件数据包括发电机（含调相机）所连接的节点号、有功与无功功率；负荷所连接的节点号、有功与无功功率；PV节点与给定电压值；平衡节点的节点号与给定电压值。 1.2.1节点功率参数 电力系统中有流入流出功率的称为功率节点，有流入功率的称发电节点，一般为各发电站、枢纽变电站等节点；有流出功率的称负荷节点。对于电力系统稳态计算来说，功率节点都用有功功率P和无功功率Q来简单表示。其等价回路如下： Q G P G P L Q L 节点功率参数的数据文件格式一般可写为：
目录 摘要 (1) 1.任务及题目要求 (2) 2.计算原理 (3) 2.1牛顿—拉夫逊法简介 (3) 2.2牛顿—拉夫逊法的几何意义 (7) 3计算步骤 (7) 4.结果分析 (9) 小结 (11) 参考文献 (12) 附录：源程序 (13) 本科生课程设计成绩评定表 (32)
摘要 电力系统的出现，使高效，无污染，使用方便，易于调控的电能得到广泛应用，推动了社会生产各个领域的变化，开创了电力时代，发生率第二次技术革命。电力系统的规模和技术水准已经成为一个国家经济发展水平的标志之一。 电力系统稳态分析包括潮流计算和静态安全分析。电力系统潮流计算是电力系统最基本的计算，也是最重要的计算。所谓潮流计算，就是已知电网的接线方式与参数及运行条件，计算电力系统稳态运行各母线电压、个支路电流与功率及网损。对于正在运行的电力系统，通过潮流计算可以判断电网母线电压、支路电流和功率是否越限，如果有越限，就应采取措施，调整运行方式。对于正在规划的电力系统，通过潮流计算，可以为选择电网供电方案和电气设备提供依据。潮流计算还可以为继电保护和自动装置定整计算、电力系统故障计算和稳定计算等提供原始数据。 在电力系统运行方式和规划方案的研究中，都需要进行潮流计算以比较运行方式或规划供电方案的可行性、可靠性和经济性。同时，为了实时监控电力系统的运行状态，也需要进行大量而快速的潮流计算。因此，潮流计算是电力系统中应用最广泛、最基本和最重要的一种电气运算。在系统规划设计和安排系统的运行方式时，采用离线潮流计算；在电力系统运行状态的实时监控中，则采用在线潮流计算。 关键词：电力系统潮流计算牛顿-拉夫逊法
课程设计报告 学生:学号： 学院: 班级: 题目: 电力系统潮流计算课程设计
课设题目及要求 一 .题目原始资料 1、系统图：两个发电厂分别通过变压器和输电线路与四个变电所相连。 2、发电厂资料： 母线1和2为发电厂高压母线，发电厂一总装机容量为（ 300MW ），母线3为机压母线，机压母线上装机容量为（ 100MW ），最大负荷和最小负荷分别为50MW 和20MW ；发电厂二总装机容量为（ 200MW ）。 3、变电所资料： (一) 变电所1、2、3、4低压母线的电压等级分别为：35KV 10KV 35KV 10KV (二) 变电所的负荷分别为： 60MW 40MW 40MW 50MW (三) 每个变电所的功率因数均为cos φ=0.85； 变电所1 变电所母线 电厂一 电厂二
(四) 变电所1和变电所3分别配有两台容量为75MVA 的变压器，短路损 耗414KW ，短路电压（%）=16.7；变电所2和变电所4分别配有两台容 量为63MVA 的变压器，短路损耗为245KW ，短路电压（%）=10.5； 4、输电线路资料： 发电厂和变电所之间的输电线路的电压等级及长度标于图中，单位长度的电阻为Ω17.0，单位长度的电抗为Ω0.402，单位长度的电纳为S -610*2.78。 二、 课程设计基本容： 1. 对给定的网络查找潮流计算所需的各元件等值参数，画出等值电路图。 2. 输入各支路数据，各节点数据利用给定的程序进行在变电所在某一负荷 情况下的潮流计算，并对计算结果进行分析。 3. 跟随变电所负荷按一定比例发生变化，进行潮流计算分析。 1) 4个变电所的负荷同时以2%的比例增大； 2) 4个变电所的负荷同时以2%的比例下降 3) 1和4号变电所的负荷同时以2%的比例下降，而2和3号变电所的 负荷同时以2%的比例上升； 4. 在不同的负荷情况下，分析潮流计算的结果，如果各母线电压不满足要 求，进行电压的调整。（变电所低压母线电压10KV 要求调整围在9.5-10.5 之间；电压35KV 要求调整围在35-36之间） 5. 轮流断开支路双回线中的一条，分析潮流的分布。（几条支路断几次） 6. 利用DDRTS 软件，进行绘制系统图进行上述各种情况潮流的分析，并进 行结果的比较。 7. 最终形成课程设计成品说明书。 三、课程设计成品基本要求： 1. 在读懂程序的基础上画出潮流计算基本流程图 2. 通过输入数据，进行潮流计算输出结果 3. 对不同的负荷变化，分析潮流分布，写出分析说明。 4. 对不同的负荷变化，进行潮流的调节控制，并说明调节控制的方法，并 列表表示调节控制的参数变化。 5. 打印利用DDRTS 进行潮流分析绘制的系统图，以及潮流分布图。
第四章潮流计算的计算机算法 第一节概述 潮流计算是电力系统最基本、最常用的计算。根据系统给定的运行条件、网络接线及元件参数，通过潮流计算可以确定各母线的电压（幅值及相角），各元件中流过的功率、整个系统的功率损耗等。潮流计算是实现电力系统安全经济发供电的必要手段和重要工作环节。因此潮流计算在电力系统的规划设计、生产运行、调度管理及科学研究中都有着广泛的应用。 电力系统潮流计算分为离线潮流计算和在线潮流计算。前者主要用于系统规划设计和安排系统的运行方式，后者则用于正在运行系统的经常监视及实时控制。本章主要讨论离线潮流计算问题，它的基本算法同样适用于在线潮流计算。 潮流计算在数学上是多元非线性方程组的求解问题，求解的方法有很多种。自从五十年代计算机应用于电力系统以来，当时求解潮流的方法是以节点导纳矩阵为基础的逐次代入法（导纳法），后来为解决导纳法的收敛性较差的问题，出现了以阻抗矩阵为基础的逐次代入法（阻抗法）。到六十年代，针对阻抗法占用计算机内存大的问题又出现了分块阻抗法及牛顿－拉夫逊（Newton－Raphson）法。Newton —Raphson法是数学上解非线形方程式的有效方法，有较好的收敛性。将N－R法用于潮流计算是以导纳矩阵为基础的，由于利用了导纳矩阵的对称性、稀疏性及节点编号顺序优化等技巧，使N－R法在收敛性、占用内存、计算速度方面的优点都超过了阻抗法，成为六十年代末期以后普遍采用的方法。同时国内外广泛研究了诸如非线形规划法、直流法、交流法等各种不同的潮流计算方法。七十年代以来，又涌现出了更新的潮流计算方法。其中有1974年由B、Stott、O、Alsac 提出的快速分解法以及1978年由岩本伸一等提出的保留非线性的高 129
附录1 使用牛顿拉夫逊法进行潮流计算的Matlab程序代码 % 牛拉法计算潮流程序 %----------------------------------------------------------------------- % B1矩阵：1、支路首端号；2、末端号；3、支路阻抗；4、支路对地电纳 % 5、支路的变比；6、支路首端处于K侧为1，1侧为0 % B2矩阵：1、该节点发电机功率；2、该节点负荷功率；3、节点电压初始值 % 4、PV节点电压V的给定值；5、节点所接的无功补偿设备的容量 % 6、节点分类标号：1为平衡节点（应为1号节点）；2为PQ节点；3为PV节点； %------------------------------------------------------------------------ clear all; format long; n=input('请输入节点数:nodes='); nl=input('请输入支路数:lines='); isb=input('请输入平衡母线节点号:balance='); pr=input('请输入误差精度:precision='); B1=input('请输入由各支路参数形成的矩阵:B1='); B2=input('请输入各节点参数形成的矩阵:B2='); Y=zeros(n);e=zeros(1,n);f=zeros(1,n);V=zeros(1,n);sida=zeros(1,n);S1=zeros(nl); %------------------------------------------------------------------ for i=1:nl %支路数 if B1(i,6)==0 %左节点处于1侧 p=B1(i,1);q=B1(i,2); else %左节点处于K侧 p=B1(i,2);q=B1(i,1); end Y(p,q)=Y(p,q)-1./(B1(i,3)*B1(i,5)); %非对角元 Y(q,p)=Y(p,q); %非对角元 Y(q,q)=Y(q,q)+1./(B1(i,3)*B1(i,5)^2)+B1(i,4); %对角元K侧 Y(p,p)=Y(p,p)+1./B1(i,3)+B1(i,4); %对角元1侧 end %求导纳矩阵 disp('导纳矩阵 Y='); disp(Y) %------------------------------------------------------------------- G=real(Y);B=imag(Y); %分解出导纳阵的实部和虚部 for i=1:n %给定各节点初始电压的实部和虚部
题目：潮流计算与matlab 教学单位电气信息学院姓名 学号 年级 专业电气工程及其自动化指导教师 职称副教授
摘要 电力系统稳态分析包括潮流计算和静态安全分析。本文主要运用的事潮流计算，潮流计算是电力网络设计与运行中最基本的运算，对电力网络的各种设计方案及各种运行方式进行潮流计算，可以得到各种电网各节点的电压，并求得网络的潮流及网络中的各元件的电力损耗，进而求得电能损耗。本位就是运用潮流计算具体分析，并有MATLAB仿真。 关键词：电力系统潮流计算 MATLAB Abstract Electric power system steady flow calculation and analysis of the static safety analysis. This paper, by means of the calculation, flow calculation is the trend of the power network design and operation of the most basic operations of electric power network, various design scheme and the operation ways to tide computation, can get all kinds of each node of the power grid voltage and seek the trend of the network and the network of the components of the power loss, and getting electric power. The standard is to use the power flow calculation and analysis, the specific have MATLAB simulation. Key words: Power system; Flow calculation; MATLAB simulation
第四章 电力系统潮流分析与计算 电力系统潮流计算是电力系统稳态运行分析与控制的基础，同时也是安全性分析、稳定性分析电磁暂态分析的基础（稳定性分析和电磁暂态分析需要首先计算初始状态，而初始状态需要进行潮流计算）。其根本任务是根据给定的运行参数，例如节点的注入功率，计算电网各个节点的电压、相角以及各个支路的有功功率和无功功率的分布及损耗。 潮流计算的本质是求解节点功率方程，系统的节点功率方程是节点电压方程乘以节点电压构成的。要想计算各个支路的功率潮流，首先根据节点的注入功率计算节点电压，即求解节点功率方程。节点功率方程是一组高维的非线性代数方程，需要借助数字迭代的计算方法来完成。简单辐射型网络和环形网络的潮流估算是以单支路的潮流计算为基础的。 本章主要介绍电力系统的节点功率方程的形成，潮流计算的数值计算方法，包括高斯迭代法、牛顿拉夫逊法以及PQ 解藕法等。介绍单电源辐射型网络和双端电源环形网络的潮流估算方法。 4-1 潮流计算方程--节点功率方程 1. 支路潮流 所谓潮流计算就是计算电力系统的功率在各个支路的分布、各个支路的功率损耗以及各个节点的电压和各个支路的电压损耗。由于电力系统可以用等值电路来模拟，从本质上说，电力系统的潮流计算首先是根据各个节点的注入功率求解电力系统各个节点的电压，当各个节点的电压相量已知时，就很容易计算出各个支路的功率损耗和功率分布。 假设支路的两个节点分别为k 和l ，支路导纳为kl y ，两个节点的电压已知，分别为k V 和l V ，如图4-1所示。 图4-1 支路功率及其分布 那么从节点k 流向节点l 的复功率为（变量上面的“－”表示复共扼）： )]([l k kl k kl k kl V V y V I V S (4-1) 从节点l 流向节点k 的复功率为： )]([k l kl l lk l lk V V y V I V S (4-2) 功率损耗为： 2)()(kl kl l k kl l k lk kl kl V y V V y V V S S S (4-3)
潮流例题：根据给定的参数或工程具体要求（如图），收集和查阅资料；学习相关软件（软件自选：本设计选择Matlab进行设计）。 2.在给定的电力网络上画出等值电路图。 3.运用计算机进行潮流计算。 4.编写设计说明书。 一、设计原理 1．牛顿-拉夫逊原理 牛顿迭代法是取x0 之后，在这个基础上，找到比x0 更接近的方程的跟，一步一步迭代，从而找到更接近方程根的近似跟。牛顿迭代法是求方程根的重要方法之一，其最大优点是在方程f(x) = 0 的单根附近具有平方收敛，而且该法还可以用来求方程的重根、复根。电力系统潮流计算，一般来说，各个母线所供负荷的功率是已知的，各个节点电压是未知的（平衡节点外）可以根据网络结构形成节点导纳矩阵，然后由节点导纳矩阵列写功率方程，由于功率方程里功率是已知的，电压的幅值和相角是未知的，这样潮流计算的问题就转化为求解非线性方程组的问题了。为了便于用迭代法解方程组，需要将上述功率方程改写成功率平衡方程，并对功率平衡方程求偏导，得出对应的雅可比矩阵，给未知节点赋电压初值，一般为额定电压，将初值带入功率平衡方程，得到功率不平衡量，这样由功率不平衡量、雅可比矩阵、节点电压不平衡量（未知的）构成了误差方程，解误差方程，得到节点电压不平衡量，节点电压加上节点电压不平衡量构成新的节点电压初值，将新的初值带入原来的功率平衡方程，并重新形成雅可比矩阵，然后计算新
的电压不平衡量，这样不断迭代，不断修正，一般迭代三到五次就能收敛。 牛顿—拉夫逊迭代法的一般步骤： （1）形成各节点导纳矩阵Y。 （2）设个节点电压的初始值U和相角初始值e 还有迭代次数初值为0。 （3）计算各个节点的功率不平衡量。 （4）根据收敛条件判断是否满足，若不满足则向下进行。 （5）计算雅可比矩阵中的各元素。 （6）修正方程式个节点电压 （7）利用新值自第（3）步开始进入下一次迭代，直至达到精度退出循环。 （8）计算平衡节点输出功率和各线路功率 2．网络节点的优化 １）静态地按最少出线支路数编号 这种方法由称为静态优化法。在编号以前。首先统计电力网络个节点的出线支路数，然后，按出线支路数有少到多的节点顺序编号。当由n 个节点的出线支路相同时，则可以按任意次序对这n 个节点进行编号。这种编号方法的根据是导纳矩阵中，出线支路数最少的节点所对应的行中非零元素也２）动态地按增加出线支路数最少编号在上述的方法中，各节点的出线支路数是按原始网络统计出来的，在编号过程中认为固定不变的，事实上，在节点消去过程中，每消去一个节点以后，与该节点相连的各节点的出线支路数将发生变化(增加，减少或保持不变)。因此，如果每消去一个节点后，立即修正尚未编号节点的出线支路数，然后选其中支路数最少的一个节点进行编号，就可以预期得到更好的效果，动态按最少出线支路数编号方法的特点就是按出线最少原则编号时考虑了消去过程中各节点出线支路数目的变动情况。 3．MATLAB编程应用 Matlab 是“Matrix Laboratory”的缩写，主要包括：一般数值分析，矩阵运算、数字信号处理、建模、系统控制、优化和图形显示等应用程序。由于使用Matlab 编程运算与人进行科学计算的思路和表达方式完全一致，所以不像学习高级语言那样难于掌握，而且编程效率和计算效率极高，还可在计算机上直接输出结果和精美的图形拷贝，所以它的确为一高效的科研助手。 二、设计内容 1.设计流程图
电力系统分析潮流计算报告
目录 一．配电网概述 (3) 1.1 配电网的分类 (3) 1.2 配电网运行的特点及要求 (3) 1.3 配电网潮流计算的意义 (4) 二．计算原理及计算流程 (4) 2.1 前推回代法计算原理 (4) 2.2 前推回代法计算流程 (7) 2.3主程序清单： (9) 2.4 输入文件清单： (11) 2.5计算结果清单： (12) 三．前推回代法计算流程图 (13) 参考文献 (14)
一．配电网概述 1.1 配电网的分类 在电力网中重要起分配电能作用的网络就称为配电网； 配电网按电压等级来分类，可分为高压配电网（35—110KV），中压配电网（6—10KV，苏州有20KV的），低压配电网（220/380V）； 在负载率较大的特大型城市，220KV电网也有配电功能。 按供电区的功能来分类，可分为城市配电网，农村配电网和工厂配电网等。 在城市电网系统中,主网是指110KV及其以上电压等级的电网，主要起连接区域高压（220KV及以上）电网的作用。 配电网是指35KV及其以下电压等级的电网，作用是给城市里各个配电站和各类用电负荷供给电源。 从投资角度看，我国与国外先进国家的发电、输电、配电投资比率差异很大，国外基本上是电网投资大于电厂投资，输电投资小于配电投资。我国刚从重发电轻供电状态中转变过来，而在供电投资中，输电投资大于配电投资。从我国城网改造之后，将逐渐从输电投资转入配电建设为主。 本文是基于前推回代法的配电网潮流分析计算的研究，研究是是以根节点为10kV的电压等级的配电网。 1.2 配电网运行的特点及要求 配电系统相对于输电系统来说，由于电压等级低、供电范围小，但与用户直接相连，是供电部门对用户服务的窗口，因而决定了配电网运行有如下特点和基本要求：
%*************************电力系统直角坐标系下的牛顿拉夫逊法潮流计算********** clear clc load E:\data\IEEE014_Node.txt Node=IEEE014_Node; weishu=size(Node); nnum=weishu(1,1); %节点总数 load E:\data\IEEE014_Branch.txt branch=IEEE014_Branch; bwei=size(branch); bnum=bwei(1,1); %支路总数 Y=(zeros(nnum)); Sj=100; %********************************节点导纳矩阵******************************* for m=1:bnum; s=branch(m,1); %首节点 e=branch(m,2); %末节点 R=branch(m,3); %支路电阻 X=branch(m,4); %支路电抗 B=branch(m,5); %支路对地电纳 k=branch(m,6); if k==0 %无变压器支路情形 Y(s,e)=-1/(R+j*X); %互导纳 Y(e,s)=Y(s,e); end if k~=0 %有变压器支路情形 Y(s,e)=-(1/((R+j*X)*k)); Y(e,s)=Y(s,e); Y(s,s)=-(1-k)/((R+j*X)*k^2); Y(e,e)=-(k-1)/((R+j*X)*k); %对地导纳 end Y(s,s)=Y(s,s)-j*B/2; Y(e,e)=Y(e,e)-j*B/2; %自导纳的计算情形 end for t=1:nnum; Y(t,t)=-sum(Y(t,:))+Node(t,12)+j*Node(t,13); %求支路自导纳 end G=real(Y); %电导 B=imag(Y); %电纳 %******************节点分类************************************* * pq=0; pv=0; blancenode=0; pqnode=zeros(1,nnum); pvnode=zeros(1,nnum); for m=1:nnum; if Node(m,2)==3 blancenode=m; %平衡节点编号 else if Node(m,2)==0 pq=pq+1; pqnode(1,pq)=m; %PQ 节点编号 else if Node(m,2)==2 pv=pv+1; pvnode(1,pv)=m; %PV 节点编号 end end end end %*****************************设置电压初值********************************** Uoriginal=zeros(1,nnum); %对各节点电压矩阵初始化 for n=1:nnum Uoriginal(1,n)=Node(n,9); %对各点电压赋初值 if Node(n,9)==0;
《电力系统潮流上机》课程设计报告 院系：电气工程学院 班级：电088班 学号： 0812002221 学生姓名：刘东昇 指导教师：张新松 设计周数：两周 日期：2010年 12 月 25 日
一、课程设计的目的与要求 目的：培养学生的电力系统潮流计算机编程能力，掌握计算机潮流计算的相关知识 要求：基本要求： 1.编写潮流计算程序； 2.在计算机上调试通过； 3.运行程序并计算出正确结果； 4.写出课程设计报告 二、设计步骤： 1.根据给定的参数或工程具体要求（如图），收集和查阅资料；学习相关软件（软件自选：本设计选择Matlab进行设计）。 2.在给定的电力网络上画出等值电路图。 3.运用计算机进行潮流计算。 4.编写设计说明书。 三、设计原理 1．牛顿-拉夫逊原理 牛顿迭代法是取x0 之后，在这个基础上，找到比x0 更接近的方程的跟，一步一步迭代，从而找到更接近方程根的近似跟。牛顿迭代法是求方程根的重要方法之一，其最大优点是在方程f(x) = 0 的单根附近具有平方收敛，而且该法还可以用来求方程的重根、复根。电力系统潮流计算，一般来说，各个母线所供负荷的功率是已知的，各个节点电压是未知的（平衡节点外）可以根据网络结构形成节点导纳矩阵，然后由节点导纳矩阵列写功率方程，由于功率方程里功率是已知的，电压的幅值和相角是未知的，这样潮流计算的问题就转化为求解非线性方程组的问题了。为了便于用迭代法解方程组，需要将上述功率方程改写成功率平衡方程，并对功率平衡方程求偏导，得出对应的雅可比矩阵，给未知节点赋电压初值，一般为
额定电压，将初值带入功率平衡方程，得到功率不平衡量，这样由功率不平衡量、雅可比矩阵、节点电压不平衡量（未知的）构成了误差方程，解误差方程，得到节点电压不平衡量，节点电压加上节点电压不平衡量构成新的节点电压初值，将新的初值带入原来的功率平衡方程，并重新形成雅可比矩阵，然后计算新的电压不平衡量，这样不断迭代，不断修正，一般迭代三到五次就能收敛。 牛顿—拉夫逊迭代法的一般步骤： （1）形成各节点导纳矩阵Y。 （2）设个节点电压的初始值U和相角初始值e 还有迭代次数初值为0。 （3）计算各个节点的功率不平衡量。 （4）根据收敛条件判断是否满足，若不满足则向下进行。 （5）计算雅可比矩阵中的各元素。 （6）修正方程式个节点电压 （7）利用新值自第（3）步开始进入下一次迭代，直至达到精度退出循环。 （8）计算平衡节点输出功率和各线路功率 2．网络节点的优化 １）静态地按最少出线支路数编号 这种方法由称为静态优化法。在编号以前。首先统计电力网络个节点的出线支路数，然后，按出线支路数有少到多的节点顺序编号。当由n 个节点的出线支路相同时，则可以按任意次序对这n 个节点进行编号。这种编号方法的根据是导纳矩阵中，出线支路数最少的节点所对应的行中非零元素也２）动态地按增加出线支路数最少编号在上述的方法中，各节点的出线支路数是按原始网络统计出来的，在编号过程中认为固定不变的，事实上，在节点消去过程中，每消去一个节点以后，与该节点相连的各节点的出线支路数将发生变化(增加，减少或保持不变)。因此，如果每消去一个节点后，立即修正尚未编号节点的出线支路数，然后选其中支路数最少的一个节点进行编号，就可以预期得到更好的效果，动态按最少出线支路数编号方法的特点就是按出线最少原则编号时考虑了消去过程中各节点出线支路数目的变动情况。 3．MATLAB编程应用 Matlab 是“Matrix Laboratory”的缩写，主要包括：一般数值分析，矩阵运算、数字信号处理、建模、系统控制、优化和图形显示等应用程序。由于使用Matlab 编程运算与人进行科学计算的思路和表达方式完全一致，所以不像学习高级语言那样难于掌握，而且编程效率和计算效率极高，还可在计算机上直接输出结果和精美的图形拷贝，所以它的确为一高效的科研助手。 四、设计内容
高等电力系统分析 （潮流计算的计算机算法）PQ分解法潮流计算（IEEE14)
目录 一、MATLAB源程序 二、对支路参数（B1）、节点参数（B2）的说明 三、带入数据，运行结果
一、MATLAB源程序 clear close all n=input('请输入节点数:n='); n1=input('请输入支路数:n1='); isb=input('请输入平衡节点号:isb='); pr=input('请输入误差精度:pr='); B1=input('请输入支路参数:B1='); B2=input('请输入节点参数:B2='); n2=input('请输入PQ节点个数:n2='); Y=zeros(n); for i=1:n1 p=B1(i,1); q=B1(i,2); Y(p,q)=Y(p,q)-1/(B1(i,3)+B1(i,4)*1j); %非对角元 Y(q,p)=Y(p,q); Y(p,p)=Y(p,p)+1/(B1(i,3)+B1(i,4)*1j)+B1(i,6)*1j; %对角元 Y(q,q)=Y(q,q)+1/(B1(i,3)+B1(i,4)*1j)+B1(i,6)*1j; end disp('导纳矩阵Y='); disp(Y) %--------------------------------------------- %---------------下面是求P,Q,V,O矩阵--------------- V=zeros(1,n);O=zeros(1,n);P=zeros(1,n);Q=zeros(1,n); G=real(Y);B=imag(Y); for i=1:n P(i)=B2(i,3); Q(i)=B2(i,4); V(i)=B2(i,5); O(i)=B2(i,6); end B3=B(1:n-1,1:n-1); %不含平衡节点，由节点导纳虚部构成 B4=B(1:n2,1:n2); %所有PQ节点 %---------------------------------------------- %---------------下面是求ΔP,ΔQ矩阵--------------- DX=0;ICT=1;Mp=1;Mq=1; while ICT~=0 m1=1;m2=1; for i=1:n
电力系统分析课程设计 设计题目9节点电力网络潮流计算 指导教师 院（系、部）电气与控制工程学院 专业班级 学号 姓名 日期
电气工程系课程设计标准评分模板
目录 1 PSASP软件简介 (1) 1.1 PSASP平台的主要功能和特点 (1) 1.2 PSASP的平台组成 (2) 2 牛顿拉夫逊潮流计算简介 (3) 2.1 牛顿—拉夫逊法概要 (3) 2.2 直角坐标下的牛顿—拉夫逊潮流计算 (5) 2.3 牛顿—拉夫逊潮流计算的方法 (6) 3 九节点系统单线图及元件数据 (7) 3.1 九节点系统单线图 (7) 3.2 系统各项元件的数据 (8) 4 潮流计算的结果 (10) 4.1 潮流计算后的单线图 (10) 4.2 潮流计算结果输出表格 (10) 5 结论 (14)
电力系统分析课程设计任务书9节点系统单线图如下： 基本数据如下：
表3 两绕组变压器数据 负荷数据
1 PSASP软件简介 “电力系统分析综合程序”（Power System Analysis Software Package,PSASP）是一套历史悠久、功能强大、使用方便的电力系统分析程序，是高度集成和开发具有我国自主知识产权的大型软件包。 基于电网基础数据库、固定模型库以及用户自定义模型库的支持，PSASP可进行电力系统（输电、供电和配电系统）的各种计算分析，目前包括十多个计算机模块，PSASP的计算功能还在不断发展、完善和扩充。 为了便于用户使用以及程序功能扩充，在PSASP7.0中设计和开发了图模一体化支持平台，应用该平台可以方便地建立电网分析的各种数据，绘制所需要的各种电网图形（单线图、地理位置接线图、厂站主接线图等）；该平台服务于PSASP 的各种计算，在此之外可以进行各种分析计算，并输出各种计算结果。 1.1PSASP平台的主要功能和特点 PSASP图模一体化支持平台的主要功能和特点可概括为： 1. 图模支持平台具备MDI多文档操作界面，是一个单线图图形绘制、元件数据录入编辑、各种计算功能、结果显示、报表和曲线输出的集成环境。用户可以方便地建立电网数据、绘制电网图形、惊醒各种分析计算。人机交互界面全部汉化，界面良好，操作方便。 2. 真正的实现了图模一体化。可边绘图边建数据，也可以在数据已知的情况下进行图形自动快速绘制；图形、数据自动对应，所见即所得。 3. 应用该平台可以绘制各种电网图形，包括单线图、地理位置接线图、厂站主接线图等。 ●所有图形独立于各种分析计算，并为各计算模块所共享； ●可在图形上进行各种计算操作，并在图上显示各种计算结果； ●同一系统可对应多套单线图，多层子图嵌套； ●单线图上可细化到厂站主接线结构；
课程设计论文 基于MATLAB的电力系统潮流计算 学院:电气工程学院 专业:电气工程及自动化 班级:电自0710班 学号:0703110304 姓名: 马银莎
内容摘要 潮流计算是电力系统最基本最常用的计算。根据系统给定的运行条件，网络接线及元件参数，通过潮流计算可以确定各母线的电压（幅值和相角），各支路流过的功率，整个系统的功率损耗。潮流计算是实现电力系统安全经济发供电的必要手段和重要工作环节。因此，潮流计算在电力系统的规划计算，生产运行，调度管理及科学计算中都有着广泛的应用。 潮流计算在数学上是多元非线性方程组的求解问题，牛顿—拉夫逊Newton-Raphson法是数学上解非线性方程组的有效方法，有较好的收敛性。运用电子计算机计算一般要完成以下几个步骤：建立数学模型，确定解算方法，制订计算流程，编制计算程序。 关键词 牛顿-拉夫逊法（Newton-Raphson）变压器及非标准变比无功调节 高斯消去法潮流计算Mtlab
一 .电力系统潮流计算的概述 在电力系统的正常运行中，随着用电负荷的变化和系统运行方式的改变，网络中的损耗也将发生变化。要严格保证所有的用户在任何时刻都有额定的电压是不可能的，因此系统运行中个节点出现电压的偏移是不可避免的。为了保证电力系统的稳定运行，要进行潮流调节。 随着电力系统及在线应用的发展，计算机网络已经形成，为电力系统的潮流计算提供了物质基础。电力系统潮流计算是电力系统分析计算中最基本的内容，也是电力系统运行及设计中必不可少的工具。根据系统给定的运行条件、网络接线及元件参数，通过潮流计算可以确定各母线电压的幅值及相角、各元件中流过的功率、整个系统的功率损耗等。潮流计算是实现电力系统安全经济发供电的必要手段和重要工作环节，因此潮流计算在电力系统的规划设计、生产运行、调度管理及科学研究中都有着广泛的应用。它的发展主要围绕这样几个方面：计算方法的收敛性、可靠性；计算速度的快速性；对计算机存储容量的要求以及计算的方便、灵活等。 常规的电力系统潮流计算中一般具有三种类型的节点：PQ 、PV 及平衡节点。一个节点有四个变量，即注入有功功率、注入无功功率，电压大小及相角。常规的潮流计算一般给定其中的二个变量：PQ 节点（注入有功功率及无功功率），PV 节点（注入有功功率及电压的大小），平衡节点（电压的大小及相角）。 1、变量的分类： 负荷消耗的有功、无功功率——1L P 、1L Q 、2L P 、2L Q 电源发出的有功、无功功率——1G P 、1G Q 、2G P 、2G Q 母线或节点的电压大小和相位——1U 、2U 、1δ、2δ 在这十二个变量中，负荷消耗的有功和无功功率无法控制，因它们取决于用户，它们就称为不可控变量或是扰动变量。电源发出的有功无功功率是可以控制的自变量，因此它们就称为控制变量。母线或节点电压的大小和相位角——是受控制变量控制的因变量。其中, 1U 、2U 主要受1G Q 、2G Q 的控制, 1δ、2δ主要受 1G P 、2G P 的控制。这四个变量就是简单系统的状态变量。 为了保证系统的正常运行必须满足以下的约束条件： 对控制变量 max min max min ;Gi Gi Gi Gi Gi Gi Q Q Q P P P <<<< 对没有电源的节点则为 0;0==Gi Gi Q P 对状态变量i U 的约束条件则是 m a x m i n i i i U U U <<
程序为计算书3-4的过程 程序可以解决开式单直网络和树状网络的计算。树状网络计算时要自己先设定好支路的起始节点和终止节点标号以及计算顺序源代码： #include
基于MATLAB的电力系统潮流计算 %简单潮流计算的小程序，相关的原始数据数据数据输入格式如下： %B1是支路参数矩阵，第一列和第二列是节点编号。节点编号由小到大编写%对于含有变压器的支路，第一列为低压侧节点编号，第二列为高压侧节点%编号，将变压器的串联阻抗置于低压侧处理。 %第三列为支路的串列阻抗参数。 %第四列为支路的对地导纳参数。 %第五烈为含变压器支路的变压器的变比 %第六列为变压器是否是否含有变压器的参数，其中“1”为含有变压器，%“0”为不含有变压器。 %B2为节点参数矩阵，其中第一列为节点注入发电功率参数；第二列为节点%负荷功率参数；第三列为节点电压参数；第六列为节点类型参数，其中 %“1”为平衡节点，“2”为PQ节点，“3”为PV节点参数。 %X为节点号和对地参数矩阵。其中第一列为节点编号，第二列为节点对地%参数。 n=input('请输入节点数:n='); n1=input('请输入支路数:n1='); isb=input('请输入平衡节点号:isb='); pr=input('请输入误差精度:pr='); B1=input('请输入支路参数:B1='); B2=input('请输入节点参数:B2='); X=input('节点号和对地参数:X='); Y=zeros(n); Times=1; %置迭代次数为初始值 %创建节点导纳矩阵 for i=1:n1 if B1(i,6)==0 %不含变压器的支路 p=B1(i,1); q=B1(i,2); Y(p,q)=Y(p,q)-1/B1(i,3); Y(q,p)=Y(p,q); Y(p,p)=Y(p,p)+1/B1(i,3)+0.5*B1(i,4); Y(q,q)=Y(q,q)+1/B1(i,3)+0.5*B1(i,4); else %含有变压器的支路 p=B1(i,1); q=B1(i,2); Y(p,q)=Y(p,q)-1/(B1(i,3)*B1(i,5)); Y(q,p)=Y(p,q); Y(p,p)=Y(p,p)+1/B1(i,3);