资源描述：
《节点导纳矩阵的建立》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、如上图所示的简单电力系统中，网络各元件参数的标幺值如下：z12=0.10+j0.40y120=y210=j0.01528z13=j0.3,k=1.1z14=0.12+j0.50y140=y410=j0.01920z24=0.08+j0.40y240=y420=j0.01413系统中节点1、2为PQ节点，节点3为PV节点，节点4为平衡节点。节点导纳矩阵的运行程序如下：clcCleardisp('网络各元件参数用标幺值表示！！！');N0=input('请输入节点数:N0=');n1=input('请输入支路数:n1=');l=input('请
2、输入PQ节点的个数=');form=1:lc(m)=input(['请输入第',num2str(m),'个PQ节点的节点号为：']);endt=input('请输入PV节点的个数=');form=1:t c(m)=input(['请输入第',num2str(m),'个PV节点的节点号为：']);endb=input('请输入平衡节点号:b=');%%由支路参数形成矩阵B1disp('各支路连接情况：') i=1;form=1:n1 symsYN p=input(['第',num2str(m),'条支路的起始节点']);
3、q=input(['第',num2str(m),'条支路的终止节点']); mn=input(['第',num2str(m),'条支路是否有变压器（请输入‘Y’或‘N’）']); y=0;k=1; ifmn=='Y'; k=input('请输入变压器变比(标幺值):'); z=input(['请输入第',num2str(m),'条支路的线路阻抗']); else z=input(['请输入第',num2str(m),'条支路的线路阻抗：']); y=input(['请输入第
4、',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);fori=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)
5、; 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)程序运行结果：
