% 定义圆波导几何尺寸和介质参数
a = 1; % 内半径
b = 2; % 外半径
epsilon_r = 2.2; % 相对介电常数
% 计算TE11模式的截止频率和波长
c = 299792458; % 真空中的光速
f_cutoff = c/(2*pi*a*sqrt(epsilon_r));
lambda_cutoff = 2*a/sqrt(epsilon_r);
disp(['TE11 mode cutoff frequency: ', num2str(f_cutoff/1e9), ' GHz']);
disp(['TE11 mode cutoff wavelength: ', num2str(lambda_cutoff*1e3), ' mm']);
% 定义计算区域的范围和分辨率
x_range = [-b b];
y_range = [-b b];
n_points = 100;
% 生成计算区域的网格
[x, y] = meshgrid(linspace(x_range(1), x_range(2), n_points), ...
linspace(y_range(1), y_range(2), n_points));
r = sqrt(x.^2+y.^2);
% 计算TE11模式的电场分量
E_theta = zeros(size(x));
for i = 1:size(x, 1)
for j = 1:size(x, 2)
if r(i, j) > a && r(i, j) < b
E_theta(i, j) = besselj(1, 2*pi*r(i, j)/lambda_cutoff)/r(i, j)*cos(