含风电场的电力系统潮流计算程序(共6页).doc

精选优质文档-倾情为你奉上%本程序的功能是用牛拉法进行含风电场的电力系统潮流计算 function s1=pf(a,B1,B2) n=a(1);%节点数 nl=a(2);%支路数 isb=a(3);%平衡节点号 pr=a(4);%误差精度for i=1:n for j=1:n G(i,j)=0; B(i,j)=0; endend%求导纳矩阵%B1为支路参数矩阵，其每一行格式为首节点号，末节点号，支路电导，支路电纳，首节点对地电纳，末节点对地电纳for i=1:nl p=B1(i,1);q=B1(i,2) G(p,q)=G(p,q)-B1(i,3); B(p,q)=B(p,q)- B1(i,4); G(q,p)=G(p,q); B(q,p)=B(p,q); G(p,p)=G(p,p)+ B1(i,3); B(p,p)=B(p,p)+B1(i,4); G(q,q)=G(q,q)+ B1(i,3);