【物理應用】基于matlab直角坐標的電力系統潮流計算【含Matlab原始碼 341期】

一、簡介

參照何仰贊電力系統分析（下）的牛拉法程式流程框圖，將程式分為節點導納矩陣的形成、判斂、不平衡量計算、雅可比矩陣計算、節點電壓修正、計算功率等幾個部分，

二、源代碼

clear
%先求節點導納矩陣
N1=30;      %節點數
L1=43;      %支路數
%N1=5;L1=7;
B1=xlsread('XIANLU2.xlsx','B2:F44');
%B1=xlsread('XIANLU.xlsx','B2:F44');
G=zeros(N1);
B=zeros(N1);
%每行存盤一條支路
%第一列存盤支路的一個端點I,變壓器所在端加個負號，勵磁支路為0
%第二列存盤另一個端點J，變壓器所在端加個負號，勵磁支路為0
%第三列存盤支路的電阻R
%第四列存盤支路電抗X
%第五列存盤支路對地電納或變壓器變比k,k在首端為正
for m1=1:L1
     I=B1(m1,1);J=B1(m1,2);R=B1(m1,3);X=B1(m1,4);k=B1(m1,5);
     if I*J==0  %勵磁支路？
        if I==0
             G(J,J)=G(J,J)+R;
             B(J,J)=B(J,J)+X;
        end
         if J==0
             G(I,I)=G(I,I)+R;
             B(I,I)=B(I,I)+X;
         end
     else
         if I*J>0 %線路支路？
            B(I,I)= B(I,I)+ k;
            B(J,J)= B(J,J)+ k;
            k=1;
         else
              if I<0  
                 t=I;
                 I=J;
                 J=t;
              end
              J=abs(J);
              if k<0
                  k=-1/k;
              end
         end
                    G(I,J)=-(1/k)*R/(R^2+X^2);
                    G(J,I)=G(I,J);
                    B(I,J)=(1/k)*X/(R^2+X^2);
                    B(J,I)=B(I,J);
                    G(I,I)=G(I,I)+R/(R^2+X^2);
                    B(I,I)=B(I,I)-X/(R^2+X^2);
                    G(J,J)=G(J,J)+1/(k^2)*R/(R^2+X^2);
                    B(J,J)=B(J,J)-1/(k^2)*X/(R^2+X^2);      
     end
end
Y=G+B*1i;
%求差.
b=xlsread('JIEDIAN2.xlsx','C2:F31');
%b=xlsread('JIEDIAN.xlsx','C2:F31');
precision=1;          %誤差
t=0;                  %迭代次數
pq=24;               %pq節點數
%pq=4;
U=b(pq+2:N1,1);
%開始牛拉法
while precision>0.00001
    P1=zeros(pq,1);
    Q=zeros(pq,1);
    P2=zeros(N1-pq-1,1);
    U2=zeros(N1-pq-1,1);
    deltpqu=zeros(2*(N1-1),1);
    deltu2=zeros(N1-pq-1,1);
    for i=1:pq
        for j=1:N1
            P1(i,1)=P1(i,1)+b(i+1,1)*(G(i+1,j)*b(j,1)-B(i+1,j)*b(j,2))+b(i+1,2)*(G(i+1,j)*b(j,2)+B(i+1,j)*b(j,1));
            Q(i,1) =Q(i,1) +b(i+1,2)*(G(i+1,j)*b(j,1)-B(i+1,j)*b(j,2))-b(i+1,1)*(G(i+1,j)*b(j,2)+B(i+1,j)*b(j,1)); 
        end
    end
    for i=1:N1-pq-1
        for j=1:N1
            P2(i,1)=P2(i,1)+b(pq+i+1,1)*(G(pq+i+1,j)*b(j,1)-B(pq+i+1,j)*b(j,2))+b(pq+i+1,2)*(G(pq+i+1,j)*b(j,2)+B(pq+i+1,j)*b(j,1));
        end
        U2(i,1)=b(pq+i+1,1)^2+b(pq+i+1,2)^2;
    end
    deltp1=b(2:pq+1,3)-P1(1:pq,1);
    deltq =b(2:pq+1,4)-Q(1:pq,1);
    deltp2=b(pq+2:N1,3)-P2(1:N1-pq-1,1);
    for i=1:N1-pq-1
        deltu2(i,1)=U(i,1)^2-U2(i,1);
    end

三、檔案

四、備注

完整代碼或者代寫添加QQ 1564658423
往期回顧>>>>>>
【物理應用】基于matlab楊氏雙縫干涉實驗【含Matlab原始碼 001期】
【物理應用】基于matlab井筒多相流【含Matlab原始碼 055期】
【物理應用】基于matlab 二維對流擴散溫度場【含Matlab原始碼 056期】
【物理應用】基于matlab內彈道【含Matlab原始碼 057期】
【物理應用】基于matlab兩個光纖光柵的新光子微波發生器的仿真分析【含Matlab原始碼 058期】
【物理應用】基于matlab計算油氣井井底壓力【含Matlab原始碼 059期】
【物理應用】基于matlab光流場【含Matlab原始碼 060期】
【物理應用】基于matlab 帶電粒子在放射狀電場和均勻磁場下的混沌運動模擬【含Matlab原始碼 061期】
【物理應用】基于matlab異構交通流【含Matlab原始碼 062期】
【物理應用】基于matlab模擬斜拋運動含GUI界面【含Matlab原始碼 063期】
【物理應用】基于matlab計算晶體結構的x射線衍射圖譜【含Matlab原始碼 064期】
【物理應用】基于matlab衍射光柵教學版【含Matlab原始碼 065期】
【物理應用】基于matlab電荷電場線和電勢【含Matlab原始碼 066期】
【物理應用】基于matlab速度場找到漩渦【含Matlab原始碼 067期】
【物理應用】基于matlab車道檢測仿真【含Matlab原始碼 068期】
【物理應用】基于matalb GUI界面的電偶極子【含Matlab原始碼 069期】
【物理應用】基于matlab車載自組織網路中路邊性能及防碰撞演算法研究【含Matlab原始碼 263期】
【物理應用】基于matlab Q學習的無線體域網路由方法【含Matlab原始碼 264期】
【物理應用】基于matlab GUI家庭室內溫濕度控制【含Matlab原始碼 326期】
【電路仿真】基于simulink三相相控變流器的設計【含Matlab原始碼 327期】
【物理應用】基于matlab的導線網坐標計算【含Matlab原始碼 328期】
【物理應用】基于matlab的功率譜估計【含Matlab原始碼 329期】
【物理應用】基于動態矩陣控制方法的熱水器出口溫度控制【含Matlab原始碼 330期】
【電路仿真】基于matlab simulink三相半波可控整流電路【含Matlab原始碼 331期】
【電路仿真】基于matlab simulink三相橋式全控整流電路【含Matlab原始碼 332期】
【電路仿真】基于simulink模擬鎖相環設計【含Matlab原始碼 339期】
【電路仿真】基于simulink全數字鎖相環設計【含Matlab原始碼 340期】