logo资料库

PQ分解法的程序.pdf

发布时间：2022-06-08 发布人：admin 分类：说明书 资料大小：0.06M 资料格式：pdf 举报 版权申诉
第1页 / 共8页
第2页 / 共8页
第3页 / 共8页
第4页 / 共8页
第5页 / 共8页
第6页 / 共8页
第7页 / 共8页
第8页 / 共8页
资料共8页，全文预览结束
    n=input('请输入节点数:n=');n1=input('请输入支路数:n1=');isb=input('请输入平衡母线节点号:isb=');pr=input('请输入误差精度:pr=');B1=input('请输入由各支路参数形成的矩阵:B1=');B2=input('请输入由各节点参数形成的矩阵:B2=');na=input('请输入PQ节点数na=');Y=zeros(n);YI=zeros(n);e=zeros(1,n);f=zeros(1,n);V=zeros(1,n);o=zeros(1,n);fori=1:n1ifB1(i,6)==0p=B1(i,1);q=B1(i,2);elsep=B1(i,2);q=B1(i,1);endY(p,q)=Y(p,q)-1./(B1(i,3)*B1(i,5));YI(p,q)=YI(p,q)-1./B1(i,3);Y(q,p)=Y(p,q);YI(q,p)=YI(p,q);Y(q,q)=Y(q,q)+1./(B1(i,3)*B1(i,5)^2)+B1(i,4)./2;YI(q,q)=YI(q,q)+1./B1(i,3);Y(p,p)=Y(p,p)+1./B1(i,3)+B1(i,4)./2;YI(p,p)=YI(p,p)+1./B1(i,3);end%求导纳矩阵disp('导纳矩阵Y=');disp(Y);G=real(Y);B=imag(YI);BI=imag(Y);fori=1:nS(i)=B2(i,1)-B2(i,2);BI(i,i)=BI(i,i)+B2(i,5);endP=real(S);Q=imag(S);fori=1:ne(i)=real(B2(i,3));f(i)=imag(B2(i,3));V(i)=B2(i,4);endfori=1:nifB2(i,6)==2V(i)=sqrt(e(i)^2+f(i)^2);o(i)=atan(f(i)./e(i));endendfori=2:n
    ifi==nB(i,i)=1./B(i,i);elseIC1=i+1;forj1=IC1:nB(i,j1)=B(i,j1)./B(i,i);endB(i,i)=1./B(i,i);fork=i+1:nforj1=i+1:nB(k,j1)=B(k,j1)-B(k,i)*B(i,j1);endendendendp=0;q=0;fori=1:nifB2(i,6)==2p=p+1;k=0;forj1=1:nifB2(j1,6)==2k=k+1;A(p,k)=BI(i,j1);endendendendfori=1:naifi==naA(i,i)=1./A(i,i);elsek=i+1;forj1=k:naA(i,j1)=A(i,j1)./A(i,i);endA(i,i)=1./A(i,i);fork=i+1:naforj1=i+1:naA(k,j1)=A(k,j1)-A(k,i)*A(i,j1);endendendendICT2=1;ICT1=0;kp=1;kq=1;K=1;DET=0;ICT3=1;whileICT2~=0|ICT3~=0ICT2=0;ICT3=0;
    fori=1:nifi~=isbC(i)=0;fork=1:nC(i)=C(i)+V(k)*(G(i,k)*cos(o(i)-o(k))+BI(i,k)*sin(o(i)-o(k)));endDP1(i)=P(i)-V(i)*C(i);DP(i)=DP1(i)./V(i);DET=abs(DP1(i));ifDET>=prICT2=ICT2+1;endendendNp(K)=ICT2;ifICT2~=0fori=2:nDP(i)=B(i,i)*DP(i);ifi~=nIC1=i+1;fork=IC1:nDP(k)=DP(k)-B(k,i)*DP(i);endelseforLZ=3:iL=i+3-LZ;IC4=L-1;forMZ=2:IC4I=IC4+2-MZ;DP(I)=DP(I)-B(I,L)*DP(L);endendendendfori=2:no(i)=o(i)-DP(i);endkq=1;L=0;fori=i:nifB2(i,6)==2C(i)=0;L=L+1;fork=1:nC(i)=C(i)+V(k)*(G(i,k)*sin(o(i)-o(k))-BI(i,k)*cos(o(i)-o(k)));end
    DQ1(i)=Q(i)-V(i)*C(i);DQ(L)=DQ1(i)./V(i);DET=abs(DQ1(i));ifDET>=prICT3=ICT3+1;endendendelsekp=0;ifkp~=0L=0;fori=1:nifB2(i,6)==2C(i)=0;L=L+1;fork=1:nC(i)=C(i)+V(k)*(G(i,k)*sin(o(i)-o(k))-BI(i,k)*cos(o(i)-o(k)));endDQ1(i)=Q(i)-V(i)*C(i);DQ(L)=DQ1(i)./V(i);DET=abs(DQ1(i));endendendendNq(K)=ICT3;ifICT3~=0L=0;fori=1:naDQ(i)=A(i,i)*DQ(i);ifi==naforLZ=2:iL=i+2-LZ;IC4=L-1;forMZ=1:IC4I=IC4+1-MZ;DQ(I)=DQ(I)-A(I,L)*DQ(L);endendelseIC1=i+1;fork=IC1:naDQ(k)=DQ(k)-A(k,i)*DQ(i);endend
    endL=0;fori=1:nifB2(i,6)==2L=L+1;V(i)=V(i)-DQ(L);endendkp=1;K=K+1;elsekq=0;ifkp~=0;K=K+1;endendfori=1:nDy(K-1,i)=V(i);endenddisp('迭代次数');disp(K);disp('没有达到精度要求的有功功率个数');disp(Np);disp('没有达到精度要求的无功功率个数');disp(Nq);fork=1:nE(k)=V(k)*cos(o(k))+V(k)*sin(o(k))*j;o(k)=o(k)*180/pi;enddisp('各节点的实际电压标么值E为(节点号从小到大排列):');disp(E);disp('各节点的电压大小V为(节点号从小到大排列):');disp(V);disp('各节点的电压相角时o为(节点号从小到大排列):');disp(o);forp=1:nC(p)=0;forq=1:nC(p)=C(p)+conj(Y(p,q))*conj(E(q));endS(p)=E(p)*C(p);enddisp('各节点的功率S为(节点号从小到大排列):');
    disp(S);disp('各条支路的首端功率Si为(顺序同您输入B1时一样):');fori=1:n1ifB1(i,6)==0p=B1(i,1);q=B1(i,2);elsep=B1(i,2);q=B1(i,1);endSi(p,q)=E(p)*(conj(E(p))*conj(B1(i,4)./2)+(conj(E(p)*B1(i,5))-conj(E(q)))*conj(1./(B1(i,3)*B1(i,5))));disp(Si(p,q));enddisp('各条支路的末端功率Sj为(顺序同您输入B1时一样):');fori=1:n1ifB1(i,6)==0p=B1(i,1);q=B1(i,2);elsep=B1(i,2);q=B1(i,1);endSj(q,p)=E(q)*(conj(E(q))*conj(B1(i,4)./2)+(conj(E(q)./B1(i,5))-conj(E(p)))*conj(1./(B1(i,3)*B1(i,5))));disp(Sj(q,p));enddisp('各条支路的功率损耗DS为(顺序同您输入B1时一样):');fori=1:n1ifB1(i,6)==0p=B1(i,1);q=B1(i,2);elsep=B1(i,2);q=B1(i,1);endDS(i)=Si(p,q)+Sj(q,p);disp(DS(i));endfori=1:KCs(i)=i;forj=1:nDy(K,j)=Dy(K-1,j);endend请输入节点数:n=5请输入支路数:n1=5请输入平衡母线节点号:isb=1请输入误差精度:pr=0.00001请输入由各支路参数形成的矩阵:B1=[120.03i01.050;230.08+0.3i0.5i10;240.1+0.35i01
    0;340.04+0.25i0.5i10;350.015i01.051]请输入由各节点参数形成的矩阵:B2=[001.051.0501;03.7+1.3i1.05002;02+1i1.05002;01.6+0.8i1.05002;501.051.0503]请输入PQ节点数na=3导纳矩阵Y=0-33.3333i0+31.7460i0000+31.7460i1.5846-35.7379i-0.8299+3.1120i-0.7547+2.6415i00-0.8299+3.1120i1.4539-66.9808i-0.6240+3.9002i0+63.4921i0-0.7547+2.6415i-0.6240+3.9002i1.3787-6.2917i0000+63.4921i00-66.6667i迭代次数6没有达到精度要求的有功功率个数444330没有达到精度要求的无功功率个数000000各节点的实际电压标么值E为(节点号从小到大排列):1.05001.0473-0.0759i1.0017+0.3149i1.0465-0.0862i0.9766+0.3856i各节点的电压大小V为(节点号从小到大排列):1.05001.05001.05001.05001.0500各节点的电压相角时o为(节点号从小到大排列):0-4.145617.4510-4.711721.5470各节点的功率S为(节点号从小到大排列):2.5302+1.8416i-3.7000-1.2812i-2.0000-3.7435i-1.6000+0.3100i5.0000+3.6788i各条支路的首端功率Si为(顺序同您输入B1时一样):2.5302+1.8416i-1.1986+0.3020i0.0288-0.0081i
    1.6729-0.2175i5.0000+3.6788i各条支路的末端功率Sj为(顺序同您输入B1时一样):-2.5302-1.5751i1.3271-0.3715i-0.0287+0.0084i-1.5713+0.3016i-5.0000-3.1545i各条支路的功率损耗DS为(顺序同您输入B1时一样):0+0.2665i0.1285-0.0695i8.1234e-005+2.8432e-004i0.1017+0.0841i0.0000+0.5243i
    分享到：
    赞 收藏

    相关推荐

    资料库

    课程资源

    共收录17145份资料，累计13个分类，关注成员有19位，主要包括：PHP,网络管理,网页制作,Java,.Net,数据库,3G/移动开发,C/C++,游戏开发,嵌入式,讲义,软件测试,专业指导

    热门标签

    PHP 网络管理 网页制作 Java .Net 数据库 3G/移动开发 C/C++ 游戏开发 嵌入式 讲义 软件测试 专业指导

    最新资料