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-10-13 17:57:34 +08:00
parent d058163923
commit 60434809c4
1 changed files with 26 additions and 18 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

44
run.m
View File

@ -73,24 +73,24 @@ while(k<=kmax && maxD> EPS)
iterPhaseASpotLoadQ=phaseASpotLoadQ+imag(SA); iterPhaseASpotLoadQ=phaseASpotLoadQ+imag(SA);
iterPhaseBSpotLoadQ=phaseBSpotLoadQ+imag(SB); iterPhaseBSpotLoadQ=phaseBSpotLoadQ+imag(SB);
iterPhaseCSpotLoadQ=phaseCSpotLoadQ+imag(SC); iterPhaseCSpotLoadQ=phaseCSpotLoadQ+imag(SC);
%全部转换为负荷电流 % %全部转换为负荷电流
CurpA=conj((iterPhaseASpotLoadP+1j*iterPhaseASpotLoadQ)./VoltpA); % CurpA=conj((iterPhaseASpotLoadP+1j*iterPhaseASpotLoadQ)./VoltpA);
CurpB=conj((iterPhaseBSpotLoadP+1j*iterPhaseBSpotLoadQ)./VoltpB); % CurpB=conj((iterPhaseBSpotLoadP+1j*iterPhaseBSpotLoadQ)./VoltpB);
CurpC=conj((iterPhaseCSpotLoadP+1j*iterPhaseCSpotLoadQ)./VoltpC); % CurpC=conj((iterPhaseCSpotLoadP+1j*iterPhaseCSpotLoadQ)./VoltpC);
%转换为序电流 % %转换为序电流
f012=Tp2f*conj([CurpA';CurpB';CurpC']); % f012=Tp2f*conj([CurpA';CurpB';CurpC']);
%把三序电流分离出来 % %把三序电流分离出来
If0=conj(f012(1,:)'); % If0=conj(f012(1,:)');
If1=conj(f012(2,:)'); % If1=conj(f012(2,:)');
If2=conj(f012(3,:)'); % If2=conj(f012(3,:)');
%试着算一下正序电流 % %试着算一下正序电流
% fsY11*V1; % fsY11*V1;
%形成负荷序电流的测量值 % %形成负荷序电流的测量值
mIf0=If0; % mIf0=If0;
mIf1=If1; % mIf1=If1;
mIf1(3)=-mIf1(2); % mIf1(3)=-mIf1(2);
mIf2=If2; % mIf2=If2;
%计算 % %计算
% fsY11=fsY11+sparse(Balance,Balance,ones(length(Balance),1),busNum,bus % fsY11=fsY11+sparse(Balance,Balance,ones(length(Balance),1),busNum,bus
% Num);%这里要置0，置1，否则是奇异的 % Num);%这里要置0，置1，否则是奇异的
%%做最小二乘法 %%做最小二乘法
@ -104,7 +104,15 @@ while(k<=kmax && maxD> EPS)
%转换为三相电压 %转换为三相电压
VoltpABC=Tp2f\conj([ Vf0'; (Vmf1.*exp(1j*Vaf1))'; Vf2']);%用Tp2f\ 代替Tf2p* VoltpABC=Tp2f\conj([ Vf0'; (Vmf1.*exp(1j*Vaf1))'; Vf2']);%用Tp2f\ 代替Tf2p*
VoltpA=conj(VoltpABC(1,:)'); VoltpA=conj(VoltpABC(1,:)');
CurpA=-conj((iterPhaseASpotLoadP+1j*iterPhaseASpotLoadQ)./VoltpA);
VoltpB=conj(VoltpABC(2,:)');
CurpB=-conj((iterPhaseBSpotLoadP+1j*iterPhaseBSpotLoadQ)./VoltpB);
VoltpC=conj(VoltpABC(3,:)');
CurpC=-conj((iterPhaseCSpotLoadP+1j*iterPhaseCSpotLoadQ)./VoltpC);
f012=Tp2f*conj([CurpA';CurpB';CurpC']);
If0=conj(f012(1,:)');
If1=conj(f012(2,:)');
If2=conj(f012(3,:)');
If0(Balance)=0; If0(Balance)=0;