哈工大电力系统分析短路及潮流计算实验上机程序文档格式.doc
- 文档编号:6986217
- 上传时间:2023-05-07
- 格式:DOC
- 页数:13
- 大小:247.50KB
哈工大电力系统分析短路及潮流计算实验上机程序文档格式.doc
《哈工大电力系统分析短路及潮流计算实验上机程序文档格式.doc》由会员分享,可在线阅读,更多相关《哈工大电力系统分析短路及潮流计算实验上机程序文档格式.doc(13页珍藏版)》请在冰点文库上搜索。
,2);
B=fscanf(fid,'
%f'
[5,A(2,1)]);
fclose(fid);
B=B'
;
%求节点导纳矩阵
Y=zeros(A(1,1));
%形成5×
5的0阵
%%%%%====================================================
%导纳阵元素计算
fora=1:
1:
A(2,1)
m=B(a,1);
n=B(a,2);
ifB(a,5)>
Y(m,m)=Y(m,m)+1/(B(a,3)+j*B(a,4));
Y(n,n)=Y(n,n)+1/(B(a,3)+j*B(a,4))/(B(a,5)^2);
Y(m,n)=-1/(B(a,3)+j*B(a,4))/B(a,5);
Y(n,m)=Y(m,n);
elseifB(a,5)<
Y(m,m)=Y(m,m)+1/(B(a,3)+j*B(a,4))-j*B(a,5);
Y(n,n)=Y(n,n)+1/(B(a,3)+j*B(a,4))-j*B(a,5);
Y(m,n)=-1/(B(a,3)+j*B(a,4));
end
end
%%%%%=======================================================
Y(2,2)=Y(2,2)+j*0.01;
Z=inv(Y);
If=1/Z(3,3);
%接金属短路,Zf=0
fori=1:
5
V(1,i)=1-Z(i,3)/Z(3,3);
%%%========================================================
%计算各支路电流
I(m,n)=(B(a,5)*V(1,m)-V(1,n))/(B(a,3)+j*B(a,4));
0
I(m,n)=(V(1,m)-V(1,n))/(B(a,3)+j*B(a,4));
%%%=======================================================
ans.txt'
w'
fprintf(fid,'
Y_matrix\n'
%输出导纳阵
forj=1:
k=Y(i,j);
re=real(k);
fprintf(fid,'
%8.4f'
re);
im=imag(k);
ifim<
fprintf(fid,'
%8.4fj\t'
im);
elseifim>
=0
+%8.4fj\t'
end
end
\n'
%%%==========================================================
%输出阻抗阵,导纳阵求逆
Z_matrix\n'
k=Z(i,j);
%%%=========================================================
%打印输出短路电流
fprintf(fid,'
If='
re=real(If);
im=imag(If);
\nV\n'
%打印输出节点电压
Note%d:
V%d='
i,i);
k=V(1,i);
%输出打印支路电流
Ibr\n'
4
k=I(i,j);
if(re~=0||im~=0)
fprintf(fid,'
Branch%d--%d:
I%d%d='
i,j,i,j);
fprintf(fid,'
%8.4fj\n'
+%8.4fj\n'
%关闭文件
附:
'
2、程序输出结果
节点导纳阵:
节点阻抗阵:
短路电流:
If=0.0001+0.1082j
节点电压:
各支路电流:
四、实验总结
这是我的第一次上机实验,感觉稍微有点难,主要还是在工具软件—C语言或者MATLAB的运用上,但是我相信,以后学习中,我会努力掌握的,这是我把理论应用于实际中的必要桥梁!
根据所给的电力系统,编制牛顿法潮流计算程序,通过计算机进行调试,最后完成一个切实可行的电力系统计算应用程序。
电力系统潮流计算的计算程序设计及编制和调试。
三、程序框图
四、实验过程
1、实验程序及说明
node4.txt'
%打开输入数据
8);
%读8个数
C=fscanf(fid,'
[3,(A(1,1)-1)]);
C=C'
B(2,5)=1/B(2,5);
%得4×
4的0阵
%%%===============================================================
%与前一实验同法求导纳阵
elseifB(a,5)<
%%%=================================================================
%节点电压赋初值
Q=zeros(2*A(1,1),1);
fori=1:
(A(1,1)-1)
ifC(i,3)>
Q(2*i-1,1)=C(i,3);
Q(2*i,1)=0;
else
Q(2*i-1,1)=1;
Q(2*i,1)=0;
Q(2*A(3,1)-1,1)=A(4,1);
Q(2*A(3,1),1)=0;
answer.txt'
====================节点电压V=============================\n'
迭代计数\t\tV1=e1+jf1\t\t\t\tV2=e2+jf2\t\t\t\tV3=e3+jf3\t\t\n'
%%%=======================求W阵===================================
W=zeros(2*(A(1,1)-1),50);
forx=1:
50%设置迭代次数为50次
A(7,1)
k=C(i,1);
p=0;
q=0;
m=0;
n=0;
forj=1:
A(1,1)
g=real(Y(k,j));
b=imag(Y(k,j));
e=Q(2*j-1,1);
f=Q(2*j,1);
p=p+g*e-b*f;
q=q+g*f+b*e;
m=m+g*e-b*f;
n=n+g*f+b*e;
end
W(2*k,x)=C(k,2)-Q(2*k-1,1)*p-Q(2*k)*q;
W(2*k-1,x)=C(k,3)-Q(2*k)*m+Q(2*k-1,1)*n;
End
%%%=========================PQ节点=================================
%%%=========================PV节点=================================
forl=1:
A(8,1)
k=C(l+A(7,1),1);
forj=1:
end
W(2*k-1,x)=(C(k,3))^2-(Q(2*k-1,1)^2+Q(2*k)^2);
%%%================================================================
%比较是否符合条件
Max=0;
2*(A(1,1)-1)
Max=max(abs(W(i,x)),Max);
ifMax<
A(5,1)
break;
%求雅克比行列式
A(1,1)-1
%===========================================
%===========================i~=j
ifi~=j
k=C(i,1);
g=real(Y(k,j));
b=imag(Y(k,j));
e=Q(2*k-1,1);
f=Q(2*k,1);
J(2*k-1,2*j-1)=-(g*e+b*f);
J(2*k-1,2*j)=b*e-g*f;
ifi~=3
J(2*k,2*j)=-J(2*k-1,2*j-1);
J(2*k,2*j-1)=J(2*k-1,2*j);
else
J(2*k,2*j-1)=0;
J(2*k,2*j)=0;
end
%===============================i=j==================
ifi==j
k=C(i,1);
p=0;
q=0;
m=0;
n=0;
g=real(Y(k,k));
b=imag(Y(k,k));
J(2*k-1,2*j-1)=-g*Q(2*i-1,1)-b*Q(2*i,1);
J(2*k-1,2*j)=b*Q(2*i-1,1)-g*Q(2*i,1);
ifk~=3
J(2*k,2*j-1)=b*Q(2*i-1,1)-g*Q(2*i,1);
J(2*k,2*j)=g*Q(2*i-1,1)+b*Q(2*i,1);
end
forr=1:
g=real(Y(k,r));
b=imag(Y(k,r));
e=Q(2*r-1,1);
f=Q(2*r,1);
p=p+g*e-b*f;
q=q+g*f+b*e;
end
J(2*k-1,2*j-1)=-p+J(2*k-1,2*j-1);
J(2*k-1,2*j)=-q+J(2*k-1,2*j);
ifi~=3
J(2*k,2*j-1)=q+J(2*k,2*j-1);
J(2*k,2*j)=-p+J(2*k,2*j);
else
J(2*k,2*j-1)=-2*Q(2*i-1,1);
J(2*k,2*j)=-2*Q(2*i,1);
end
end
%%%==============================================================
%解修正方程得修正量Q0
fori=1:
2:
J0(i+1,:
)=J(i,:
J0(i,:
)=J(i+1,:
J=J0;
Q0=-inv(J)*W(:
x);
fori=1:
6
Q(i,1)=Q(i,1)+Q0(i,1);
End
\t%d\t\t'
k=C(i,1);
ifQ(2*k,1)>
%8.4f+%8.4fj\t\t'
Q(2*k-1,1),Q(2*k,1));
else
%8.4f%8.4fj\t\t'
end
%%%==================================平衡点功率======================
k=A(3,1);
v=0;
j=sqrt(-1);
forb=1:
m=conj(Y(k,b));
p=Q(2*b-1,1)+Q(2*b,1)*j;
n=conj(p);
v=v+m*n;
Wp=(Q(2*k-1,1)+j*Q(2*k,1))*v;
%%%===========================节点电压V======================
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 哈工大 电力系统 分析 短路 潮流 计算 实验 上机 程序