牛拉潮流程序_第1页
牛拉潮流程序_第2页
牛拉潮流程序_第3页
牛拉潮流程序_第4页
牛拉潮流程序_第5页
已阅读5页,还剩5页未读 继续免费阅读

下载本文档

版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领

文档简介

1、cleaar;clc;%+输入带带有变压压器的支支路矩阵阵中各节节点对应应各变比比% fuuncttionn%=NODDE,BBrannch=OppDF_; %打开矩矩阵 (tesst.mm)文件件Nodee=NOODE;N=Noode(:,11); %节节点号Typee=Noode(:,22); %节点点类型BR=BBrannch %将将支路信信息保存存在BRR中%K=Brrancch(:,6); %支支路变压压器变比比,0代代表没有有变压器器n=leengtth(NN); %节节点数nbr=lenngthh(K); %支路路数Totaal_oof_BBus11=siize(NODDE);%

2、取节点点矩阵的的行和列列Totaal_oof_BBus=Tottal_of_Buss1(11,1)%bbus矩矩阵的行行数即 节点点数Totaal_oof_BBrannch11=siize(Braanchh);%取支路路矩阵的的行和列列Totaal_oof_BBrannch=Tottal_of_Braanchh1(11,1);% 支路路braanchh矩阵行行数即 支路路数Z=zeeross(Tootall_off_Buus1); %将节点点排序重重新存储储节点信信息%定定义为 节点数数的方阵阵formmat shoortb=1; %排排序标志志位pq=00; %PPQ节点点标志位位pv=00

3、; %PPV节点点标志位位ph=00; %平衡节节点标志志位%-按照PPQ,PPV,平平衡节点点的次序序排序各各种节点点%-统计计PQ节节点数 0代表表是pqq节点for a=11:Tootall_off_Buus if NNODEE(a,2) = 0 ZZ(b,:)=NODDE(aa,:); bb=b+1; ppq=ppq+11; endend%-统计计PV节节点数22代表ppv节点点for a=11:Tootall_off_Buus if NNODEE(a,2) = 2 ZZ(b,:)=NODDE(aa,:); bb=b+1; ppv=ppv+11; endend%-统计计平衡节节点数 3

4、代代表平衡衡节点for a=11:Tootall_off_Buus if NNODEE(a,2) = 3 Z(bb,:)=NOODE(a,:); b=bb+1; ph=ph+1; endendZZ2=ZZ;%将节节点进行行重新排排序%mm=zzeroos(nn,1);for i=11:n mm(i,11)=ii;endZ1(:,1)=mmm(:,1);Brannch11=zeeross(nbbr,22);for i=11:nif ZZ(i,1)=Z11(i)for j=11:nbbrif Braanchh(j,1)=Z(i,11) BBrannch11(j,1)=Z1(i);endif BBr

5、annch(j,22)=Z(ii,1) BBrannch11(j,2)=Z1(i);endendelseefor j=11:nbbrif BBrannch(j,11)=Z(ii,1) BBrannch11(j,1)=Z(ii,1);endif BBrannch(j,22)=Z(ii,1) BBrannch11(j,2)=Z(ii,1);endendendendBrannch(:,11)=BBrannch11(:,1);Brannch(:,22)=BBrannch11(:,2);Z(:,1)=Z1(:,11);j=sqqrt(-1);%-矩矩阵已经经完成按按照PQQ,PVV,平衡衡节点的的顺序排

6、排列起来来-YSNOODE=Z; %保存排排序后的的原始节节点数据据%=Y=zeeross(n,n);%求互互导纳%for i=11:nfor t=11:nbbrif (Braanchh(t,1)=i|Brrancch(tt,2)=ii)& Brrancch(tt,6)=00 %非变压压器支路路% YY(Brrancch(tt,1),Brrancch(tt,2)=-1/(Braanchh(t,3)+j*BBrannch(t,44);Y(Brrancch(tt,2),Brrancch(tt,1)=YY(Brrancch(tt,1),Brrancch(tt,2);elseeif (Braanchh

7、(t,1)=i|Brrancch(tt,2)=ii)&Braanchh(t,6)=0 %变压压器支路路% Y(BBrannch(t,11),BBrannch(t,22)=(-11/(jj*Brrancch(tt,4)/Braanchh(t,6); Y(BBrannch(t,22),BBrannch(t,11)=Y(BBrannch(t,11),BBrannch(t,22);endendendend%求自导导纳%for i=11:nfor t=11:nbbrif (Braanchh(t,1)=i|Brrancch(tt,2)=ii)& Brrancch(tt,6)=00 %非变压压器支路路%Y(

8、i,i)=Y(ii,i)+1/(Brrancch(tt,3)+j*Braanchh(t,4)+(11/2)*j*Braanchh(t,5);elseeif BBrannch(t,11)=i&Braanchh(t,6)=0 %变压压器支路路且ii为首节节点% YY(i,i)=Y(ii,i)+1/(j*Braanchh(t,4);elseeif Braanchh(t,2)=i&Brrancch(tt,6)=00%_变压压器支路路且ii为末节节点% YY(i,i)=Y(ii,i)+(11/(jj*Brrancch(tt,4)/(Brrancch(tt,6)*Brrancch(tt,6);endend

9、endendend%若若有并联联电容器器组,则则自导纳纳要加上上并联电电容器的的导纳%for i=11:nif NNODEE(i,13)=00 Y(i,ii)=YY(i,i)+j*NNODEE(i,13)endendYn=leengtth(NN);G=reeal(Y); %实实部,即即电导B=immag(Y); %虚虚部,即即电纳%给给定初始始的电压压值与相相位值%U_fiirstt=Z(:,33); %初初始电压压幅值phasse_ffirsst=ZZ(:,4); %初初始相位位值e=U_firrst.*coos(pphasse_ffirsst);f=U_firrst.*siin(pphas

10、se_ffirsst);%计算算Dellta_P初始始功率量量%P=Z(:,55); %节节点负荷荷有功分分量Q=Z(:,66); %节节点负荷荷无功分分量PG=ZZ(:,7); %发发电机发发出的有有功QG=ZZ(:,8); %发发电机发发出的无无功U0=ZZ(:,9); %节节点电压压都的初初始值Deltta_PP=zeeross(1,n-11);for i=11:n-1for j=11:n DDeltta_PP(i)=Deeltaa_P(i)-e(ii)*(G(ii,j)*e(j)-B(ii,j)*f(j)+f(j)*(G(i,jj)*ff(j)+B(i,jj)*ee(j);endend

11、for i=11:n-1 Dellta_P(ii)=DDeltta_PP(i)-(PP(i)-PGG(i);endDeltta_PP%计计算Deeltaa_Q初初始功率率量%m=0;for i=11:n;if TTypee(i)=22; %计算算PV节节点的个个数 m=m+11; %mm代表ppv节点点个数endendDeltta_QQ=zeeross(1,n-mm-1);for i=11:n-m-11for j=11:n Deeltaa_Q(i)=Dellta_Q(ii)-ff(i)*(GG(i,j)*e(jj)-BB(i,j)*f(jj)+e(ii)*(G(ii,j)*f(j)+B(ii,

12、j)*e(j);endendfor i=11:n-m-11 Deeltaa_Q(i)=Dellta_Q(ii)-(Q(ii)-QQG(ii);endDeltta_QQDeltta_VV=zeeross(1,m);for i=11:mfor j=11:nif TTypee(j)=22 DDeltta_VV(i)=U00(j)2-(e(j)2+ff(i)2);endendendDeltta_VVnum=0;dispp(第,numm2sttr(nnum),次次时的DDeltta总的的失配量量为:) %-进进入循环环体判断断是否满满足条件件-%-先先算出最最大值,作为判判断是否否收敛的的依据-%DEL

13、=Deeltaa_P Dellta_Q; %Deeltaa_P Dellta_Q_%MAX =maax(aabs(DELL); MAXThetta_ffirsst=zzeroos(11,n);U_f=U_ffirsst;Deltta_FF_E11=TThetta_ffirsst(11:n-1) U_ff(1:n-mm-1);Deltta_FF=Deeltaa_F_E1;Deltta_CCor=Dellta_F_EE1; %_DDeltta_tthe Dellta_u_%dispp(第一次次最大失失配量误误差:,nuum2sstr(MAXX)%-循环环判断-%if MMAX1e-0044 % 判

14、断依依据 diisp(-下面开开始下一一次迭代代过程!-)endwhille MMAX1e-0044 nnum=numm+1;%形成雅雅克比矩矩阵%-先先求非对对角元素素-(HH)-%Hik=zerros(n-11,n-1);for i=11:n-1for k=11:n-1if ii=kk ttheiik=TThetta_ffirsst(ii)-TThetta_ffirsst(kk); HHik(i,kk)=-U_ffirsst(ii)*UU_fiirstt(k)*(GG(i,k)*sinn(thheikk)-BB(i,k)*coss(thheikk);endendend%-再求求对角元元素-

15、(HH)-%for i=11:n-1for k=11:nif ii=kk thheikk=Thhetaa_fiirstt(i)-Thhetaa_fiirstt(k); Hiik(ii,i)=Hiik(ii,i)+U_firrst(k)*(G(i,kk)*ssin(theeik)-B(i,kk)*ccos(theeik);endend HHik(i,ii)=UU_fiirstt(i)*Hiik(ii,i);endHik%-先求求非对角角元素-NN-%Nik=zerros(n-11,n-m-11);for i=11:n-1for k=11:n-m-11if ii=kk ttheiik=TThett

16、a_ffirsst(ii)-TThetta_ffirsst(kk); NNik(i,kk)=-U_ffirsst(ii)*UU_fiirstt(k)*(GG(i,k)*coss(thheikk)+BB(i,k)*sinn(thheikk);endendend%-再求求对角元元素-%for i=11:n-m-11for k=11:nif ii=kk thheikk=Thhetaa_fiirstt(i)-Thhetaa_fiirstt(k); Niik(ii,i)=Niik(ii,i)+U_firrst(k)*(G(i,kk)*ccos(theeik)+B(i,kk)*ssin(theeik);e

17、ndend NNik(i,ii)=-U_ffirsst(ii)*NNik(i,ii)-22*U_firrst(i)*U_ffirsst(ii)*GG(i,i);endNik%-先求非非对角元元素-(M)-%Mik=zerros(n-mm-1,n-11);for i=11:n-m-11for k=11:n-1if ii=kk ttheiik=TThetta_ffirsst(ii)-TThetta_ffirsst(kk); Mikk(i,k)=U_ffirsst(ii)*UU_fiirstt(k)*(GG(i,k)*coss(thheikk)+BB(i,k)*sinn(thheikk);enden

18、dend%-再求求对角元元素-%for i=11:n-m-11for k=11:nif ii=kk thheikk=Thhetaa_fiirstt(i)-Thhetaa_fiirstt(k); Miik(ii,i)=Miik(ii,i)+U_firrst(k)*(G(i,kk)*ccos(theeik)+B(i,kk)*ssin(theeik);endend MMik(i,ii)=-U_ffirsst(ii)*MMik(i,ii);endMik%-先求非非对角元元素-(L)-%Lik=zerros(n-mm-1,n-mm-1);for i=11:n-m-11for k=11:n-m-11if

19、ii=kk ttheiik=TThetta_ffirsst(ii)-TThetta_ffirsst(kk); Likk(i,k)=-U_firrst(i)*U_ffirsst(kk)*(G(ii,k)*siin(ttheiik)-B(ii,k)*coos(ttheiik);endendend%-再求求对角元元素-%for i=11:n-m-11for k=11:nif ii=kk thheikk=Thhetaa_fiirstt(i)-Thhetaa_fiirstt(k); Liik(ii,i)=Liik(ii,i)+U_firrst(k)*(G(i,kk)*ssin(theeik)-B(i,k

20、k)*ccos(theeik);endend LLik(i,ii)=-U_ffirsst(ii)*LLik(i,ii)+22*U_firrst(i)*U_ffirsst(ii)*BB(i,i);endLik%-至至此雅可可比矩阵阵已经形形成-%-开开始构造造Deeltaa_f;Dellta_ekacbb=HHik Nikk;Miik LLik;kacbb %雅克比比矩阵%-修正正各个量量,包括e,f,PP,Q,U22(重要要!)-% DEL=DELL;Deltta_FF_E=(-11*innv(kkacbb)*DELL;Deltta_FF=Deeltaa_F_E;Deltta_CCor=Del

21、lta_F+DDeltta_CCor;Thetta_ffirsst(11,1:n-11)=DDeltta_CCor(1,11:n-1); Thetta_ffirsst(11,n)=0;%初始始相角的的修正%Thetta_ffirsst=TThetta_ffirsst;Thetta_ffirsst %修正正后的角角度值%Deltta_CC=Deeltaa_Coor;U_fiirstt(1:n-mm-1,1)=Dellta_C(nn:2*n-mm-2,1); U_fiirstt %修修正后的的电压值值%e=U_firrst.*coos(TThetta_ffirsst);f=U_firrst.*si

22、in(TThetta_ffirsst);e f%-计算算修正DDeltta_PP%Deltta_PP=zeeross(1,n-11);for i=11:n-1for k=11:n DDeltta_PP(i)=Deeltaa_P(i)-e(ii,1)*(GG(i,k)*e(kk,1)-B(i,kk)*ff(k,1)-f(i,11)*(G(ii,k)*f(k,11)+BB(i,k)*e(kk,1);endendfor i=11:n-1 Dellta_P(ii)=DDeltta_PP(i)-(PP(i,1)-PG(i,11);endDeltta_PP%-DDeltta_PP-计算算完成-%-计算算D

23、ellta_Q-%Deltta_QQ=zeeross(1,n-mm-1);for i=11:n-m-11for k=11:n Deeltaa_Q(i)=Dellta_Q(ii)-ff(i)*(GG(i,k)*e(kk)-BB(i,k)*f(kk)+e(ii)*(G(ii,k)*f(k)+B(ii,k)*e(k);endendfor i=11:n-m-11 Deeltaa_Q(i)=Dellta_Q(ii)-(Q(ii)-QQG(ii);endDeltta_QQDEL=Deeltaa_P Dellta_Q;dispp(第 ,nuum2sstr(numm), 次时时的Deeltaa总的失失配量为为

24、:)% DEEL%-继续续判断最最大值 MAXX =mmax(abss(DEEL); Theeta_firrst=Theeta_firrst;end%求求平衡节节点的有有功功率率和无功功功率%Ps0=0;i=n;for t=11:n theeij=Theeta_firrst(i)-Theeta_firrst(t); Ps00=Pss0+UU_fiirstt(t)*(GG(i,t)*coss(thheijj)+BB(i,t)*sinn(thheijj);endPs0=U_ffirsst(ii)*PPs0;Z(i,7)=Ps00;Qs0=0;i=n;for t=11:n theeij=Theeta

25、_firrst(i)-Theeta_firrst(t); Qs00=Qss0+UU_fiirstt(t)*(GG(i,t)*sinn(thheijj)-BB(i,t)*coss(thheijj);endQs0=U_ffirsst(ii)*QQs0;Z(i,8)=Qs00;%计计算PVV节点的的无功功功率%Qv0=zerros(1,nn);for i=nn-m:n-11for t=11:n thheijj=Thhetaa_fiirstt(i)-Thhetaa_fiirstt(t); Qvv0(ii)=QQv0(i)+U_ffirsst(tt)*(G(ii,t)*siin(ttheiij)-B(ii,t)*coos(ttheiij);endQv0(i)=U_ffirsst(ii)*QQv0(i);endfor i=nn-m:n-11 Z(ii,8)=Qvv0(ii);endZj=sqqrt(-1);for i=11:n Anggle(i)=Theeta_firrst(i)*1800/pii; P_iin(ii)=ZZ(i

温馨提示

  • 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

最新文档

评论

0/150

提交评论