function  [x, y]= dbphase2xy(db,phase)
          phase=unwrap(phase,175);                %  8/25/02 rab
          x = (10.^(db/20)).*cos( (pi/180)*phase );
          y = (10.^(db/20)).*sin( (pi/180)*phase );