1.完整项目描述和程序获取
>面包多安全交易平台:https://mbd.pub/o/bread/Y56bm51w
>如果链接失效，可以直接打开本站店铺搜索相关店铺：
>如果链接失效，程序调试报错或者项目合作也可以加微信或者QQ联系。
2.部分仿真图预览
3.算法概述
在电力工程中，“潮流”还特指电网各处电压（包括幅值与相角）、有功功率、无功功率等的分布。潮流的分布是运行调度单位和维修部门所必须知道的事项。而潮流计算，是指给定电网中一些参数、已知值和未知值中假设的初始值，通过重复迭代，最终求出潮流分布的精确值，常用方法有牛顿-拉夫逊法和PQ分解法。
4.部分源码
............................................................
Y=G+B*1i;
%求差.
b=xlsread('TEST2.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
.....................................................
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
deltu2;
for i=1:pq
deltpqu(2*i-1,1)=deltp1(i);
deltpqu(2*i,1)=deltq(i);
end
for i=1:N1-pq-1
deltpqu(2*pq+2*i-1,1)=deltp2(i);
deltpqu(2*pq+2*i,1)=deltu2(i);
end
deltpqu; %delta
%求雅克比
aii=zeros(N1);bii=zeros(N1);
H1=zeros(N1);N11=zeros(N1);
for i=1:N1
for j=1:N1
aii(i,1)=aii(i,1)+G(i,j)*b(j,1)-B(i,j)*b(j,2);
bii(i,1)=bii(i,1)+G(i,j)*b(j,2)+B(i,j)*b(j,1);
end
end
%H、N
...................................................
Jacobi;
%修正
Correction=Jacobi^-1*deltpqu;
for i=1:N1-1
b(i+1,1)=b(i+1,1)+Correction(2*i);
b(i+1,2)=b(i+1,2)+Correction(2*i-1);
end
precision=max(abs(deltpqu));
%precision=max(abs(Correction));
t=t+1;
b;
end
%平衡节点功率
pp=0;
for j=1:N1
pp=pp+conj(Y(1,j))*(b(j,1)-b(j,2)*1i);
end
b(1,3)=real(b(1,1)*pp);
b(1,4)=imag(b(1,1)*pp);
%各线路功率
S=zeros(N1); y=zeros(N1); y0=zeros(N1);
for i=1:N1
for j=1:N1
if i~=j
y(i,j)=-Y(i,j);
end
end
end
y;
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
y0(I,J)=k*1i; y0(J,I)=k*1i;
end
if I*J<0
if I<0
t=I;
I=J;
J=t;
end
J=abs(J);
if k<0
k=-1/k;
end
y0(I,J)=(R/(R^2+X^2)-X/(R^2+X^2)*1i)*(k-1)/k;
y0(J,I)=(R/(R^2+X^2)-X/(R^2+X^2)*1i)*(1-k)/k^2;
end
end
.......................................
A343
