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

电力系统分析作业——电网节点导纳矩阵的计算机形成编程软件：matlab R2025b程序说明：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;end[m,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经手算校验，程序结果准确。