19 lines
398 B
Matlab
19 lines
398 B
Matlab
txt=dlmread('E:\书籍\论文\配电网状态估计\负荷.txt');
|
|
txt=txt./10000;
|
|
x=1:length(txt);
|
|
plot(x,txt,'k');
|
|
hold on
|
|
plot(x,txt*0.8,'k--');
|
|
plot(x,txt*1.2,'k--');
|
|
legend('Load profile','Lower/Upper bound')
|
|
|
|
ylabel('Load /kW')
|
|
xlabel('Time /h')
|
|
k=1;
|
|
for i=0:length(x)-1
|
|
x1(k)=x(length(x)-i);
|
|
Z2(k)=txt(length(x)-i)*1.2;
|
|
k=k+1;
|
|
end
|
|
fill([x x1],[txt'*0.8 Z2],[240, 240, 240]./255,'edgealpha',0)
|