diff --git a/lineWithConfig.m b/lineWithConfig.m
index 6ae23e9..2761b36 100644
--- a/lineWithConfig.m
+++ b/lineWithConfig.m
@@ -16,17 +16,26 @@ phase3X=lineZ(entry+5:entry+7,:);
 phase3B2=lineZ(entry+9:entry+11,:);
 phase3Z=phase3R+1j*phase3X;
 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(3,1)=phase3Y(2,1);
 phase3Y(3,2)=phase3Y(3,1);
 phase3Y(2,3)=phase3Y(1,3);
 phase3Y(2,2)=phase3Y(1,1);
 phase3Y(3,3)=phase3Y(1,1);
-%标幺值
-baseS=1;
-baseY=baseS/(baseHighU^2);
-phase3Y=phase3Y/baseY;
 %进行Fortescue变换
 %直接写结果
 a=exp(1j*2*pi/3);
@@ -39,8 +48,10 @@ Tf2p=inv(Tp2f);
 %     1 a^2 a;
 %     1 a a^2];
 % fs3Y=Tp2f*phase3Y*Tf2p;
-Zl=phase3Y(1,1);
-Zm=phase3Y(2,1);
+% Zl=phase3Y(1,1);
+% Zm=phase3Y(2,1);
+Zl=diagEle;
+Zm=offDiagEle;
 fs3Y=3/3*diag([Zl+2*Zm,Zl-Zm,Zl-Zm]);%这里不是1/3，IEEE上的公式错了。
 fs3Y(abs(fs3Y)<1e-5)=0;
 fs3Y=sparse(fs3Y);