使用 UPFC 计算电力系统网络潮流（Matlab代码实现）

💥💥💞💞欢迎来到本博客❤️❤️💥💥

🏆博主优势：🌞🌞🌞博客内容尽量做到思维缜密，逻辑清晰，为了方便读者。

⛳️座右铭：行百里者，半于九十。

📋📋📋本文目录如下：🎁🎁🎁

目录

💥1 概述

📚2 运行结果

🎉3 参考文献

🌈4 Matlab代码实现

💥1 概述

本文为电力系统网络执行潮流计算，同时集成了统一潮流控制器（UPFC）的功能。UPFC是一种灵活的交流输电系统（FACTS）设备，可以调节输电线路参数，如阻抗和相位角，以控制功率流和电压曲线。

关键步骤：

初始化：该代码使用与发电机设置、负载需求、输电线路和 UPFC 相关的数据初始化电力系统网络。
潮流方程的制定：潮流方程是根据基尔霍夫电流定律（KCL）和基尔霍夫电压定律（KVL）建立的，考虑到UPFC的存在来控制特定区域的潮流。
UPFC控制：UPFC调节功率流和电压曲线，以提高系统稳定性。它动态调整传输线的阻抗和相位角。
迭代求解：采用迭代方法利用UPFC模型求解非线性潮流方程。该代码使用雅可比矩阵迭代更新电压幅度和角度（VM 和 DL），直到实现收敛。
收敛检查：通过比较连续迭代的结果来验证收敛性。潮流计算将继续进行，直到满足指定的收敛标准。

📚2 运行结果


itn =1itn =2itn =3itn =4THE SOLUTION IS ------------------------minsing =0.2198Pmk =-0.4000Qmk =-0.2500VMcr =0.0489DLcrD =-129.6673VMvr =0.9708DLvrD =-6.4230Pbb =-9.6091e-08Pcr =5.1958e-04Pvr =-5.1968e-04Qcr =0.0230Qvr =-0.2837VM =1.05001.04501.01001.05001.01001.05000.99931.01711.00000.99391.00871.03301.00541.01441.00641.00850.99320.98940.98250.98460.98200.98300.98940.97540.98440.96610.99901.00100.97860.96691.0019DLD =0-1.7871-5.5002-4.1493-6.4175-6.3774-6.1804-3.5264-6.4199-8.3436-4.2953-7.8566-5.1138-8.8013-8.8296-8.2868-8.5811-9.4083-9.5479-9.3068-8.8566-8.8489-9.2266-9.3918-9.5838-10.0326-9.4077-5.6022-10.6993-11.6280-6.3644P =0.99190.5830-0.10000.2000-0.44200.2000-0.2280-0.0240-0.0000-0.0580-0.0760-0.11200.0000-0.0620-0.0820-0.0350-0.0900-0.0320-0.0950-0.0220-0.1750-0.0000-0.0320-0.08700.0000-0.03500.00000.0000-0.0240-0.10600.0000Q =-0.03510.34900.15830.25640.00580.1297-0.1090-0.01200.0000-0.0200-0.0160-0.07500.0000-0.0160-0.0250-0.0180-0.0580-0.0090-0.0340-0.0070-0.11200.0000-0.0160-0.06700.0000-0.02300.00000.0000-0.0090-0.0190-0.0000Ploss =0.0579Qloss =0.2191Lj =0.02370.01680.06150.09090.01980.05480.01820.07990.08640.07840.09410.10730.11480.11020.10660.10580.10450.11800.11310.13330.10040.02370.13280.15470.0594t =0.1190>>

部分代码：

JM3;
JM4=zeros(nb-g+1,nb-g+1);
for im4=1:nb-g
for jm4=1:nb-g
JM4(im4,jm4)=JM4(im4,jm4)+J4(im4,jm4);
end
end
JM4;

JM1(ku-1,ku-1)=JM1(ku-1,ku-1)+PkDLk;
JM1(ku-1,nb)=JM1(ku-1,nb)+PkDLCR;
JM1(ku-1,nb+1)=JM1(ku-1,nb+1)+PkDLVR;
JM1(mu-1,mu-1)=JM1(mu-1,mu-1)+PmDLm;
JM1(mu-1,nb)=JM1(mu-1,nb)+PmDLCR;
JM1(nb,ku-1)=JM1(nb,ku-1)+PmkDLk;
JM1(nb,mu-1)=JM1(nb,mu-1)+PmkDLm;
JM1(nb,nb)=JM1(nb,nb)+PmkDLCR;
JM1(nb+1,ku-1)=JM1(nb+1,ku-1)+PBBDLk;
JM1(nb+1,mu-1)=JM1(nb+1,mu-1)+PBBDLm;
JM1(nb+1,nb)=JM1(nb+1,nb)+PBBDLCR;
JM1(nb+1,nb+1)=JM1(nb+1,nb+1)+PBBDLVR;

JM2(ku-1,ku-g)=JM2(ku-1,ku-g)+PkVk;
JM2(ku-1,nb-g+1)=JM2(ku-1,nb-g+1)+PkVCR;
JM2(mu-1,mu-g)=JM2(mu-1,mu-g)+PmVm;
JM2(mu-1,nb-g+1)=JM2(mu-1,nb-g+1)+PmVCR;
JM2(nb,ku-g)=JM2(nb,ku-g)+PmkVk;
JM2(nb,mu-g)=JM2(nb,mu-g)+PmkVm;
JM2(nb,nb-g+1)=JM2(nb,nb-g+1)+PmkVCR;
JM2(nb+1,ku-g)=JM2(nb+1,ku-g)+PBBVk;
JM2(nb+1,mu-g)=JM2(nb+1,mu-g)+PBBVm;
JM2(nb+1,nb-g+1)=JM2(nb+1,nb-g+1)+PBBVCR;

JM3(ku-g,ku-1)=JM3(ku-g,ku-1)+QkDLk;
JM3(ku-g,nb)=JM3(ku-g,nb)+QkDLCR;
JM3(ku-g,nb+1)=JM3(ku-g,nb+1)+QkDLVR;
JM3(mu-g,mu-1)=JM3(mu-g,mu-1)+QmDLm;
JM3(mu-g,nb)=JM3(mu-g,nb)+QmDLCR;
JM3(nb-g+1,ku-1)=JM3(nb-g+1,ku-1)+QmkDLk;
JM3(nb-g+1,mu-1)=JM3(nb-g+1,mu-1)+QmkDLm;
JM3(nb-g+1,nb)=JM3(nb-g+1,nb)+QmkDLCR;

JM4(ku-g,ku-g)=JM4(ku-g,ku-g)+QkVk;
JM4(ku-g,nb-g+1)=JM4(ku-g,nb-g+1)+QkVCR;
JM4(mu-g,mu-g)=JM4(mu-g,mu-g)+QmVm;
JM4(mu-g,nb-g+1)=JM4(mu-g,nb-g+1)+QmVCR;
JM4(nb-g+1,ku-g)=JM4(nb-g+1,ku-g)+QmkVk;
JM4(nb-g+1,mu-g)=JM4(nb-g+1,mu-g)+QmkVm;
JM4(nb-g+1,nb-g+1)=JM4(nb-g+1,nb-g+1)+QmkVCR;