电磁场与仿真软件(15)

最新推荐文章于 2024-08-11 23:46:27 发布

Walt Lu

最新推荐文章于 2024-08-11 23:46:27 发布

阅读量316

点赞数 1

文章标签：射频工程

版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。

本文链接：https://blog.csdn.net/weixin_46649493/article/details/126290970

版权

上面2篇都是理论计算,接下来举个实例来说明:

有一列垂直入射介质的电磁波, 介质的相对介电常数为 $\varepsilon _{r}=9$ , 相对磁导率 $\mu _{r}=1$ 即(n=3)

我们按照之前推导的公式:

$\Omega =\begin{bmatrix} 0 &0 &\tilde{k}_{x} \tilde{k}_{y}/\varepsilon_{r} &\mu_{r}-\tilde{k}_{x}^{2} \\ 0 &0 &\tilde{k}_{y}^{2}-\mu_{r} &-\tilde{k}_{x} \tilde{k}_{y}/\varepsilon_{r} \\ \tilde{k}_{x} \tilde{k}_{y}/\mu_{r}&\varepsilon_{r}-\tilde{k}_{x}^{2}/\mu_{r} &0 &0 \\ \tilde{k}_{y}^{2}-\varepsilon_{r} &-\tilde{k}_{x} \tilde{k}_{y}/\mu_{r} &0 &0 \end{bmatrix}$

我们架设的是沿z轴传播的电磁波,k只有在z轴方向上有分量,即 $k_{x},k_{y}=0$

所以得到:

用matlab可以算出:

还是用教程上的图清楚一点:

按照上一篇所写的,代入公式得到:

$\frac{\mathrm{d} }{\mathrm{d} {z}'}\begin{bmatrix} E_{x}\\ E_{y}\\ \vec{\tilde{H}}_{x}\\ \vec{\tilde{H}}_{y} \end{bmatrix}$ = $\begin{bmatrix} -j0.32 &0.32j &0 &0 \\ 0&0 &j0.32 &-j0.32 \\ 0&0 &0.95 &0.95 \\ 0.95&0.95 &0 & 0 \end{bmatrix}$ $\begin{bmatrix} e^{j3.0{z}'} &0 &0 &0 \\ 0&e^{-j3.0{z}'}&0 &0 \\ 0&0 &e^{j3.0{z}'} &0 \\ 0&0 &0 & e^{-j3.0{z}'} \end{bmatrix}$ c

W仅包含了E和H之间的幅度关系, $\lambda$ 描述了这些模式如何在z方向上累积相位,如下图.

接下来看下这个模型的解,中间层的电磁场解

和之前类似的步骤, 先直接写出中间层的微分方程,再代入边界条件.

其中 ${z}'_{i}=k_{0}z_{i}$ 归一化的阻抗

写出在 $L_{i-1},L_{i}$ 边界的边界条件:

$\Psi _{i-1}(k_{0}L_{i-1})=\Psi _{i}(0)$ 即:

$W_{i-1}e^{\lambda_{i-1}k_{0}L_{i-1} }c_{i-1}=W_{i}c_{i}$

写出在 $L,L_{i+1}$ 边界的边界条件上的边界条件:

$\Psi _{i}(k_{0}L_{i})=\Psi _{i+1}(0)$ 即:

$W_{i}e^{\lambda_{i}k_{0}L_{i} }c_{i}=W_{i+1}c_{i+1}$ ---a

两个边界条件看起来很像,我们把a式两边同时乘以 $W_{i+1}^{-1}$ 得到:

$c_{i+1}=W_{i+1}^{-1}W_{i}e^{\lambda_{i}k_{0}L_{i} }c_{i}=T_{i}*c_{i}$

然后如果有多级相连,如下图:

我们只要知道初始入射条件 $c_{0}$ ，然后再通过矩阵相乘: ${T}'_{total}=T5*T4*T3*T2*T1$ 就可以得到经过这些介电物质之后电磁波的解了.矩阵计算非常复杂,模型建立好之后,交由计算机去计算就可以了. 当然我们也需要把 $c_{0}$ 统一写成矩阵形式,

还有补充一下, ${T_{total}}'$ 的乘顺序看起来式反的,实际上这样才是对的,比较直观的理解是 $c_{0}$ 在最后面,

和 $c_{0}$ 最先相乘的一定是T1.当然也可以从传输矩阵的定义去理解~~

到这里,我们还需要入射到第一个边界时产生的反射,这样才是完成方程的解

这里我们去第一个入射边界的边界条件:

$W_{ref}c_{ref}=W_{1}c_{1} ==> c_{1} =W_{1}^{-1}W_{ref}c_{ref}$

我们可以得到完成的传输矩阵:

$T_{global}={T}'_{total}W_{1}^{-1}W_{ref}$

但传输矩阵法有一定局限性,这个会在下一篇中再写,这边先放一张原教材中的图

关注

1
点赞
踩
1

收藏

觉得还不错? 一键收藏
打赏
0
评论
复制链接

分享到 QQ

分享到新浪微博

扫一扫

评论

被折叠的条评论为什么被折叠?

到【灌水乐园】发言

查看更多评论

添加红包

成就一亿技术人!

hope_wisdom

发出的红包

打赏作者

Walt Lu 你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20

扫码支付：¥1

获取中

扫码支付

您的余额不足，请更换扫码支付或充值

实付元

使用余额支付

点击重新获取

扫码支付

钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。