《IEEE33节点matlab潮流程序》由会员分享，可在线阅读，更多相关《IEEE33节点matlab潮流程序（2页珍藏版）》请在装配图网上搜索。
1、Bus=1,0,0;2,100, 60; 3,90,40; 4,120,80; 5,60,30; 6,60,20; 7, 200, 100 ; 8,200,100 ;9, 60, 20;10.60, 20; 11,45,30;12.60, 35;13.60, 35; 14,120,80;15.60, 10;16.60, 20;17.60, 20;18.90, 40;19.90, 40;20.90, 40;21.90, 40;22.90, 40;23.90, 40;24.420, 200;25.420, 200 ;26.60, 25; 27,60 ,25;28.60, 20; 29,120, 7
2、0; 30,200,600 ; 31,150, 70; 32,210, 100 ; 33,60,40;2,2,3,0.4930,0.2511;3,3,4,0.3660,0.1864;4,4,5,0.3811,0.1941;5,5,6,0.8190,0.7070;6,6,7,0.1872,0.6188;7,7,8,0.7144,0.2351;8,8,9,1.0300,0.7400;9,9,10,1.0440,0.7400;Branch=1 ,1 ,2 ,0.0922,0.0407;10,10,11,0.1966,0.065;11,11,12,0.3744,0.1238;12,12,13,1.46
3、80,1.1550;13,13,14,0.5416,0.7129;14,14,15,0.5910,0.5260;15,15,16,0.7463,0.5450;16,16,17,1.2890,1.7210;17,17,18,0.7320,0.5740;18.2, 19,0.1640,0.1565;19,19,20,1.5042,1.3554;20,20,21,0.4095,0.4784;21,21,22,0.7089,0.9373;22.3, 23,0.4512,0.3083;23,23,24,0.8980,0.7091;24,24,25,0.8960,0.7011;25,6, 26,0.203
4、0,0.1034;26,26,27,0.2842,0.1447;27,27,28,1.0590,0.9337;28,28,29,0.8042,0.7006;29,29,30,0.5075,0.2585;30,30,31,0.9744,0.9630;31,31,32,0.3105,0.3619;32,32,33,0.3410,0.5302;busnum,row=size(Bus); branchnum,row=size(Branch); soubus=Branch(:,2);mobus=Branch(:,3);Vbus=ones(busnum,1);Vbus(:,1)=12.66;Vbus1=V
5、bus;Ploss=zeros(busnum,1);Qloss=zeros(busnum,1); e=1; k=0; Branch1=Branch; n=1;while isempty(Branch1) m=1;s,row=size(Branch1);while s0 t=find(Branch1(:,2)= Branch1(s,3); if isempty(t)T1(n,:)= Branch1(s,:); n=n+1; elseT2(m,:)= Branch1(s,:);m=m+1; end; s=s-1; end; Branch1=T2; T2=; end; while e1.0e-05
6、P=zeros(busnum,1);Q=zeros(busnum,1); for s=1:branchnum i=T1(s,2); j=T1(s,3);R=T1(s,4);X=T1(s,5);Pload=Bus(j,2);Qload=Bus(j,3);II=(Pload+P(j)人2+(Qload+Q(j)人2)/(Vbus(j)人2*1000);Ploss(i,j)=II*R;Qloss(i,j)=II*X;P(i,j)=Pload+Ploss(i,j)+P(j);Q(i,j)=Qload+Ploss(i,j)+Q(j);P(i)=P(i)+P(i,j);Q(i)=Q(i)+Q(i,j); end;for s=branchnum:-1:1i=T1(s,2);j=T1(s,3);R=T1(s,4);X=T1(s,5);Vbus(j) = (Vbus(i)-(P(i,j)*R+Q(i,j)*X)/(Vbus(i)*1000)A2+(P(i,j)*X-Q(i,j )*R)/(Vbus(i)*1000)A2;Vbus(j)=sqrt(Vbus(j);end;e=max(abs(Vbus1-Vbus);Vbus1=Vbus;k=k+1;end
IEEE33节点matlab潮流程序
