电力系统短路电流计算-matlab 程序
%电力系统极坐标下的牛顿-拉夫逊法潮流计算
disp('电力系统极坐标下的牛顿-拉夫逊法潮流计算:');
clear
n=input('请输入结点数：n=');
n1=input('请输入 PV 结点数：n1=');
n2=input('请输入 PQ 结点数：n2=');
isb=input('请输入平衡结点：isb=');
pr=input('请输入精确度：pr=');
K=input('请输入变比矩阵看:K=');
C=input('请输入支路阻抗矩阵：C=');
y=input('请输入支路导纳矩阵：y=');
U=input('请输入结点电压矩阵：U=');
S=input('请输入各结点的功率：S=');
Z=zeros(1,n);N=zeros(n1+n2,n2);L=zeros(n2,n2);QT1=zeros(1,n1+n2);
for m=1:n
for R=1:n
C(m,m)=C(m,m)+y(m,R);
if K(m,R)~=0
C(m,m)=C(m,m)+1/((K(m,R)*C(m,R))/(K(m,R)-1));
C(R,R)=C(R,R)+1/((K(m,R)^2*C(m,R))/(1-K(m,R)));
C(m,R)=C(m,R)*K(m,R);
C(R,m)=C(m,R);
end
end
end
for m=1:n
for R=1:n
if m~=R
Z(m)=Z(m)+1/C(m,R);
end
end
end
for m=1:n
for R=1:n
if m==R
Y(m,m)=C(m,m)+Z(m);
else
Y(m,R)=-1/C(m,R);
end
end
end
disp('结点导纳矩阵:');
disp(Y);
disp('迭代中的雅克比矩阵:');
1
资源评论
