潮流计算的计算机算法

《潮流计算的计算机算法》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、高等电力系统分析（潮流计算的计算机算法）PQ分解法潮流计算（IEEE14)目录一、MATLAB源程序二、对支路参数（B1）、节点参数（B2）的说明三、带入数据，运行结果一、MATLAB源程序clearclosealln=input('请输入节点数:n=');n1=input('请输入支路数:n1=');isb=input('请输入平衡节点号:isb=');pr=input('请输入误差精度:pr=');B1=input('请输入支路参数:B1=');B2=input('请输入节点参数:B2=');n2=input('请输入PQ节点个数:n2=');Y=zeros(n);fori=

2、1:n1p=B1(i,1);q=B1(i,2);Y(p,q)=Y(p,q)-1/(B1(i,3)+B1(i,4)*1j);%非对角元Y(q,p)=Y(p,q);Y(p,p)=Y(p,p)+1/(B1(i,3)+B1(i,4)*1j)+B1(i,6)*1j;%对角元Y(q,q)=Y(q,q)+1/(B1(i,3)+B1(i,4)*1j)+B1(i,6)*1j;enddisp('导纳矩阵Y=');disp(Y)%---------------------------------------------%---------------下面是求P,Q,V,O矩阵-------------

3、--V=zeros(1,n);O=zeros(1,n);P=zeros(1,n);Q=zeros(1,n);G=real(Y);B=imag(Y);fori=1:nP(i)=B2(i,3);Q(i)=B2(i,4);V(i)=B2(i,5);O(i)=B2(i,6);endB3=B(1:n-1,1:n-1);%不含平衡节点，由节点导纳虚部构成B4=B(1:n2,1:n2);%所有PQ节点%----------------------------------------------%---------------下面是求ΔP,ΔQ矩阵---------------DX=0;ICT=

4、1;Mp=1;Mq=1;whileICT~=0m1=1;m2=1;fori=1:nifi~=isbC(i)=0;D(i)=0;forj1=1:nC(i)=C(i)+V(i)*V(j1)*(G(i,j1)*cos(O(i)-O(j1))+B(i,j1)*sin(O(i)-O(j1)));D(i)=D(i)+V(i)*V(j1)*(G(i,j1)*sin(O(i)-O(j1))-B(i,j1)*cos(O(i)-O(j1)));endDP(m1)=P(i)-C(i);m1=m1+1;ifB2(i,2)==1DQ(m2)=Q(i)-D(i);m2=m2+1;endendendm1=m1

5、-1;%所有节点数m2=m2-1;%PQ节点数DPQ=[DP';DQ'];%求DP,DQV1=V(:,1:m1);V2=diag(V1);V3=inv(V2);%对V矩阵求逆H=V3*DP';%ΔP/VK=-inv(B3)*H;%-ΔP/V/B3deltO=V3*K;%Δ角=-ΔP/V/V/B3max1=max(abs(DP));fori=1:m1ifmax1

6、=-ΔQ/V/Bmax2=max(abs(DQ));fori=1:m2ifmax2

7、迭代次数为：');disp(DX);fori=1:nE(i)=V(i)*cos(O(i))+1j*V(i)*sin(O(i));o(i)=180*angle(E(i))/pi;enddisp('---------------------------------------------');disp('修正后各节点电压标么值为（节点号从小到大排列）:');disp(V);disp('---------------------------------------------');di