《节点导纳矩阵的建立.doc》由会员分享，可在线阅读，更多相关《节点导纳矩阵的建立.doc（3页珍藏版）》请在金锄头文库上搜索。
1、如上图所示的简单电力系统中，网络各元件参数的标幺值如下：z12=0.10+j0.40 y120=y210=j0.01528 z13=j0.3,k=1.1 z14=0.12+j0.50y140=y410=j0.01920 z24=0.08+j0.40 y240=y420=j0.01413系统中节点1、2为PQ节点，节点3为PV节点，节点4为平衡节点。节点导纳矩阵的运行程序如下：clcCleardisp(网络各元件参数用标幺值表示！);N0=input(请输入节点数:N0=); n1=input(请输入支路数:n1=); l=input(请输入PQ节点的个数=); for m=1:lc(m)=in
2、put(请输入第,num2str(m),个PQ节点的节点号为：); endt=input(请输入PV节点的个数=); for m=1:t c(m)=input(请输入第,num2str(m),个PV节点的节点号为：); endb=input(请输入平衡节点号:b=); %由支路参数形成矩阵B1disp(各支路连接情况：)i=1; for m=1:n1 syms Y N p=input(第,num2str(m),条支路的起始节点); q=input(第,num2str(m),条支路的终止节点); mn=input(第,num2str(m),条支路是否有变压器（请输入Y或N）); y=0;k=1;
3、 if mn=Y; k=input(请输入变压器变比(标幺值):); z=input(请输入第,num2str(m),条支路的线路阻抗); else z=input(请输入第,num2str(m),条支路的线路阻抗：); y=input(请输入第,num2str(m),条支路线路的对地阻抗：); end B1(i,1)=p;B1(i,2)=q;B1(i,3)=z;B1(i,4)=y;B1(i,5)=1/k; i=i+1; enddisp(由支路参数形成的矩阵B1) B1%求节点导纳矩阵Y=zeros(N0); e=zeros(1,N0); f=zeros(1,N0); for i=1:n1 p=B1(i,1);q=B1(i,2);Y(p,q)=Y(p,q)-1./(B1(i,3)*B1(i,5); Y(q,p)=Y(p,q); Y(q,q)=Y(q,q)+1./(B1(i,3)*B1(i,5)2)+B1(i,4)./2; Y(p,p)=Y(p,p)+1./B1(i,3)+B1(i,4)./2; Enddisp(导纳矩阵 Y=); disp(Y)程序运行结果：
