请教各位大神，关于潮流计算不收敛的问题，自己看了好几遍，没看出哪有问题，请求各路大神有时间的话，看一下，指导一下。牛拉法直角坐标，存在问题：迭代时，DW和DU值不停变大。

weixin_45105277

于 2020-03-07 14:58:33 发布

阅读量4.3k

点赞数

文章标签： matlab 算法

本文链接：https://blog.csdn.net/weixin_45105277/article/details/104714858

版权

博主在进行潮流计算时遇到不收敛的困扰，DW和DU值在迭代过程中持续增大，针对牛拉法直角坐标系统，多次检查未发现错误，期望专业人员提供指导。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

Ybus=[1.25-5.5i,-0.5+3i,-0.75+2.5i;-0.5+3i,1.3-7i,-0.8+4i;-0.75+2.5i,-0.8+4i,1.55-6.5i];
sbus1=[2+1i;0.5;0];
bus1=[1,1;2,2;3,3];
n=3;
JKB=zeros(2*(n-1),2*(n-1));
H=zeros(n-1,n-1);
J=zeros(n-1,n-1);
N=zeros(n-1,n-1);
L=zeros(n-1,n-1);
R=zeros(n-1,n-1);
S=zeros(n-1,n-1);
aii=zeros(n-1,1);
bii=zeros(n-1,1);
G=real(Ybus);
B=imag(Ybus);
P=real(sbus1);%给定有功功率
Q=imag(sbus1);%给定无功功率
V1=ones(3,1);
e=real(V1);%节点初始电压实部
f=imag(V1);
DW=zeros(2*(n-1),1);
DU=zeros(2*(n-1),1);
maxdata=1;Maxdata=0.00001;
T=1;%注入功率
while maxdata>Maxdata;
    for i=1:n-1%n是平衡节点
        P1=0;Q1=0;U1=0;
        if bus1(i,2)==1%当PQ节点时
            for j=1:n
                P1=P1+e(i)*(G(i,j)*e(j)-B(i,j)*f(j))+f(i)*(G(i,j)*f(j)+B(i,j)*e(j));%功率方程
                Q1=Q1+f(i)*(G(i,j)*e(j)-B(i,j)*f(j))-e(i)*(