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)