电力系统分析作业——电网节点导纳矩阵的计算机形成

电力系统分析作业——电网节点导纳矩阵的计算机形成编程软件：matlabR2010b程序说明：1.如果已经输入i-j支路的信息，则不可再输入j-i支路的信息。2.变压器支路的第一个节点编号默认为变压器一次侧，即变压器的等值电路中的阻抗归算侧，亦即变压器非标准变比的1:k中的‘1’。3.标幺值等值电路中，如果变比为1:1,则默认为线路，因此，变压器的非标准变比不可以是1:1。5.如果变压器支路也有导纳B不为零，则说明此导纳就是励磁导纳，与线路的导纳B/2不同含义，只算作变压器原边的自导纳。4.由于程序执行的是复数运算，所以即使实部为零时，也会输出实部‘0’。程序代码:a=load('data.txt');%从’data.txt’中读入数据[m,n]=size(a);w=1i;u=1;while(u=m)hnode=a(u,1);enode=a(u,2);z=a(u,3)+a(u,4)*w;b=a(u,5)*w;k=a(u,6);y(hnode,enode)=-1/(k*z);y(enode,hnode)=-1/(k*z);y(hnode,hnode)=y(hnode,hnode)+1/(k*z)+(k-1)/(k*z);y(enode,enode)=y(enode,enode)+1/(k*z)+(1-k)/(k*k*z);if(abs(k-1)0.0001)%如果为线路y(hnode,hnode)=y(hnode,hnode)+b;y(enode,enode)=y(enode,enode)+b;endif(abs(k-1)0.0001)%如果为变压器y(hnode,hnode)=y(hnode,hnode)-b;endu=u+1;end[m,n]=size(y);disp(‘Y=’);disp(y(1:m,1:n));clear;算例输入数据：首端编号末端编号电阻电抗电纳/2变比230.080.300.2514200.01501.055300.0301.05120.040.250.251130.10.3501输出数据：Y=1.3787-6.2917i-0.6240+3.9002i-0.7547+2.6415i00-0.6240+3.9002i1.4539-66.9808i-0.8299+3.1120i0+63.4921i0-0.7547+2.6415i-0.8299+3.1120i1.5846-35.7379i00+31.7460i00+63.4921i00-66.6667i0000+31.7460i00-33.3333i经手算校验，程序结果准确。