%% 画电压，为了写论文用
load('case1V.mat');
CV1=Volt;
load('case2V.mat');
CV2=Volt;
load('case3V.mat');
CV3=Volt;
load('case4V.mat');
CV4=Volt;%负荷曲线偏差大
load('caseR.mat');
CVR=Volt;
load('caseM.mat');
CVM=Volt;
subplot(3,2,1)
hist(CVR)
xlabel('电压/p.u');
ylabel('数量/个');
title('真实值电压分布')
subplot(3,2,2)
hist(CVM)
xlabel('电压/p.u');
ylabel('数量/个');
title('测量值电压分布')
subplot(3,2,3)
hist(CV1)
xlabel('电压/p.u');
ylabel('数量/个');
title('100%可知情形电压分布')
subplot(3,2,4)
hist(CV2)
xlabel('电压/p.u');
ylabel('数量/个');
title('50%可知情形电压分布')
subplot(3,2,5)
hist(CV3)
xlabel('电压/p.u');
ylabel('数量/个');
title('始端覆盖情形电压分布')
subplot(3,2,6)
hist(CV4)
xlabel('电压/p.u');
ylabel('数量/个');
title('典型负荷曲线偏差大情形电压分布')

% figure()
% hold on
% 
% plot(CVR)
% title('真实值')
% 
% plot(CVM)
% title('测量值')
% 
% plot(CV1)
% title('Case 1')
% 
% plot(CV2)
% title('Case 2')
% 
% plot(CV3)
% title('Case 3')