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

传真(Fax)编程

开发平台:

Matlab

  1. function rcs = cylinder (r, h, lambda)
  2. % This program compute monostatic RCS for a cylinder
  3. % Circular symmetry is assumed.
  4. % Plot of RCS versus aspect angle is produced
  5. %clear all
  6. index = 0;
  7. eps =0.00001;
  8. % Enter dimensions of cylinder
  9. %r = .125;
  10. %h = 1.;
  11. % Compute wavelength
  12. %lambda = 3.0e+8 /9.5e+9;
  13. % Compute RCS from zero aspect to broadside
  14. for theta = 0.0:.1:90-.5
  15.    index = index +1;
  16.    theta = theta * pi /180.;
  17.    rcs(index) = (lambda * r * sin(theta) / ...
  18.       (8. * pi * (cos(theta))^2)) + eps;
  19. end
  20. % Compute RCS for broadside specular
  21. theta = pi/2;
  22. index = index +1;
  23. rcs(index) = (2. * pi * h^2 * r / lambda )+ eps;
  24. % Compute RCS from 90 to 180 degrees
  25. for theta = 90+.5:.1:180.
  26.    index = index + 1;
  27.    theta = theta * pi / 180.;
  28.    rcs(index) = ( lambda * r * sin(theta) / ...
  29.       (8. * pi * (cos(theta))^2)) + eps;
  30. end
  31. % Plot results
  32. delta= 180/(index-1)
  33. angle = 0:delta:180;
  34. plot(angle,10*log10(rcs),'k');
  35. grid;
  36. xlabel ('Aspect angle - degrees');
  37. ylabel ('RCS - dBsm');
  38. title  ('Frequency = 9.5 GHz');