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

传真(Fax)编程

开发平台:

Matlab

  1. clear all
  2. indes = 0;
  3. eps =0.00001;
  4. a1 =.125;
  5. h = 1.;
  6. lambda = 3.0e+8 /9.5e+9;
  7. lambda = 0.00861;
  8. index = 0;
  9. for theta = 0.0:.1:90-.1
  10.    index = index +1;
  11.    theta = theta * pi /180.;
  12.    rcs(index) = (lambda * a1 * sin(theta) / ...
  13.       (8 * pi * (cos(theta))^2)) + eps;
  14. end
  15. theta*180/pi
  16. theta = pi/2;
  17. index = index +1
  18. rcs(index) = (2 * pi * h^2 * a1 / lambda )+ eps;
  19. for theta = 90+.1:.1:180.
  20.    index = index + 1;
  21.    theta = theta * pi / 180.;
  22.    rcs(index) = ( lambda * a1 * sin(theta) / ...
  23.       (8 * pi * (cos(theta))^2)) + eps;
  24. end
  25. %%%%%%%%%%%
  26. r = a1;
  27. index = 0;
  28. for aspect_deg = 0.:.1:180
  29.    index = index +1;
  30.    aspect = (pi /180.) * aspect_deg; 
  31. % Compute RCS using Eq. (2.37)
  32.    if (aspect == 0 | aspect == pi)
  33.        rcs_po(index) = (4.0 * pi^3 * r^4 / lambda^2) + eps;
  34.        rcs_mu(index) = rcs_po(1);
  35.     else
  36.        x = (4. * pi * r / lambda) * sin(aspect);
  37.        val1 = 4. * pi^3 * r^4 / lambda^2;
  38.        val2 = 2. * besselj(1,x) / x;
  39.        rcs_po(index) = val1 * (val2 * cos(aspect))^2 + eps;
  40.     end
  41.  end
  42. rcs_t =(rcs_po + rcs);
  43. %%%%%%%%%%%%%
  44. angle = 0:.1:180;
  45. plot(angle,10*log10(rcs_t(1:1801)));
  46. grid