stateestimateionyalmip-lu9-.../plotError.m

40 lines
720 B
Matlab

function [ output_args ] = plotError( SEVolt,SEVAngle,rVolt,rVAngle,mVolt,PD,PD0,rPD,QD,QD0,rQD,PDi )
subplot(4,1,1);
x=1:length(SEVolt);
plot(x,SEVolt);
hold on
plot(x,rVolt,'r-');
plot(x,mVolt,'g-');
title('电压')
legend('SE','Actual','Measurement')
%画负荷的标志
% t=sparse(length(rVolt),1);
% t(PDi)=.01;
% bar(1:length(t),t);
subplot(4,1,2);
plot(x,SEVAngle);
hold on
plot(x,rVAngle,'r-');
title('相角')
legend('SE','Actual')
x=1:length(PD);
%% PD
subplot(4,1,3);
plot(x,PD);
hold on
plot(x,rPD,'r-');
plot(x,PD0,'g-');
legend('SE','Actual','Measurement')
title('有功负荷')
%% QD
subplot(4,1,4);
plot(x,QD);
hold on
plot(x,rQD,'r-');
plot(x,QD0,'g-');
legend('SE','Actual','Measurement')
title('无功负荷')
end