facat
/
stateestimation-ipm
1
0
Fork 0
Code Issues Pull Requests Releases Wiki Activity

加入了结果曲线的比较 

Signed-off-by: facat <dugg@21cn.com>
This commit is contained in:
facat 2013-04-14 22:24:19 +08:00
parent 059c5c2074
commit 7d4f318b58
3 changed files with 20 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

6
plotAndComparison.asv Normal file
View File

@ -0,0 +1,6 @@
function [ output_args ] = plotAndComparison( rVolt,rVAngel,SEVolt,SEVAngel )
x=1:length(rVolt);
plot(x,rVolt)
end

12
plotAndComparison.m Normal file
View File

@ -0,0 +1,12 @@
function plotAndComparison( rVolt,rVAngel,SEVolt,SEVAngel )
x=1:length(rVolt);
subplot(2,1,1);
plot(x,rVolt,'r-');
hold on
plot(x,SEVolt);
subplot(2,1,2);
plot(x,rVAngel,'r-');
hold on
plot(x,SEVAngel);
end

3
run.m
View File

@ -10,7 +10,7 @@ addpath('.\Powerflow')
% 电压 相角 % 电压 相角
%% %%
%% 开始生成量测量 %% 开始生成量测量
sigma=0.05;% ±ê×¼²î sigma=0.01;% ±ê×¼²î
%% 电压 %% 电压
%电压幅值 %电压幅值
rVolt=Volt; %幅值 rVolt=Volt; %幅值
@ -117,3 +117,4 @@ fprintf('Ŀ
fprintf('相对误差\n'); fprintf('相对误差\n');
(abs(rVolt-double(SEVolt)))./(rVolt) (abs(rVolt-double(SEVolt)))./(rVolt)
MaxDeviation(rVolt,SEVolt,rVAngel,SEVAngel) MaxDeviation(rVolt,SEVolt,rVAngel,SEVAngel)
plotAndComparison( rVolt,rVAngel,SEVolt,SEVAngel )