matlab潮流计算毕业设计(编辑修改稿)

matlab潮流计算毕业设计(编辑修改稿)内容摘要：

i=1 j=1 while i=^m+1 if j=n i=i+1 elseif i=j J(i,j)=[(dot G(i,:)*edot B(i,:)*f)G(i,i)*e(i)B(i,i)*f(i) (dot G(i,:)*f+dot B(i,:)*e)+B(i,i)*e(i)G(i,i)*f(i)。 dot G(i,:)*f+dot B(i,:)*e+B(i,i)*e(i)G(i,i)*f(i) (dot G(i,:)*edot B(i,:)*f)+G(i,i)*e(i)+B(i,i)*f(i)]。 j=j+1 else J(i,j)=[(G(i,j)*e(i)+B(i,j)f(i)) B(i,j)*e(i)G(i,j)*f(i)。 B(i,j)*e(i)G(i,j)*f(i) G(i,j)*e(i)+B(i,j)f(i)] j=j+1 end end while i=^n if j=n i=i+1 elseif i=j J(i,j)=[(dot G(i,:)*edot B(i,:)*f)G(i,i)*e(i)B(i,i)*f(i) (dot G(i,:)*f+dot B(i,:)*e)+B(i,i)*e(i)G(i,i)*f(i)。 2*e(i) 2*f(i)]。 j=j+1 else J(i,j)=[(G(i,j)*e(i)+B(i,j)f(i)) B(i,j)*e(i)G(i,j)*f(i)。 0 0]。 j=j+1 end end 计算过程及数据分析 本文用 MATLAB 结合牛顿 拉夫逊法，牛顿 拉夫逊法是求解非线性方程有效的迭代计算方法，在牛顿 拉夫逊法的每一次迭代过程中，非线性问题通过线性化逐步近似，可提高其精确度。 潮流计算中复杂矩阵的输入问题可通过创建 M 文件来解决 ； MATLAB 称为矩阵实验室 ， 其能进行潮流计算中的各种矩阵运算 ， 程 序的编写因 MATLAB 提供了许多功能函数而变得简单易行。 并且 MATLAB 提供了可视化技术，使图形和数据联系起来。 南华大学 电气工程 学院毕业设计 第 15 页 共 63 页 以 IEEE6BUS 标准试验系统的潮流计算进行仿真 ,具体计算中系统节点分三类 :(1)PQ 节点 ,即节点的有功功率 P 和无功功率 Q 已知 ,节点电压 (V,δ )未知。 (2)PV节点 ,即节点的有功功率 P 和电压幅值 V 已知 ,节点的无功功率 Q 和电压的相位δ未知。 (3)平衡节点 ,即节点的电压幅值 V和相位δ已知 ,节点的有功功率 P和无功功率 Q未知 .图 1 所示的算例中 ,节点总数为 6(n=6)。 PQ 节点有 4 个 (m=4),如节 点 1~4。 PV 节点有 1个 (nm1=1),如节点 5。 平衡节点有 1个 ,如节点 PQ节点 ,再PV 节点 ,最后平衡节点的顺序进行编号 ,即 :1,2,„ ,m 为 PQ 节点。 m+1,m+2,„ ,n1 为PV节点。 n 为平衡节点 .参数如表 所示 表 系统参数表 Line Number Bus From Number To R X Tap Ratio 1 1 2 0.000 0.300 0.025 2 1 4 0.097 0.407 3 1 6 0.123 0.518 4 2 5 0.828 0.640 5 3 5 0.723 1.050 6 4 3 0.000 0.133 1.100 7 4 6 0.080 0.370 供电网络化简并作出等值电路 如图 所示 网络参数： 南华大学 电气工程 学院毕业设计 第 16 页 共 63 页 0||,||0||||||||)0(6)0(5)0(4)0(3)0(2)0(1)0(6)0(5)0(4)0(3)0(2)0(1VVVVVV 负荷参数 , ,, 544 332211   PQP QPQPQP 图 IEE－ 6系统图 节点导纳矩阵数据如 表 表 节点导纳矩阵 0 + 0 + 0 + 0 + 0 0 + 0 0 0 0 + + 0 + 0 0 + 0 + 0 + + 0 0 + 0 0 + 0 用公式 ： ),2,1()()()()(1 11 1 mieBfGefBeGfeBfGffBeGePPPPnjnjjijjiijijijiisiisinjnjjijjiijijijiisiisi     计算 iP , iQ 南华大学 电气工程 学院毕业设计 第 17 页 共 63 页 用公式 ： )1,2,1()()()(222222 1 1    nmmifeVVVVeBfGffBeGePPPPiiisiisinjnj jijjiijijijiisiisi计算 iP , iV 形成雅克比矩阵 ：  13211,12,11,11,332311,222211,112111321*nnnnnnnnn VVVVJJJJJJJJJJJJWWWW 解之可得 121 ,  nVVV 由 11121211 ,   nnnn fVeVfVeV 可得第一次迭代误差 111111111111 ,   nnnnnn fffeeefffeee 继续以上计算，直到满足条件。 计算全线功率和平衡节点功率。 MATLAB 潮流 计算程序 潮流计算程序如下： clear clc Z12=。 R12=real(Z12)。 X12=imag(Z12)。 R14=。 X14=。 R16=。 X16=。 R25=。 X25=。 R35=。 X35=。 R43=0。 X43=。 R46=。 X46=。 K1=。 K2=。 y12=1/(R12+X12*j)。 y14=1/(R14+X14*j)。 y16=1/(R16+X16*j)。 y25=1/(R25+X25*j)。 y35=1/(R35+X35*j)。 y43=1/(R43+X43*j)。 y46=1/(R46+X46*j)。 Y(1,1)=y16+y14+y12/(K1*K1)。 Y(2,1)=y12/K1。 Y(1,2)=Y(2,1)。 Y(4,1)=y14。 Y(1,4)=Y(4,1)。 Y(6,1)=y16。 Y(1,6)=Y(6,1)。 Y(2,2)=y12+y25。 Y(5,2)=y25。 Y(2,5)=Y(5,2)。 Y(3,3)=y35+y43。 Y(4,3)=y43/K2。 Y(3,4)=Y(4,3)。 Y(5,3)=y35。 Y(3,5)=Y(5,3)。 Y(4,4)=y14+y46+y43/(K2*K2)。 Y(6,4)=y46。 Y(4,6)=Y(6,4)。 Y(5,5)=y35+y25。 Y(6,6)=y16+y46。 G=real(Y)。 B=imag(Y)。 S1=input(39。 节点 1功率 =39。 ) S2=input(39。 节点 2功率 =39。 ) S3=input(39。 节点 3功率 =39。 ) S4=input(39。 节点 4功率 =39。 ) P(5)=input(39。 节点 5有功 =39。 ) e(5)=input(39。 节点 5电压 =39。 ) 南华大学 电气工程 学院毕业设计 第 18 页 共 63 页 e(6)=input(39。 节点 6电压 =39。 ) x=input(39。 误差允许范围 39。 ) P(1)=。 Q(1)=。 P(2)=。 Q(2)=。 P(3)=。 Q(3)=。 P(4)=。 Q(4)=。 P(5)=。 m=4。 n=6。 V1=1。 e(1)=real(V1)。 f(1)=imag(V1)。 e(2)=1。 e(3)=1。 e(4)=1。 f(2)=0。 f(3)=0。 f(4)=0。 f(5)=0。 f(6)=0。 v(n1)=e(n1)。 k=0 while 35 h=1。 while h~=m+1 P1(h)=P(h)e(h)*(dot(G(h,:),e)dot(B(h,:),f))f(h)*(dot(G(h,:),f)+dot(B(h,:),e))。 Q1(h)=Q(h)f(h)*(dot(G(h,:),e)dot(B(h,:),f))+e(h)*(dot(G(h,:),f)+dot(B(h,:),e))。 W(2*h1)=P1(h)。 W(2*h)=Q1(h)。 h=h+1。 end while h~=n P1(h)=P(h)e(h)*(dot(G(h,:),e)dot(B(h,:),f))f(h)*(dot(G(h,:),f)+dot(B(h,:),e))。 V1(h)=v(h)*v(h)(e(h)*e(h)+f(h)*f(h))。 W(2*h1)=P1(h)。 W(2*h)=V1(h)。 h=h+1。 end plot(k,P1(1),39。 *39。 ) hold on plot(k,Q1(1),39。 +39。 ) hold on Xlabel(39。 k39。 ) Ylabel(39。 P(*) Q(+)39。 ) if abs(W(1))x|abs(W(2))x|abs(W(3))x|abs(W(4))x|abs(W(5))x|abs(W(6))x|abs(W(7))x|abs(W(8))x|abs(W(9))x|abs(W(10))x i=1。 y=1 while i~=m+1 if y==n i=i+1 y=1 elseif i==y J(2*i1,2*y1)=(dot(G(i,:),e)dot(B(i,:),f))G(i,i)*e(i)B(i,i)*f(i) J(2*i1,2*y)=(dot(G(i,:),f)+dot(B(i,:),e))+B(i,i)*e(i)G(i,i)*f(i) J(2*i,2*y1)=dot(G(i,:),f)+dot(B(i,:),e)+B(i,i)*e(i)G(i,i)*f(i) J(2*i,2*y)=(dot(G(i,:),e)dot(B(i,:),f))+G(i,i)*e(i)+B(i,i)*f(i) y=y+1 else J(2*i1,2*y1)=(G(i,y)*e(i)+B(i,y)*f(i)) J(2*i1,2*y)=B(i,y)*e(i)G(i,y)*f(i) J(2*i,2*y1)=B(i,y)*e(i)G(i,y)*f(i) J(2*i,2*y)=G(i,y)*e(i)+B(i,y)*f(i) y=y+1 南华大学 电气工程 学院毕业设计 第 19 页 共 63 页 end end while i~=n if y==n i=i+1 y=1 elseif i==y J(2*i1,2*y1)=(dot(G(i,:),e)dot(B(i,:),f))G(i,i)*e(i)B(i,i)*f(i) J(2*i1,2*y)=(dot(G(i,:),f)+dot(B(i,:),e))+B(i,i)*e(i)G(i,i)*f(i) J(2*i,2*y1)=2*e(i) J(2*i,2*y)=2*f(i) y=y+1 else J(2*i1,2*y1)=(G(i,y)*e(i)+B(i,y)*f(i)) J(2*i1,2*y)=B(i,y)*e(i)G。