function [U,I]= BackforwardSweep(Dr,S,U0) % Dr-直接求解矩阵,S-节电功率需求
Uo-根节点电压
% U-节点电压 I-节点电流
%设置收敛精度
e=1e-5;;%由3改成5更满足实际工程
%设置最大迭代次数
max_iter=50;%由30改成50,更满足实际工程
I=conj(S./U0);
Uf=U0;
for ii=1:max_iter
U=U0-Dr*I;
et=abs(sqrt(Uf.*conj(Uf))-sqrt(U.*conj(U)));
if sum(et)<e
break;
end
I=conj(S./U);
Uf=U;
end
I=conj(S./U);
end
评论0