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

引入符号分解，加快计算速度。 

Signed-off-by: facat <facat@ipso.laptop>
This commit is contained in:
facat 2014-04-23 00:39:13 +08:00
parent 732b61b2a9
commit 8163082770
2 changed files with 12 additions and 10 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

7
dataRead.m
View File

@ -2,9 +2,10 @@ function [ fsY0, fsY1, fsY2,phaseASpotLoadP,phaseBSpotLoadP,phaseCSpotLoadP ...
phaseASpotLoadQ,phaseBSpotLoadQ,phaseCSpotLoadQ,setIJ,nodeNum,Balance,phaseABCY] = dataRead(lineZ,dataFile ) phaseASpotLoadQ,phaseBSpotLoadQ,phaseCSpotLoadQ,setIJ,nodeNum,Balance,phaseABCY] = dataRead(lineZ,dataFile )
data=dlmread(dataFile); data=dlmread(dataFile);
baseHighU=data(1,2); baseHighU=data(1,2);
fsY0=0; fsY0=sparse(0);
fsY1=0; fsY1=sparse(0);
fsY2=0; fsY2=sparse(0);
zeroEntries=find(data(:,1)==0); zeroEntries=find(data(:,1)==0);
lines=data(zeroEntries(1)+1:zeroEntries(2)-1,:); lines=data(zeroEntries(1)+1:zeroEntries(2)-1,:);
[setIJ,nodeNum]=numberNode(lines); [setIJ,nodeNum]=numberNode(lines);

15
run.m
View File

@ -12,11 +12,6 @@ Tf2p=inv(Tp2f);
fsY1amp=abs(fsY1); fsY1amp=abs(fsY1);
[r,c,fsY1ang]=find(fsY1); [r,c,fsY1ang]=find(fsY1);
fsY1ang=angle(fsY1ang); fsY1ang=angle(fsY1ang);
%%²Î¿¼µçѹ
refV0=1;
refV1=0;
refV2=0;
%refV1=1;
Pabc=phaseASpotLoadP+phaseBSpotLoadP+phaseCSpotLoadP; Pabc=phaseASpotLoadP+phaseBSpotLoadP+phaseCSpotLoadP;
Qabc=phaseASpotLoadQ+phaseBSpotLoadQ+phaseCSpotLoadQ; Qabc=phaseASpotLoadQ+phaseBSpotLoadQ+phaseCSpotLoadQ;
busNum=length(phaseASpotLoadP); busNum=length(phaseASpotLoadP);
@ -51,6 +46,9 @@ fsY2=fsY2+sparse(Balance,Balance,ones(length(Balance),1),busNum,busNum);
fsY0(Balance,:)=0; fsY0(Balance,:)=0;
fsY0(:,Balance)=0; fsY0(:,Balance)=0;
fsY0=fsY0+sparse(Balance,Balance,ones(length(Balance),1),busNum,busNum); fsY0=fsY0+sparse(Balance,Balance,ones(length(Balance),1),busNum,busNum);
%%LU分解
[fsY0L,fsY0U,fsY0P,fsY0Q,fsY0R]=lu(fsY0);
[fsY2L,fsY2U,fsY2P,fsY2Q,fsY2R]=lu(fsY2);
tic tic
while(k<=kmax && maxD> EPS) while(k<=kmax && maxD> EPS)
k=k+1; k=k+1;
@ -75,8 +73,11 @@ while(k<=kmax && maxD> EPS)
If2=conj(f012(3,:)'); If2=conj(f012(3,:)');
If0(Balance)=0; If0(Balance)=0;
If2(Balance)=0; If2(Balance)=0;
Vf0=fsY0\If0; %Vf0=fsY0\If0;
Vf2=fsY2\If2; %[fsY0L,fsY0U,fsY0P,fsY0Q,fsY0R]
Vf0=fsY0Q*(fsY0U\(fsY0L\(fsY0P*(fsY0R\If0))));
%Vf2=fsY2\If2;
Vf2=fsY2Q*(fsY2U\(fsY2L\(fsY2P*(fsY2R\If2))));
end end
fprintf('Fortiscue·¨¼ÆËãʱ¼ä %f\n',toc); fprintf('Fortiscue·¨¼ÆËãʱ¼ä %f\n',toc);
Vf1=Vmf1.*exp(1j*Vaf1); Vf1=Vmf1.*exp(1j*Vaf1);