众文网1.如何用matlab分析潮流计算
原发布者：gzACpN3C7ilwf
南华大学电气工程学院毕业设计前言电力系统潮流计算是研究电力系统稳态运行情况的一种计算，它根据给定的运行条件及系统接线情况确定整个电力系统各部分的运行状态：各母线的电压，各元件中流过的功率，系统的功率损耗等等。在电力系统规划的设计和现有电力系统运行方式的研究中，都需要利用潮流计算来定量地分析比较供电方案或运行方式的合理性。可靠性和经济性。此外，电力系统潮流计算也是计算系统动态稳定和静态稳定的基础。所以潮流计算是研究电力系统的一种很重要和很基础的计算。随着科学技术的发展，电力系统变得越来越复杂，电气工程师掌握一种好的能对电力系统进行仿真的软件是学习和研究的需要。文章简要介绍了MATLAB发展历史、组成和强大的功能，并用简单例子分别就编程和仿真两方面分析了MATIAB软件在电力系统研究中的具体应用。采取等效电路法，能对特殊、复杂地电力系统进行高效仿真研究，因此，掌握编程和仿真是学好MATLAB的基础。与众多专门的电力系统仿真软件相比，MATLAB软件具有易学、功能强大和开放性好，是电力系统仿真研究的有力工具。第1页共60页南华大学电气工程学院毕业设计1电力系统的基本概念电力系统：发电机把机械能转化为电能，电能经变压器和电力线路输送并分配到用户，在那里经电动机、电炉和电灯等设备又将电能转化为机械能、热能和光能等。这些生产、变换、输送、分配、消费电能的发电机、变压器、变换器、电力线
2.MATLAB的电力系统潮流计算程序
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 % 配电网前推回代潮流计算程序 % 使用IEEE 33节点配电系统作为算例，可实现弱环网情况下的潮流计算 countnum=0; BranchData = [ 1 2 0.0922 0.0470; 2 3 0.4930 0.2511; 3 4 0.3660 0.1864; 4 5 0.3811 0.1941; 5 6 0.8190 0.7070; 6 7 0.1872 0.6188; 7 8 0.7114 0.2351; 8 9 1.0300 0.7400; 9 10 1.0440 0.7400; 10 11 0.1966 0.0650; 11 12 0.3744 0.1238; 12 13 1.4680 1.1550; 13 14 0.5416 0.7129; 14 15 0.5910 0.5260; 15 16 0.7463 0.5450; 16 17 1.2890 1.7210; 17 18 0.7320 0.5740; 2 19 0.1640 0.1565; 19 20 1.5042 1.3554; 20 21 0.4095 0.4784; 21 22 0.7089 0.9373; 3 23 0.4512 0.3083; 23 24 0.8980 0.7091; 24 25 0.8960 0.7011; 6 26 0.2030 0.1034; 26 27 0.2842 0.1447; 27 28 1.0590 0.9337; 28 29 0.8042 0.7006; 29 30 0.5075 0.2585; 30 31 0.9744 0.9630; 31 32 0.3105 0.3619; 32 33 0.3410 0.5302； ]； % 支路，阻抗 NodeData = [ 2 100.00 60.00; 3 90.00 40.00; 4 120.00 80.00; 5 60.00 30.00; 6 60.00 20.00; 7 200.00 100.00; 8 200.00 100.00; 9 60.00 20.00; 10 60.00 20.00; 11 45.00 30.00; 12 60.00 35.00; 13 60.00 35.00; 14 120.00 80.00; 15 60.00 10.00; 16 60.00 20.00; 17 60.00 20.00; 18 90.00 40.00; 19 90.00 40.00; 20 90.00 40.00; 21 90.00 40.00; 22 90.00 40.00; 23 90.00 50.00; 24 420.00 200.00; 25 420.00 200.00; 26 60.00 25.00; 27 60.00 25.00; 28 60.00 20.00; 29 120.00 70.00; 30 200.00 600.00; 31 150.00 70.00; 32 210.00 100.00; 33 60.00 40.00； ]； % 节点，负荷 UB = 12.66； % 电压基准 kV SB = 10； % 功率基准 MVA ZB = UB^2/SB； % 阻抗基准 ohm BranchData(:,[3,4]) = BranchData(:,[3,4]) / ZB； % 阻抗标幺化 NodeData(:,[2,3]) = NodeData(:,[2,3]) / SB / 1000；% 功率标幺化 NN = 33； % 节点数 A0 = zeros(NN); for n = 1:NN-1 A0(BranchData(n,1),BranchData(n,2)) = 1; end % 形成 A0 AssociatedMatrix=0; for n=2:NN-1 AssociatedMatrix(n,n)=1; temp=BranchData(n-1,1); AssociatedMatrix(n,1:n-1)=AssociatedMatrix(temp,1:n-1); end A0T = A0'； % 形成 A0 的转置 S = [0;-NodeData(:,2) - i*NodeData(:,3)]； % 形成 S ZL = [0;BranchData(:,3) + i*BranchData(:,4)]； % 形成 ZL V = ones(NN,1); V(1) = 1； % 各个节点电压赋初值 IL(NN,1) = -conj(S(NN) / V(NN)）； % 最末支路电流赋初值 Delta = 1； % 收敛判据赋初值 TempV = V； % 赋初值，用于记忆上次迭代结果 while Delta > 1e-8 countnum=countnum+1; IN = conj(S ./ V)； % 节点注入电流 for n = 1:NN-1 IL(NN-n) = A0(NN-n,NN-n+1:end) * IL(NN-n+1:end) - IN(NN-n); end % 电流回代过程 for n = 2:NN V(n) = A0T(n,1:n-1) * V(1:n-1) - ZL(n) * IL(n); end % 电压前推过程 Delta = max(abs(V-TempV)）； % 更新收敛判据 TempV = V； % 记忆迭代结果 end Vangle(:,1)=abs(V); Vangle(:,2)=angle(V)/3.1415*180; for i=1:NN-1 st=BranchData(i,1); en=BranchData(i,2); Sij(i,1)=V(st)*conj((V(st)-V(en))/ZL(i+1)); Sji(i,1)=V(en)*conj((V(en)-V(st))/ZL(i+1)); end 希望你能用到。
3.MATLAB的电力系统潮流计算程序
% 配电网前推回代潮流计算程序 % 使用IEEE 33节点配电系统作为算例，可实现弱环网情况下的潮流计算 countnum=0; BranchData = [ 1 2 0.0922 0.0470; 2 3 0.4930 0.2511; 3 4 0.3660 0.1864; 4 5 0.3811 0.1941; 5 6 0.8190 0.7070; 6 7 0.1872 0.6188; 7 8 0.7114 0.2351; 8 9 1.0300 0.7400; 9 10 1.0440 0.7400; 10 11 0.1966 0.0650; 11 12 0.3744 0.1238; 12 13 1.4680 1.1550; 13 14 0.5416 0.7129; 14 15 0.5910 0.5260; 15 16 0.7463 0.5450; 16 17 1.2890 1.7210; 17 18 0.7320 0.5740; 2 19 0.1640 0.1565; 19 20 1.5042 1.3554; 20 21 0.4095 0.4784; 21 22 0.7089 0.9373; 3 23 0.4512 0.3083; 23 24 0.8980 0.7091; 24 25 0.8960 0.7011; 6 26 0.2030 0.1034; 26 27 0.2842 0.1447; 27 28 1.0590 0.9337; 28 29 0.8042 0.7006; 29 30 0.5075 0.2585; 30 31 0.9744 0.9630; 31 32 0.3105 0.3619; 32 33 0.3410 0.5302； ]； % 支路，阻抗 NodeData = [ 2 100.00 60.00; 3 90.00 40.00; 4 120.00 80.00; 5 60.00 30.00; 6 60.00 20.00; 7 200.00 100.00; 8 200.00 100.00; 9 60.00 20.00; 10 60.00 20.00; 11 45.00 30.00; 12 60.00 35.00; 13 60.00 35.00; 14 120.00 80.00; 15 60.00 10.00; 16 60.00 20.00; 17 60.00 20.00; 18 90.00 40.00; 19 90.00 40.00; 20 90.00 40.00; 21 90.00 40.00; 22 90.00 40.00; 23 90.00 50.00; 24 420.00 200.00; 25 420.00 200.00; 26 60.00 25.00; 27 60.00 25.00; 28 60.00 20.00; 29 120.00 70.00; 30 200.00 600.00; 31 150.00 70.00; 32 210.00 100.00; 33 60.00 40.00； ]； % 节点，负荷 UB = 12.66； % 电压基准 kV SB = 10； % 功率基准 MVA ZB = UB^2/SB； % 阻抗基准 ohm BranchData(:,[3,4]) = BranchData(:,[3,4]) / ZB； % 阻抗标幺化 NodeData(:,[2,3]) = NodeData(:,[2,3]) / SB / 1000；% 功率标幺化 NN = 33； % 节点数 A0 = zeros(NN); forn = 1:NN-1 A0(BranchData(n,1),BranchData(n,2)) = 1; end % 形成 A0 AssociatedMatrix=0; forn=2:NN-1 AssociatedMatrix(n,n)=1; temp=BranchData(n-1,1); AssociatedMatrix(n,1:n-1)=AssociatedMatrix(temp,1:n-1); end A0T = A0'； % 形成 A0 的转置 S = [0;-NodeData(:,2) - i*NodeData(:,3)]； % 形成 S ZL = [0;BranchData(:,3) + i*BranchData(:,4)]； % 形成 ZL V = ones(NN,1); V(1) = 1； % 各个节点电压赋初值 IL(NN,1) = -conj(S(NN) / V(NN)）； % 最末支路电流赋初值 Delta = 1； % 收敛判据赋初值 TempV = V； % 赋初值，用于记忆上次迭代结果 whileDelta > 1e-8 countnum=countnum+1; IN = conj(S ./ V)； % 节点注入电流 forn = 1:NN-1 IL(NN-n) = A0(NN-n,NN-n+1:end) * IL(NN-n+1:end) - IN(NN-n); end % 电流回代过程 forn = 2:NN V(n) = A0T(n,1:n-1) * V(1:n-1) - ZL(n) * IL(n); end % 电压前推过程 Delta = max(abs(V-TempV)）； % 更新收敛判据 TempV = V； % 记忆迭代结果 end Vangle(:,1)=abs(V); Vangle(:,2)=angle(V)/3.1415*180; fori=1:NN-1 st=BranchData(i,1); en=BranchData(i,2); Sij(i,1)=V(st)*conj((V(st)-V(en))/ZL(i+1)); Sji(i,1)=V(en)*conj((V(en)-V(st))/ZL(i+1)); end 希望你能用到。
4.求帮忙分析电力系统潮流计算的MATLAB程序
% 配电网前推回代潮流计算程序 % 使用IEEE 33节点配电系统作为算例，可实现弱环网情况下的潮流计算countnum=0;BranchData = [ 1 2 0.0922 0.0470; 2 3 0.4930 0.2511; 3 4 0.3660 0.1864; 4 5 0.3811 0.1941; 5 6 0.8190 0.7070; 6 7 0.1872 0.6188; 7 8 0.7114 0.2351; 8 9 1.0300 0.7400; 9 10 1.0440 0.7400; 10 11 0.1966 0.0650; 11 12 0.3744 0.1238; 12 13 1.4680 1.1550; 13 14 0.5416 0.7129; 14 15 0.5910 0.5260; 15 16 0.7463 0.5450; 16 17 1.2890 1.7210; 17 18 0.7320 0.5740; 2 19 0.1640 0.1565; 19 20 1.5042 1.3554; 20 21 0.4095 0.4784; 21 22 0.7089 0.9373; 3 23 0.4512 0.3083; 23 24 0.8980 0.7091; 24 25 0.8960 0.7011; 6 26 0.2030 0.1034; 26 27 0.2842 0.1447; 27 28 1.0590 0.9337; 28 29 0.8042 0.7006; 29 30 0.5075 0.2585; 30 31 0.9744 0.9630; 31 32 0.3105 0.3619; 32 33 0.3410 0.5302；]； % 支路，阻抗NodeData = [ 2 100.00 60.00; 3 90.00 40.00; 4 120.00 80.00; 5 60.00 30.00; 6 60.00 20.00; 7 200.00 100.00; 8 200.00 100.00; 9 60.00 20.00; 10 60.00 20.00; 11 45.00 30.00; 12 60.00 35.00; 13 60.00 35.00; 14 120.00 80.00; 15 60.00 10.00; 16 60.00 20.00; 17 60.00 20.00; 18 90.00 40.00; 19 90.00 40.00; 20 90.00 40.00; 21 90.00 40.00; 22 90.00 40.00; 23 90.00 50.00; 24 420.00 200.00; 25 420.00 200.00; 26 60.00 25.00; 27 60.00 25.00; 28 60.00 20.00; 29 120.00 70.00; 30 200.00 600.00; 31 150.00 70.00; 32 210.00 100.00; 33 60.00 40.00；]； % 节点，负荷UB = 12.66； % 电压基准 kVSB = 10； % 功率基准 MVAZB = UB^2/SB； % 阻抗基准 ohmBranchData(:,[3,4]) = BranchData(:,[3,4]) / ZB； % 阻抗标幺化NodeData(:,[2,3]) = NodeData(:,[2,3]) / SB / 1000；% 功率标幺化NN = 33； % 节点数A0 = zeros(NN);for n = 1:NN-1 A0(BranchData(n,1),BranchData(n,2)) = 1;end % 形成 A0AssociatedMatrix=0;for n=2:NN-1 AssociatedMatrix(n,n)=1; temp=BranchData(n-1,1); AssociatedMatrix(n,1:n-1)=AssociatedMatrix(temp,1:n-1);endA0T = A0'； % 形成 A0 的转置S = [0;-NodeData(:,2) - i*NodeData(:,3)]； % 形成 SZL = [0;BranchData(:,3) + i*BranchData(:,4)]； % 形成 ZLV = ones(NN,1);V(1) = 1； % 各个节点电压赋初值IL(NN,1) = -conj(S(NN) / V(NN)）； % 最末支路电流赋初值Delta = 1； % 收敛判据赋初值TempV = V； % 赋初值，用于记忆上次迭代结果while Delta > 1e-8 countnum=countnum+1; IN = conj(S ./ V)； % 节点注入电流 for n = 1:NN-1 IL(NN-n) = A0(NN-n,NN-n+1:end) * IL(NN-n+1:end) - IN(NN-n); end % 电流回代过程 for n = 2:NN V(n) = A0T(n,1:n-1) * V(1:n-1) - ZL(n) * IL(n); end % 电压前推过程 Delta = max(abs(V-TempV)）； % 更新收敛判据 TempV = V； % 记忆迭代结果endVangle(:,1)=abs(V);Vangle(:,2)=angle(V)/3.1415*180;for i=1:NN-1 st=BranchData(i,1); en=BranchData(i,2); Sij(i,1)=V(st)*conj((V(st)-V(en))/ZL(i+1)); Sji(i,1)=V(en)*conj((V(en)-V(st))/ZL(i+1));end 希望你能用到。
5.求帮忙分析电力系统潮流计算的MATLAB程序
% 配电网前推回代潮流计算程序 % 使用IEEE 33节点配电系统作为算例，可实现弱环网情况下的潮流计算countnum=0;BranchData = [ 1 2 0.0922 0.0470; 2 3 0.4930 0.2511; 3 4 0.3660 0.1864; 4 5 0.3811 0.1941; 5 6 0.8190 0.7070; 6 7 0.1872 0.6188; 7 8 0.7114 0.2351; 8 9 1.0300 0.7400; 9 10 1.0440 0.7400; 10 11 0.1966 0.0650; 11 12 0.3744 0.1238; 12 13 1.4680 1.1550; 13 14 0.5416 0.7129; 14 15 0.5910 0.5260; 15 16 0.7463 0.5450; 16 17 1.2890 1.7210; 17 18 0.7320 0.5740; 2 19 0.1640 0.1565; 19 20 1.5042 1.3554; 20 21 0.4095 0.4784; 21 22 0.7089 0.9373; 3 23 0.4512 0.3083; 23 24 0.8980 0.7091; 24 25 0.8960 0.7011; 6 26 0.2030 0.1034; 26 27 0.2842 0.1447; 27 28 1.0590 0.9337; 28 29 0.8042 0.7006; 29 30 0.5075 0.2585; 30 31 0.9744 0.9630; 31 32 0.3105 0.3619; 32 33 0.3410 0.5302；]； % 支路，阻抗NodeData = [ 2 100.00 60.00; 3 90.00 40.00; 4 120.00 80.00; 5 60.00 30.00; 6 60.00 20.00; 7 200.00 100.00; 8 200.00 100.00; 9 60.00 20.00; 10 60.00 20.00; 11 45.00 30.00; 12 60.00 35.00; 13 60.00 35.00; 14 120.00 80.00; 15 60.00 10.00; 16 60.00 20.00; 17 60.00 20.00; 18 90.00 40.00; 19 90.00 40.00; 20 90.00 40.00; 21 90.00 40.00; 22 90.00 40.00; 23 90.00 50.00; 24 420.00 200.00; 25 420.00 200.00; 26 60.00 25.00; 27 60.00 25.00; 28 60.00 20.00; 29 120.00 70.00; 30 200.00 600.00; 31 150.00 70.00; 32 210.00 100.00; 33 60.00 40.00；]； % 节点，负荷UB = 12.66； % 电压基准 kVSB = 10； % 功率基准 MVAZB = UB^2/SB； % 阻抗基准 ohmBranchData(:,[3,4]) = BranchData(:,[3,4]) / ZB； % 阻抗标幺化NodeData(:,[2,3]) = NodeData(:,[2,3]) / SB / 1000；% 功率标幺化NN = 33； % 节点数A0 = zeros(NN);for n = 1:NN-1 A0(BranchData(n,1),BranchData(n,2)) = 1;end % 形成 A0AssociatedMatrix=0;for n=2:NN-1 AssociatedMatrix(n,n)=1; temp=BranchData(n-1,1); AssociatedMatrix(n,1:n-1)=AssociatedMatrix(temp,1:n-1);endA0T = A0'； % 形成 A0 的转置S = [0;-NodeData(:,2) - i*NodeData(:,3)]； % 形成 SZL = [0;BranchData(:,3) + i*BranchData(:,4)]； % 形成 ZLV = ones(NN,1);V(1) = 1； % 各个节点电压赋初值IL(NN,1) = -conj(S(NN) / V(NN)）； % 最末支路电流赋初值Delta = 1； % 收敛判据赋初值TempV = V； % 赋初值，用于记忆上次迭代结果while Delta > 1e-8 countnum=countnum+1; IN = conj(S ./ V)； % 节点注入电流 for n = 1:NN-1 IL(NN-n) = A0(NN-n,NN-n+1:end) * IL(NN-n+1:end) - IN(NN-n); end % 电流回代过程 for n = 2:NN V(n) = A0T(n,1:n-1) * V(1:n-1) - ZL(n) * IL(n); end % 电压前推过程 Delta = max(abs(V-TempV)）； % 更新收敛判据 TempV = V； % 记忆迭代结果endVangle(:,1)=abs(V);Vangle(:,2)=angle(V)/3.1415*180;for i=1:NN-1 st=BranchData(i,1); en=BranchData(i,2); Sij(i,1)=V(st)*conj((V(st)-V(en))/ZL(i+1)); Sji(i,1)=V(en)*conj((V(en)-V(st))/ZL(i+1));end 希望你能用到。
6.如何用matlab编出电网潮流计算程序
% 配电网前推回代潮流计算程序
% 使用IEEE 33节点配电系统作为算例，可实现弱环网情况下的潮流计算
countnum=0;
BranchData = [
1 2 0.0922 0.0470;
2 3 0.4930 0.2511;
3 4 0.3660 0.1864;
4 5 0.3811 0.1941;
5 6 0.8190 0.7070;
6 7 0.1872 0.6188;
7 8 0.7114 0.2351;
8 9 1.0300 0.7400;
7.电力系统潮流计算的潮流计算的发展趋势
通过几十年的发展，潮流算法日趋成熟。
近几年，对潮流算法的研究仍然是如何改善传统的潮流算法，即高斯-塞德尔法、牛顿法和快速解耦法。牛顿法，由于其在求解非线性潮流方程时采用的是逐次线性化的方法，为了进一步提高算法的收敛性和计算速度，人们考虑采用将泰勒级数的高阶项或非线性项也考虑进来，于是产生了二阶潮流算法。
后来又提出了根据直角坐标形式的潮流方程是一个二次代数方程的特点，提出了采用直角坐标的保留非线性快速潮流算法。对于保留非线性算法典型论文有：1.文献[保留非线性的电力系统概率潮流计算]提出了它在电力系统概率潮流计算中的应用。
该文献提出了一种新的概率潮流计算方法，它保留了潮流方程的非线性，又利用了P-Q解耦方法，因而数学模型精度较高，且保留了P-Q解耦的优点，有利于大电网的随机潮流计算，用提出的方法对一个典型的系统进行了计算，其数值用MonteCarlo随机模拟作了验证，得到了满意的结果。2.文献[基于系统分割的保留非线性的快速P-Q解耦潮流计算法]分析研究了保留非线性的P-Q解耦快速潮流计算法。
该文献提出了一种新的状态估计算法，既保留了量测方程非线性又利用了快速P-Q分解方法，因此数学模型精度高且保留了快速P-Q分解的优点，提高了状态估计的计算精度和速度.采用系统分割方法将大系统分割为多个小系统，分别对每个小系统进行状态估计，然后对各小系统的状态估计结果进行协调，得到整个系统具有同一参考节点的状态估计结果，这样可大大提高状态估计的计算速度，有利于进行大电网的状态估计.在18节点系统上进行的数字仿真实验验证了该方法的有效性。岩本伸一等提出了一种保留非线性的快速潮流计算法，但用的是直角坐标系，因而没法利用P-Q解耦。
为了更有利于大电网的潮流计算，将此原理推广用于P-Q解耦。这样，既利用了保留非线性的快速算法，在迭代中使用常数雅可比矩阵，又保留了P-Q解耦的优点。
对于一些病态系统，应用非线性潮流计算方法往往会造成计算过程的振荡或者不收敛，从数学上讲，非线性的潮流计算方程组本来就是无解的。这样，人们提出来了将潮流方程构造成一个函数，求此函数的最小值问题，称之为非线性规划潮流的计算方法。
优点是原理上保证了计算过程永远不会发散。如果将数学规划原理和牛顿潮流算法有机结合一起就是最优乘子法。
另外，为了优化系统的运行，从所有以上的可行潮流解中挑选出满足一定指标要求的一个最佳方案就是最优潮流问题。最优潮流是一种同时考虑经济性和安全性的电力网络分析优化问题。
OPF 在电力系统的安全运行、经济调度、可靠性分析、能量管理以及电力定价等方面得到了广泛的应用。最优潮流方面的典型论文有：1.文献[电力系统最优潮流新算法的研究]以NCP 方法为基础，提出了一种新的求解最优潮流算法——投影渐近半光滑牛顿型算法。
该文献以NCP方法为基础，提出了一种新的求解OPF算法——投影渐近半光滑牛顿型算法。针对电力系统的特点，本文的研究工作如下： 1.建立了与OPF问题的KKT系统等价的带界约束的半光滑方程系统。
与已有的NCP方法相比，新的模型由于无需考虑界约束对应的对偶变量（乘子变量），降低了问题的维数，从而适用于解大规模的电力系统问题。 2.基于建立的新模型，本文提出了一类新的Newton型算法，该算法一方面保持界约束的相容性，另一方面有较好的全局与局部超线性收敛性，同时，算法结构简单，易于实现。
3.考虑到电力系统固有的弱耦合特性，受传统解耦最优潮流方法的启示，在所提出的新Newton型方法的基础上，本文又设计了一类分解方法。新方法基于解耦——校正的策略实现算法，不仅充分利用了系统的弱耦合特性，同时保证分解算法在理论上的收敛性。
4.根据所提出的两种算法，用标准的IEEE电力测试系统进行数值实验，并与已有的其他方法进行比较。结果显示新算法具有良好的收敛性和计算效果，在电力系统的规划与运行方面将有广阔的应用前景。
2.文献[基于可信域内点法的最优潮流问题研究]介绍了OPF内点法具有收敛性强、多项式时间复杂性等优点，是极具潜力的优秀算法之一。电力系统不断发展，使得OPF算法跻身于极其困难、非凸的大规模非线性规划行列。
可信域和线性搜索方法是保证最优化算法全局收敛性能的两类技术，将内点法和可信域、线性搜索方法有机结合，构造新的优化算法，是数学规划领域的研究热点。此方面的典型文献有：1.文献[电力市场环境下基于最优潮流的输电容量充裕度研究]首先以最优潮流为工具，选取系统中的关键线路作为系统输电容量充裕度的研究对象，从电网运行的安全性、可靠性的角度系统地研究了输电线路稳定限额对输电容量充裕度的影响，指出稳定限额因子与影子价格的乘积可直接反应出稳定限额水平的经济价值，同时也可以较好的指示出系统运行相对安全、经济的稳定限额水平区间。
2.文献[电力市场环境下基于最优潮流的节点实时电价和购电份额研究]为了为配电公司最优购电模型提供价格参考依据，以发电成本最小为目标函数，考虑电力需求价格弹性的影响，建立了实时电价。