电力系统短路电流计算matlab程序.pdf
版权申诉
5星 · 超过95%的资源 86 浏览量
2021-12-20
09:46:20
上传
评论 1
收藏 6KB PDF 举报 
%电力系统极坐标下的牛顿 -拉夫逊法潮流计算
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('迭代中的雅克比矩阵 :');
liuyeping111
- 粉丝: 1
- 资源: 4万+
下载权益
C知道特权
VIP文章
课程特权
VIP享7折，此内容立减5.97元 开通VIP
相关推荐
资源上传下载、课程学习等过程中有任何疑问或建议，欢迎提出宝贵意见哦~我们会及时处理！
点击此处反馈
评论2
最新资源