%b1是支路引數矩陣，第一列和第二列是起始節點編號和終點節點編號
%第三列、第四列、第五列、第六列、第七列、第八列分別為：支路電阻、電抗、電導、電納、變壓器變比、是否有變壓器（1為有、0為無）。
%b2為節點引數矩陣，其中第一列到第六列為節點編號；為節點型別；注入有功、注入無功、電壓幅值、電壓相位。
%「1」為平衡節點，「2」為pq節點，「3」為pv節點引數。
n=input('請輸入節點數:n=');
n1=input('請輸入支路數:n1=');
isb=input('請輸入平衡節點號:isb=');
pr=input('請輸入誤差精度:pr=');
b1=input('請輸入支路引數:b1=');
b2=input('請輸入節點引數:b2=');
y=zeros(n
n=1%建立節點導納矩陣
for i=1:n1
if b1(i,8)==0
p=b1(i,1);
q=b1(i,2);
y(p,q)=y(p,q)-1/(b1(i,3)+b1(i,4))-b1(i,5);
y(q,p)=y(p,q);
y(p,p)=y(p,p)+1/(b1(i,3)+b1(i,4))+0.5*b1(i,6)+b1(i,5);
y(q,q)=y(q,q)+1/(b1(i,3)+b1(i,4))+0.5*b1(i,6)+b1(i,5);
else
p=b1(i,1);
q=b1(i,2);
y(p,q)=y(p,q)-1/((b1(i,3)+b1(i,4))*b1(i,7)) -b1(i,5);
y(q,p)=y(p,q);
y(p,p)=y(p,p)+1/(b1(i,3)+b1(i,4))+b1(i,5);
y(q,q)=y(q,q)+1/(b1(i,7)^2*(b1(i,3)+b1(i,4)))+b1(i,5);
endendy
g=real(y);
b=imag(y);
pris=zeros(2*n-2,1);
imbs=zeros(2*n-2,1);
%建立pris，用於儲存初始功率引數
h=0;
j=0;
for i=1:n
if i~=isb&b2(i,2)==1
h=h+1;
for j=1:n
pris(2*h-1,1)=pris(2*h-1,1)+b2(i,5)*(g(i,j)*b2(j,5)-b(i,j)*b2(j,6))+b2(i,6)*(g(i,j)*b2(j,6)+b(i,j)*b2(j,5
pris(2*h,1)=pris(2*h,1)+b2(i,6)*(g(i,j)*b2(j,5)-b(i,j)*b2(j,6))-b2(i,5)*(g(i,j)*b2(j,6)+b(i,j)*b2(j,5
endendendfor i=1:n
if i~=isb&b2(i,2)==2
h=h+1;
for j=1:n
pris(2*h-1,1)=pris(2*h-1,1)+ b2(i,5)*(g(i,j)*b2(j,5)-b(i,j)*b2(j,6))+b2(i,6)*(g(i,j)*b2(j,6)+b(i,j)*b2(j,5));
pris(2*h,1)=pris(2*h,1)+ b2(i,6)*(g(i,j)*b2(j,5)-b(i,j)*b2(j,6))-b2(i,5)*(g(i,j)*b2(j,6)+b(i,j)*b2(j,5))
endendendpris
u3=zeros(n-h-1,2);%u3儲存pv節點的初始電壓
t=0;
for i=1:n
if b2(i,2)==2
t=t+1;
u3(t,1)=b2(i,5);
u3(t,2)=b2(i,6);
endendu3
%imbs於儲存有功功率、無功功率和電壓幅值的不平衡量
h=0;
for i=1:n
if i~=isb&b2(i,2)==1
h=h+1;
imbs(2*h-1,1)=b2(i,3)-pris(2*h-1,1);
imbs(2*h,1)=b2(i,4)-pris(2*h,1);
endendt=0;
for i=1:n
if i~=isb&b2(i,2)==2
h=h+1;
t=t+1;
imbs(2*h-1,1)=b2(i,3)-pris(2*h-1,1);
imbs(2*h,1)=u3(t,1)^2+u3(t,2)^2-b2(i,5)^2-b2(i,6)^2;
endendimbs
i=zeros(n-1,1);%i，儲存節點電流引數
h=0;
for i=1:n
if i~=isb
h=h+1;
i(h,1)=(pris(2*h-1,1)-pris(2*h,1)*sqrt(-1))/conj(b2(i,5)+b2(i,6)*sqrt(-1
endendi
jacbi=zeros(2*n-2);%jacbi(雅可比矩陣)
h=0;
k=0;
for i=1:n
if b2(i,2)==1
h=h+1;
for j=1:n
if j~=isb
k=k+1;
if i==j
jacbi(2*h-1,2*k-1)=-b(i,j)*b2(i,5)+g(i,j)*b2(i,6)+imag(i(h,1
jacbi(2*h-1,2*k)=g(i,j)*b2(i,5)+b(i,j)*b2(i,6)+real(i(h,1));
jacbi(2*h,2*k-1)=-jacbi(2*h-1,2*k)+2*real(i(h,1));
jacbi(2*h,2*k)=jacbi(2*h-1,2*k-1)-2*imag(i(h,1));
else
jacbi(2*h-1,2*k-1)=-b(i,j)*b2(i,5)+g(i,j)*b2(i,6);
jacbi(2*h-1,2*k)=g(i,j)*b2(i,5)+b(i,j)*b2(i,6);
jacbi(2*h,2*k-1)=-jacbi(2*h-1,2*k);
jacbi(2*h,2*k)=jacbi(2*h-1,2*k-1);
endif k==(n-1)
k=0;
endend
endendendk=0;
for i=1:n
if b2(i,2)==2
h=h+1;
for j=1:n
if j~=isb
k=k+1;
if i==j
jacbi(2*h-1,2*k-1)= -b(i,j)*b2(i,5)+g(i,j)*b2(i,6)+imag(i(h,1));
jacbi(2*h-1,2*k)= g(i,j)*b2(i,5)+b(i,j)*b2(i,6)+real(i(h,1));
jacbi(2*h,2*k-1)=2*b2(i,6);
jacbi(2*h,2*k)=2*b2(i,5);
else
jacbi(2*h-1,2*k-1)= -b(i,j)*b2(i,5)+g(i,j)*b2(i,6);
jacbi(2*h-1,2*k)= g(i,j)*b2(i,5)+b(i,j)*b2(i,6);
jacbi(2*h,2*k-1)=0;
jacbi(2*h,2*k)=0;
endif k==(n-1)
k=0;
endend
endend
潮流計算說明
用matlab計算 電子數字計算機的應用 電力系統計算 西安交通大學 清華大學 浙江大學 湖南大學 成都工學院 水利電力部電網排程研究所合編,水利電力出版社，1978。計算第151頁例3 1。節點電壓計算結果在第156頁例表3 5中。功率計算結果在第168頁例圖3 4中。1.使用雅克比迭代法進行潮流...
潮流計算方法
本人乙個星期以來的汗水，弄清楚了計算機演算法計算潮流的基礎，如果有什麼不懂的可以發資訊到郵箱 接下來開始弄潮流的優化問題，吼吼！電力系統的潮流計算的計算機演算法 以matlab為環境 這裡理論不做過多介紹，推薦一本專門講解電力系統分析的計算機演算法的書籍電力系統分析的計算機演算法 邱曉燕 劉天琪編著...
什麼是潮流計算
潮流計算，電力學名詞，指在給定電力系統網路拓撲 元件引數和發電 負荷參量條件下，計算有功功率 無功功率及電壓在電力網中的分布。潮流計算是根據給定的電網結構 引數和發電機 負荷等元件的執行條件，確定電力系統各部分穩態執行狀態引數的計算。通常給定的執行條件有系統中各電源和負荷點的功率 樞紐點電壓 平衡點...