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-21 22:10:03 +08:00
parent f196750c69
commit ee14c207aa
1 changed files with 17 additions and 6 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

23
lineWithConfig.m
View File

@ -16,17 +16,26 @@ phase3X=lineZ(entry+5:entry+7,:);
phase3B2=lineZ(entry+9:entry+11,:); phase3B2=lineZ(entry+9:entry+11,:);
phase3Z=phase3R+1j*phase3X; phase3Z=phase3R+1j*phase3X;
phase3Y=inv(phase3Z); phase3Y=inv(phase3Z);
%标幺值
baseS=1;
baseY=baseS/(baseHighU^2);
phase3Y=phase3Y/baseY;
diagEle=diag(phase3Y);%对角元素
diagEle=sum(diagEle)/3;
offDiagEle=phase3Y-diag(diag(phase3Y));
offDiagEle=sum(sum(offDiagEle))/6;
%让对角和非对角相等 %让对角和非对角相等
phase3Y(1,1)=diagEle;
phase3Y(1,2)=offDiagEle;
phase3Y(2,1)=offDiagEle;
phase3Y(3,1)=offDiagEle;
phase3Y(1,3)=offDiagEle;
phase3Y(1,3)=phase3Y(1,2); phase3Y(1,3)=phase3Y(1,2);
phase3Y(3,1)=phase3Y(2,1); phase3Y(3,1)=phase3Y(2,1);
phase3Y(3,2)=phase3Y(3,1); phase3Y(3,2)=phase3Y(3,1);
phase3Y(2,3)=phase3Y(1,3); phase3Y(2,3)=phase3Y(1,3);
phase3Y(2,2)=phase3Y(1,1); phase3Y(2,2)=phase3Y(1,1);
phase3Y(3,3)=phase3Y(1,1); phase3Y(3,3)=phase3Y(1,1);
%±êçÛÖµ
baseS=1;
baseY=baseS/(baseHighU^2);
phase3Y=phase3Y/baseY;
%进行Fortescue变换 %进行Fortescue变换
%直接写结果 %直接写结果
a=exp(1j*2*pi/3); a=exp(1j*2*pi/3);
@ -39,8 +48,10 @@ Tf2p=inv(Tp2f);
% 1 a^2 a; % 1 a^2 a;
% 1 a a^2]; % 1 a a^2];
% fs3Y=Tp2f*phase3Y*Tf2p; % fs3Y=Tp2f*phase3Y*Tf2p;
Zl=phase3Y(1,1); % Zl=phase3Y(1,1);
Zm=phase3Y(2,1); % Zm=phase3Y(2,1);
Zl=diagEle;
Zm=offDiagEle;
fs3Y=3/3*diag([Zl+2*Zm,Zl-Zm,Zl-Zm]);%这里不是1/3，IEEE上的公式错了。 fs3Y=3/3*diag([Zl+2*Zm,Zl-Zm,Zl-Zm]);%这里不是1/3，IEEE上的公式错了。
fs3Y(abs(fs3Y)<1e-5)=0; fs3Y(abs(fs3Y)<1e-5)=0;
fs3Y=sparse(fs3Y); fs3Y=sparse(fs3Y);