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

传真(Fax)编程

开发平台:

Matlab

  1. % Use this program to reproduce Fig. 13.12aa
  2. %This program computes the back-scattered RCS for an ellipsoid.
  3. % The angle phi is fixed to three values 0, 45, and 90 degrees
  4. % The angle theta is varied from 0-180 deg.
  5. % A plot of RCS versus theta is generated
  6. % Last modified on July 16, 2003
  7. clear all;
  8. % ===   Input parameters   ===
  9. a = .15;            % 15 cm
  10. b = .20;            % 20 cm
  11. c = .95 ;           % 95 cm
  12. % ===   End of Input parameters   ===
  13. as = num2str(a);
  14. bs = num2str(b);
  15. cs = num2str(c);
  16. eps = 0.00001;
  17. dtr = pi/180;
  18. for q = 1:3
  19.     if q == 1 
  20.         phir = 0;       % the first value of the angle phi
  21.     elseif q == 2
  22.         phir = pi/4;    % the second value of the angle phi
  23.     elseif q == 3
  24.         phir = pi/2;    % the third value of the angle phi
  25.     end   
  26.     sin_phi_s = sin(phir)^2;
  27.     cos_phi_s = cos(phir)^2;
  28.     % Generate aspect angle vector
  29.     theta = 0.:.05:180;
  30.     thetar = theta * dtr;
  31.     if(a ~= b & a ~= c)
  32.         rcs(q,:) = (pi * a^2 * b^2 * c^2) ./ (a^2 * cos_phi_s .* (sin(thetar).^2) + ...
  33.             b^2 * sin_phi_s .* (sin(thetar).^2) + ...
  34.             c^2 .* (cos(thetar).^2)).^2 ;
  35.     elseif(a == b & a ~= c)
  36.         rcs(q,:) = (pi * b^4 * c^2) ./ ( b^2 .* (sin(thetar).^2) + ...
  37.             c^2 .* (cos(thetar).^2)).^2 ;
  38.     elseif (a == b & a ==c)
  39.         rcs(q,:) = pi * c^2;
  40.     end
  41. end
  42. rcs_db = 10.0 * log10(rcs);
  43. figure (1);
  44. plot(theta,rcs_db(1,:),'b',theta,rcs_db(2,:),'r:',theta,rcs_db(3,:),'g--','linewidth',1.5);
  45. xlabel ('Aspect angle, Theta [Degrees]');
  46. ylabel ('RCS - dBsm');
  47. title (['Ellipsoid with (a,b,c) = (', [as],', ', [bs],', ', [cs], ')  meter'])
  48. legend ('phi = 0^o','phi = 45^o','phi = 90^o')
  49. grid;