牛顿拉夫逊法潮流计算matlab程序

所需积分/C币:34 2019-04-30 13:08:52 7KB PDF
2
收藏 收藏
举报

MATLAB实现电力系统分析里的牛顿—拉夫逊法计算潮流分布的过程代码以及注释,根据华科版电力系统分析课本上的过程编程的
%设雅可比矩阵为J,以下语句用来实现雅可比矩阵中对角线上元素HNJL的排列 for m=1: n1 (2*m-12“m-1)=H(mm) JJ(2*m-1,2m)=N(mm) J(2m,2m-1)=J(m,m J(2m,2m)=(m,m); %以下语句用于实现雅可比矩阵非对角线上元素的排列 1:N1 if else H(m, n)=-B(m, n)*E(m)+G(m, n)*F(m) N(m, n)=G(m, n) *E(m)+B(m,n)*F(m) J(m, n=-B(m, n)*F(m)-G(m, n)*E(m); L(m, n )=G(m, n)*F(m)-B(m, n)*E(m) JJ(2m-1,2”n-1)=H(mn) 1(2m-12n)=N(mn); J(2*m,2*n-1)=J(m,n) J(2m,2"n)=L(m,n); end %设由P,Q的改变量组成的8×1矩阵为PQ,由E,F的改变量组成的8×1短阵为dU for m=1: N1 PQ(2*m-1)=dP(m);PQ(2*m)=dQ(m); end dU=inv小J)PQ precision=max(abs(dU)); for n=1: N1 F(n)=F(n)+dU(2*n-1) E(n=E(n)+dU(2*n) for n=1: n1+1 U(n)=E(n)+(F(n); nd k-1. dU=dUPQ, U %计算S(5),也就是题目中的S1,即平衡节点功率 for m=1: n1+1 I(m)=YB(5, m)*U(m) S(5=U(5*sum(conj() %设网络总损耗为Ss,计算输电效率 efficiency for m=1: n1+1 So(m=S(m): P(m)=real(s(m)) end Ss=sum(so efficiency=(abs(P(3)+P(4)+P(2))(P(5)+(P(1))*100 %计算线功率S(m,n),与课本中各元素的相对位置有所不同 for m=1: n1+1 fon=1:N1+1 S(m, n=U(m) "(conj(U(m)-conj(U(n))*conj (-YB(m, n)); end end S

...展开详情
试读 5P 牛顿拉夫逊法潮流计算matlab程序
立即下载 低至0.43元/次 身份认证VIP会员低至7折
抢沙发
一个资源只可评论一次,评论内容不能少于5个字
  • 分享学徒

    成功上传1个资源即可获取
关注 私信 TA的资源
上传资源赚积分or赚钱
    最新推荐
    牛顿拉夫逊法潮流计算matlab程序 34积分/C币 立即下载
    1/5
    牛顿拉夫逊法潮流计算matlab程序第1页

    试读结束, 可继续读1页

    34积分/C币 立即下载 >