资源描述：
《哈工大电力系统分析短路及潮流计算实验上机程序.doc》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、上机实验实验一节点导纳阵生成和短路电流计算实验二简单系统的牛顿法潮流计算姓名：班级：141班学号：实验说明：本程序的电路结构来自翁增银、何仰赞主编的《电力系统分析》的例题实验一节点导纳阵生成和短路电流计算一、实验目的根据所给的电力系统，编制短路电流计算程序，通过计算机进行调试，最后完成一个切实可行的电力系统计算应用程序。通过自己设计电力系统计算程序加深对电力系统分析的理解，同时加强计算机实际应用能力。二、实验内容1、编写数据输入、形成导纳阵程序2、电力系统短路计算实用公式的计算程序及编制和调试三、实验过程1、程序代码及说明%清屏clear%读数据fid=fopen(
2、'node5.txt','r');A=fscanf(fid,'%d',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)>0Y(m,m)=Y(m,m)+1/(B(a,3)+j*B(a,4));Y(n,n)=Y(n,n)+1/(B(a,3)+j*
3、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)<0Y(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));Y(n,m)=Y(m,n);endend%%%%%=======================================================Y(2,2)=Y(2
4、,2)+j*0.01;Z=inv(Y);If=1/Z(3,3);%接金属短路，Zf=0fori=1:1:5V(1,i)=1-Z(i,3)/Z(3,3);end%%%========================================================%计算各支路电流fora=1:1:A(2,1)m=B(a,1);n=B(a,2);ifB(a,5)>0I(m,n)=(B(a,5)*V(1,m)-V(1,n))/(B(a,3)+j*B(a,4));elseifB(a,5)<0I(m,n)=(V(1,m)-V(1,n))/(B(a,3)+j*B(
5、a,4));endend%%%=======================================================fid=fopen('ans.txt','w');fprintf(fid,'Y_matrix');%%%=======================================================%输出导纳阵fori=1:1:5forj=1:1:5k=Y(i,j);re=real(k);fprintf(fid,'%8.4f',re);im=imag(k);ifim<0fprintf(fid,'%8.4fj
6、t',im);elseifim>=0fprintf(fid,'+%8.4fjt',im);endendfprintf(fid,'');end%%%==========================================================%%%==========================================================%输出阻抗阵，导纳阵求逆fprintf(fid,'Z_matrix');fori=1:1:5forj=1:1:5k=Z(i,j);re=real(k);fprintf(fid,'
7、%8.4f',re);im=imag(k);ifim<0fprintf(fid,'%8.4fjt',im);elseifim>=0fprintf(fid,'+%8.4fjt',im);endendfprintf(fid,'');end%%%=========================================================%%%=========================================================%打印输出短路电流fprintf(fid,'If=');re=real(If);fpr
