I want to plot the attenuation coefficient of circular waveguide using zeros of Bessel function, can someone kindly suggest how to do it for multimodes not for single mode??

7 views (last 30 days)
u=1.84
nu=linspace(0.1e12,1e12,100);
lambda=3e8./nu;
n=1
ko=(2*pi)./lambda;
T=4e-3;
zo=-i*n.*ko.*T;
m=besselj(1,u)./besselj(0,u);
zTE=-i.*(m./u).*zo;
rm=(zTE-zo)./(zTE+zo);
rm1=-(besselj(0,u)+i*besselj(1,u))./(besselj(0,u)-i*besselj(1,u));
Rm=abs(rm).^2;
% Rm1=abs(rm1)
thetaz=asin((u)./(n.*ko.*T));
m1=(2./u).*(zTE.*zo)./(zTE.^2-zo.^2)
m2=-(2./u).*(besselj(1,u)*besselj(0,u))./(besselj(1,u).^2+besselj(0,u).^2)
alpha=(1-Rm)./(2.*T.*cot(thetaz)+(2*T*cot(thetaz).*real(m2)))
plot(nu,(alpha))

Answers (1)

Epsilon
Epsilon on 18 Sep 2024
Edited: Epsilon on 18 Sep 2024
Hi Deepak,
To plot the attenuation coefficient of circular waveguides using the zeros of Bessel function for multiple modes, start with identifying the zeros of the Bessel function to determine the cutoff frequencies for each mode. Then use a loop to calculate the attenuation coefficient for each mode separately.
Here is an exemplar code for your reference:
% Parameters
nu = linspace(0.1e12, 1e12, 100); % Frequency range
lambda = 3e8 ./ nu; % Wavelength
n = 1; % Refractive index
ko = (2 * pi) ./ lambda; % Wave number
T = 4e-3; % Thickness
modes = [1.84, 3.05, 4.2]; % Zeros of the Bessel function for different modes
% Preallocate for efficiency
alpha_modes = zeros(length(modes), length(nu));
% Loop over each mode
for mode_index = 1:length(modes)
u = modes(mode_index); % Current mode zero
zo = -1i * n .* ko .* T;
m = besselj(1, u) ./ besselj(0, u);
zTE = -1i .* (m ./ u) .* zo;
rm = (zTE - zo) ./ (zTE + zo);
Rm = abs(rm).^2;
thetaz = asin((u) ./ (n .* ko .* T));
m2 = -(2 ./ u) .* (besselj(1, u) * besselj(0, u)) ./ (besselj(1, u).^2 + besselj(0, u).^2);
alpha = (1 - Rm) ./ (2 .* T .* cot(thetaz) + (2 * T * cot(thetaz) .* real(m2)));
% Store the result for this mode
alpha_modes(mode_index, :) = alpha;
end
% Plot the results
figure;
hold on;
colors = lines(length(modes));
for mode_index = 1:length(modes)
plot(nu, alpha_modes(mode_index, :), 'DisplayName', ['Mode ' num2str(mode_index)], 'Color', colors(mode_index, :));
end
hold off;
xlabel('Frequency (Hz)');
ylabel('Attenuation Coefficient');
title('Attenuation Coefficient for Multiple Modes');
legend show;
Hope this helps!

Products


Release

R2022b

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!