哈工大电力系统分析短路及潮流计算实验上机程序文档格式.doc

1、,2）;B=fscanf（fid,%f,5,A（2,1）;fclose（fid）;B=B;%求节点导纳矩阵Y=zeros（A（1,1）; % 形成55的0阵%=%导纳阵元素计算for a=1:1:A（2,1） m=B（a,1）; n=B（a,2）; if B（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）; elseif B（a,5） Y（m,m）=Y（m,m）+1/（B（a,3）

2、+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）; endend%=Y（2,2）=Y（2,2）+j*0.01;Z=inv（Y）;If=1/Z（3,3）; %接金属短路，Zf=0for i=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.txtwfprin

3、tf（fid,Y_matrixn%输出导纳阵 for j=1: k=Y（i,j）; re=real（k）; fprintf（fid,%8.4f,re）; im=imag（k）; if im=0+%8.4fjt end endn%=%输出阻抗阵，导纳阵求逆Z_matrixn k=Z（i,j）;%=%打印输出短路电流 fprintf（fid,If= re=real（If）; im=imag（If）;nVn%打印输出节点电压Note%d:V%d=,i,i）; k=V（1,i）;%输出打印支路电流Ibrn4 k=I（i,j）; if（re=0|im=0） fprintf（fid,Branch%d-%d

4、:I%d%d=,i,j,i,j）; fprintf（fid,%8.4fjn+%8.4fjn %关闭文件附：2、程序输出结果节点导纳阵： 节点阻抗阵： 短路电流：If= 0.0001+ 0.1082j节点电压：各支路电流：四、实验总结这是我的第一次上机实验，感觉稍微有点难，主要还是在工具软件C语言或者MATLAB的运用上，但是我相信，以后学习中，我会努力掌握的，这是我把理论应用于实际中的必要桥梁！根据所给的电力系统，编制牛顿法潮流计算程序，通过计算机进行调试，最后完成一个切实可行的电力系统计算应用程序。电力系统潮流计算的计算程序设计及编制和调试。三、程序框图四、实验过程1、实验程序及说明node

5、4.txt %打开输入数据,8）; %读8个数C=fscanf（fid,3,（A（1,1）-1）;C=CB（2,5）=1/B（2,5）; %得44的0阵%=%与前一实验同法求导纳阵 else if B（a,5） 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迭代计数tt V1=e1+jf1tttt V2=e2+jf2ttttV3=e3+jf3ttn%=求W阵= W=zeros（2*（A（1,1）-1）,

6、50）; for x=1:50 %设置迭代次数为50次A（7,1） k=C（i,1）; p=0; q=0; m=0; n=0; for j=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节点=for l=1:A（8,1

7、） k=C（l+A（7,1）,1）; for j=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）;if Max%8.4f+%8.4fjtt ,Q（2*k-1,1）,Q（2*k,1）; else%8.4f%8.4fjtt end%=平衡点功率= k=A（3,1）; v=0; j=sqrt（-1）; for b=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=