1
0

pv2params.m 541 B

1234567891011121314151617
  1. function pv2params (pv)
  2. display(sprintf('mic\t sens\t dir\t fbr\t az\t el\n'));
  3. for i = 1:size(pv,1)
  4. p = pv(i,1);
  5. v = pv(i,2:4);
  6. s = ( p + norm(v) )/ 2;
  7. d = norm(v) / ( norm(v) + p );
  8. fbr = (p+norm(v)) / (p-norm(v));
  9. fbr_sign = sign(fbr);
  10. fbr_dB = 20 * log10( abs(fbr) );
  11. o = v / norm(v);
  12. az_el_rad = [ atan2( o(2), o(1) ) asin( o(3) ) ];
  13. az_el_deg = 180/pi * az_el_rad;
  14. display(sprintf('%5g\t %5g\t %5g\t %5g\t %5g\t %5g\n', ...
  15. i, s, d, fbr_dB, az_el_deg(1), az_el_deg(2)));
  16. end
  17. end