RF高频腔设计（3）

我们从能量守恒的角度来通俗地理解波导。
波从波源A点产生，如果是在自由空间传输，则波是以球面的形式扩散，要传输到1000米后的B点，则球面积变为了$4\ast \pi \ast r^2$，假设接收能量的入口面积为100平方厘米，则接收到的能量只有$\frac{1}{10^{-8}}$。

而换一种方式，我们通过一根管道来传输，如果不考虑介质的损耗，则有可能实现接近100%的能量传输。

但是波在管道传输是来回反射传播的，不同的波束之间有相互的干扰，如果波峰恰好遇到波峰，则波会增强，如果每次都如此，则波会一直传播下去。如果波峰遇到波谷，则波会湮灭，则无法传播下去。因此只有特定模式的波可以在波导中传播。
以上是我个人对波导的通俗的理解。

3.3 圆波导

之前的方波导的构建是从两个平面波的叠加开始的。
现在考虑无穷个平面波的叠加，$\frac{1}{2\pi }\int \mathrm{d}\alpha$，并在圆柱坐标中分析：
即
之前的
$$\begin{gathered} E_y\propto \left(\cos \left(\omega t-\frac{\omega }{c}(z\cos \alpha +x\sin \alpha )\right)+\cos \left(\omega t-\frac{\omega }{c}(z\cos \alpha -x\sin \alpha )\right)\right) \\ =2\cos \left(\frac{\omega }{c}\sin \alpha x\right)\cos \left(\omega t-\frac{\omega }{c}\cos \alpha z\right)=2\cos (k_{\perp }x)\cos (\omega t-k_{z}z) \end{gathered}$$
变为
$\frac{1}{2\pi }\int 2\cos \left(\frac{\omega }{c}\sin \alpha x\right)\cos \left(\omega t-\frac{\omega }{c}\cos \alpha z\right)d\alpha$
代入圆柱坐标：
$x=\rho \sin \vartheta ,z=\rho \cos \vartheta$
得到：
$$\begin{array}{r}\frac{1}{2\pi }{\int }_0^{2\pi }\cos (k\rho (\cos \alpha \cos \vartheta +\sin \alpha \sin \vartheta ))\mathrm{d}\alpha =\frac{1}{2\pi }{\int }_0^{2\pi }\cos \left(k\rho \right(\cos (\alpha -\vartheta )\left)\right)\mathrm{d}(\alpha -\vartheta )=J_0(k\rho )\end{array}$$
这是贝塞尔函数 𝐽₀ 的积分表示。如下图，我们现在得到了一个由贝塞尔函数 𝐽₀ 描述的径向驻波模式。由于贝塞尔函数 𝐽₀ 的零点，存在一些径向位置，电场始终为零。在这些半径处，可以插入理想导电的边界条件，而不会扰乱这种场分布，这样我们就可以构建一个圆形波导！

TM模式，截止频率${\omega }_c$与波导半径的关系为：
$\frac{{\omega }_{\mathrm{c}}}{c}=\frac{{\chi }_{m,n}}{a}$
其中${\chi }_{m,n}$是m阶贝塞尔函数的第n个零点。
TE模式，截止频率${\omega }_c$与波导半径的关系为：
$\frac{{\omega }_{\mathrm{c}}}{c}=\frac{{\chi }_{m,n}^{\prime }}{a}$
其中${\chi }_{m,n}^{\prime }$是m阶贝塞尔函数导数的第n个零点。

我们可以将电磁波理解为一只老鼠，波长为老鼠的大小，去钻一个管道，频率越小，则波长越大，则老鼠越大，到了截止频率，则老鼠就卡住进不去了。

这也是为什么在很多电脑和电子射频的金属外壳上有很多小孔的原因，可以增加散热，同时因为孔很小，可以挡住电磁波的泄漏，起到电磁屏蔽的作用。

我们可以看到在半径$\alpha$固定的情况下，模式的截止频率和贝塞尔函数或贝塞尔函数的导数的零点成正比。
圆波导的最低阶模式为TE11，此时${\chi }_{1,1}^{\prime }=1.84118$
最低阶的TM模式为TM01，此时${\chi }_{0,1}=2.40483$

3.4 矩形波导和圆波导

通过横向特征函数 𝑇 可以将波导推广到具有任意横截面的情况。
横向的波动方程为：
$\Delta T+{\left(\frac{{\omega }_{\mathrm{c}}}{c}\right)}^{2}T=0$
符号 △ 表示 Laplace 算子，也称为拉普拉斯算子，通常用于描述空间中的二阶偏微分。
${\nabla }^2=\frac{{\partial }^2}{\partial x^2}+\frac{{\partial }^2}{\partial y^2}+\frac{{\partial }^2}{\partial z^2}$

矩形波导TE的横向特征函数为：
$$T_{mn}=\frac{1}{\pi }\sqrt{\frac{ab{\epsilon }_m{\epsilon }_n}{(mb{)}^2+(na{)}^2}}\cos \left(\frac{m\pi }{a}x\right)\cos \left(\frac{n\pi }{b}y\right)$$
矩形波导TM的横向特征函数为：
$$T_{mn}=\frac{2}{\pi }\sqrt{\frac{ab}{(mb{)}^2+(na{)}^2}}\sin \left(\frac{m\pi }{a}x\right)\sin \left(\frac{n\pi }{b}y\right)$$
圆波导的TE横向特征函数为：
$$T_{mn}=\sqrt{\frac{{\epsilon }_m}{\pi ({\chi }_{mn}^{\prime 2}-m^2)}}\frac{J_m\left({\chi }_{mn}^{\prime }\frac{\rho }{a}\right)}{J_m({\chi }_{mn}^{\prime })}\left\{\begin{array}{c}\cos (m\phi )\\ \sin (m\phi )\end{array}\right\}$$
圆波导的TM横向特征函数为：
T m n = ε m π J m ( χ m n ρ a ) J m − 1 ( χ m n ) { sin ⁡ ( m φ ) } T_{mn}=\sqrt{\frac{\varepsilon_{m}}{\pi}}\frac{J_{m}\left(\chi_{mn}\frac{\rho}{a}\right)}{J_{m-1}(\chi_{mn})}\left\{\sin(m\varphi)\right\} Tmn​=πεm​​ ​Jm−1​(χmn​)Jm​(χmn​aρ​)​{sin(mφ)}
每种模式的纵向横电场和磁场分布由函数 𝑈(𝑧) 和 𝐼(𝑧) 描述。
𝑈(𝑧) 和 𝐼(𝑧) 是从纵向波动方程中导出的，此外还可以导出特征阻抗 𝑍₀
$$\frac{\mathrm{d}U(z)}{\mathrm{d}z}+j{k}_z{Z}_{0}I(z)=0,\frac{\mathrm{d}I(z)}{\mathrm{d}z}+\frac{j{k}_z}{Z_0}U(z)=0$$

---

推导过程如下：
欧姆定律：
$\mathbf{J}=\sigma \mathbf{E}$
理想波导中忽略电流密度，则有：
$0=\sigma \mathbf{E}$
理想波导中E和H沿着波导轴向方向传播，因此它们只依赖于z和t。
Maxwell’s 梯度方程：
$\nabla \times \mathbf{E}=-\mu \frac{\partial \mathbf{H}}{\partial t}$
$\nabla \times \mathbf{H}=\epsilon \frac{\partial \mathbf{E}}{\partial t}$
可变换为：
$\frac{\partial E_z}{\partial t}=-\mu \frac{\partial H_z}{\partial z}$
$\frac{\partial H_z}{\partial t}=\epsilon \frac{\partial E_z}{\partial z}$
将磁场和电场表示为电压和电流的形式有：

1. 电场$E_z=-\frac{\mathrm{d}U}{\mathrm{d}z}$
2. 磁场$H_z=Z_{0}I$
   其中$Z_0=\sqrt{\frac{\mu }{\epsilon }}$是波导的特征阻抗。
   将其代入上面的方程有：
   $$-\frac{\partial }{\partial t}\left(-\frac{\mathrm{d}U}{\mathrm{d}z}\right)=-\mu \frac{\partial }{\partial z}(Z_{0}I)\Rightarrow \frac{\mathrm{d}U}{\mathrm{d}z}=j{k}_z{Z}_{0}I$$
   $$\frac{\partial }{\partial t}(Z_{0}I)=\epsilon \frac{\partial }{\partial z}\left(-\frac{\mathrm{d}U}{\mathrm{d}z}\right)\Rightarrow \frac{\mathrm{d}I}{\mathrm{d}z}=\frac{j{k}_z}{Z_0}U$$

---