潮流计算课题.docx

一、系统结构图： 二、网络参数：  = 1 \* Arabic 1、支路参数： 网络类型支路 编号电压 等级 （kV）装机容量 (MW)导线的技术参数（Ω）（Ω） （） 环 网1-2220 10013.6125.567.851-38.321130.552.243-510.2128.874.992-38.5105.428.361-47.579129.651.454-513.84125.312.78 辐 射 网1-2—— ——24——2-3363-4482-5125-6442、节点参数： 节点类型节点编号发电功率（MW）负荷视在功率 环 网10未知(平衡节点)21000（PV节点）3015+9.4i4027+6i5035.5+25.5i 辐 射 网10未知(平衡节点)204+2i306+3.2i403+1.44i504+3.2i602+1.1i三、潮流计算流程图: 四、matlab程序： clear; Un=input(请输入Un:); %输入所需的额定电压 PQ=[ %节点电压 有功 无功 Un 0 0 Un 4 2 Un 6 3.2 Un 3 1.44 Un 4 3.2 Un 2 1.1 ]; FT=[ %首端 末端 4 3 3 2 6 5 5 2 2 1 ]; RX=[ % R X 4 8 3 6 4 4 1 2 2 4 ]; NN=size(PQ,1); %节点数 NB=size(FT,1); %支路数数 V=PQ(:,1); %V初始电压相量 maxd=1 k=1 while maxd0.0001 k=k+1; PQ2=PQ; %每一次迭代各节点的注入有功和无功相同 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,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 单位Mvar end 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); %计算支路电压损耗的横分量dv2 V2(kt)=sqrt((V(kf)-dv1)^2+dv2^2); %计算支路末端电压/kV angle(kt)=angle(kf)+atand(dv2/(V(kf)-dv1)); %计算支路 end maxd=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=