CST微带传输线 - 本征模 - Beta相位常数,色散图, 相速度, Pierce耦合阻抗

关于Beta之前的文章解释了很多,这期说说Alphaα 是衰减常数(attenuation constant),表示波损耗,和S21插损相关,但这几个量很多人还是搞不清楚。

首先,S21和插损Insertion Loss严格上讲是不一样的,但是非常相近,所以很多时候通用就通用了。初学者对Alpha的理解往往比较混乱,单位是什么?有没有dB?计算需要一条传输线线还是两条?你就说头大不?甚至本征模也可以计算α ,听说过么?

先澄清下文中用到的定义:

α 的单位是Np/m, 乘上8.685就变成了dB/m,广义上的α 就是某个数每单位距离。

S21(dB) = 20log(S21(linear))

L, L1, L2 是传输线长度,单位都是米。

直接上表:CST中计算α 的8个方法:

方法

求解器

传输线数量

CST后处理方式

公式

大概数值或量级

1

T, F

1

S21(dB)+手动

α= -S21/L

dB/m

0.1dB/L

>0

2

T, F

1

S21(linear)+手动

α= S21^(-2/L)

1/m

0.9^(-2/L)

>1

3

T, F

S21(dB)+手动

α= (S21- S21)/(L1-L2)

dB/m

0.1dB/L

>0

4

T, F

S21(linear)+手动

α= (S21/ S21)^(-2/(L1-L2))

1/m

0.9^(-2/L)

>1

5

T, F

1

S11+S21+自动

Power Attenuation = 

[(1-|S11|^2)/|S21|^2]^(0.001/L)

1/mm

(可转换为1/m, dB/m)

1.1^(0.001/L)

>1

6

T, F

1

S11+S21+自动

Gamma=

1/L*cosh^-1(A)

Np/m

>0

7

E

1 (单元)

Q, Vg, ω + 自动

α= ω/(Vg*2*Q)

Np/m

(可转换为1/m, dB/m)

>0

8

T, F

1

Macro

Results->

1D Results->  Calculate ABCD – Zc – Beta - Alpha

过时不用

 我们准备两条有铜损和介质损的传输线,依次验证这些计算方法:

             

40mm 

60mm 

 

方法1

这种方法比较简单,假设S21代表整条传输线的衰减,那么在计算好S21后,可直接用后处理计算,不同长度α 相同:

方法2

这个和方法1一样,只不过是线性计算式:

方法3

这个适用于只有两个不同长度的S21相减才能表征传输线的损耗的情况。

 方法4

同方法3一样,只不过是线性式:

方法5

这个方法是直接用后处理S-parameters -> Extract TL Properties from S-parameters自动给出传输线属性结果,其中结果之一就是Power Attenuation功率衰减:

所用公式更完整,考虑到了S11回波损耗,不再是简单的S21了,这个就是插损。单位是1/mm,也可手动转化成1/m:

方法6

其实后处理S-parameters -> Extract TL Properties fromS-parameters自动给出AlphaBeta, 所以Alpha就是Gamma的实部:

但是这个单位是Np/m, 可以换成dB/m进行比较:

 以上这些方法的前提都是S参数要计算准确,所以传输线模式要控制好。下面的方法7比较特殊,只需要传输线单元仿本征模就可以了,有时候更容易一些:

 

方法7

公式如下:

需要群速度,角频率,品质因子Q

与本征模仿真慢波结构单元基本相同,传输方向是周期结构。

 扫描角为phase

 

求解器不计算损耗,留给后处理:

先添加本征模后处理,获得Q

运行慢波结构macro

参数扫描开始:

仿真结束后,可看Q值品质因子:

 将Q值对phase变成一条线,也将模式频率对phase变成一条线:

获得Q对频率的曲线:

 再提取一条频率对频率的曲线:

 将这三个量组合成Alpha

与其他方法获得的Alpha基本相同:

 方法8

有些老版用户可能知道有个宏也可以计算AlphaBeta,代码比较老,不推荐用。

最后一个问题,dB/m与1/m 直接的换算:

α (dB/m) = 10log(α (1/m))

小结:

  1.  可用S21或俩S21的差计算Alpha;
  2.  可用插损计算Alpha; 后处理自动得到3) 可用本征模计算Alpha; 有时比S21更方便容易, 比如复杂周期性的慢波结构。 

以下是使用MATLAB绘制费米能级-本征费米能级的像的代码和步骤: 1.首先,我们需要定义一些变量,例如温度,费米能级,本征费米能级等。这些变量的值可以从引用中获取。 2.接下来,我们需要使用MATLAB的plot函数绘制费米能级和本征费米能级的曲线。我们可以使用不同的颜色和线型来区分它们。 3.最后,我们需要添加一些标签和标题,以使像更加清晰和易于理解。 ```matlab % 定义变量 T = 300; % 温度 k = 8.617e-5; % 玻尔兹曼常数 Nc = 2.8e19; % 晶格中的电子数 Nv = 1.04e19; % 晶格中的空穴数 Eg = 1.12; % 硅的带隙宽度 ni = sqrt(Nc * Nv) * exp(-Eg / (2 * k * T)); % 本征载流子密度 Ef = -0.05; % 费米能级 Ec = Ef + Eg / 2; % 导带底部 Ev = Ef - Eg / 2; % 价带顶部 % 绘制费米能级和本征费米能级的曲线 E = linspace(-0.5, 0.5, 1000); % 能量范围 f = 1 ./ (1 + exp((E - Ef) / (k * T))); % 费米分布函数 f0 = 1 ./ (1 + exp((E - Ec) / (k * T))); % 导带底部的费米分布函数 f1 = 1 ./ (1 + exp((E - Ev) / (k * T))); % 价带顶部的费米分布函数 plot(E, f, 'r-', 'LineWidth', 2); % 绘制费米能级曲线 hold on; plot(E, f0, 'b--', 'LineWidth', 2); % 绘制本征费米能级曲线 plot(E, f1, 'b--', 'LineWidth', 2); ylim([0 1.2]); % 设置y轴范围 legend('费米能级', '本征费米能级'); % 添加例 xlabel('能量 (eV)'); % 添加x轴标签 ylabel('费米分布函数'); % 添加y轴标签 title('费米能级-本征费米能级像'); % 添加标题 ```
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

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

余额充值