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

电力系统分析作业电网节点导纳矩阵的计算机形成编程软件：matlab R2010b程序说明：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; end u=u+1;endm,n=size(y);disp(Y=);disp(y(1:m,1:n);clear;算例输入数据：首端编号 末端编号 电阻 电抗 电纳/2 变比 2 3 0.08 0.30 0.25 1 4 2 0 0.015 0 1.05 5 3 0 0.03 0 1.05 1 2 0.04 0.25 0.25 1 1 3 0.1 0.35 0 1输出数据：Y=1.3787 - 6.2917i -0.6240 + 3.9002i -0.7547 + 2.6415i 0 0 -0.6240 + 3.9002i 1.4539 -66.9808i -0.8299 + 3.1120i 0 +63.4921i 0 -0.7547 + 2.6415i -0.8299 + 3.1120i 1.5846 -35.7379i 0 0 +31.7460i0 0 +63.4921i 0 0 -66.6667i 0 0 0 0 +31.7460i 0 0 -33.3333i经手算校验，程序结果准确。