哈工大电力系统分析短路及潮流计算实验上机程序文档格式.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======================
升级会员 