题号:72编程计算潮流(电力系统分析)节点数:3 支路数:3 计算精度:0.00010支路1:0.0300+j0.09001┠—————□—————┨2支路2:0.0200+j0.09002┠—————□—————┨3支路3:0.0300+

来源:学生作业帮助网 编辑:六六作业网 时间:2024/12/23 21:55:33
题号:72编程计算潮流(电力系统分析)节点数:3支路数:3计算精度:0.00010支路1:0.0300+j0.09001┠—————□—————┨2支路2:0.0200+j0.09002┠—————□

题号:72编程计算潮流(电力系统分析)节点数:3 支路数:3 计算精度:0.00010支路1:0.0300+j0.09001┠—————□—————┨2支路2:0.0200+j0.09002┠—————□—————┨3支路3:0.0300+
题号:72编程计算潮流(电力系统分析)
节点数:3 支路数:3 计算精度:0.00010
支路1:0.0300+j0.0900
1┠—————□—————┨2
支路2:0.0200+j0.0900
2┠—————□—————┨3
支路3:0.0300+j0.0900
3┠—————□—————┨1
节点1:PQ节点,S(1)=-0.5000-j0.2000
节点2:PQ节点,S(2)=-0.6000-j0.2500
节点3:平衡节点,U(3)=1.0000∠0.0000

题号:72编程计算潮流(电力系统分析)节点数:3 支路数:3 计算精度:0.00010支路1:0.0300+j0.09001┠—————□—————┨2支路2:0.0200+j0.09002┠—————□—————┨3支路3:0.0300+

n=5;

nl=5;

isb=1;

pr=0.00001; 

B1=[1 2 0.03i 0 1.05 0;2 3 0.08+0.3i 0.5i 1 0;2 4 0.1+0.35i 0 1 0;3 4 0.04+0.25i 0.5i 1 0;3 5 0.015i 0 1.05 1];

B2=[0 0 1.05 1.05 0 1;0 3.7+1.3i 1.05 0 0 2;0 2+1i 1.05 0 0 2;0 1.6+0.8i 1.05 0 0 2;5 0 1.05 1.05 0 3];

X=[1 0;2 0;3 0;4 0;5 0];

na=3;

Y=zeros(n);YI=zeros(n);e=zeros(1,n);f=zeros(1,n);V=zeros(1,n);O=zeros(1,n);

for i=1:n

    if X(i,2)~=0;

        p=X(i,1);

        Y(p,p)=1./X(i,2);

    end

end

for i=1:nl

    if B1(i,6)==0

        p=B1(i,1);q=B1(i,2);

    else p=B1(i,2);q=B1(i,1);

    end

    Y(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

G=real(Y);B=imag(YI);BI=imag(Y);

for i=1:n

    S(i)=B2(i,1)-B2(i,2);

    BI(i,i)=BI(i,i)+B2(i,5);

end

P=real(S);Q=imag(S);

for i=1:n

    e(i)=real(B2(i,3));

    f(i)=imag(B2(i,3));

    V(i)=B2(i,4);

end

for i=1:n

    if B2(i,6)==2

        V(i)=sqrt(e(i)^2+f(i)^2);

        O(i)=atan(f(i)./e(i));

    end

end

for i=2:n

    if i==n

        B(i,i)=1./B(i,i);

    else IC1=i+1;

        for j1=IC1:n

            B(i,j1)=B(i,j1)./B(i,i);

        end

        B(i,i)=1./B(i,i);

        for k=i+1:n

            for j1=i+1:n

                B(k,j1)=B(k,j1)-B(k,i)*B(i,j1);

            end

        end

    end

end

p=0;q=0;

for i=1:n

    if B2(i,6)==2

      p=p+1;k=0;

        for j1=1:n

            if B2(j1,6)==2

                k=k+1;

                A(p,k)=BI(i,j1);

            end

        end

    end

end

for i=1:na

    if i==na

        A(i,i)=1./A(i,i);

    else k=i+1;

        for j1=k:na

            A(i,j1)=A(i,j1)./A(i,i);

        end

        A(i,i)=1./A(i,i);

        for k=i+1:na

            for j1=i+1:na

                A(k,j1)=A(k,j1)-A(k,i)*A(i,j1);

            end

        end

    end

end

ICT2=1;ICT1=0;kp=1;kq=1;K=1;DET=0;ICT3=1;

while ICT2~=0|ICT3~=0

ICT2=0;ICT3=0;

for i=1:n

    if i~=isb

        C(i)=0;

        for k=1:n

            C(i)=C(i)+V(k)*(G(i,k)*cos(O(i)-O(k))+BI(i,k)*sin(O(i)-O(k)));

        end

        DP1(i)=P(i)-V(i)*C(i);

        DP(i)=DP1(i)./V(i);

        DET=abs(DP1(i));

        if DET>=pr

            ICT2=ICT2+1;

        end

    end

end

Np(K)=ICT2;

   if ICT2~=0

   for i=2:n

       DP(i)=B(i,i)*DP(i);

       if i~=n

           IC1=i+1;

           for k=IC1:n

               DP(k)=DP(k)-B(k,i)*DP(i);

           end

       else

           for LZ=3:i

               L=i+3-LZ;

               IC4=L-1;

               for MZ=2:IC4

                   I=IC4+2-MZ;

                   DP(I)=DP(I)-B(I,L)*DP(L);

               end

           end

       end

   end

       for i=2:n

           O(i)=O(i)-DP(i);

       end

           kq=1;L=0;

           for i=1:n

               if B2(i,6)==2

                   C(i)=0;L=L+1;

                   for k=1:n

                       C(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));

                   if DET>=pr

                       ICT3=ICT3+1;

                   end

               end

           end

   else kp=0;

       if kq~=0;

           L=0;

           for i=1:n

               if B2(i,6)==2

                   C(i)=0;L=L+1;

                   for k=1:n

                       C(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));

               end

           end

       end

   end

   Nq(K)=ICT3;

   if ICT3~=0

       L=0;

       for i=1:na

           DQ(i)=A(i,i)*DQ(i);

           if i==na

               for LZ=2:i

                   L=i+2-LZ;

                   IC4=L-1;

                   for MZ=1:IC4

                       I=IC4+1-MZ;

                       DQ(I)=DQ(I)-A(I,L)*DQ(L);

                   end

               end

           else

               IC1=i+1;

               for k=IC1:na

                   DQ(k)=DQ(k)-A(k,i)*DQ(i);

               end

           end

       end

       L=0;

       for i=1:n

           if B2(i,6)==2

               L=L+1;

               V(i)=V(i)-DQ(L);

           end

       end

       kp=1;

       K=K+1;

   else

       kq=0;

       if kp~=0

           K=K+1;

       end

   end

   for i=1:n

       Dy(K-1,i)=V(i);

   end

end

disp('迭代次数')

disp(K);

disp('每次没有达到精度要求的有功功率个数为');

disp(Np);

disp('每次没有达到精度要求的无功功率个数为');

disp(Nq);

for k=1:n

    E(k)=V(k)*cos(O(k))+V(k)*sin(O(k))*j;

    O(k)=O(k)*180./pi;

end

disp('各节点的电压标么值E为');

disp(E);

disp('各节点的电压V大小');

disp(V);

disp('各节点的电压相角O');

disp(O);

for p=1:n

    C(p)=0;

    for q=1:n

        C(p)=C(p)+conj(Y(p,q))*conj(E(q));

    end

    S(p)=E(p)*C(p);

end

disp('各节点的功率为');

disp(S);

disp('各条支路的首端功率为');

for i=1:nl

    if B1(i,6)==0

        p=B1(i,1);q=B1(i,2);

    else p=B1(i,2);q=B1(i,1);

    end

    Si(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));

end

disp('各条支路的末端功率为');

for i=1:nl

    if B1(i,6)==0

        p=B1(i,1);q=B1(i,2);

    else p=B1(i,2);q=B1(i,1);

    end

    Sj(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));

end

disp('各条支路的功率损耗为');

for i=1:nl

    if B1(i,6)==0

        p=B1(i,1);q=B1(i,2);

    else p=B1(i,2);q=B1(i,1);

    end

    DS(i)=Si(p,q)+Sj(q,p);

    disp(DS(i));

end

for i=1:K

    Cs(i)=i;

    for j=1:n

        Dy(K,j)=Dy(K-1,j);

    end

end

disp('每次迭代后各节点的电压值如图所示');

plot(Cs,Dy)

xlabel('迭代次数')

ylabel('电压')

title('电压迭代次数曲线');

题号:72编程计算潮流(电力系统分析)节点数:3 支路数:3 计算精度:0.00010支路1:0.0300+j0.09001┠—————□—————┨2支路2:0.0200+j0.09002┠—————□—————┨3支路3:0.0300+ 题号:2编程计算潮流(电力系统分析)题号:2节点数:4 支路数:4 计算精度:0.00010支路 1:0.0200+j0.08001┠—————□—————┨3支路 2:0.0400+j0.12001┠—————□—————┨4支路 急求:电力系统分析计算 电力系统分析 电力系统分析 牛顿拉夫逊法潮流计算公式请大神指教,怎么能算出来这个公式Qi,-Pi,-Qi 是从哪里出来的, 牛顿拉夫逊法过程以及其难点和重点?电力系统分析潮流计算中,牛顿拉夫逊法过程的难点和重点是什么?能否简要说明一下其过程, 电力系统分析,简单系统的潮流分析,有功功率无功功率分布是怎样的?在不计网络功率损耗的情况下,下图所示网络各段电路中( C )A、仅有有功功率; B、仅有无功功率; 电路 电力系统分析哪个简单 什么是网损率,怎样计算呢?网损率的定义和计算公式是什么?电力系统分析 在电力系统分析计算中,相比于有名值,标幺制是如何简化计算的?举例说明 电力系统分析和计算中,功率和阻抗分别是用一相还是用三相? 用matlab进行牛顿拉夫逊法的潮流计算的编程, 请问电力系统分析是什么性质的分析? 哪里有电力系统分析的课件 电气工程基础 电力系统分析 这两门课有什么联系? 电力系统分析,额定容量公式分析 电力潮流控制的手段? 跪求《电力系统分析》(华科第三版)习题答案如题[email protected]