Facebook
From Sharp Hornbill, 5 Years ago, written in Matlab M-file.
Embed
Download Paste or View Raw
Hits: 279
  1. re = (4/3)*6371000;
  2. hr = 5:5:100;
  3. ht = 500
  4. lambda = 0.03
  5. f = 10000
  6. r = 10000
  7. fi = pi
  8. p = (2/sqrt(3))*sqrt(re.*(ht+hr)+((r.^2)/4))
  9. ksi = asin((2*re.*r.*(ht-hr))/p.^3)
  10. r1 = (r/2) - (p.*sin(ksi./3))
  11. r2 = r - r1
  12. phi1 = r1/re
  13. phi2 = r2/re
  14. R1 = sqrt((hr.^2)+4*re.*(re+hr).*(sin(phi1./2)).^2)
  15. R2 = sqrt((ht.^2)+4*re.*(re+ht).*(sin(phi2./2)).^2)
  16. Rd = sqrt(((ht-hr).^2)+4*(re+ht).*(re+hr).*(sin((phi1+phi2)/2)).^2)
  17. psi = asin((ht./R1)-(R1./(2*re)))
  18. deltaR = (4*R1.*R2.*(sin(psi)).^2)/(R1+R2+Rd)
  19. eps = 15 + ((i * 0.15)/(2*pi.*f))
  20. arg1 = eps-((cos(psi)).^2);
  21. arg2 = sqrt(arg1);
  22. arg3 = sin(psi);
  23. arg4 = eps.*arg3;
  24. rv = (arg4-arg2)./(arg4+arg2)
  25. rh = (arg3-arg2)./(arg3+arg2);
  26. ro = rv/(exp(i*fi))
  27. deltaPhi = ((2*pi)./lambda).*deltaR
  28. alpha = deltaPhi + fi
  29. F = abs(1+ro.*exp(i*alpha))
  30. %figure(1)
  31. plot(hr, F)
  32. axis( [0 400 0 4])
  33. grid on;
  34. xlabel('Wysokość nadajnika');
  35. ylabel('Współczynnik F');