%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
%  delwye.m - performs delta-wye or wye-delta conversion for
%             balanced or unbalanced load.
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
clear; disp(' ');
Bal=input('  Balanced Load? (1=Yes, 2=No) '); disp(' ');
Dir=input('  Conversion Direction? (1=D->Y, 2=Y->D) '); disp(' ');

if Dir==1
   disp('  DELTA to WYE CONVERSION'); disp(' ');
   if Bal==1
      magZab=input('  magZab = ');
      magZbc=magZab; magZca=magZab;
      angZab=input('  angZab = ')*pi/180;
      angZbc=angZab; angZca=angZab;
   elseif Bal==2
      magZab=input('  magZab = ');
      angZab=input('  angZab = ')*pi/180;
      magZbc=input('  magZbc = ');
      angZbc=input('  angZbc = ')*pi/180;
      magZca=input('  magZca = ');
      angZca=input('  angZca = ')*pi/180;
   end
   Zab=magZab*exp(j*angZab); Zbc=magZbc*exp(j*angZbc);
   Zca=magZca*exp(j*angZca); den=Zab+Zbc+Zca;
   Za=Zab*Zca/den; Zb=Zbc*Zab/den; Zc=Zca*Zbc/den;
   disp(' '); disp(' WYE IMPEDANCES in POLAR FORM');
   disp([abs(Za) angle(Za)*180/pi; abs(Zb) angle(Zb)*180/pi; ...
      abs(Zc) angle(Zc)*180/pi]);
elseif Dir==2
   disp('  WYE to DELTA CONVERSION'); disp(' ');
   if Bal==1
      magZa=input('  magZa = ');
      magZb=magZa; magZc=magZa;
      angZa=input('  angZa = ')*pi/180;
      angZb=angZa; angZc=angZa;
   elseif Bal==2
      magZa=input('  magZa = ');
      angZa=input('  angZa = ')*pi/180;
      magZb=input('  magZb = ');
      angZb=input('  angZb = ')*pi/180;
      magZc=input('  magZc = ');
      angZc=input('  angZc = ')*pi/180;
   end
   Za=magZa*exp(j*angZa); Zb=magZb*exp(j*angZb);
   Zc=magZc*exp(j*angZc); num=Za*Zb+Zb*Zc+Zc*Za;
   Zab=num/Zc; Zbc=num/Za; Zca=num/Zb;
   disp(' '); disp(' DELTA IMPEDANCES in POLAR FORM');
   disp([abs(Zab) angle(Zab)*180/pi; ...
         abs(Zbc) angle(Zbc)*180/pi; ...
         abs(Zca) angle(Zca)*180/pi]);
end