《基于MATLAB电力系统PQ时尚计算程序设计.doc》由会员分享，可免费在线阅读全文，更多与《基于MATLAB电力系统PQ时尚计算程序设计(最终版)》相关文档资源请在帮帮文库(www.woc88.com)搜索数亿文档库存。
1、=B(i,j)B(i,i)；wwwwoccomendB(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)；wwwwoccomendendendendICT=；ICT=；k=,kq=；K=；DET=；ICT=；whileICT~=|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。
2、i； i； i；]请输入节点号及其对地阻抗形成的矩阵：X=[；；；]wwwwoccom请输入PQ节点数na=如图所示，输入框图InutBlockDiagram输入所需数据后，单击操作，结果如下:结果:每次迭代次数不符合精度要求的有功功率数为每个节点的电压标准值E(从小到大排节点号)为:i ii iwwwwoccom各节点的电压U大小(从小到大排节点号)为:各节点的电压相角O(从小到大排节点号)为:各节点的功率S(从小到大排节点号)为: iiii i各支路的一端功率Si(顺序与您输入B时相同)为: i i i i i各支路的终端功率Sj(顺序与您输入B时相同)为:iiiii各支路的功率损失DS(顺序与您输入B时相同)为: ii i i i如图所示，点击计算后的软件操作：wwwwoccom图仿真图FigureSimulationdiagram手动计算后，得到号节点的功率分别为： iiii i举例验算结论：通过手动运算，发现运算结果与手动结果运算结果一样，故软件计算正确。wwwwoccom通过学习趋势计算原始总结。
3.电厂（或发电机）有时可以根据其他原则进行选择。例如，为了提高计算的收敛性，可以选择大量出线或靠近电网中心的电厂母线作为平衡节点。以上三类节点个运行参数P、Q、U、θ已知量是两个，待求量也是两个，只是类型不同。在趋势计算中，任何复杂的电力系统都可以归因于以下部件（参数）：发电机（注入电流或功率）、负载（注入负电流或功率）、输电线路支路（电阻、电抗）、变压器支路（电阻、电抗、变比）、母线对地支路（电阻、导电）、线路对地支路（一般为线路充电容器导电）。简单的网络集中了上述各种组件，如图所示(a)所示。必须指出的是，如果只研究稳态下的趋势而不涉及临时过程的计算，则不需要发电机和负载的阻抗参数，只需给发电机和负载注入功率和电流，wwwwoccom并规定发电机和负载的注入功率和电流正确，负载正确。(a)电网结构图用于潮流计算(a)Theflowcalculationgridchart无源线性网络(可用导纳矩阵或阻抗矩阵表示)IIIIggPjQ?PjQ?RRPjQ??PjQ??RRPjQj?RRR(b)等值网络的趋势计算(b)thetrendtoc。
4、ifB(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),xlable(‘迭代次数’)，ylable(电压)，title(‘电压迭代次数曲线’)；设计验证现在以一个例子来说明软件的使用方法:如所示，电力系统接线()n为节点数、nl为支路数、isb平衡母线节点号(固定为)r误差精度。()输入由支路参数形成的矩阵B矩阵B每行由以下参数组成:①支路的第一端号P；②末端号Q，且PltQ；③支路的阻抗(R jX)；④支路对地容抗；⑤支路的变比K；wwwwoccom⑥将标志转换为哪一侧(如果支路首端P处于高压侧，请输入，否则请输入)。()输入由支路参数形成的矩阵B矩阵B每行由以下参数组成:①支路的第一端号P；②末端号Q，且PltQ；③支路的阻抗(R jX)；④支路对地容抗；⑤支路的变比K；wwwwoccom⑥将标志转换为哪一侧(如果支路首端P处于高压侧，请输入，否则请输入)。①②③④⑤P?j:U?j?jjjjj?j?j???j:j?j?U图电路图Figureschematic()请输入每个节点参数形成的矩阵B矩阵B。
5、(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=:iwwwwoccomL=i LZ；IC=L；forMZ=:ICI=IC MZ；DP(I)=DP(I)B(I,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=:nwwwwoccomifB(i,)==C(i)==[；。
6.参考一些参考资料和文献，了解传统趋势的死算，全面了解一些新的趋势计算技术。趋势计算是电力系统的基本计算。根据给定的运行条件和系统接线条件，确定系统的运行状态、母线的电压、各部件流动的功率、系统的功率损耗等，比较运行模式的合理性、可靠性和经济性。本文对电力系统的趋势计算方法进行了分析和比较PQ分解法具有计算速度快、内存少等优点。本设计首先进行PQ分解法的理论研究，然后应用MATLAB软件编制PQ通过趋势程序，并通过电力系统实例验证编程程序，得出结论。wwwwoccom；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))；endendendendNq(K)=ICT；ifICT~=L=；fori=:naDQ(i)=A(i,i)*DQ(i)；ifi==naforLZ=:iL=i LZ；IC=L；forMZ=:ICI=IC MZ；。
7、Q(I)=DQ(I)A(I,L)*DQ(L)；endendelseIC=i ；fork=IC:naDQ(k)=DQ(k)A(k,i)*DQ(i)；endendwwwwoccomendL=；fori=:nifB(i,)==L=L ；V(i)=V(i)DQ(L)；endendk=；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)；wwwwoccomdis(每个节点的电压U大小(节点号从小到大):‘)为；dis(U)；dis(每个节点的电压相角O(从小到大排节点号):‘)为；dis(O)；for=:nC()=；forq=:nC()=C() conj(Y(,。
8、))*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=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,))))；wwwwoccomdis(Sj(q,))；enddis(各支路的功率损失DS(顺序与输入B时相同):‘)为；fori=:nl。
9、B(i,)；q=B(i,)；else=B(i,)；q=B(i,)；endwwwwoccomY(,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,)；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。
10、MATLAB广泛应用于自动控制、数学运算、信号分析、计算机技术、图像信号处理、财务分析、航天工业和生物医学工程等领域。由于MATLAB语言功能强大，人际界面友好，编程效率高，作图功能强大智能，编程句简洁灵活，表达和操作能力强。本程序的功能是程序清单%PQ趋势计算采用分解法n=inu(‘请输入节点数：n=’)；nl==inut(‘请输入支路数：nl==’)；isb==inut(‘请输入平衡母线节点号：isb=’)；r==inut(‘请输入误差精度:r=’)；B==inut(‘请输入由支路参数形成的矩阵:B=’)；B==inut(‘请输入由节点参数形成的矩阵:B=’)；X==inut(‘请输入由节点号及其对地阻抗形成的矩阵: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,)===。
11、是由下列参数构成的：①节点所接发电机的功率SG；②节点负荷的功率SL；③节点电压的初始值；④PU节点电压U的给定值；⑤节点所接的无功补偿设备的容量；⑥节点分类标号igl。输入数据：请输入节点数：n=请输入支路数：nl=请输入平衡母线节点号：isb=请输入误差精度：r=请输入由支路参数形成的矩阵：B=[i；+ii；+i；+ii；i]请输入由节点参数形成的矩阵：及年由岩本伸一等提出的保留非线性的高速潮流计算法，其中快速分解法从年开始已在国内使用，并习惯称之为PQ分解法。由于PQ分解法在计算速度上大大超过了NR法，不但能应用于离线潮流计算，而且也能应用于在线潮流计算。潮流计算的基本方程节点的分类用一般的电路理论求解网络方程，目的是给出电压源(或电流源)研究网络内的电流(或电压)分布，作为基础的方程式，一般用线性代数方程式表示。然而在电力系统中，给出发电机或负荷连接母线上电压或电流(都是向量)的情况是很少的，一般是给出发电机母线上发电机的有功功率(P)和母线电压的幅值(U)，给出负荷母线上负荷消耗的有功功率(P)和无功功率(Q)。主要目的是由这些已知量去求电力系统内的各种电气量 。
12、。所以，根据电力系统中各节点性质的不同，很自然的把节点分成三种类型。PQ节点对这一类节点，事先给定的是节点功率(P、Q)，待求的未知量是节点电wwwwoccom压向量(U、θ)，所以叫“PQ节点”。通常变电所母线都是PQ节点，当某些发电机的输出功率P、Q给定时，也作为PQ节点。PQ节点上的发电机称之为PQ机(或PQ给定型发电机)。在潮流计算中，系统大部分节点属于PQ节点。PU节点这类节点给出的参数是该节点的有功功率P及电压幅值U，待求量为该节点的无功功率Q及电压向量的相角θ。这类节点在运行中往往要有一定可调节的无功电源，用于维持给定的电压值。通常选择有一定无功功率储备的发电机母线或者变电所有无功补偿设备的母线作PU节点处理。PU节点上的发电机称之为PU机。平衡节点在潮流计算中，这类节点一般只设一个。对该节点，给定其电压值，并在计算中取该节点电压向量的方向作为参考轴，相当于给定该点电压向量的角度为零。也就是说，对平衡节点给定的运行参数是U和θ，因此又称为Uθ节点，而待求量是该节点的P，Q，整个系统的功率平衡由这一节点承担。关于平衡节点的选择，一般选择系统中担任调频调压的某