《基于MATLAB电力系统PQ潮流计算程序设计.doc》由会员分享，可免费在线阅读全文，更多与《基于MATLAB电力系统PQ潮流计算程序设计》相关文档资源请在帮帮文库(www.woc88.com)数亿文档库存里搜索。

1、=:nC()=C()+conj(Y(,q))*conj(E(q))；endS()=E()*C()；enddis(‘各节点的功率S(节点号从小到大排)：‘)为；dis(S)；dis(‘各条支路的首端功率Si为(顺序同您输入B时的一样)：‘)为；fori=:nlifB(i,)===B(i,)；q=B(i,)；else=B(i,)；q=B(i,)；endSi(,q)=E()*(conj(E())*conj(B(i,))+(conj(E()*B(i,))conj(E(q)))*conj((B(i,)*B(i,))))；dis(Si(,q))；enddis(‘各条支路的末端功率Sj(顺序同您输入B时的一样)：‘)为；fori=:nlifB(i,)===B(i,)；q。

2、ndendk=；K=K+；elsekq=；ifk~=K=K+；endendfori=:nDy(K,i)=V(i)；endenddis(‘迭代次数‘)；dis(K)；dis(‘每次没有达到精度要求的有功功率个数为‘)；dis(NP)；dis(‘每次没有达到精度要求的无功功率个数为‘)；dis(Nq)；fork=:nE(k)=V(k)*cos(O(K))+V(k)*sin(O(k))*j；O(k)=O(k)*i；enddis(‘各节点的电压标么值E为(节点号从小到大排)：‘)；dis(E)；dis(‘各节点的电压U大小(节点号从小到大排)：‘)为；dis(U)；dis(‘各节点的电压相角O(节点号从小到大排)：‘)为；dis(O)；for=:nC()=；for。

3、的PQ分解法潮流计算程序MATLAB已广泛应用于自动控制、数学运算、信号分析、计算机技术、图像信号处理、财务分析、航天工业和生物医学工程等领域。由于MATLAB语言功能强大、人际界面友好、编程效率高、强大而智能化的作图功能，且具有编程语句简洁、灵活、表达和运算能力强等显著特点。程序清单%本程序的功能是用PQ分解法进行潮流计算n=inut(‘请输入节点数：n=’)；nl==inut(‘请输入支路数：nl==’)；isb==inut(‘请输入平衡母线节点号：isb=’)；r==inut(‘请输入误差精度:r=’)；B==inut(‘请输入由支路参数形成的矩阵:B=’)；B==inut(‘请输入由节点参数形成的矩阵:B=’)；X==inut(‘请输入由节点号及其。

4、对地阻抗形成的矩阵:X=’)；na==inut(‘请输入PQ节点号:na=’)；Y==zeros(n)；YI=zeros(n)；e=zeros(,n)；f=zeros(,n)；V==zeros(,n)；O=zeros(,n)；fori=:nifX(i,)~=；=x(i,)；Y(,)=X(i,)；endendfori=:nlifB(i,)===B(i,)；q=B(i,)；else=B(i,)；q=B(i,)；endY(,q)=Y(,q)(B(i,)*B(i,))；YI(,q)=YI(,q)B(i,)；Y(q,)=Y(,q)；YI(q,)=YI(,q)；Y(q,q)=Y(q,q)+(B(i,)*B(i,)~)+B(i,)；YI(q,q)=YI(q,q)+B(i。

5、，得至号节点功率分别为：+iiii+i举例验算结论：通过手动运算，发现运算结果与手动结果运算结果一样，故软件计算正确。总结通过学习潮流计算原理参阅一些参考资料及文献，对传统潮流死算有所理解，对一些新的潮流计算技术有了整体的认识。潮流计算是电力系统的一项基本计算，它根据给定的运行条件及系统接线情况确定系统的运行状态、母线的电压，各元件中流过的功率，系统的功率损耗等，进而比较运行方式的合理性、可靠性和经济性。本文分析和比较了电力系统潮流计算的各种方法，其中PQ分解法具有计算速度快、所占内存少等优势被广泛采用。本设计首先进行了PQ分解法的理论研究，然后应用MATLAB软件编制PQ分解法潮流程序，并通过一个电力系统实例对所编程序进行验证，得出结论。=:nC(i)=。

6、输入框图InutBlockDiagram在输入需要的数据后,点击运算,出现结果如下：结果：迭代次数每次没有达到精度要求的有功功率个数为每次没有达到精度要求的无功功率个数为各节点的电压标么值E(节点号从小到大排)为：i+ii+i各节点的电压U大小(节点号从小到大排)为：各节点的电压相角O(节点号从小到大排)为：各节点的功率S(节点号从小到大排)为：+iiii+i各条支路的首端功率Si(顺序同您输入B时的一样)为：+i+i+i+i+i各条支路的末端功率Sj(顺序同您输入B时的一样)为：iiiii各条支路的功率损耗DS(顺序同您输入B时的一样)为：+ii+i+i+i点击计算后软件运行情况如图所示：图仿真图FigureSimulationdiagram经过手动计算。

7、L)*DP(L)；endendendendfori=:nO(i)=O(i)DP(i)；endkq=；L=；fori=:nifB(i,)==C(i)=；L=L+；fork=:nC(i)=C(i)+V(k)*(G(i,k)*sin(O(i)O(k))+BI(i,k)*cos(O(i)O(k)))；endDQ(i)=Q(i)V(i)*C(i)；DQ(L)=DQ(i)V(i)；DET=abs(DQ(i))；ifDETgt=rICT=ICT+；endendendendelsek=；ifkq~=；L=；fori=:nifB(i,)==C(i)=；L=L+；for+i；+i；+i；]请输入由节点号及其对地阻抗形成的矩阵：X=[；；；；]请输入PQ节点数na=输入如图所示。

8、,i)；endB(i,i)=B(i,i)；fork=i+:nforj=i+:nB(k,j)=B(k,j)B(k,i)*B(i,j)；endendendend=；q=；fori=:nifB(i,)===+；k=；forj=:nifB(i,)==k=k+；A(,k)=BI(i,j)；endendendendfori=:naifi==naA(i,i)=A(i,i)；elsek=i+；forj=k:naA(i,j)=A(i,j)A(i,i)；endA(i,i)=A(i,i)；fork=i+:naforj=i+:naA(k,j)=A(k,j)A(k,i)*A(i,j)；endendendendICT=；ICT=；k=,kq=；K=；DET=；ICT=；whileIC。

9、(i)+V(k)*(G(i,k)*sin(O(i)O(k))BI(i,k)*cos(O(i)O(k)))；endDQ(i)=Q(i)V(i)*C(i)；DQ(L)=DQ(i)V(i)；DET=abs(DQ(i))；endendendendNq(K)=ICT；ifICT~=L=；fori=:naDQ(i)=A(i,i)*DQ(i)；ifi==naforLZ=:iL=i+LZ；IC=L；forMZ=:ICI=IC+MZ；DQ(I)=DQ(I)A(I,L)*DQ(L)；endendelseIC=i+；fork=IC:naDQ(k)=DQ(k)A(k,i)*DQ(i)；endendendL=；fori=:nifB(i,)==L=L+；V(i)=V(i)DQ(L)；。

10、B(i,)；else=B(i,)；q=B(i,)；endSj(q,)=E(q)*(conj(E(q))*conj(B(i,))+(conj(E(q)*B(i,))conj(E()))*conj((B(i,)*B(i,))))；dis(Sj(q,))；enddis(‘各条支路的功率损耗DS(顺序同您输入B时的一样)：‘)为；fori=:nlifB(i,)===B(i,)；q=B(i,)；else=B(i,)；q=B(i,)；endDS(i)=Si(,q)+Sj(q,)；dis(DS(i))；endfori=:KCs(i)=i；forj=:nDy(K,j)=Dy(K,j)；endenddis(‘以下是每次迭代后各节点的电压值(如图所示)‘)；lot(Cs,Dy。

11、~=|ICT~=ICT=；ICT=；fori=:nifi~=isbC(i)=；fork=:nC(i)=C(i)+V(k)*(G(i,k)*cos(O(i)O(k))+BI(i,k)*sin(O(i)O(k)))；endDP(i)=P(i)V(i)*C(i)；DP(i)=DP(i)V(i)；DET=abs(DP(i))；ifDETgt=rICT=ICT+；endendendN(K)=ICT；ifICT~=fori=:nDP(i)=B(i,i)*DP(i)；ifi~=nIC=i+；fork=ic:nDP(k)=DP(k)B(k,i)*DP(i)；endelseforLZ=:iL=i+LZ；IC=L；forMZ=:ICI=IC+MZ；DP(I)=DP(I)B(I 。

12、)；Y(,)=Y(,)+B(i,)+B(i,)；YI(,)=YI(,)+B(i,)；end%求导纳矩阵G=real(Y)；B=imag(YI)；BI=imag(Y)；fori=:nS(i)=B(i,)B(i,)；BI(i,i)=BI(i,i)+B(i,)；endP=real(S)；Q=imag(S)；fori=:ne(i)=real(B(i,))；f(i)=imag(B(i,))；V(i)=B(i,)；endfori=:nifB(i,)==V(i)=sqrt(e(i)~+f(i)~)；O(i)=atan(f(i)e(i))；endendfori=:nifi==nB(i,i)=B(i,i)；elseIC=i+；forj=IC:nB(i,j)=B(i,j)B(