facat
/
distributionnetwork-power2current
1
0
Fork 0
Code Issues Pull Requests Releases Wiki Activity

改了迭代的初值 

Signed-off-by: dugg@lab-desk <dugg@lab-desk>
This commit is contained in:
dugg@lab-desk 2014-12-10 15:14:11 +08:00
parent 0b8015e9d0
commit 7d05bcee3a
3 changed files with 20 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
IPMLoop.m
View File

@ -2,8 +2,8 @@ function [ V1r,V1i,I1r,I1i,isConverged ] = IPMLoop(V1measurement,wV1r,wV1i,I1mea
%把每个序的循环写在这个函数中。其实也就是内点法循环。
V1r=1*ones(busNum,1);
V1i=0*ones(busNum,1);
I1r=0.0*ones(length(Loadi),1);%注入电流，相当于放电机电流，与负荷电流负荷相反
I1i=0.0*ones(length(Loadi),1);%注入电流，相当于放电机电流，与负荷电流负荷相反
I1r=0.1*ones(length(Loadi),1);%注入电流，相当于放电机电流，与负荷电流负荷相反
I1i=0.1*ones(length(Loadi),1);%注入电流，相当于放电机电流，与负荷电流负荷相反
KK=0;
plotGap=zeros(1,60);
%初始化

0
feeder13/data1.txt → feeder13/data1三相不平衡.txt
View File

18
run.m
View File

@ -269,15 +269,33 @@ for II=1:3
if II==1
fprintf('ÕýÐò\n');
wV1r(1:2)=1e10;
wV1i(1:2)=1e10;
wI1r(1:2)=1e10;
wI1i(1:2)=1e10;
tic
[ V1r,V1i,I1r,I1i,isConverged1 ]=IPMLoop(V1measurement,wV1r,wV1i,I1measurement,wI1r,wI1i,BalI1r,BalI1i,busNum,Loadi,fsY11,Balance,1 );
toc
end
if II==2
fprintf('¸ºÐò\n');
wV2r(1:2)=1e10;
wV2i(1:2)=1e10;
wI2r(1:2)=1e10;
wI2i(1:2)=1e10;
tic
[ V2r,V2i,I2r,I2i,isConverged2 ]=IPMLoop(V2measurement,wV2r,wV2i,I2measurement,wI2r,wI2i,BalI2r,BalI2i,busNum,Loadi,fsY22,Balance,0 );
toc
end
if II==3
fprintf('ÁãÐò\n');
wV0r(1:2)=1e10;
wV0i(1:2)=1e10;
wI0r(1:2)=1e10;
wI0i(1:2)=1e10;
tic
[ V0r,V0i,I0r,I0i,isConverged0 ]=IPMLoop(V0measurement,wV0r,wV0i,I0measurement,wI0r,wI0i,BalI0r,BalI0i,busNum,Loadi,fsY00,Balance,0 );
toc
end
end