电力系统分析作业——电网节点导纳矩阵的计算机形成(共2页).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; y(enode,enode)=y(enode,enode)+b; end 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经手算校验，程序结果准确。专心-专注-专业
