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

传真(Fax)编程

开发平台:

Matlab

  1. % This program calculates the back-scattered RCS for a perfectly
  2. % conducting sphere using Eq.(13.28), and produce plots similar to Fig.10.8 
  3. % Spherical Bessel functions are computed using series approximation and recursion.
  4. clear all
  5. eps   = 0.00001;
  6. index = 0;
  7. % kr limits are [0.05 - 15] ===> 300 points
  8. for kr = 0.05:0.05:15
  9.    index = index + 1;
  10.    sphere_rcs   = 0. + 0.*i;
  11.    f1    = 0. + 1.*i;
  12.    f2    = 1. + 0.*i;
  13.    m     = 1.;
  14.    n     = 0.;
  15.    q     = -1.;
  16.    % initially set del to huge value
  17.    del =100000+100000*i;
  18.    while(abs(del) > eps)
  19.       q   = -q;
  20.       n   = n + 1;
  21.       m   = m + 2;
  22.       del = (2.*n-1) * f2 / kr-f1;
  23.       f1  = f2;
  24.       f2  = del;
  25.       del = q * m /(f2 * (kr * f1 - n * f2));
  26.       sphere_rcs = sphere_rcs + del;
  27.    end
  28.    rcs(index)   = abs(sphere_rcs);
  29.    sphere_rcsdb(index) = 10. * log10(rcs(index));
  30.    end
  31. figure(1);
  32. n=0.05:.05:15;
  33. plot (n,rcs,'k');
  34. set (gca,'xtick',[1 2 3 4 5 6 7 8 9 10 11 12 13 14 15]);
  35. %xlabel ('Sphere circumference in wavelengths');
  36. %ylabel ('Normalized sphere RCS');
  37. grid;
  38. figure (2);
  39. plot (n,sphere_rcsdb,'k');
  40. set (gca,'xtick',[1 2 3 4 5 6 7 8 9 10 11 12 13 14 15]);
  41. xlabel ('Sphere circumference in wavelengths');
  42. ylabel ('Normalized sphere RCS - dB');
  43. grid;
  44. figure (3);
  45. semilogx (n,sphere_rcsdb,'k');
  46. xlabel ('Sphere circumference in wavelengths');
  47. ylabel ('Normalized sphere RCS - dB');