matlab fidp,MATLAB潮流计算仿真

《MATLAB潮流计算仿真》由会员分享，可在线阅读，更多相关《MATLAB潮流计算仿真(8页珍藏版)》请在装配图网上搜索。

1、附录A MATLAB程序本程序的功能是用牛顿拉夫逊法进行潮流计算 B1矩阵1、支路首端号；2、末端号；3、支路阻抗；4、支路对地电纳 5、支路的变比；6、支路首端处于K侧为1，1侧为0 B2矩阵1、该节点发电机功率；2、该节点负荷功率；3、节点电压初始值 4、PV节点电压V的给定值；5、节点所接的无功补偿设备的容量 6、节点分类标号clear;n10;请输入节点数n;nl11;请输入支路数nl;isb1;请输入平衡母线节点号isb;pr0.00001;请输入误差精度pr;B11 2 1.755e-24.155e-2i 0.26i 1 0; 1 4 3.159e-27.479e-2i 0.121。

2、5i 1 0; 1 6 3.159e-27.479e-2i 0.1215i 1 0; 2 3 3.68e-30.11135i 0 0.909 1; 4 5 3.68e-30.11135i 0 0.909 1; 4 6 2.808e-26.648e-2i 0.108i 1 0; 6 7 3.0865e-30.0833i 0 0.909 1; 6 8 3.159e-27.479e-2i 0.1215i 1 0; 6 10 2.457e-25.817e-2i 0.0945i 1 0; 8 9 3.0865e-30.0833i 0 0.909 1; 8 10 2.808e-26.648e-2i 0.1。

3、08i 1 0;请输入由支路参数形成的矩阵 B1;B20 0 1.05 1.05 0 1; 0 0 1 0 0 2; 0 0.60.3718i 1 0 0 2; 0 0 1 0 0 2; 0 0.40.247i 1 0 0 2; 0 0 1 0 0 2; 0 0.350.2169i 1 0 0 2; 0 0 1 0 0 2; 0 0.50.3099i 1 0 0 2; 0.8 0 1.05 1.05 0 3;请输入各节点参数形成的矩阵 B2;Yzerosn;ezeros1,n;fzeros1,n;Vzeros1,n;sidazeros1,n;S1zerosnl;修改部分ym0;SB100;UB。

4、220;ym您输入的参数是标么值(若不是则输入一个不为零的数值);if ym0 SB请输入功率基准值SB; UB请输入电压基准值UB; YBSB./UB./UB; BB1B1; BB2B2; for i1nl B1i,3B1i,3*YB; B1i,4B1i,4./YB; end dispB1矩阵B1; dispB1 for i1n B2i,1B2i,1./SB; B2i,2B2i,2./SB; B2i,3B2i,3./UB; B2i,4B2i,4./UB; B2i,5B2i,5./SB; end dispB2矩阵B2; dispB2end -for i1nl 支路数 if B1i,60 左节。

5、点处于低压侧 pB1i,1;qB1i,2; else pB1i,2;qB1i,1; end Yp,qYp,q-1./B1i,3*B1i,5; 非对角元 Yq,pYp,q; Yq,qYq,q1./B1i,3*B1i,52B1i,4./2; 对角元K侧 Yp,pYp,p1./B1i,3B1i,4./2; 对角元1侧 end求导纳矩阵disp导纳矩阵 Y;dispYGrealY;BimagY; 分解出导纳阵的实部和虚部 for i1n 给定各节点初始电压的实部和虚部 eirealB2i,3; fiimagB2i,3; ViB2i,4; PV节点电压给定模值 endfor i1n 给定各节点注入功率 。

6、SiB2i,1-B2i,2; i节点注入功率SG-SL Bi,iBi,iB2i,5; i节点无功补偿量 endPrealS;QimagS;ICT10;IT21;N02*n;NN01;a0;while IT20 IT20;aa1; for i1n if iisb 非平衡节点 Ci0;Di0; for j11n CiCiGi,j1*ej1-Bi,j1*fj1;Gij*ej-Bij*fj DiDiGi,j1*fj1Bi,j1*ej1;Gij*fjBij*ej end P1Ci*eifi*Di;节点功率P计算eiGij*ej-Bij*fjfiGij*fjBij*ej Q1Ci*fi-ei*Di;节点功。

7、率Q计算fiGij*ej-Bij*fj-eiGij*fjBij*ej求P,Q V2ei2fi2; 电压模平方 以下针对非PV节点来求取功率差及Jacobi矩阵元素 if B2i,63 非PV节点 DPPi-P1; 节点有功功率差 DQQi-Q1; 节点无功功率差 以上为除平衡节点外其它节点的功率计算 求取Jacobi矩阵 for j11n if j1isb dP/de-dQ/df X2Bi,j1*ei-Gi,j1*fi; dP/dfdQ/de X3X2; X2dp/df X3dQ/de X4-X1; X1dP/de X4dQ/df p2*i-1;q2*j1-1; Jp,qX3。

8、;Jp,NDQ;mp1; Jm,qX1;Jm,NDP;qq1; Jp,qX4;Jm,qX2; elseif j1i dP/de X2-DiBi,i*ei-Gi,i*fi; dP/df X3DiBi,i*ei-Gi,i*fi; dQ/de X4-CiGi,i*eiBi,i*fi; dQ/df p2*i-1;q2*j1-1;Jp,qX3;Jp,NDQ;扩展列Q mp1; Jm,qX1;qq1;Jp,qX4;Jm,NDP;扩展列P Jm,qX2; end end else 下面是针对PV节点来求取Jacobi矩阵的元素 DPPi-P1; PV节点有功误差 DVVi2-V2; PV节点电压误。

9、差 for j11n if j1isb dP/de X2Bi,j1*ei-Gi,j1*fi; dP/df X50;X60; p2*i-1;q2*j1-1;Jp,qX5;Jp,NDV; mp1; Jm,qX1;Jm,NDP;qq1;Jp,qX6; Jm,qX2; elseif j1i dP/de X2-DiBi,i*ei-Gi,i*fi; dP/df X5-2*ei; X6-2*fi; p2*i-1;q2*j1-1;Jp,qX5;Jp,NDV; mp1; Jm,qX1;Jm,NDP;qq1;Jp,qX6; Jm,qX2; end end end end end 以上为求雅可比矩阵的各个元素 。

10、 for k3N0 N02*n (从第三行开始，第一、二行是平衡节点) k1k1;N1N; NN01 即 N2*n1扩展列P、Q for k2k1N1 扩展列P、Q Jk,k2Jk,k2./Jk,k; 非对角元规格化 end Jk,k1; 对角元规格化 if k3 不是第三行 k4k-1; for k33k4 用k3行从第三行开始到当前行前的k4行消去 for k2k1N1 k3行后各行下三角元素 Jk3,k2Jk3,k2-Jk3,k*Jk,k2;消去运算 end Jk3,k0; end if kN0 break; end for k3k1N0 for k2k1N1 Jk3,。

11、k2Jk3,k2-Jk3,k*Jk,k2;消去运算 end Jk3,k0; end else for k3k1N0 for k2k1N1 Jk3,k2Jk3,k2-Jk3,k*Jk,k2;消去运算 end Jk3,k0; end end end 上面是用线性变换方式将Jacobi矩阵化成单位矩阵 for k32N0-1 Lk1./2; eLeL-Jk,N; 修改节点电压实部 k1k1; fLfL-Jk1,N; 修改节点电压虚部 end 修改节点电压- for k3N0 DETabsJk,N; if DETpr 电压偏差量是否满足要求 IT2IT21; 不满足要求的节点数加1 end end I。

12、CT2aIT2; ICT1ICT11;end用高斯消去法解w-J*Vdisp迭代次数;dispICT1;disp没有达到精度要求的个数;dispICT2;for k1n Vksqrtek2fk2; sidakatanfk./ek*180./pi; Ekekfk*j;end 计算各输出量 disp各节点的实际电压标幺值E为节点号从小到大排列;dispE;EEE*UB;dispEE;disp-;disp各节点的电压大小V为节点号从小到大排列;dispV;VVV*UB;dispVV;disp-;disp各节点的电压相角sida为节点号从小到大排列;dispsida;for p1n Cp0; for 。

13、q1n CpCpconjYp,q*conjEq; end SpEp*Cp;enddisp各节点的功率S为节点号从小到大排列;dispS;disp;SSS*SB;dispSS;disp-;disp各条支路的首端功率Si为顺序同您输入B1时一致;for i1nl pB1i,1;qB1i,2; if B1i,60 Sip,qEp*conjEp*conjB1i,4./2conjEp*B1i,5-conjEq*conj1./B1i,3*B1i,5; SiziSip,q; else Sip,qEp*conjEp*conjB1i,4./2conjEp./B1i,5-conjEq*conj1./B1i,3*B。

14、1i,5; SiziSip,q; end dispSip,q; SSip,qSip,q*SB; ZFS,num2strp,num2strq,,num2strSSip,q; dispZF; dispSSip,q; disp-;enddisp各条支路的末端功率Sj为顺序同您输入B1时一致;for i1nl pB1i,1;qB1i,2; if B1i,60 Sjq,pEq*conjEq*conjB1i,4./2conjEq./B1i,5-conjEp*conj1./B1i,3*B1i,5; SjyiSjq,p; else Sjq,pEq*conjEq*conjB1i,4./2conjEq*B1i,5。

15、-conjEp*conj1./B1i,3*B1i,5; SjyiSjq,p; end dispSjq,p; SSjq,pSjq,p*SB; ZFS,num2strq,num2strp,,num2strSSjq,p; dispZF; dispSSjq,p; disp-;enddisp各条支路的功率损耗DS为顺序同您输入B1时一致;for i1nl pB1i,1;qB1i,2; DSiSip,qSjq,p; dispDSi; DDSiDSi*SB; ZFDS,num2strp,num2strq,,num2strDDSi; dispZF; dispDDSi; disp-;endfigure1;subplot2,2,1;plotV;xlabel节点号;ylabel电压标幺值;grid on;subplot2,2,2;plotsida;xlabel节点号;ylabel电压角度;grid on;subplot2,2,3;barrealS;ylabel节点注入有功;grid on;subplot2,2,4;barSiz;ylabel支路首端无功;grid on; 。