facat
/
pes2014-debug2-laplace
1
0
Fork 0
Code Issues Pull Requests Releases Wiki Activity

用200PD0和200QD0作为20点叠加了噪音的数据

This commit is contained in:
facat 2013-03-06 22:16:56 +08:00
parent 520d749263
commit 3b2edd87e0
6 changed files with 27 additions and 11 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

BIN
200PD0.mat Normal file
View File

Binary file not shown.

BIN
200QD0.mat Normal file
View File

Binary file not shown.

14
NormalDistributionCurve.m Normal file
View File

@ -0,0 +1,14 @@
%%
% 画正态分布曲线
%%
x=-4:0.1:4;
y=gaussmf(x,[1 0]);
plot(x,y)
hold on
% line([0 0],[ 0 1])
% set(gca,'xtick',[],'xticklabel',[])
% set(gca,'ytick',[],'yticklabel',[])
axis([-4 4 0 1.2])
axis normal
xlabel('gaussmf, P=[2 5]')

12
OPF.m
View File

@ -53,10 +53,10 @@ for II=1:1
%% %%
Precision=Precision/1; Precision=Precision/1;
%% ¼ÓÎó²î %% ¼ÓÎó²î
PD0(Loadi)=PD0(Loadi).*(1+normrnd(0,0.5,length(Loadi),1)); PD0(Loadi)=PD0(Loadi).*(1+normrnd(0,0.05,length(Loadi),1));
QD0(Loadi)=QD0(Loadi).*(1+normrnd(0,0.5,length(Loadi),1)); QD0(Loadi)=QD0(Loadi).*(1+normrnd(0,0.05,length(Loadi),1));
load('0PD0.mat'); load('200PD0.mat');
load('0QD0.mat'); load('200QD0.mat');
% PG0(PGi)=PG0(PGi).*(1+normrnd(0,0.01,length(PGi),1)); % PG0(PGi)=PG0(PGi).*(1+normrnd(0,0.01,length(PGi),1));
% QG0(PVi)=QG0(PVi).*(1+normrnd(0,0.01,length(PVi),1)); % QG0(PVi)=QG0(PVi).*(1+normrnd(0,0.01,length(PVi),1));
%% Á±äѹÆ÷ÈÝÁ¿ %% Á±äѹÆ÷ÈÝÁ¿
@ -132,12 +132,14 @@ PD(Loadi)=thesis.MeanPD();
QD(Loadi)=thesis.MeanQD(); QD(Loadi)=thesis.MeanQD();
PG(Balance)=thesis.MeanPG(); PG(Balance)=thesis.MeanPG();
QG(PVi)=thesis.MeanQG(); QG(PVi)=thesis.MeanQG();
thesis.MaxDeviation(RealPG(Balance),RealQG(PVi),RealPD(Loadi),RealQD(Loadi))
thesis.StatDeviation(RealPG(Balance),RealQG(PVi),RealPD(Loadi),RealQD(Loadi))
thesis.MaxDeviation(PG0(Balance),QG0(PVi),PD0(Loadi),QD0(Loadi)) thesis.MaxDeviation(PG0(Balance),QG0(PVi),PD0(Loadi),QD0(Loadi))
thesis.StatDeviation(PG0(Balance),QG0(PVi),PD0(Loadi),QD0(Loadi)) thesis.StatDeviation(PG0(Balance),QG0(PVi),PD0(Loadi),QD0(Loadi))
% thesis.SqareDeviation(RealPG(Balance(1)),RealQG(PVi(1)),RealPD(Loadi),RealQD(Loadi)) % thesis.SqareDeviation(RealPG(Balance(1)),RealQG(PVi(1)),RealPD(Loadi),RealQD(Loadi))
thesis.SqareDeviation(PG0(Balance(1)),QG0(PVi(1)),PD0(Loadi),QD0(Loadi)) thesis.SqareDeviation(PG0(Balance(1)),QG0(PVi(1)),PD0(Loadi),QD0(Loadi))
thesis.PercentOfPass(PG0(Balance(1)),QG0(PVi(1)),PD0(Loadi),QD0(Loadi)) thesis.PercentOfPass(PG0(Balance(1)),QG0(PVi(1)),PD0(Loadi),QD0(Loadi))
thesis.MaxBranchDeviation(Linei,Linej,Liner,Linex,Lineb,Transfori,Transforj,Transforr,Transforx,Branchi,Branchg,Branchb,Transfork0,Volt0,UAngel0,Volt,UAngel,'E:\ËãÀý\Áø½ð¢ñ926_21671693_2012-09-06\newFIle0.txt',PD0,QD0) thesis.MaxBranchDeviation(Linei,Linej,Liner,Linex,Lineb,Transfori,Transforj,Transforr,Transforx,Branchi,Branchg,Branchb,Transfork0,Volt0,UAngel0,Volt,UAngel,'E:\ËãÀý\Áø½ð¢ñ926_21671693_2012-09-06\newFIle20.txt',PD0,QD0)
thesis.StatBranchDeviation(Linei,Linej,Liner,Linex,Lineb,Transfori,Transforj,Transforr,Transforx,Branchi,Branchg,Branchb,Transfork0,Volt0,UAngel0,Volt,UAngel) thesis.StatBranchDeviation(Linei,Linej,Liner,Linex,Lineb,Transfori,Transforj,Transforr,Transforx,Branchi,Branchg,Branchb,Transfork0,Volt0,UAngel0,Volt,UAngel)
toc toc

4
OPF_Init.asv
View File

@ -20,9 +20,9 @@ wQG=0;
%randInt=randperm(size(Loadi,1)); %randInt=randperm(size(Loadi,1));
%randPDind=randInt(1:10); %randPDind=randInt(1:10);
randPDind=0; randPDind=0;
wPD=1/.5^2*zeros(size(Loadi,1),1); wPD=1/.05^2*zeros(size(Loadi,1),1);
wPD(1:2:end)=0; wPD(1:2:end)=0;
wQD=1/.5^2*ones(size(Loadi,1),1); wQD=1/.05^2*ones(size(Loadi,1),1);
wQD(1:2:end)=0; wQD(1:2:end)=0;
%wD(randPDind)=0;%一些负荷不约束 %wD(randPDind)=0;%一些负荷不约束
%wD(7)=0; %wD(7)=0;

8
OPF_Init.m
View File

@ -20,10 +20,10 @@ wQG=0;
%randInt=randperm(size(Loadi,1)); %randInt=randperm(size(Loadi,1));
%randPDind=randInt(1:10); %randPDind=randInt(1:10);
randPDind=0; randPDind=0;
wPD=1/.5^2*ones(size(Loadi,1),1); wPD=1/.05^2*zeros(size(Loadi,1),1);
% wPD(1:2:end)=0; wPD(1:2:end)=0;
wQD=1/.5^2*ones(size(Loadi,1),1); wQD=1/.05^2*zeros(size(Loadi,1),1);
% wQD(1:2:end)=0; wQD(1:2:end)=0;
%wD(randPDind)=0;%一些负荷不约束 %wD(randPDind)=0;%一些负荷不约束
%wD(7)=0; %wD(7)=0;
% wD(11)=0; % wD(11)=0;