一种基于场路协同的太赫兹功率放大器电路设计方法与流程

本发明涉及一种基于场路协同的太赫兹功率放大器电路设计方法，属于太赫兹电路设计技术领域。

背景技术：

太赫兹的频率范围为100ghz-10thz，位于微波与红外辐射之间，具有许多独特的优势。比如，与微波毫米波相比，太赫兹频率很高，具有空间分辨率高的优点；与红外辐射相比，太赫兹受恶劣天气影响较小，可以有效地穿透沙尘和烟雾。因此太赫兹在探测防撞系统、安全检查、雷达成像、无损检测、宽带通信、射电天文以及生物医疗具有广阔的应用前景。

但是太赫兹电路的设计依然面临诸多问题，其中之一就是需要进行电路全版图仿真优化设计，这样会带来设计周期长，成本高的问题。对于微波毫米波电路，采用传统的场路协同仿真方法设计的电路与全版图仿真结果较为一致，无需进行进一步的优化设计。而太赫兹电路由于频率很高，分布参数效应异常明显，较小的尺寸变化就会影响设计结果，因此，往往需要进行电路全版图迭代优化设计。但是进行全版图迭代优化设计往往需要很长的仿真时间，仿真优化效率低下，带来极高的仿真代价。

技术实现要素：

本发明的目的是提出一种基于场路协同的太赫兹功率放大器电路设计方法，在进行太赫兹电路全版图仿真设计时，通过引入重心插值算法，建立了电路版图中某一元件的参数化模型。通过所建立的参数化模型可快速得到元件不同设计变量取值下的等效端口网络模型，缩短了版图的仿真时间，达到了提高版图设计效率的目的。

本发明提出的基于场路协同的太赫兹功率放大器电路设计方法，包括以下步骤：

(1)设定待设计太赫兹功率放大器的增益g；

(2)利用场路协同仿真方法，对太赫兹功率放大器电路全版图，进行仿真计算得到太赫兹功率放大器电路的增益g0，将步骤(1)的增益g与仿真计算得到的增益g0比较，若g0≥g，则判定太赫兹功率放大器电路已达标，完成设计；若g0＜g，则判定太赫兹功率放大器电路未达标，进行步骤(3)；

(3)建立太赫兹功率放大器电路的元件参数化模型，包括以下步骤：

(3-1)从太赫兹功率放大器电路的所有元件中任意选取一个元件r，设定该元件r的设计变量为l，并设定设计变量l的取值范围为[l1,ln]，[l1,ln]中的插值节点值为(l1,l2,….,ln-1,ln)，设定太赫兹功率放大器元件r的工作频带为[f1,fm],[f1,fm]中所要计算的离散频点值为(f1,f2,….,fm-1,fm)；

(3-2)采用有限元方法，计算出太赫兹功率放大器版图中元件r在插值节点值(l1,l2,….,ln-1,ln)和离散频点值(f1,f2,….,fm-1,fm)处的散射参数值s(fi,lv)，其中fi为(3-1)中第i个离散频点值，i＝1,...,m，lv为(3-1)中第v个插值节点值，v＝1,...,n；

(3-3)利用重心插值方法，建立元件r的设计变量l在频点(f1,f2,….,fm-1,fm)处的参数化模型如下：

其中，ω是权重系数，π表示连乘，v＝1,...,n，i＝1,...,m，

设计变量l的取值范围为[l1,ln],且l的取值不等于插值节点值，即l≠l1,l2,….,ln-1,ln；

(4)从设计变量l的取值范围[l1,ln]中取出最小值l1，令lj＝l1，并判断lj是否是插值节点值，若是，则由步骤(3-2)直接获得相应的散射参数值s(fi,lj)，若不是，则求解步骤(3)的元件参数化模型，得到并将s(fi,lj)转换为等效端口网络模型，模型为touchstonefile格式；

(5)用步骤(4)的等效端口网络模型代替太赫兹功率放大器全版图中的元件r，利用场路协同仿真方法，对太赫兹功率放大器全版图进行仿真计算，得到太赫兹功率放大器的增益g1，将g1与待设计太赫兹功率放大器的增益g比较，若g1≥g，则判定太赫兹功率放大器电路已达标，完成设计，若g1＜g，则判定太赫兹功率放大器电路未达标，进行步骤(6)；

(6)使lj＝lj+d，其中d＝(ln-l1)/(10n)，n为步骤(3-1)中插值节点的个数，重复步骤(4)和步骤(5)，经过10n次迭代，遍历设计变量l取值范围[l1,ln]中的最小值l1到最大值ln，lj具有10n个值，通过10n次重复步骤(4)和步骤(5)的迭代过程，得到10n个太赫兹功率放大器的增益值gs(s＝1,2,...,10n)，取10n个太赫兹功率放大器的增益值gs中的最大值max(gs)，将max(gs)与步骤(1)的增益g进行比较，若max(gs)≥g，则完成设计；若max(gs)＜g，则重新从太赫兹功率放大器电路中选取元件r以外的其它元件，重复步骤(3)-步骤(6)。

本发明提出的基于场路协同的太赫兹功率放大器电路设计方法，其优点是：本发明利用重心插值算法可建立版图元件关于设计变量的参数化模型，可快速得到版图元件在设计变量多个取值下的等效端口网络模型，缩短了基于场路协同仿真的全版图设计时间，有效提高了设计效率。

附图说明

图1为本发明方法的流程框图。

图2为本发明方法涉及的太赫兹功率放大器的全版图。

图3为设计前的太赫兹功率放大器增益仿真曲线。

图4为用本方法设计后的太赫兹功率放大器增益仿真曲线。

具体实施方式

本发明提出的基于场路协同的太赫兹功率放大器电路设计方法，其流程框图如图1所示，包括以下步骤：

(1)设定待设计太赫兹功率放大器的增益g；

(2)利用场路协同仿真方法，对太赫兹功率放大器电路全版图，如图2所示，进行仿真计算得到太赫兹功率放大器电路的增益g0，将步骤(1)的增益g与仿真计算得到的增益g0比较，若g0≥g，则判定太赫兹功率放大器电路已达标，完成设计；若g0＜g，则判定太赫兹功率放大器电路未达标，进行步骤(3)；

(3)建立太赫兹功率放大器电路的元件参数化模型，包括以下步骤：

(3-1)从太赫兹功率放大器电路的所有元件中任意选取一个元件r，太赫兹功率放大器电路中的元件可以是传输线长度ltl、传输线宽度wtl、变压器线圈的半径r和叉指电容器长度lcap等，设定该元件r的设计变量为l，并设定设计变量l的取值范围为[l1,ln]，[l1,ln]中的插值节点值为(l1,l2,….,ln-1,ln)，设定太赫兹功率放大器元件r的工作频带为[f1,fm],[f1,fm]中所要计算的离散频点值为(f1,f2,….,fm-1,fm)；

(3-2)采用有限元方法，计算出太赫兹功率放大器版图中元件r在插值节点值(l1,l2,….,ln-1,ln)和离散频点值(f1,f2,….,fm-1,fm)处的散射参数值s(fi,lv)，其中fi为(3-1)中第i个离散频点值，i＝1,...,m，lv为(3-1)中第v个插值节点值，v＝1,...,n；

(3-3)利用重心插值方法，建立元件r的设计变量l在频点(f1,f2,….,fm-1,fm)处的参数化模型如下：

其中，ω是权重系数，π表示连乘，v＝1,...,n，i＝1,...,m

设计变量l的取值范围为[l1,ln],且l的取值不等于插值节点值，即l≠l1,l2,….,ln-1,ln。

(4)从设计变量l的取值范围[l1,ln]中取出最小值l1，令lj＝l1，并判断lj是否是插值节点值，若是，则由步骤(3-2)直接获得相应的散射参数值s(fi,lj)，若不是，则求解步骤(3)的元件参数化模型，得到并将s(fi,lj)转换为等效端口网络模型，模型为touchstonefile格式，具体转换方法为已有技术，可参见touchstonefileformatspecification,eia/ibisopenforum,2002http://vhdl.org/ibis/connector/touchstone_spec11.pdf。

(5)用步骤(4)的等效端口网络模型代替太赫兹功率放大器全版图中的元件r，利用场路协同仿真方法，对太赫兹功率放大器全版图进行仿真计算，得到太赫兹功率放大器的增益g1，将g1与待设计太赫兹功率放大器的增益g比较，若g1≥g，则判定太赫兹功率放大器电路已达标，完成设计，若g1＜g，则判定太赫兹功率放大器电路未达标，进行步骤(6)；

(6)令lj＝lj+d，其中d＝(ln-l1)/(10n)，n为步骤(3-1)中插值节点的个数。经过10n次迭代，lj的值就从设计变量l取值范围[l1,ln]中的最小值l1变化到了最大值ln。此时，lj具有10n个值，通过迭代重复10n次步骤(4)和步骤(5)的过程，会计算出10n个太赫兹功率放大器的增益值gs(s＝1,2,...,10n)，取其中的最大值max(gs)与步骤(1)的增益g进行比较，若max(gs)≥g，则完成设计，记录下此时对应的lj的值；若max(gs)＜g，则重新从太赫兹功率放大器电路中选取元件r以外的其它元件，重复步骤(3)-步骤(6)。

以下介绍本发明方法的一个实施例：

(1)设定待设计太赫兹功率放大器在频点130ghz处的增益g＝10db，并设定该太赫兹功率放大器电路中的传输线长度l、传输线宽度w、变压器线圈的半径r和叉指电容器长度lcap；

(2)利用场路协同仿真方法，对太赫兹功率放大器电路全版图，如图2所示，进行仿真计算得到太赫兹功率放大器电路工作在频点130ghz处的增益为g0＝9.1db，仿真结果如图3所示。图3为设计前太赫兹功率放大器增益仿真曲线。从图3中可以看出，太赫兹功率放大器增益曲线的峰值点处于约126ghz频点处，并非130ghz频点处。将步骤(1)的增益g＝10db与仿真计算得到的增益g0＝9.1db比较。由于g0＜g，则判定太赫兹功率放大器电路未达标，进行步骤(3)；

(3)建立太赫兹功率放大器电路的元件参数化模型，包括以下步骤：

(3-1)从太赫兹功率放大器电路的所有元件中选取一段传输线r，设定该传输线r的设计变量为传输线长度l，并设定传输线长度l的取值范围为[5,14]um，[5,14]um中的4个插值节点值为(5,8,11,14)um，设定太赫兹功率放大器电路中传输线r的工作频带为[90,150]ghz,[90,150]ghz中所要计算的61个离散频点值为(90,91,….,149,150)ghz，均匀间隔1ghz；

(3-2)采用有限元方法，计算出太赫兹功率放大器版图中传输线r在插值节点值(5,8,11,14)um和离散频点值(90,91,….,149,150)ghz处的散射参数值s(fi,lv)，其中fi为(3-1)中第i个离散频点值，i＝1,...,61，lv为(3-1)中第v个插值节点值，v＝1,...,4；

(3-3)利用重心插值方法，建立传输线r关于长度l在频点(90,91,….,149,150)ghz处的参数化模型如下：

其中，ω是权重系数，符号π表示连乘

fi为(3-1)中第i个离散频点值，i＝1,...,61

lv为(3-1)中第v个插值节点值，v＝1,...,4

传输线长度l的取值范围为[5,14]um,且l的取值不等于插值节点值，即l≠(5,8,11,14)um。

(4)从传输线长度l的取值范围[5,14]中取最小值lj＝5um，经过判断，此时lj＝5um属于插值节点值，则直接从步骤(3-2)中获得s(fi,5)，若不属于插值节点值，则求解步骤(3)的元件参数化模型，得到然后将s(fi,lj)转换为等效端口网络模型，模型为touchstonefile格式，具体转换方法为已有技术，可参见touchstonefileformatspecification,eia/ibisopenforum,2002http://vhdl.org/ibis/connector/touchstone_spec11.pdf

(5)用步骤(4)的等效端口网络模型代替太赫兹功率放大器全版图中的传输线r，利用场路协同仿真方法，对太赫兹功率放大器全版图进行仿真计算，得到太赫兹功率放大器的增益g1＝8.9db，将g1与待设计太赫兹功率放大器的增益g比较，由于8.9db<10db，则判定太赫兹功率放大器电路未达标，进行步骤(6)；

(6)令lj＝lj+d，其中d＝0.225um。经过40次迭代，lj的值就从设计变量l取值范围[5,14]um中的最小值5um变化到了最大值14um。此时，lj具有40个值，将其中lj≠(8,11,14)um的值都代入公式计算，对于s(fi,8)，s(fi,11)和s(fi,14)可以通过步骤(3-2)得到。这样就可以得到40组s(fi,5.225)，s(fi,5.45)，…,s(fi,13.775),s(fi,14)，其中fi＝(90,91,….,149,150)ghz。然后转换得到40个touchstonefile格式的等效端口网络模型。重复步骤(5)，计算出40个太赫兹功率放大器的增益值gs(s＝1,2,...,40)，得到最大值max(gs)＝11.32db和对应lj的值13.1um。将图2中太赫兹功率放大器全版图中的传输线长度l设定为13.1um，采用场路协同仿真方法，得到设计后太赫兹功率放大器增益仿真曲线，如图4所示。从图中可以看出，在中心频点130ghz处，太赫兹功率放大器增益为11.32db，与步骤(1)的增益g＝10db进行比较，11.32db>10db，完成设计。

以上实施例验证了本发明的正确性和有效性。