《潮流计算报告》由会员分享，可在线阅读，更多相关《潮流计算报告（33页珍藏版）》请在金锄头文库上搜索。
1、一、系统结构图：二、网络参数：1、支路参数：网络类型支路编号电压等级（kV）装机容量(MW)导线的技术参数（）（）（）环网1-222010013.6125.567.851-38.321130.552.243-510.2128.874.992-38.5105.428.361-47.579129.651.454-513.84125.312.78辐射网1-2242-3363-4482-5125-6442、节点参数：节点类型节点编号发电功率（MW）负荷视在功率环网10未知(平衡节点)21000（PV节点）3015+9.4i4027+6i5035.5+25.5i辐射网10未知(平衡节点)204+2i30
2、6+3.2i403+1.44i504+3.2i602+1.1i三、潮流计算流程图:四、matlab程序：clear;Un=input(请输入Un:); %输入所需的额定电压PQ=%节点电压 有功 无功Un 0 0Un 4 2 Un 6 3.2Un 3 1.44Un 4 3.2Un 2 1.1 ;FT=%首端 末端 4 33 26 55 22 1;RX=% R X 4 83 64 41 2 2 4; NN=size(PQ,1); %节点数NB=size(FT,1); %支路数数V=PQ(:,1); %V初始电压相量maxd=1k=1while maxd0.0001 k=k+1; PQ2=PQ;
3、%每一次迭代各节点的注入有功和无功相同 PL=0.0;for i=1:NB kf=FT(i,1); %前推始节点号 kt=FT(i,2); %前推终节点号 x=(PQ2(kf,2)2+PQ2(kf,3)2)/V(kf)/V(kf);%计算沿线电流 /平方A losss(i,1)=RX(i,1)*x; %计算线路有功损耗 /MW losss(i,2)=RX(i,2)*x; %计算线路无功损耗/MW PQ1(i,1)=PQ2(kf,2)+RX(i,1)*x; %计算支路首端有功/MW RX(i,1)*R PQ1(i,2)=PQ2(kf,3)+RX(i,2)*x; %计算沿支路的无功/MW RX(i
4、,2)*X PQ2(kt,2)= PQ2(kt,2)+PQ1(i,1); %用PQ1去修正支路末端节点的有功P 单位MW PQ2(kt,3)= PQ2(kt,3)+PQ1(i,2); %用PQ1去修正支路末端节点的有功Q 单位Mvarend angle(1)=0.0;for i=NB:-1:1 kf=FT(i,2); %回代始节点号 kt=FT(i,1); %回代终节点号 dv1=(PQ1(i,1)*RX(i,1)+PQ1(i,2)*RX(i,2)/V(kf); %计算支路电压损耗的纵分量dv1 dv2=(PQ1(i,1)*RX(i,2)-PQ1(i,2)*RX(i,1)/V(kf); %计算
5、支路电压损耗的横分量dv2 V2(kt)=sqrt(V(kf)-dv1)2+dv22); %计算支路末端电压/kV angle(kt)=angle(kf)+atand(dv2/(V(kf)-dv1); %计算支路endmaxd=abs(V2(2)-V(2);V2(1)=V(1);for i=3:1:NN if abs(V2(i)-V(i)maxd; maxd=abs(V2(i)-V(i); end end fullloss(1,1)=0;%计算线路总损耗 fullloss(1,2)=0; finalPQ=max(PQ1); for i=1:NB fullloss(1,1)=fullloss(1
6、,1)+losss(i,1); fullloss(1,2)=fullloss(1,2)+losss(i,2); enddisp(辐射网迭代次数：)kdisp(辐射网系统电压差精度：) maxddisp(辐射网系统末端节点有功和无功：)finalPQ %潮流分布 即支路首端潮流MVAdisp(辐射网系统总功率损耗：)fullloss %线路总损耗MVAdisp(辐射网系统各支路功率损耗：)losss %各支路损耗MVAdisp(辐射网系统各节点电压幅值：)V=V2 %节点电压模计算结果kVdisp(辐射网系统各节点电压相角：)angle %节点电压角度计算结果单位度endclcdisp(辐射网迭
7、代次数：)kdisp(辐射网系统电压差精度：) maxddisp(辐射网系统末端节点有功和无功/MVA：)FinPQ=finalPQ(1,1)+finalPQ(1,2)*j %潮流分布 即支路首端潮流MVAdisp(辐射网系统总功率损耗/MVA：)Fulloss=fullloss (1,1)+fullloss(1,2)*j %线路总损耗MVAdisp(辐射网系统各支路功率损耗/MVA：)for(a=1:5)LOSS=losss (a,1)+losss(a,2)*j %各支路损耗MVAenddisp(辐射网系统各节点电压幅值/KV：)V=V2 %节点电压模计算结果kVdisp(辐射网系统各节点电
8、压相角：)angle %节点电压角度计算结果单位度n=5; %input(节点数);nl=6; %input(支路数);isb=1; %input(平衡母线节点号);pr=0.000001; %input(误差精度：pr=);B1=1,2,13.6+125.5i,0.00006785i,1,0;1,3,8.321+130.5i,0.00005224i,1,0;3,5,10.2+128.8i,0.00007499i,1,0;2,3,8.5+105.4i,0.00002836i,1,0;1,4,7.579+129.6i,0.00005145i,1,0;4,5,13.84+125.31i,0.000
9、0278i,1,0; %input(由支路参数形成的矩阵);B2=-FinPQ,0,Un,0,0,1;100,0,Un,Un,0,3;0,15+9.4i,Un,0,0,2;0,27+6i,Un,0,0,2;0,35.5+25.5i,Un,0,0,2; %input(各节点参数形成的矩阵);Y=zeros(n);e=zeros(1,n);f=zeros(1,n);V=zeros(1,n);sida=zeros(1,n);S1=zeros(nl);%对各矩阵置零%修改部分ym=1;SB=100;UB=Un; %定义视在功率和电压基值if ym=0 %若不是标幺值 YB=SB./UB./UB; %定
10、义导纳标幺值 BB1=B1; BB2=B2; for i=1:nl B1(i,3)=B1(i,3)*YB; %切换为阻抗标幺值 B1(i,4)=B1(i,4)./YB; %切换为导纳标幺值 end disp(支路矩阵B1=); sparseB1=sparse(B1); disp(sparseB1) %输出标幺值稀疏矩阵B1 disp(-); for i=1:n B2(i,1)=B2(i,1)./SB; %切换为视在功率标幺值 B2(i,2)=B2(i,2)./SB; %切换为视在功率标幺值 B2(i,3)=B2(i,3)./UB; %切换为电压标幺值 B2(i,4)=B2(i,4)./UB; %切换为电压标幺值 B2(i,5)=B2(i,5)./SB; %切换为视在功率标幺值 end disp(节点矩阵B2=); sparseB2=sparse(B2); disp(sparseB2) %输出标幺值稀疏矩阵B2enddisp(-);% % %-for i=1:nl %支路数 if B1(i,6)=0 %左节点处于1侧 p=B1(i,1);q=B1(i,2); else p=B1(i,2);q=B1(i,1); %使左节点处于低压侧 end Y(p,q)=Y(p,q)-1./(B1(i,3)*B1(i,5); %求解非对角元导纳 Y(q,p)=Y(p,q); %对角元两侧对称
