电力系统分析大作业matlab三机九节点潮流计算报告.docx

1、电力系统分析大作业matlab三机九节点潮流计算报告电力系统分析大作业一、设计题目本次设计题目选自课本第五章例5-8,美国西部联合电网WSCC系统的简化三机九节点系统，例题中已经给出了潮流结果，计算结果可以与之对照。取=0.00001。二、计算步骤第一步，为了方便编程，修改节点的序号，将平衡节点放在最后。如下图：第二步，这样得出的系统参数如下表所示：第三步，形成节点导纳矩阵。第四步，设定初值：；，。第五步，计算失配功率=0，=-1.25，=-0.9，=0，=-1，=0，=1.63，=0.85；=0.8614，=-0.2590，=-0.0420，=0.6275，=-0.1710，=0.7101。

2、显然，。第六步，形成雅克比矩阵（阶数为1414）第七步，解修正方程，得到：-0.0371，-0.0668，-0.0628，0.0732，0.0191，0.0422，0.1726，0.0908；0.0334，0.0084，0.0223，0.0372，0.0266，0.0400。从而-0.0371，-0.0668，-0.0628，0.0732，0.0191，0.0422，0.1726，0.0908；1.0334，1.0084，1.0223，1.0372，1.0266，1.0400。然后转入下一次迭代。经三次迭代后。迭代过程中失配功率的变化情况如下表：k0123=P10-0.01060.00010.

3、0000000845=P2-1.250.03790.00050.0000000896=P3-0.90.04390.00050.0000000778=P40-0.0421-0.0012-0.0000003421=P5-10.0610.00090.0000001845=P60-0.0269-0.0007-0.0000001631=P71.63-0.0579-0.0004-0.0000000278=P80.85-0.0336-0.0002-0.0000000103=Q10.8614-0.0501-0.0004-0.0000000561=Q2-0.259-0.0714-0.0012-0.0000002

4、774=Q3-0.042-0.0424-0.0006-0.0000001236=Q40.6275-0.1875-0.0021-0.0000003279=Q5-0.171-0.0241-0.0004-0.0000000805=Q60.7101-0.0828-0.0007-0.0000000799max1.630.0610.00090.0000001845迭代过程中节点电压变化情况如下表：kU1U2U3U4U5U6011111111.03341.00841.02231.03721.02661.040021.02590.99581.01281.02591.01601.032431.02580.995

5、61.01271.02581.01591.0324迭代收敛后各节点的电压和功率：kUPQ11.0258-2.21680.00000.000020.9956-3.9888-1.2500-0.500031.0127-3.6874-0.9000-0.300041.02583.71970.00000.000051.01590.7275-1.0000-0.350061.03241.96670.00000.000071.02509.28001.63000.066581.02504.66480.8500-0.108691.04000.00000.71640.2705同课本上给出的潮流相比较，结果完全一致，证

6、明计算过程与程序编写正确。最后得出迭代收敛后各支路的功率和功率损耗：ijPijQijIijPjiQjiIjiPLQL120.40940.22890.4572-0.4068-0.38690.56390.0026-0.1579130.30700.01030.2995-0.3054-0.16540.34300.0017-0.155124-0.8432-0.11310.85450.8662-0.08380.84840.0230-0.196936-0.5946-0.13460.60200.6082-0.18070.61460.0135-0.3153450.7638-0.00800.7447-0.7590

7、-0.10700.75460.0048-0.115056-0.2410-0.24300.33680.24180.03120.23620.0009-0.2118910.71640.27050.7363-0.7164-0.23920.73630.00000.0312741.63000.06651.5916-1.63000.09181.59160.00000.1583860.8500-0.10860.8360-0.85000.14960.83600.00000.0410三、源程序及注释由于计算流程比较简单，所以编写程序过程中没有采用模块化的形式，直接按顺序一步步进行。disp(【节点数:】);n1=

8、xlsread(input.xls,A3:A3)%节点数disp(【支路数:】);n=xlsread(input.xls,B3:B3)%支路数disp(【精度:】);Accuracy=xlsread(input.xls,B4:B4)%精度branch=xlsread(input.xls,E4:K12);node=xlsread(input.xls,M4:S12);Data_B1=branch;%支路参数Data_B2=node;%节点参数T1=zeros(n,2);T2=zeros(n1,3);i=sqrt(-1);formatshortforj=1:nT1(j,1)=Data_B1(j,3)

9、+Data_B1(j,4)*1i;T1(j,2)=Data_B1(j,5)*1i;endforj=1:n1T2(j,1)=Data_B2(j,1)+Data_B2(j,2)*1i;T2(j,2)=Data_B2(j,3)+Data_B2(j,4)*1i;endB1=zeros(n,6);B2=zeros(n1,5);forj=1:nB1(j,1)=Data_B1(j,1);B1(j,2)=Data_B1(j,2);B1(j,3)=T1(j,1);B1(j,4)=T1(j,2);B1(j,5)=Data_B1(j,6);B1(j,6)=Data_B1(j,7);endforj=1:n1B2(j,

10、1)=T2(j,1);B2(j,2)=T2(j,2);B2(j,3)=Data_B2(j,5);B2(j,4)=Data_B2(j,6);B2(j,5)=Data_B2(j,7);enddisp(【支路参数矩阵:】);B1%显示支路参数矩阵disp(【节点参数矩阵:】);B2%显示节点参数矩阵%以上为从excel中导入初值的程序Y=zeros(n1);fori=1:nifB1(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

11、(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);Y(q,q)=Y(q,q)+1/(B1(i,5)2*B1(i,3);endenddisp(【导纳矩阵:】);Y%显示导纳矩阵m=0;fori=1:n1ifB2(i,5)=2m=m+1;endendm%PQ节点个数l=0;fori=1:n1ifB2(i,5)=1l=l+1;endendl%PV节点个数Mismatch_po

12、wer=zeros(l+m*2,1);fori=1:n1-1Pj=0;forj=1:n1Pj=Pj+(B2(i,3)*B2(j,3)*(real(Y(i,j)*cos(B2(i,4)-B2(j,4)+imag(Y(i,j)*sin(B2(i,4)-B2(j,4);endMismatch_power(i,1)=real(B2(i,1)-real(B2(i,2)-Pj;endfork=n1:(l+m*2)Qj=0;forj=1:n1Qj=Qj+B2(k-n1+1),3)*B2(j,3)*(real(Y(k-n1+1),j)*sin(B2(k-n1+1),4)-B2(j,4)-imag(Y(k-n1

13、+1),j)*cos(B2(k-n1+1),4)-B2(j,4);endMismatch_power(k,1)=imag(B2(k-n1+1),1)-imag(B2(k-n1+1),2)-Qj;end%Mismatch_power%计算失配功率times=0;while(max(Mismatch_power)Accuracy)fori=1:(n1-1)Pj=0;forj=1:n1Pj=Pj+B2(i,3)*B2(j,3)*(real(Y(i,j)*cos(B2(i,4)-B2(j,4)+imag(Y(i,j)*sin(B2(i,4)-B2(j,4);endMismatch_power(i,1)

14、=real(B2(i,1)-real(B2(i,2)-Pj;endfork=n1:(l+m*2)Qj=0;forj=1:n1Qj=Qj+B2(k-n1+1),3)*B2(j,3)*(real(Y(k-n1+1),j)*sin(B2(k-n1+1),4)-B2(j,4)-imag(Y(k-n1+1),j)*cos(B2(k-n1+1),4)-B2(j,4);endMismatch_power(k,1)=imag(B2(k-n1+1),1)-imag(B2(k-n1+1),2)-Qj;enddisp(【当前迭代次数:】);timesdisp(【失配功率:】);Mismatch_powerJacob

15、ian=zeros(l+m*2);%雅克比矩阵7*7%Hfori=1:(n1-1)forj=1:(n1-1)ifi=jP_H=0;fork=1:n1P_H=P_H+B2(i,3)*B2(k,3)*(real(Y(i,k)*sin(B2(i,4)-B2(k,4)-imag(Y(i,k)*cos(B2(i,4)-B2(k,4);endJacobian(i,i)=P_H-B2(i,3)*B2(i,3)*(0-imag(Y(i,i);elseJacobian(i,j)=0-B2(i,3)*B2(j,3)*(real(Y(i,j)*sin(B2(i,4)-B2(j,4)-imag(Y(i,j)*cos(

16、B2(i,4)-B2(j,4);endendend%Nfori=1:(n1-1)forj=1:mifi=jP_N=0;fork=1:n1P_N=P_N+B2(k,3)*(real(Y(i,k)*cos(B2(i,4)-B2(k,4)+imag(Y(i,k)*sin(B2(i,4)-B2(k,4);endJacobian(i,n1-1+i)=0-B2(i,3)*real(Y(i,i)-P_N;elseJacobian(i,n1-1+j)=0-B2(i,3)*(real(Y(i,j)*cos(B2(i,4)-B2(j,4)+imag(Y(i,j)*sin(B2(i,4)-B2(j,4);enden

17、dend%Kfori=1:mforj=1:(n1-1)ifi=jP_K=0;fork=1:n1P_K=P_K+B2(i,3)*B2(k,3)*(real(Y(i,k)*cos(B2(i,4)-B2(k,4)+imag(Y(i,k)*sin(B2(i,4)-B2(k,4);endJacobian(n1-1+i,i)=0+B2(i,3)*B2(i,3)*real(Y(i,i)-P_K;elseJacobian(n1-1+i,j)=B2(i,3)*B2(j,3)*(real(Y(i,j)*cos(B2(i,4)-B2(j,4)+imag(Y(i,j)*sin(B2(i,4)-B2(j,4);ende

18、ndend%Lfori=1:mforj=1:mifi=jP_L=0;fork=1:n1P_L=P_L+B2(k,3)*(real(Y(i,k)*sin(B2(i,4)-B2(k,4)-imag(Y(i,k)*cos(B2(i,4)-B2(k,4);endJacobian(n1-1+i,n1-1+i)=0-P_L+B2(i,3)*imag(Y(i,i);elseJacobian(n1-1+i,n1-1+j)=0-B2(i,3)*(real(Y(i,j)*sin(B2(i,4)-B2(j,4)-imag(Y(i,j)*cos(B2(i,4)-B2(j,4);endendendS=zeros(l+m

19、*2,1);%初始化电压角度变化量S=inv(Jacobian)*(0-Mismatch_power);%求解修正方程S=(Jacobian)(0-Mismatch_power);%求解修正方程fori=1:(n1-1)%角度初值加变化量B2(i,4)=B2(i,4)+S(i,1);endfori=1:m%电压初值加变化量B2(i,3)=B2(i,3)+S(n1-1+i,1);enddisp(【雅克比矩阵:】);Jacobian%显示雅克比矩阵%S=inv(Jacobian)times=times+1;endtimes=times-1;disp(【共计迭代次数:】);times%显示迭代次数U

20、_It=zeros(n1,1);%初始化电压向量fori=1:n1U_It(i,1)=B2(i,3)*cos(B2(i,4)+B2(i,3)*sin(B2(i,4)*1j;endangle_It=zeros(n1,1);%将电压角度的弧度值转为角度值fori=1:n1angle_It(i,1)=B2(i,4)*180/pi;endNode_S_It=U_It.*(conj(Y)*conj(U_It);%求解节点功率disp(【迭代收敛后各节点的电压幅值:】);Node_U_It=abs(U_It)%显示迭代收敛后各节点的电压幅值disp(【迭代收敛后各节点的电压角度:】);angle_It%显

21、示迭代收敛后各节点的电压角度disp(【迭代收敛后各节点的功率:】);Node_S_It%显示迭代收敛后各节点的功率Branch_It=zeros(n,10);fori=1:n;ifB1(i,6)=0;%不带变压器支路m=B1(i,1);%得到支路号n=B1(i,2);Branch_It(i,1)=m;%显示支路号Branch_It(i,2)=n;a=U_It(m,1)*(conj(U_It(m,1)*conj(B1(i,4)*0.5+(conj(U_It(m,1)-conj(U_It(n,1)/conj(B1(i,3);Branch_It(i,3)=real(a);%显示PijBranch_

22、It(i,4)=imag(a);%显示Qijb=U_It(m,1)*B1(i,4)*0.5+(U_It(m,1)-U_It(n,1)/B1(i,3);Branch_It(i,5)=sqrt(real(b)2+imag(b)2);%显示Iijc=U_It(n,1)*(conj(U_It(n,1)*conj(B1(i,4)*0.5+(conj(U_It(n,1)-conj(U_It(m,1)/conj(B1(i,3);Branch_It(i,6)=real(c);%显示PjiBranch_It(i,7)=imag(c);%显示Qjid=U_It(n,1)*B1(i,4)*0.5+(U_It(n,1

23、)-U_It(m,1)/B1(i,3);Branch_It(i,8)=sqrt(real(d)2+imag(d)2);%显示Ijie=a+c;Branch_It(i,9)=real(e);%显示线路损耗有功分量Branch_It(i,10)=imag(e);%显示线路损耗无功分量else%带变压器支路(同以上内容)m=B1(i,1);n=B1(i,2);Branch_It(i,1)=m;Branch_It(i,2)=n;a=U_It(m,1)*(conj(U_It(m,1)/conj(B1(i,3)-conj(U_It(n,1)*conj(1/(B1(i,5)*B1(i,3);Branch_I

24、t(i,3)=real(a);Branch_It(i,4)=imag(a);b=U_It(m,1)*(B1(i,5)-1)/B1(i,3)/B1(i,5)+(U_It(m,1)-U_It(n,1)/(B1(i,5)*B1(i,3);Branch_It(i,5)=sqrt(real(b)2+imag(b)2);c=U_It(n,1)*(conj(U_It(n,1)/(conj(B1(i,5)*B1(i,5)*B1(i,3)-conj(U_It(m,1)*conj(1/(B1(i,5)*B1(i,3);Branch_It(i,6)=real(c);Branch_It(i,7)=imag(c);d=

25、U_It(n,1)*(1-B1(i,5)/B1(i,5)/B1(i,5)/B1(i,3)+(U_It(n,1)-U_It(m,1)/B1(i,5)/B1(i,3);Branch_It(i,8)=sqrt(real(d)2+imag(d)2);e=a+c;Branch_It(i,9)=real(e);Branch_It(i,10)=imag(e);endenddisp(【迭代收敛后各支路的功率和功率损耗:】);Branch_It%显示迭代收敛后各支路的功率和功率损耗%向Excel表中输出数据%Node_S_It_Real=real(Node_S_It);%Node_S_It_imag=imag(

26、Node_S_It);%xlswrite(output.xls,Node_U_It,1,B3);%xlswrite(output.xls,angle_It,1,C3);%xlswrite(output.xls,Node_S_It_Real,1,D3);%xlswrite(output.xls,Node_S_It_imag,1,E3);%xlswrite(output.xls,Branch_It,1,G3);程序中还有将数据从Excel表格中读入输出的xlsread和xlswrite功能，直接将数据输入到Excel表格中，可以省略将数据写在程序中或者一一输入的步骤，适用于任何节点的电力系统潮流计算。四、程序运行结果五、手算结果（第一次迭代）六、总结通过采用计算机和手算进行潮流计算，我对潮流计算的计算过程和MATLAB软件的使用有了更深层次的了解。我们已经将如此复杂的问题通过矩阵这样的方式得以简化，但仍然有庞大的计算过程是难以通过手算方式进行解决的。同时也深深感叹于电力系统的庞大而精细，为自己能在以后为之付出感到期待。