用matlab写出来

以下是一个用Matlab计算矩形波导TM模式截止波长和截止频率的示例代码：

% 定义矩形波导尺寸和模式阶数
a = 0.02; % 宽度
b = 0.01; % 高度
m = 1; % TE模式阶数
n = 1; % TM模式阶数

% 计算TE模式截止波长和截止频率
lambda_c_TE = 2*a / (m + 1);
f_c_TE = 3e8 / lambda_c_TE; % 光速为3e8 m/s

% 计算TM模式截止波长和截止频率
lambda_c_TM = 2*b / (n + 1);
f_c_TM = 3e8 / lambda_c_TM;

% 输出结果
disp(['TE mode cutoff wavelength: ', num2str(lambda_c_TE), ' m']);
disp(['TE mode cutoff frequency: ', num2str(f_c_TE), ' Hz']);
disp(['TM mode cutoff wavelength: ', num2str(lambda_c_TM), ' m']);
disp(['TM mode cutoff frequency: ', num2str(f_c_TM), ' Hz']);

这个Matlab代码中，首先定义了矩形波导的宽度a和高度b，以及TE模式和TM模式的阶数m和n。然后，分别根据TE模式和TM模式的截止波长公式，计算出对应的截止波长和截止频率。最后，用disp函数输出结果。

需要注意的是，这个示例代码只是一个简单的例子，实际应用中需要根据具体情况进行计算。例如，如果需要考虑波导的损耗和非理想性等因素，需要采用更复杂的计算模型。