RCS Calculation Cylinder From Formula

Hello,

I am trying to formulate a RCS formula using MATLAB. The formula is given in the screenshot.

My code is given below.

N=4000;
sum=0;

for n=(-4000):N
    
    k=2;
    c=physconst('lightspeed');
    f=2.52e12;
    lambda=c/f;
    k0=(2*pi)/(lambda);
 	a = [0.005, 0.010, 0.015, 0.020, 0.025, 0.030, 0.035, 0.040];
    Z=k0*a;
    teta=90;
    teta=deg2rad(teta);
    H=besselh(n,k,Z);
    B=besselj(n,Z);
    
    sum=sum+(B./H)*exp(1i*n*teta);
    
end
    
    summ=abs(sum).^2;
    RCS=10*log((4./(k0*a*pi)).*summ);
    plot(a,RCS)
    xlabel('a(mm)')
    ylabel('RCS')
    grid on

Can you help me with this because I am getting wrong results.

Thank you in advence.

Answer 1

[Tarunbir Gambhir]

A better approach to solve a series formula is to use the symsum function. I suggest you create an anonymous function in terms of 'theta' and use symsum for the series.

[Deniz Sezgin]

Thank you very much.