dugg@lab-desk
parent
3231d81b99
commit
e7144d70ae
|
|
@ -0,0 +1,73 @@
|
|||
function createfigure(ymatrix1)
|
||||
%CREATEFIGURE(YMATRIX1)
|
||||
% YMATRIX1: bar matrix data
|
||||
|
||||
% Auto-generated by MATLAB on 11-Feb-2015 20:23:32
|
||||
|
||||
% Create figure
|
||||
figure1 = figure('Name','最大不平衡量','Color',[1 1 1]);
|
||||
|
||||
% Create axes
|
||||
axes1 = axes('Parent',figure1,'XTickLabel',{'情形A','情形B','情形C'},...
|
||||
'XTick',[1 2 3],...
|
||||
'Position',[0.13 0.11 0.29654028436019 0.391061571125265],...
|
||||
'FontName','Times New Roman');
|
||||
hold(axes1,'all');
|
||||
|
||||
% Create multiple lines using matrix input to bar
|
||||
bar1 = bar(ymatrix1,'BarLayout','stacked');
|
||||
set(bar1(1),...
|
||||
'FaceColor',[0.831372559070587 0.815686285495758 0.7843137383461],...
|
||||
'DisplayName','阶段1');
|
||||
set(bar1(2),...
|
||||
'FaceColor',[0.941176474094391 0.941176474094391 0.941176474094391],...
|
||||
'DisplayName','fdsfsd');
|
||||
|
||||
% Create xlabel
|
||||
xlabel('情形','FontSize',12,'FontName','黑体');
|
||||
|
||||
% Create ylabel
|
||||
ylabel('时间/ms','FontSize',12,'FontName','黑体');
|
||||
|
||||
% Create textbox
|
||||
annotation(figure1,'textbox',...
|
||||
[0.166919431279624 0.371781316348195 0.0888625592417062 0.059447983014862],...
|
||||
'String',{'29.68ms'},...
|
||||
'FontName','Times New Roman',...
|
||||
'LineStyle','none');
|
||||
|
||||
% Create textbox
|
||||
annotation(figure1,'textbox',...
|
||||
[0.16706161137441 0.262269639065817 0.0888625592417062 0.059447983014862],...
|
||||
'String',{'70.16ms'},...
|
||||
'FontName','Times New Roman',...
|
||||
'LineStyle','none');
|
||||
|
||||
% Create textbox
|
||||
annotation(figure1,'textbox',...
|
||||
[0.243791469194317 0.353904458598726 0.0888625592417062 0.059447983014862],...
|
||||
'String',{'31.66ms'},...
|
||||
'FontName','Times New Roman',...
|
||||
'LineStyle','none');
|
||||
|
||||
% Create textbox
|
||||
annotation(figure1,'textbox',...
|
||||
[0.243933649289103 0.244392781316348 0.0888625592417062 0.059447983014862],...
|
||||
'String',{'68.48ms'},...
|
||||
'FontName','Times New Roman',...
|
||||
'LineStyle','none');
|
||||
|
||||
% Create textbox
|
||||
annotation(figure1,'textbox',...
|
||||
[0.314881516587683 0.343288747346072 0.0888625592417062 0.059447983014862],...
|
||||
'String',{'30.08ms'},...
|
||||
'FontName','Times New Roman',...
|
||||
'LineStyle','none');
|
||||
|
||||
% Create textbox
|
||||
annotation(figure1,'textbox',...
|
||||
[0.315023696682469 0.233777070063694 0.0888625592417062 0.059447983014862],...
|
||||
'String',{'65.16ms'},...
|
||||
'FontName','Times New Roman',...
|
||||
'LineStyle','none');
|
||||
|
||||
|
|
@ -0,0 +1,40 @@
|
|||
function createfigure(ymatrix1)
|
||||
%CREATEFIGURE(YMATRIX1)
|
||||
% YMATRIX1: bar matrix data
|
||||
|
||||
% Auto-generated by MATLAB on 11-Feb-2015 21:08:10
|
||||
|
||||
% Create figure
|
||||
figure1 = figure('Color',[1 1 1]);
|
||||
|
||||
% Create axes
|
||||
axes1 = axes('Parent',figure1,...
|
||||
'XTick',[0 2 4 6 8 10 12 14 16 18 20 22 24 26 28 30 32 34 35],...
|
||||
'Position',[0.127153700189753 0.112232142857143 0.430721062618596 0.385535714285714],...
|
||||
'FontName','Times New Roman');
|
||||
% Uncomment the following line to preserve the X-limits of the axes
|
||||
xlim(axes1,[1 35]);
|
||||
% Uncomment the following line to preserve the Y-limits of the axes
|
||||
ylim(axes1,[0 0.0021]);
|
||||
hold(axes1,'all');
|
||||
|
||||
% Create multiple lines using matrix input to bar
|
||||
bar1 = bar(ymatrix1,'BarWidth',1,'Parent',axes1);
|
||||
set(bar1(1),'FaceColor',[0 0 0],'DisplayName','电压幅值');
|
||||
set(bar1(2),...
|
||||
'FaceColor',[0.941176474094391 0.941176474094391 0.941176474094391],...
|
||||
'DisplayName','电压相角');
|
||||
|
||||
% Create xlabel
|
||||
xlabel('节点号','FontSize',12,'FontName','黑体');
|
||||
|
||||
% Create ylabel
|
||||
ylabel('最大绝对偏差/(p.u，弧度)','FontSize',12,'FontName','黑体');
|
||||
|
||||
% Create legend
|
||||
legend1 = legend(axes1,'show');
|
||||
set(legend1,...
|
||||
'Position',[0.367488931056294 0.487351190476188 0.115749525616698 0.09375],...
|
||||
'FontSize',12,...
|
||||
'FontName','黑体');
|
||||
|
||||
|
|
@ -0,0 +1,20 @@
|
|||
function [ ] = DrawLoadProfile( )
|
||||
figure();
|
||||
data=dlmread('LoadProfile.txt');
|
||||
hold on;
|
||||
y1=data(:,3);
|
||||
y2=data(:,4);
|
||||
x=1:96;
|
||||
k=1;
|
||||
for i=0:96-1
|
||||
x1(k)=x(length(x)-i);
|
||||
Z2(k)=y2(length(x)-i);
|
||||
k=k+1;
|
||||
end
|
||||
h=patch([x x1],[y1' Z2],[230 230 230]./255)
|
||||
set(h,'EdgeColor','None');
|
||||
plot(1:96,data(:,2),'k');
|
||||
plot(1:96,data(:,3),'--k');
|
||||
plot(1:96,data(:,4),'--k');
|
||||
end
|
||||
|
||||
|
|
@ -0,0 +1,9 @@
|
|||
function MaxBoundErrorFigure()
|
||||
figure();
|
||||
data=dlmread('maxBoundError.txt');
|
||||
% data(1,:)=0;
|
||||
maxV=data(:,1);
|
||||
maxA=data(:,2);
|
||||
% bar([maxV';maxA']');
|
||||
CreateMaxErrorFigure([maxV';maxA']');
|
||||
end
|
||||
|
|
@ -0,0 +1,9 @@
|
|||
function MaxErrorFigure()
|
||||
figure();
|
||||
data=dlmread('maxError.txt');
|
||||
data(1,:)=0;
|
||||
maxV=data(:,1);
|
||||
maxA=data(:,2);
|
||||
% bar([maxV';maxA']');
|
||||
CreateMaxErrorFigure([maxV';maxA']');
|
||||
end
|
||||
8
OPF.m
8
OPF.m
|
|
@ -393,12 +393,14 @@ figure('Name','
|
|||
% maxDismatchPQ = [0.3123e-10 0.1497e-10 0.7351e-10; 0.6854e-10 0.1973e-10 0.5824e-10];
|
||||
% bar(maxDismatchPQ);
|
||||
%% źĆËăĘąźä
|
||||
calTime=[87 61; 211 87; 108 59;];
|
||||
calTime=[70.16 29.68; 68.48 31.661; 65.156 30.08;];
|
||||
bar(calTime,'stacked');
|
||||
% figure();
|
||||
% DeviationFigure(1:33,[CaseAREV;CaseBREV;CaseCREV],[CaseAREA;CaseBREA;CaseCREA],[CaseAREP,CaseBREP,CaseCREP],[CaseAREQ,CaseBREQ,CaseCREQ]);%,[CaseAREA;CaseBREA;CaseCREA],[CaseAREV;CaseBREV;CaseCREV],[CaseAREV;CaseBREV;CaseCREV]);
|
||||
% VoltBar(voltBarCaseAX,voltBarCaseAY,voltBarCaseBX,voltBarCaseBY,voltBarCaseCX,voltBarCaseCY);
|
||||
% AngelBar(angelBarCaseAX,angelBarCaseBY,angelBarCaseCX,angelBarCaseAY,angelBarCaseBX,angelBarCaseCY);
|
||||
% AngelBar(PDBarCaseAX,PDBarCaseAY,PDBarCaseBX,PDBarCaseBY,PDBarCaseCX,PDBarCaseCY);
|
||||
QDBar(QDBarCaseAX,QDBarCaseAY,QDBarCaseBX,QDBarCaseBY,QDBarCaseCX,QDBarCaseCY);
|
||||
|
||||
% QDBar(QDBarCaseAX,QDBarCaseAY,QDBarCaseBX,QDBarCaseBY,QDBarCaseCX,QDBarCaseCY);
|
||||
% MaxErrorFigure()
|
||||
% MaxBoundErrorFigure();
|
||||
DrawLoadProfile();
|
||||
|
|
@ -0,0 +1,49 @@
|
|||
function createfigure(VertexNormals1, YData1, XData1, Vertices1, Faces1, X1, YMatrix1)
|
||||
%CREATEFIGURE(VERTEXNORMALS1,YDATA1,XDATA1,VERTICES1,FACES1,X1,YMATRIX1)
|
||||
% VERTEXNORMALS1: patch vertexnormals
|
||||
% YDATA1: patch ydata
|
||||
% XDATA1: patch xdata
|
||||
% VERTICES1: patch vertices
|
||||
% FACES1: patch faces
|
||||
% X1: vector of x data
|
||||
% YMATRIX1: matrix of y data
|
||||
|
||||
% Auto-generated by MATLAB on 11-Feb-2015 22:03:00
|
||||
|
||||
% Create figure
|
||||
figure1 = figure('Color',[1 1 1]);
|
||||
|
||||
% Create axes
|
||||
axes1 = axes('Parent',figure1,...
|
||||
'XTickLabel',{'6:00','12:00','16:00','24:00','0:00'},...
|
||||
'XTick',[24 48 72 96 1],...
|
||||
'Position',[0.13 0.11 0.373795066413662 0.585744680851064],...
|
||||
'FontName','Times New Roman');
|
||||
% Uncomment the following line to preserve the X-limits of the axes
|
||||
xlim(axes1,[1 96]);
|
||||
hold(axes1,'all');
|
||||
|
||||
% Create patch
|
||||
patch('Parent',axes1,'VertexNormals',VertexNormals1,'YData',YData1,...
|
||||
'XData',XData1,...
|
||||
'Vertices',Vertices1,...
|
||||
'Faces',Faces1,...
|
||||
'FaceColor',[0.901960784313726 0.901960784313726 0.901960784313726],...
|
||||
'EdgeColor','none');
|
||||
|
||||
% Create multiple lines using matrix input to plot
|
||||
plot1 = plot(X1,YMatrix1,'Parent',axes1,'LineStyle','--','Color',[0 0 0]);
|
||||
set(plot1(2),'DisplayName','边界');
|
||||
set(plot1(3),'DisplayName','负荷曲线','LineStyle','-');
|
||||
|
||||
% Create xlabel
|
||||
xlabel('时刻/时','FontSize',12,'FontName','黑体');
|
||||
|
||||
% Create ylabel
|
||||
ylabel('负荷/kW','FontSize',12,'FontName','黑体');
|
||||
|
||||
% Create legend
|
||||
legend1 = legend(axes1,'show');
|
||||
set(legend1,...
|
||||
'Position',[0.155913978494625 0.517021276595745 0.108159392789374 0.168794326241135]);
|
||||
|
||||
|
|
@ -0,0 +1,96 @@
|
|||
57697.32 576.9732 663.51918 490.42722
|
||||
55270.921 552.70921 635.6155915 469.8028285
|
||||
53275.549 532.75549 612.6688135 452.8421665
|
||||
51147.685 511.47685 588.1983775 434.7553225
|
||||
49132.866 491.32866 565.027959 417.629361
|
||||
47594.538 475.94538 547.337187 404.553573
|
||||
45822.028 458.22028 526.953322 389.487238
|
||||
44106.141 441.06141 507.2206215 374.9021985
|
||||
42717.381 427.17381 491.2498815 363.0977385
|
||||
42185.436 421.85436 485.132514 358.576206
|
||||
40721.922 407.21922 468.302103 346.136337
|
||||
40139.269 401.39269 461.6015935 341.1837865
|
||||
38901.576 389.01576 447.368124 330.663396
|
||||
37535.238 375.35238 431.655237 319.049523
|
||||
37234.754 372.34754 428.199671 316.495409
|
||||
36277.407 362.77407 417.1901805 308.3579595
|
||||
35705.188 357.05188 410.609662 303.494098
|
||||
35364.886 353.64886 406.696189 300.601531
|
||||
34644.62 346.4462 398.41313 294.47927
|
||||
34811.66 348.1166 400.33409 295.89911
|
||||
34473.493 344.73493 396.4451695 293.0246905
|
||||
35018.258 350.18258 402.709967 297.655193
|
||||
35782.508 357.82508 411.498842 304.151318
|
||||
35520.4 355.204 408.4846 301.9234
|
||||
37126.94 371.2694 426.95981 315.57899
|
||||
37852.994 378.52994 435.309431 321.750449
|
||||
39096.988 390.96988 449.615362 332.324398
|
||||
40615.023 406.15023 467.0727645 345.2276955
|
||||
42599.93 425.9993 489.899195 362.099405
|
||||
44509.964 445.09964 511.864586 378.334694
|
||||
49173.813 491.73813 565.4988495 417.9774105
|
||||
53864.107 538.64107 619.4372305 457.8449095
|
||||
60260.101 602.60101 692.9911615 512.2108585
|
||||
65195.266 651.95266 749.745559 554.159761
|
||||
71841.983 718.41983 826.1828045 610.6568555
|
||||
77802.345 778.02345 894.7269675 661.3199325
|
||||
88879.392 888.79392 1022.113008 755.474832
|
||||
100781.676 1007.81676 1158.989274 856.644246
|
||||
111517.636 1115.17636 1282.452814 947.899906
|
||||
119678.504 1196.78504 1376.302796 1017.267284
|
||||
126802.705 1268.02705 1458.231108 1077.822993
|
||||
131988.024 1319.88024 1517.862276 1121.898204
|
||||
134124.319 1341.24319 1542.429669 1140.056712
|
||||
137238.631 1372.38631 1578.244257 1166.528364
|
||||
138993.205 1389.93205 1598.421858 1181.442243
|
||||
140492.797 1404.92797 1615.667166 1194.188775
|
||||
141981.796 1419.81796 1632.790654 1206.845266
|
||||
142679.116 1426.79116 1640.809834 1212.772486
|
||||
143682.364 1436.82364 1652.347186 1221.300094
|
||||
144526.194 1445.26194 1662.051231 1228.472649
|
||||
146089.742 1460.89742 1680.032033 1241.762807
|
||||
146519.133 1465.19133 1684.97003 1245.412631
|
||||
144834.953 1448.34953 1665.60196 1231.097101
|
||||
145784.897 1457.84897 1676.526316 1239.171625
|
||||
144503.807 1445.03807 1661.793781 1228.28236
|
||||
145438.77 1454.3877 1672.545855 1236.229545
|
||||
145203.708 1452.03708 1669.842642 1234.231518
|
||||
143780.956 1437.80956 1653.480994 1222.138126
|
||||
142645.898 1426.45898 1640.427827 1212.490133
|
||||
143466.557 1434.66557 1649.865406 1219.465735
|
||||
141922.25 1419.2225 1632.105875 1206.339125
|
||||
145841.261 1458.41261 1677.174502 1239.650719
|
||||
144923.155 1449.23155 1666.616283 1231.846818
|
||||
145430.575 1454.30575 1672.451613 1236.159888
|
||||
145334.611 1453.34611 1671.348027 1235.344194
|
||||
146149.143 1461.49143 1680.715145 1242.267716
|
||||
146482.33 1464.8233 1684.546795 1245.099805
|
||||
147876.05 1478.7605 1700.574575 1256.946425
|
||||
149817.698 1498.17698 1722.903527 1273.450433
|
||||
150415.094 1504.15094 1729.773581 1278.528299
|
||||
154442.119 1544.42119 1776.084369 1312.758012
|
||||
158010.804 1580.10804 1817.124246 1343.091834
|
||||
160774.574 1607.74574 1848.907601 1366.583879
|
||||
162786.258 1627.86258 1872.041967 1383.683193
|
||||
161608.648 1616.08648 1858.499452 1373.673508
|
||||
159454.941 1594.54941 1833.731822 1355.366999
|
||||
158027.837 1580.27837 1817.320126 1343.236615
|
||||
157666.495 1576.66495 1813.164693 1340.165208
|
||||
156460.303 1564.60303 1799.293485 1329.912576
|
||||
155215.52 1552.1552 1784.97848 1319.33192
|
||||
152007.853 1520.07853 1748.09031 1292.066751
|
||||
149345.228 1493.45228 1717.470122 1269.434438
|
||||
147650.189 1476.50189 1697.977174 1255.026607
|
||||
145332.662 1453.32662 1671.325613 1235.327627
|
||||
142143.342 1421.43342 1634.648433 1208.218407
|
||||
136555.706 1365.55706 1570.390619 1160.723501
|
||||
134474.784 1344.74784 1546.460016 1143.035664
|
||||
131788.941 1317.88941 1515.572822 1120.205999
|
||||
125144.813 1251.44813 1439.16535 1063.730911
|
||||
116612.327 1166.12327 1341.041761 991.2047795
|
||||
105557.928 1055.57928 1213.916172 897.242388
|
||||
90041.36 900.4136 1035.47564 765.35156
|
||||
79091.837 790.91837 909.5561255 672.2806145
|
||||
68164.237 681.64237 783.8887255 579.3960145
|
||||
61305.02 613.0502 705.00773 521.09267
|
||||
55158.221 551.58221 634.3195415 468.8448785
|
||||
|
|
@ -0,0 +1,4 @@
|
|||
0.010149012 0.000145532
|
||||
0.010613753 0.001269963
|
||||
0.010686824 0.001667111
|
||||
0.01098127 0.002552894
|
||||
|
|
@ -0,0 +1,33 @@
|
|||
0.001940128 0.001207832
|
||||
0.002036094 0.00080888
|
||||
0.002099042 0.000927452
|
||||
0.002424885 0.000935623
|
||||
0.001893664 0.001181456
|
||||
0.002018829 0.001257517
|
||||
0.00210676 0.000743724
|
||||
0.002013494 0.000936731
|
||||
0.002013003 0.001095999
|
||||
0.001842267 0.000931757
|
||||
0.002001811 0.000859772
|
||||
0.001866919 0.000876312
|
||||
0.001834921 0.00075514
|
||||
0.001937562 0.000923904
|
||||
0.00189722 0.000822962
|
||||
0.001933966 0.000901072
|
||||
0.001920252 0.000994934
|
||||
0.001747249 0.000857745
|
||||
0.002008556 0.000823106
|
||||
0.002093043 0.000753621
|
||||
0.001971762 0.000883769
|
||||
0.002987941 0.00129286
|
||||
0.002436167 0.000997165
|
||||
0.002037594 0.001015502
|
||||
0.002066403 0.000954605
|
||||
0.00276576 0.001295141
|
||||
0.001881581 0.001290146
|
||||
0.00193955 0.001351513
|
||||
0.002069665 0.001408808
|
||||
0.001936768 0.001467117
|
||||
0.001780607 0.001485602
|
||||
0.001861682 0.00143783
|
||||
0.001955221 0.001473342
|
||||
5
subOPF.m
5
subOPF.m
|
|
@ -124,6 +124,7 @@ eps=1;
|
|||
|
||||
% 第一遍，算连续的值
|
||||
fprintf('第1次迭代，算连续量。\n');
|
||||
tic;
|
||||
while(abs(Gap)>Precision*10)
|
||||
if KK>kmax
|
||||
break;
|
||||
|
|
@ -167,6 +168,7 @@ while(abs(Gap)>Precision*10)
|
|||
Gap=(Init_L*Init_Z'-Init_U*Init_W');
|
||||
KK=KK+1;
|
||||
end
|
||||
fprintf('第一遍计算时间:%f\n',toc)
|
||||
% 第二遍，算离散
|
||||
% [~,~,Init_Z,Init_W,Init_L,Init_U,~,~,~,RestraintCount,~,~,~,~,~,~,~,~,~]=OPF_Init(Busnum,Balance,PG,QG,Volt,GenU,GenL,PVi,PGi,PVQU,PVQL,RealPD,RealQD,QD,PD);
|
||||
% Gap=(Init_L*Init_Z'-Init_U*Init_W');
|
||||
|
|
@ -180,7 +182,7 @@ fprintf('\n');
|
|||
% Init_L=1*sparse(ones(1,RestraintCount));
|
||||
% Init_U=1*sparse(ones(1,RestraintCount));
|
||||
% Init_Y=sparse(1,2*Busnum);%与学姐一致
|
||||
|
||||
tic;
|
||||
fprintf('第2次迭代，算离散量。\n');
|
||||
while(abs(Gap)>Precision*10)
|
||||
if KK>kmax
|
||||
|
|
@ -228,6 +230,7 @@ while(abs(Gap)>Precision*10)
|
|||
Gap=(Init_L*Init_Z'-Init_U*Init_W');
|
||||
KK=KK+1;
|
||||
end
|
||||
fprintf('第二遍计算时间:%f\n',toc)
|
||||
%% 计算最大不平衡量
|
||||
AngleIJ=sparse(r,c,UAngel(r)-UAngel(c)-Angle',Busnum,Busnum);
|
||||
dP=PD+diag(Volt)*Y.*cos(AngleIJ)*Volt';
|
||||
|
|
|
|||
Loading…
Reference in New Issue