facat/nopgqg_lineloss
1
0
Fork 0
Code Issues Pull Requests Releases Wiki Activity

1.用这个参数

Browse Source 
2.电压分布的函数单独提出来

Signed-off-by: dmy@lab <dmy@lab.lab>
This commit is contained in:
dmy@lab
2015-04-21 16:03:11 +08:00
parent f301fe4953
commit b3410ff033
5 changed files with 15 additions and 6 deletions
Download Patch File Download Diff File Expand all files Collapse all files

3
Lineloss.m
View File

@@ -35,5 +35,6 @@ ss=Volt(Branchi)'.^2.*(Branchg)
%% %%
dispTransloss=[Transi Transj real(deltTransS)*100 imag(deltTransS)*100]; dispTransloss=[Transi Transj real(deltTransS)*100 imag(deltTransS)*100];
dispTransloss=sortrows(dispTransloss,-3); dispTransloss=sortrows(dispTransloss,-3);
full(dispTransloss); dispLineloss=full(dispLineloss);
dispTransloss=full(dispTransloss);
end end

1
MainOPF.m
View File

@@ -20,6 +20,7 @@ while 1
% totallossSum=totallossSum+totalloss; % totallossSum=totallossSum+totalloss;
[totalloss,JSM]=OPF('E:/算例/新民Ⅰ906_2729823_2012-09-06/newFIle20.txt'); [totalloss,JSM]=OPF('E:/算例/新民Ⅰ906_2729823_2012-09-06/newFIle20.txt');
totallossSum=totallossSum+totalloss; totallossSum=totallossSum+totalloss;
totallossSum*4
JSM_Sum=JSM_Sum+JSM; JSM_Sum=JSM_Sum+JSM;
loopN=loopN+1; loopN=loopN+1;
% totallossSum=totallossSum*4; % totallossSum=totallossSum*4;

3
OPF.m
View File

@@ -166,7 +166,8 @@ setIneq=0;
fprintf('总的损耗为%f(MW 有名值)\n',full(totalLoss)); fprintf('总的损耗为%f(MW 有名值)\n',full(totalLoss));
fprintf('线损率为 %f\n',full(totalLoss/sum(PG))); fprintf('线损率为 %f\n',full(totalLoss/sum(PG)));
%% 计算各线损 %% 计算各线损
%Lineloss(Linei,Linej,Liner,Linex,Lineb,Transfori,Transforj,Transforr,Transforx,Transfork0,Volt,UAngel); % Lineloss(Linei,Linej,Liner,Linex,Lineb,Transfori,Transforj,Transforr,Transforx,Transfork0,Volt,UAngel);
% Lineloss(Linei,Linej,Liner,Linex,Lineb,Transfori,Transforj,Transforr,Transforx,Branchi,Branchg,Branchb,k0,Volt,Angle)
thesis=thesis.AddPDQDPGQG(PD(Loadi),QD(Loadi),PG(Balance),QG(PVi)); thesis=thesis.AddPDQDPGQG(PD(Loadi),QD(Loadi),PG(Balance),QG(PVi));
end end

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=zeros(size(Loadi,1),1); wPD=ones(size(Loadi,1),1);
wPD(1:2:end)=0; % wPD(1:2:end)=0;
wQD=zeros(size(Loadi,1),1); wQD=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;
% wD(11)=0; % wD(11)=0;

6
VoltDistribution.m Normal file
View File

@@ -0,0 +1,6 @@
function [ output_args ] = VoltDistribution( CV4 )
hist(CV4)
xlabel('电压/p.u');
ylabel('数量/个');
end