《基于MATLAB牛顿拉夫逊法进行潮流计算》由会员分享，可在线阅读，更多相关《基于MATLAB牛顿拉夫逊法进行潮流计算（10页珍藏版）》请在金锄头文库上搜索。
1、电力系统分析N-R 迭代法计算潮流分布学 院： 信息与控制工程学院专 业： 电气工程及其自动化班 级： 电气 09-1 班 姓 名： 朱守文 学 号： 09053129流程图如下：NR迭代：在Matlab中设计的程序如下n=input(请输入节点数:n=);nl=input(请输 入支路数:nl=);isb=input(请输入平衡母线节点号:isb=);pr=input(请输 入误差精度 :pr=);B1=input(请输 入由各支路参数形成的矩 阵:B1=);B2=input(请输 入各节点参数形成的矩 阵:B2=);Y=zeros(n); e=zeros(1,n);f=zeros(1,n)
2、;V=zeros(1,n); 输入原始数据形成节点导纳矩阵按公式计算雅可比矩阵各元素计算平衡节点功率及全部线路功率输出给定节点电压初值 (0)(),iiefk用公式计算 ()()2(),kkiiiPQV及()()()max| |?kiii解修正方程式，求 ()(),kiief(1)()()(1)()(),kkkkiiiiiief1k是否O=zeros(1,n);S1=zeros(nl);for i=1:nlif B1(i,6)=0p=B1(i,1);q=B1(i,2);else p=B1(i,2);q=B1(i,1);endY(p,q)=Y(p,q)-1./(B1(i,3)*B1(i,5);Y
3、(q,p)=Y(p,q);Y(q,q)=Y(q,q)+1./(B1(i,3)*B1(i,5)2)+B1(i,4)./2;Y(p,p)=Y(p,p)+1./B1(i,3)+B1(i,4)./2;end%求 导纳矩 阵disp(导纳矩阵Y=);disp(Y);G=real(Y);B=imag(Y);for i=1:ne(i)=real(B2(i,3);f(i)=imag(B2(i,3);V(i)=B2(i,4);endfor i=1:nS(i)=B2(i,1)-B2(i,2);B(i,i)=B(i,i)+B2(i,5);endP=real(S);Q=imag(S);ICT1=0;IT2=1;N0=
4、2*n;N=N0+1;a=0;while IT2=0IT2=0;a=a+1;for i=1:nif i=isb C(i)=0;D(i)=0;for j1=1:nC(i)= C(i)+G(i,j1)*e(j1)-B(i,j1)*f(j1);D(i)= D(i)+G(i,j1)*f(j1)+B(i,j1)*e(j1);endP1=C(i)*e(i)+f(i)*D(i);Q1=f(i)*C(i)-D(i)*e(i);V2=e(i)2+f(i)2;if B2(i,6)=3DP=P(i)-P1;DQ=Q(i)-Q1;for j1=1:nif j1=isb&j1=iX1=-G(i,j1)*e(i)-B(i
5、,j1)*f(i);X2=B(i,j1)*e(i)-G(i,j1)*f(i); X3=X2;X4=-X1;p=2*i-1;q=2*j1-1;J(p,q)=X3;J(p,N)=DQ;m=p+1;J(m,q)=X1;J(m,N)=DP;q=q+1;J(p,q)=X4;J(m,q)=X2;elseif j1=i&j1=isbX1=-C(i)-G(i,i)*e(i)-B(i,i)*f(i);X2=-D(i)+B(i,i)*e(i)-G(i,i)*f(i);X3=D(i)+B(i,i)*e(i)-G(i,i)*f(i);X4=-C(i)+G(i,i)*e(i)+B(i,i)*f(i);p=2*i-1;=
6、2*j1-1;J(p,q)=X3;J(p,N)=DQ;m=p+1;J(m,q)=X1;J(m,N)=DP;q=q+1;J(p,q)=X4;J(m,q)=X2;endendelseDP=P(i)-P1;DV=V(i)2-V2;for j1=1:nif j1=isb&j1=iX1=-G(i,j1)*e(i)-B(i,j1)*f(i);X2=B(i,j1)*e(i)-G(i,j1)*f(i);X5=0;X6=0;p=2*i-1;q=2*j1-1;J(p,q)=X5;J(p,N)=DV; m=p+1;J(m,q)=X1;J(m,N)=DP;q=q+1;J(p,q)=X6;J(m,q)=X2;elsei
7、f j1=i&j1=isbX1=-C(i)-G(i,i)*e(i)-B(i,i)*f(i);X2=-D(i)+B(i,i)*e(i)-G(i,i)*f(i);X5=-2*e(i);X6=-2*f(i);p=2*i-1;q=2*j1-1;J(p,q)=X5;J(p,N)=DV;m=p+1;J(m,q)=X1;J(m,N)=DP;q=q+1;J(p,q)=X6;J(m,q)=X2;endendendend end %求雅可比矩 阵for k=3:N0 k1=k+1;N1=N; for k2=k1:N1J(k,k2)=J(k,k2)./J(k,k);endJ(k,k)=1;if k=3;k4=k-1
8、;for k3=3:k4for k2=k1:N1J(k3,k2)= J(k3,k2)-J(k3,k)*J(k,k2);endJ(k3,k)=0;endif k=N0,break;endfor k3=k1:N0for k2=k1:N1J(k3,k2)=J(k3,k2)-J(k3,k)*J(k,k2);endJ(k3,k)=0;endelsefor k3=k1:N0for k2=k1:N1J(k3,k2)= J(k3,k2)-J(k3,k)*J(k,k2);endJ(k3,k)=0;endendendfor k=3:2:N0-1 L=(k+1)./2;e(L)=e(L)-J(k,N);k1=k+1
9、;f(L)=f(L)-J(k1,N);endfor k=3:N0DET=abs(J(k,N);if DET=prIT2=IT2+1;endendICT2(a)=IT2; ICT1=ICT1+1;End %用高斯消去法解 “w=-J*V”disp(迭代次数);disp(ICT1);disp(没有达到精度要求的个数);disp(ICT2);for k=1:nV(k)=sqrt(e(k)2+f(k)2);sita(k)=atan(f(k)./e(k)*180/pi;E(k)=e(k)+f(k)*j;enddisp(各节点的实际电压标么值E为( 节点号从小到大排列 ):);disp(E);disp(各
10、节点的电压大小V为( 节点号从小到大排列 ):);disp(V);disp(各节点的电压相角时为( 节点号从小到大排列 ):);disp(sita);for p=1:nC(p)=0;for q=1:nC(p)=C(p)+conj(Y(p,q)*conj(E(q);endS(p)=E(p)*C(p);enddisp(各节点的功率S为( 节点号从小到大排列 ):);disp(S);disp(各条支路的首端功率Si为(顺序同您输 入B1时一样):);for i=1:nlif B1(i,6)=0p=B1(i,1);q=B1(i,2);else p=B1(i,2);q=B1(i,1);endSi(p,q
11、)=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);disp(Si(p,q);enddisp (各条支路的末端功率Sj为(顺序同您输入B1时一样):);for i=1:nlif B1(i,6)=0p=B1(i,1);q=B1(i,2);else p=B1(i,2);q=B1(i,1);endSj(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(
12、i,5);disp(Sj(q,p);enddisp(各条支路的功率损耗DS为(顺序同您输 入B1时一样): );for i=1:nlif B1(i,6)=0p=B1(i,1);q=B1(i,2);else p=B1(i,2);q=B1(i,1);endDS(i)=Si(p,q)+Sj(q,p);disp(DS(i);end程序中B1矩阵的每行有以下参数构成：某支路的首端号P;某支路末端号Q，且PQ;支路的阻抗（R+jX） ；支路的对地容抗；支路的变比K；折算到哪一侧的标志（如果支路的首端P处处于高压侧则输入“1”否则输入“0” ） 。程序中 B2 矩阵的每行有以下参数构成：节点所接发电机的功率
13、；节点负荷的功率；节点电压的初始值；PV 节点电压的给定值；节点所接的无功补偿设备的容量；节点分类标号：1平衡节点；2PQ 节点；3PV 节点实验程序运行结果如下：请输入节点数:n=4请输入支路数:nl=4请输入平衡母线节点号:isb=4请输入误差精度:pr=0.00001请输入由各支路参数形成的矩阵:B1=1 2 0.1+0.4i 0.01528i 1 1;1 3 0.3i 0 1.1 0;1 4 0.12+0.5i 0.0192i 1 0;2 4 0.08+0.4i 0.01413i 1 0请输入各节点参数形成的矩阵:B2=0 0.3+0.18i 1 0 0 2;0 0.55+0.13i
14、1 0 0 2;0.5 0 1.1 1.1 0 3;0 0 1.05 1.05 0 1导纳矩阵 Y=1.0421 - 7.5601i -0.5882 + 2.3529i 0 + 3.0303i -0.4539 + 1.8911i-0.5882 + 2.3529i 1.0690 - 4.7421i 0 -0.4808 + 2.4038i0 + 3.0303i 0 0 - 2.7548i 0 -0.4539 + 1.8911i -0.4808 + 2.4038i 0 0.9346 - 4.2783i迭代次数4没有达到精度要求的个数3 3 2 0 0各节点的实际电压标么值 E 为(节点号从小到大排列
15、):1.0000 0.9645 - 0.1028i 1.0876 + 0.1650i 1.0500 各节点的电压大小 V 为(节点号从小到大排列):1.0000 0.9700 1.1000 1.0500各节点的电压相角时 为(节点号从小到大排列 ):0 -6.0857 8.6269 0各节点的功率 S 为(节点号从小到大排列):-0.2598 - 0.0510i -0.5500 - 0.1300i 0.5000 + 0.0377i 0.3265 + 0.2448i 各条支路的首端功率 Si 为( 顺序同您输入 B1 时一样):-0.2559 - 0.0024i-0.5000 + 0.0377i-0.0227 - 0.1042i-0.2941 - 0.1276i各条支路的末端功率 Sj 为( 顺序同您输入 B1 时一样):0.2628
