clutter_rcs.m
上传用户:szahd2008
上传日期:2020-09-25
资源大小:1275k
文件大小:2k
源码类别:

传真(Fax)编程

开发平台:

Matlab

  1. function [sigmaC,CNR] = clutter_rcs(sigma0, thetaE, thetaA, SL, range, hr, ht, pt, f0,...
  2.     b,  t0, f, l,ant_id)
  3. % This function calculates the clutter RCS and the CNR for a ground based radar.
  4. clight = 3.e8; % speed of light in meters per second
  5. lambda = clight /f0;
  6. thetaA_deg = thetaA;
  7. thetaE_deg = thetaE;
  8. thetaA = thetaA_deg * pi /180; % antenna azimuth beamwidth in radians
  9. thetaE = thetaE_deg * pi /180.; % antenna elevation beamwidth in radians
  10. re = 6371000; % earth radius in meter
  11. rh = sqrt(8.0*hr*re/3.); % range to horizon in meters
  12. SLv = 10.0^(SL/10); % radar rms sidelobes in volts
  13. sigma0v = 10.0^(sigma0/10); % clutter backscatter coefficient 
  14. tau = 1/b; % pulse width  
  15. deltar = clight * tau / 2.; % range resolution for unmodulated pulse
  16. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  17. range_m = 1000 .* range;  % range in meters
  18. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  19. thetar = asin(hr ./ range_m);
  20. thetae = asin((ht-hr) ./ range_m);
  21. propag_atten = 1. + ((range_m ./ rh).^4); % propagation attenuation due to round earth
  22. Rg = range_m .* cos(thetar);
  23. deltaRg = deltar .* cos(thetar);
  24. theta_sum = thetae + thetar;
  25. % use sinc^2 antenna pattern when ant_id=1
  26. % use Gaussian antenna pattern when ant_id=2
  27. if(ant_id ==1) % use sinc^2 antenna pattern
  28.     ant_arg = (2.78 * theta_sum ) ./ (pi*thetaE);
  29.     gain = (sinc(ant_arg)).^2;
  30. else
  31.     gain = exp(-2.776 .*(theta_sum./thetaE).^2);
  32. end
  33. % compute sigmac
  34. sigmac = (sigma0v .* Rg .* deltaRg) .* (pi * SLv * SLv + thetaA .* gain.^2) ./ propag_atten;
  35. sigmaC = 10*log10(sigmac);
  36. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  37. figure(1)
  38. plot(range, sigmaC)
  39. grid
  40. xlabel('Slant Range in Km')
  41. ylabel('Clutter RCS in dBsm')
  42. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  43. % Calculate CNR
  44. pt = pt * 1000;
  45. g = 26000 / (thetaA_deg*thetaE_deg); % antenna gain
  46. F = 10.^(f/10); % noise figure is 6 dB
  47. Lt = 10.^(l/10); % total radar losses 13 dB
  48. k = 1.38e-23; % Boltzman抯 constant
  49. T0 = t0; % noise temperature 290K
  50. argnumC = 10*log10(pt*g*g*lambda*lambda*tau .* sigmac);
  51. argdem = 10*log10(((4*pi)^3)*k*T0*Lt*F .*(range_m).^4);
  52. CNR = argnumC - argdem;
  53. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  54. figure(2)
  55. plot(range, CNR,'r')
  56. grid
  57. xlabel('Slant Range in Km')
  58. ylabel('CNR in dB')