1、);A=fscanf(fid,%d,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)
2、=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); 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)
3、;%=ans.txtwfprintf(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.4fjtn%=%输出阻抗阵,导纳阵求逆Z_matrixn k=Z(i,j);%=%打印输出短路电流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)Branch%d-%d:I%d%d=,i,j,i,j);%
4、8.4fjn+%8.4fjn %关闭文件附:2、程序输出结果节点导纳阵: 节点阻抗阵: 短路电流:If= 0.0001+ 0.1082j 节点电压:各支路电流:四、实验总结这是我的第一次上机实验,感觉稍微有点难,主要还是在工具软件C语言或者MATLAB的运用上,但是我相信,以后学习中,我会努力掌握的,这是我把理论应用于实际中的必要桥梁!根据所给的电力系统,编制牛顿法潮流计算程序,通过计算机进行调试,最后完成一个切实可行的电力系统计算应用程序。电力系统潮流计算的计算程序设计及编制和调试。三、程序框图四、实验过程1、实验程序及说明node4.txt %打开输入数据,8); %读8个数C=fscan
5、f(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*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),50); for x=1:50 %设置迭代次数为50次A(7,1) k=C(i,1)
6、; p=0; q=0; m=0; n=0;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; 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) k=C(l+A(7,1),1); W(2*k-1,x)=(C(k,3)2-(Q(2*k-1,1)2+Q(2
7、*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);%8.4f%8.4fjtt%=平衡点功率= 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=%fid=fopen(节点电压Vnnode%d:%8.4f+%8.4fjn,k,Q(2*k-1,1),Q(2*k,1);%8.4f%8.4fjn平衡点功率P+jQ=%8.4f+%8.4fj,real(Wp),imag(Wp);附注:node4.txt %输入数据五、实验总结 本次实验是我把应用理论知识的重要实践,经过实验,我有两点感想,首先,作为一名工科学生,应该能熟练运用C语言和MATLAB等工具,其次,理论如果不用于实践,就永远不知道理论是用来干什么的,学到头一直是满脑子的浆糊,所以,以我的切身经历建议,把这门实验放在跟课程平行的时间上进行,这样不仅有利于实验开展,也有利于学生更加深刻地学习!