一种GaNHEMT晶体管小信号模型建模方法

一种gan hemt晶体管小信号模型建模方法
技术领域
1.本发明涉及一种小信号模型建模方法，尤其是涉及一种gan hemt晶体管小信号模型建模方法。


背景技术：

2.半导体器件gan hemt晶体管具有禁带宽、耐高温、击穿电场高、二维电子气浓度大、跨导大、截止频率高、噪声低和开关速度快等优点，在集成电路和微波射频领域应用极为广泛。而器件模型是半导体器件特性表征和集成电路设计的必需途径，在集成电路设计中起着至关重要的作用，模型的精确度将直接影响集成电路设计的准确性。
3.模型参数提取对于建立晶体管小信号模型至关重要，它依赖于精确的测量和可靠的参数提取方法等。现有的参数提取方法基于晶体管的小信号等效电路模型(即小信号模型)，可以分为直接提取法、数值优化法和混合提取法(直接提取和数值优化相结合的方法)。直接提取法被认为是最简单且比较可靠的提取方法，但它需要确切知晓器件的工艺特性来假设器件内元件间的倍数关从而提取元件参数，不适合应用在不同结构的器件的提参建模中，局限性较高。数值优化法虽然可以应用在不同结构的器件的提参建模中，通用性较强，然而这种方法对器件内元件初值的选取要求较高，可能导致无意义的元件参数或者陷入局部最优，参数提取结果可靠性不高。混合提取法结合了直接提取法和数值优化法两者的优点，不但可以应用在不同结构的器件的提参建模中，通用性较强，而且可以得到可靠的元件参数值，但是其参数提取的过程却复杂且耗时。


技术实现要素：

4.本发明所要解决的技术问题是提供一种能够应用在不同结构的器件中，通用性较强，参数提取结果可靠性较高，且参数提取过程简单省时的gan hemt晶体管小信号模型建模方法。
5.本发明解决上述技术问题所采用的技术方案为：一种gan hemt晶体管的小信号模型的建模方法，包括以下步骤：
6.步骤s1：提取最优寄生电容值：
7.步骤s101：在ads射频仿真软件中构建gan hemt晶体管的小信号等效电路模型，该小信号等效电路模型包括寄生元件和本征元件，所述的寄生元件包括寄生电容c1、c2、c3、c4、c5、c6，寄生电感l1、l2、l3以及寄生电阻r1、r2、r3；所述的本征元件包括本征电容c7、c8、c9，本征电阻r4、r5，本征电导g1、g2、g3以及本征压控电流源，所述的本征压控电流源的输出电流i
ds
＝v
gsgm
e-jωτ
，其中，gm为ganhemt晶体管的跨导，τ为gan hemt晶体管栅极电压的时延，v
gs
为流经本征电容c7两端的电压，ω为gan hemt晶体管的角频率，j为虚数符号，e为自然对数的底；寄生电容c1的一端、寄生电容c3的一端和寄生电感l1的一端连接且其连接端等效为gan hemt晶体管的栅极，寄生电容c1的另一端、寄生电感l3的一端和寄生电容c5的一端连接且其连接端等效为gan hemt晶体管的源极，寄生电容c2的一端、寄生电容c4的
一端、寄生电感l1的另一端和寄生电阻r1的一端连接，寄生电容c2的另一端、寄生电感l3的另一端、寄生电阻r3的一端和寄生电容c6连接，寄生电容c3的另一端、寄生电容c5的另一端和寄生电感l2的一端连接且其连接端等效为gan hemt晶体管的漏极，寄生电容c4的另一端、寄生电阻r2的一端、寄生电感l2的另一端和寄生电容c6的另一端连接，寄生电阻r1的另一端、本征电容c7的一端、本征电容c8的一端、本征电导g1的一端和本征电导g2的一端连接，寄生电阻r2的另一端、本征电容c9的一端、本征电导g2的另一端、本征电导g3的一端、本征电阻r5的一端和本征压控电流源的正极连接，寄生电阻r3的另一端、本征电阻r4的一端、本征电导g1的另一端、本征电容c9的另一端、本征电导g3的另一端和本征压控电流源的负极连接，本征电容c7的另一端与本征电阻r4的另一端连接，本征电容c8的另一端与本征电阻r5的另一端连接；将寄生电容c1的电容值记为c
pg
，寄生电容c2的电容值记为c
gsi
,寄生电容c3的电容值记为c
pgd
,寄生电容c4的电容值记为c
gdi
,寄生电容c5的电容值记为c
pd
,寄生电容c6的电容值记为c
dsi
；寄生电感l1的电感值记为lg，寄生电感l2的电感值记为ld，寄生电感l3的电感值记为ls；寄生电阻r1的电阻值记为rg，寄生电阻r2的电阻值记为rd，寄生电阻r3的电阻值记为rs；本征电容c7的电容值记为c
gs
，本征电容c8的电容值记为c
gd
，本征电容c9的电容值记为c
ds
，本征电阻r4的电阻值记为ri，本征电阻r5的电阻值记为r
gd
，本征电导g1的电导值记为g
gsf
，本征电导g2的电导值记为g
gdf
，本征电导g3的电导值记为g
ds
；
8.步骤s102：在冷夹断偏置状态下测试gan hemt晶体管的s参数，得到冷夹断偏置状态下每个频率点处gan hemt晶体管的s参数，其中将频率点的总数量记为n，第k个频率点记为fk，k＝1,2,...,n，将测试得到的冷夹断偏置状态下第k个频率点处gan hemt晶体管的s参数记为sk，sk为2
×
2的矩阵，sk中每个元素均为复数；
9.对gan hemt晶体管的小信号等效电路模型进行简化得到第一个简化电路模型，第一个简化电路模型只包含寄生电容c1、c2、c3、c4、c5、c6和本征电容c7、c8、c9，此时，寄生电容c1的一端、寄生电容c2的一端、寄生电容c5的一端、寄生电容c6的一端、本征电容c7的一端和本征电容c9的一端连接且其连接端等效为gan hemt晶体管的源极，寄生电容c1的另一端、寄生电容c2的另一端、寄生电容c3的一端、寄生电容c4的一端、本征电容c7的另一端和本征电容c8的一端连接且其连接端等效为gan hemt晶体管的栅极，寄生电容c5的另一端、寄生电容c6的另一端、寄生电容c3的另一端、寄生电容c4的另一端、本征电容c8的另一端和本征电容c9的另一端连接且其连接端等效为gan hemt晶体管的漏极；将第一个简化电路模型建立在ads射频仿真软件中；
10.将寄生电容c1、寄生电容c2和本征电容c7的电容值之和作为冷夹断偏置状态下gan hemt晶体管的栅-源分支总电容c
ga
，将寄生电容c3、寄生电容c4和本征电容c8的电容值之和作为冷夹断偏置状态下gan hemt晶体管的栅-漏分支总电容c
gdc
，将寄生电容c5、寄生电容c6和本征电容c9的电容值之和作为冷夹断偏置状态下ganhemt晶体管的漏-源分支总电容c
db
，得到式(1)至式(3)：
11.c
ga
＝c
pg
c
gsi
c
gs
ꢀꢀꢀ
(1)
12.c
gdc
＝c
pgd
c
gd
c
gdi
ꢀꢀꢀ
(2)
13.c
db
＝c
pd
c
dsi
c
ds
ꢀꢀꢀ
(3)
14.步骤s103：将测试得到的冷夹断偏置状态下每个频率点处gan hemt晶体管的s参数分别转化为y参数，得到冷夹断偏置状态下每个频率点处gan hemt晶体管的y参数，将冷
夹断偏置状态下第k个频率点处gan hemt晶体管的y参数记为yk，yk为2
×
2的矩阵，yk中每个元素均为复数，将yk的第1行第1列元素记为y
11,k
，第1行第2列元素记为y
12,k
，第2行第1列元素记为y
21,k
，第2行第2列元素记为y
22,k
；此时，栅-源分支总电容c
ga
、栅-漏分支总电容c
gdc
、漏-源分支总电容c
db
与冷夹断偏置状态下gan hemt晶体管的y参数关系式为：
15.im(y
11,k
) im(y
12,k
)＝c
ga,k
ωkꢀꢀꢀ
(4)
16.im(y
12,k
)＝-c
gdc,k
ωkꢀꢀꢀ
(5)
17.im(y
22,k
) im(y
12,k
)＝c
db,k
ωkꢀꢀꢀ
(6)
18.其中，im()表示复数的虚部，ωk＝2*π*fk，ωk为第k个频率点fk对应的角频率；im(y
ij,k
)表示y
ij,k
的虚部值，i＝1,2，j＝1,2；c
ga,k
表示第k个频率点处c
ga
的取值，c
gdc,k
表示第k个频率点处c
gdc
的取值，c
db,k
表示第k个频率点处c
db
的取值；
19.根据式(4)至式(6)拟合冷夹断偏置状态下第k个频率点处gan hemt晶体管的y参数中各元素虚部与角频率ωk的线性关系，即拟合im(y
11,k
) im(y
12,k
)与ωk、im(y
12,k
)与ωk、im(y
22,k
) im(y
12,k
)与ωk之间的线性关系，该拟合过程在matlab软件中实现，通过matlab软件直接得到im(y
11,k
) im(y
12,k
)与ωk、im(y
12,k
)与ωk、im(y
22,k
) im(y
12,k
)与ωk这三个拟合线性直线以及这三个拟合线性直线的斜率，将im(y
11,k
) im(y
12,k
)与ωk的拟合线性直线的斜率作为栅-源分支总电容c
ga
的取值、im(y
12,k
)与ωk的拟合线性直线的斜率作为栅-漏分支总电容c
gdc
的取值、im(y
22,k
) im(y
12,k
)与ωk的拟合线性直线的斜率作为漏-源分支总电容c
db
的取值，将拟合得到的c
ga
、c
gdc
、c
db
的值分别作为matlab软件中最优寄生电容数值循环程序中栅-源分支总电容、栅-漏分支总电容、漏-源分支总电容对应的搜索上限值，并将最优寄生电容数值循环程序中栅-源分支总电容、栅-漏分支总电容、漏-源分支总电容对应的搜索下限值设为0,即0≤(c
pg
c
gsi
c
gs
)≤c
ga
，0≤(c
pd
c
dsi
c
ds
)≤c
db
,0≤(c
pgd
c
gd
c
gdi
)≤c
gdc
为寄生电容的优化阈值空间；在冷夹断偏置状态下，由gan hemt器件的沟道耗尽层物理特性，假设c
gs
＝c
gd
，将c
gs
＝c
gd
作为最优寄生电容数值循环程序的约束条件；
20.步骤s104：在matlab软件中执行最优寄生电容数值循环程序，最优寄生电容数值循环程序具体流程为：首先在优化空间内产生6个随机分配的寄生电容值作为c
pg
、c
gsi
、c
pgd
、c
gdi
、c
pd
、c
dsi
的电容值，并对这6个电容值进行稳定性分析，即判断此时0≤(c
pg
c
gsi
c
gs
)≤c
ga
，0≤(c
pd
c
dsi
c
ds
)≤c
db
，0≤(c
pgd
c
gd
c
gdi
)≤c
gdc
这三个条件是否同时成立；若上述三个条件不能同时成立，则根据metropolis准则重新产生6个寄生电容数值作为c
pg
、c
gsi
、c
pgd
、c
gdi
、c
pd
、c
dsi
的电容值，若上述三个条件同时成立，则接受这6个寄生电容数值；然后采用这6个寄生电容数值对第一个简化电路模型进行仿真，得到第一个简化电路模型在每个频率点处的仿真s参数，将第一个简化电路模型在第k个频率点处的仿真s参数记为s1k，s1k为2
×
2的矩阵，s1k中每个元素均为复数；采用式(7)计算测试得到的冷夹断偏置状态下gan hemt晶体管的s参数和第一个简化电路模型的仿真s参数之间的总误差ε：
[0021][0022]
其中，为sk中第i行第j列的元素，为s1k中第i行第j列的元素，| |2表示取模后再平方，max(| |)2表示取模平方后再取最大值；
[0023]
然后判断总误差ε是否达到最小值(总误差最小值为0)，如果达到，则最优寄生电容数值循环程序结束，将本次产生的6个寄生电容数值作为最优的一组寄生电容数值；若总误差未达到最小值，则再次根据metropolis准则再次产生6个寄生电容数值进行下一次循环，周而复始，直至总误差最小；
[0024]
步骤s2：提取寄生电感的电感值和寄生电阻的电阻值：
[0025]
步骤s201：在弱前向偏置状态下测试每个频率点处gan hemt晶体管的s参数，得到弱前向偏置状态下各个频率点处gan hemt晶体管的s参数，将弱前向偏置状态下第k个频率点处gan hemt晶体管的s参数记为s2k，s2k为2
×
2的矩阵，s2k中每个元素均为复数，在弱前向偏置状态下，对gan hemt晶体管的小信号等效电路模型进行简化并增加本征电容cg、cs、cd后得到第二个简化电路模型，第二个简化电路模型只包含寄生电感l1、l2、l3，寄生电阻r1、r2、r3和本征电容cg、cs、cd，在弱前向偏置条件下，gan hemt晶体管的栅极电流较高，本征电容cg、cs、cd的电容值非常小，忽略不计；寄生电感l1的一端等效为gan hemt晶体管的栅极，寄生电感l1的另一端和寄生电阻r1的一端连接，寄生电阻r1的另一端和本征电容cg的一端连接，本征电容cg的另一端、本征电容cd的一端和本征电容cs的一端连接，本征电容cs的另一端和寄生电阻r3的一端连接，寄生电阻r3的另一端和寄生电感l3的一端连接，寄生电感l3的另一端等效为gan hemt晶体管的源极，本征电容cd的另一端和寄生电阻r2的一端连接，寄生电阻r2的另一端和寄生电感l2的一端连接，寄生电感l2的另一端等效为gan hemt晶体管的漏极；
[0026]
步骤s202：先从弱前向偏置状态下测试的gan hemt晶体管的s参数剥离步骤s1提取的一组最优寄生电容值(即c
pg
、c
gsi
、c
pgd
、c
gdi
、c
pd
、c
dsi
的最优电容值)，然后转换为z参数，得到弱前向偏置状态下每个频率点处剥离了一组最优寄生电容值的ganhemt晶体管的z参数，将弱前向偏置状态下第k个频率点处剥离了一组最优寄生电容值的gan hemt晶体管的z参数记为zk，zk为2
×
2的矩阵，zk中每个元素均为复数，re()表示复数的实部，im()表示复数的虚部，将zk的第1行第1列元素记为z
11,k
，第1行第2列元素记为z
12,k
，第2行第1列元素记为z
21,k
，第2行第2列元素记为z
22,k
；弱前向偏置状态下剥离了一组最优寄生电容值的z参数与寄生电感l1、l2、l3的电感值和寄生电阻r1、r2、r3的电阻值的关系表示为：
[0027]z11,k
＝r
g,k
r
s,k
jωk(l
g,k
l
s,k
)
ꢀꢀꢀ
(8)
[0028]z12,k
＝r
s,k
jωkl
s,k
ꢀꢀꢀ
(9)
[0029]z22,k
＝r
d,k
r
s,k
jωk(l
d,k
l
s,k
)
ꢀꢀꢀ
(10)
[0030]
其中，r
g,k
表示第k个频率点处rg的取值，r
s,k
表示第k个频率点处rs的取值，r
d,k
表示第k个频率点处rd的取值，l
g,k
表示第k个频率点处lg的取值，l
s,k
表示第k个频率点处ls的取值，l
d,k
表示第k个频率点处ld的取值；
[0031]
公式(8)、(9)、(10)中，在每个频率点处，有等式re(z
11,k
)＝r
g,k
r
s,k
、re(z
12,k
)＝r
s,k
和re(z
22,k
)＝r
d,k
r
s,k
成立，根据上述三个等式可求解出每个频率点处的r
g,k
，r
s,k
和r
d,k
值，将r
g,1
至r
g,n
的平均值作为rg，将r
s,1
至r
s,n
的平均值作为rs，将r
d,1
至r
d,n
的平均值作为rd，进一步降低测试不确定性造成的误差，寄生电阻的电阻值计算公式为：
[0032]
[0033][0034][0035]
其中，re(z
ij,k
)表示zk的第i行第j列元素的实部值，i＝1,2，j＝1,2；
[0036]
步骤s203:求解出寄生电阻后，在式(8)、(9)、(10)两端同时乘以角频率ωk，根据弱前向偏置状态下gan hemt晶体管剥离了最优寄生电容值的z参数虚部拟合ωkim(z
ij
)与ω
k2
的线性关系，即拟合ωk*im(z
12,k
)与ω
k2
、ωk*im(z
11,k
)-ωk*im(z
12,k
)与ω
k2
、ωk*im(z
22,k
)-ωk*im(z
12,k
)与ω
k2
的线性关系，该拟合借助matlab软件实现，直接得到ωk*im(z
12,k
)与ω
k2
、ωk*im(z
11,k
)-ωk*im(z
12,k
)与ω
k2
、ωk*im(z
22,k
)-ωk*im(z
12,k
)与ω
k2
这三个拟合线性直线的斜率，分别做为lg、ls和ld的值，用公式表示为：
[0037]
ω
k im(z
12,k
)＝l
s,k
ω
k2
ꢀꢀꢀ
(14)
[0038]
ω
k im(z
11,k
)-ω
k im(z
12,k
)＝l
g,k
ω
k2
ꢀꢀꢀ
(15)
[0039]
ω
k im(z
22,k
)-ωkim(z
12,k
)＝l
d,k
ω
k2
ꢀꢀꢀ
(16)
[0040]
其中，im(z
ij,k
)表示zk的第i行第j列元素的虚部值，i＝1,2，j＝1,2。
[0041]
步骤s3：提取本征本征电容的电容值、本征电阻的电阻值、本征电导的电导值，具体过程为：
[0042]
步骤s301：在热偏置状态下测试gan hemt晶体管的s参数，得到在热偏置状态下每个频率点处gan hemt晶体管的s参数，将热偏置状态下每个频率点处gan hemt晶体管的s参数剥离寄生电容、寄生电感和寄生电阻，得到剥离了寄生电容、寄生电感和寄生电阻的s参数，将剥离了寄生电容、寄生电感和寄生电阻的s参数转换为y参数，得到热偏置状态下剥离了寄生电容、寄生电感和寄生电阻的gan hemt晶体管每个频率点处的y参数，将热偏置状态下剥离了寄生电容、寄生电感和寄生电阻的gan hemt晶体管第k个频率点处的y参数记为为2
×
2的矩阵，中每个元素均为复数，将的第1行第1列元素记为第1行第2列元素记为第2行第1列元素记为第2行第2列元素记为此时本征电导g1、g2在第k个频率点的值g
gsf,k
、g
gdf,k
用表示为：
[0043][0044][0045]
针对每个频率点对应的本征电导g
gsf,k
、g
gdf,k
，取所有频率点下的该本征电导均值为该热偏置状态下本征电导的终值，由此得到热偏置状态下本征电导g
gsf
、g
gdf
的值；
[0046]
步骤s302：从中剥离本征电导g
gsf
、g
gdf
，得到剥离了本征电导g
gsf
、g
gdf
的y参数，记为为2
×
2的矩阵，的第1行第1列元素记为第1行第2列元素记为第2行第1列元素记为第2行第2列元素记为则剩余8个本征元件在第k个频率点的计算公式为：
[0047][0048][0049][0050][0051][0052][0053][0054][0055]
其中，re()表示复数的实部，im()表示复数的虚部，abs()为取绝对值，c
gs,k
表示第k个频率点处c
gs
的取值，c
gd,k
表示第k个频率点处c
gd
的取值，c
ds,k
表示第k个频率点处c
ds
的取值，g
m,k
表示第k个频率点处gm的取值，g
ds,k
表示第k个频率点处g
ds
的取值，r
i,k
表示第k个频率点处ri的取值，r
gd,k
表示第k个频率点处r
gd
的取值，τk表示第k个频率点处τ的取值；针对每个频率点处本征元件的参数值，取该本征元件所有频率点处取值的均值为该本征元件在热偏置状态下的参数终值，由此得到热偏置状态下每个本征元件的参数终值；
[0056]
步骤s4：将求解的一组最优寄生电容值、寄生电感、寄生电阻和本征元件的值均代入到ads射频仿真软件建立的小信号等效电路模型中，从而得到gan hemt晶体管完整的小信号模型。
[0057]
与现有技术相比，本发明的优点在于在构建gan hemt晶体管小信号模型的过程中，对混合提取法做进一步的改进，仅优化c
pg
、c
gsi
、c
pgd
、c
gdi
、c
pd
、c
dsi
这6个寄生电容值，优化过程简单省时，大大减少了多参数优化导致的误差累积问题；并根据器件的物理特性做出合理的假设(即步骤s103中，在冷夹断偏置状态下，根据gan hemt器件的沟道耗尽层物理特性，假设c
gs
＝c
gd
，将c
gs
＝c
gd
作为最优寄生电容数值循环程序的约束条件)，参数提取结果可靠性较高，使得到的小信号模型参数值具有正确的物理意义，能够准确反映器件的工作状态，能够用于不同结构的器件，可移植性和通用性较高。
附图说明
[0058]
图1为本发明的流程示意图；
[0059]
图2为本发明的matlab寄生电容优化程序流程图；
[0060]
图3为本发明的gan hemt器件小信号等效电路模型示意图；
[0061]
图4为本发明的gan hemt器件在冷夹断状态下等效电路示意图；
[0062]
图5为本发明的gan hemt器件在弱前向偏置状态下等效电路示意图；
[0063]
图6为本发明建立的小信号模型仿真数据与测试数据对比。
具体实施方式
[0064]
以下结合附图实施例对本发明作进一步详细描述。
[0065]
实施例：如图1所示，一种gan hemt晶体管的小信号模型的建模方法，包括以下步骤：
[0066]
步骤s1：提取最优寄生电容值：
[0067]
步骤s101：在ads射频仿真软件中构建gan hemt晶体管的小信号等效电路模型，如图3所示，该小信号等效电路模型包括寄生元件和本征元件，寄生元件包括寄生电容c1、c2、c3、c4、c5、c6，寄生电感l1、l2、l3以及寄生电阻r1、r2、r3；本征元件包括本征电容c7、c8、c9，本征电阻r4、r5，本征电导g1、g2、g3以及本征压控电流源，本征压控电流源的输出电流i
ds
＝v
gsgm
e-jωτ
，其中，gm为gan hemt晶体管的跨导(表示gan hemt晶体管的栅极电压v
gs
对漏极电流i
ds
的控制能力)，τ为gan hemt晶体管栅极电压的时延(gan hemt晶体管的栅极电压的到达时间)，v
gs
为流经本征电容c7两端的电压(也即gan hemt晶体管的栅极电压)，ω为gan hemt晶体管的角频率，j为虚数符号，e为自然对数的底；寄生电容c1的一端、寄生电容c3的一端和寄生电感l1的一端连接且其连接端等效为gan hemt晶体管的栅极，寄生电容c1的另一端、寄生电感l3的一端和寄生电容c5的一端连接且其连接端等效为gan hemt晶体管的源极，寄生电容c2的一端、寄生电容c4的一端、寄生电感l1的另一端和寄生电阻r1的一端连接，寄生电容c2的另一端、寄生电感l3的另一端、寄生电阻r3的一端和寄生电容c6连接，寄生电容c3的另一端、寄生电容c5的另一端和寄生电感l2的一端连接且其连接端等效为gan hemt晶体管的漏极，寄生电容c4的另一端、寄生电阻r2的一端、寄生电感l2的另一端和寄生电容c6的另一端连接，寄生电阻r1的另一端、本征电容c7的一端、本征电容c8的一端、本征电导g1的一端和本征电导g2的一端连接，寄生电阻r2的另一端、本征电容c9的一端、本征电导g2的另一端、本征电导g3的一端、本征电阻r5的一端和本征压控电流源的正极连接，寄生电阻r3的另一端、本征电阻r4的一端、本征电导g1的另一端、本征电容c9的另一端、本征电导g3的另一端和本征压控电流源的负极连接，本征电容c7的另一端与本征电阻r4的另一端连接，本征电容c8的另一端与本征电阻r5的另一端连接；将寄生电容c1的电容值记为c
pg
，寄生电容c2的电容值记为c
gsi
,寄生电容c3的电容值记为c
pgd
,寄生电容c4的电容值记为c
gdi
,寄生电容c5的电容值记为c
pd
,寄生电容c6的电容值记为c
dsi
；寄生电感l1的电感值记为lg，寄生电感l2的电感值记为ld，寄生电感l3的电感值记为ls；寄生电阻r1的电阻值记为rg，寄生电阻r2的电阻值记为rd，寄生电阻r3的电阻值记为rs；本征电容c7的电容值记为c
gs
，本征电容c8的电容值记为c
gd
，本征电容c9的电容值记为c
ds
，本征电阻r4的电阻值记为ri，本征电阻r5的电阻值记为r
gd
，本征电导g1的电导值记为g
gsf
，本征电导g2的电导值记为g
gdf
，本征
电导g3的电导值记为g
ds
；
[0068]
步骤s102：在冷夹断偏置状态下测试gan hemt晶体管的s参数，得到冷夹断偏置状态下每个频率点处gan hemt晶体管的s参数，其中将频率点的总数量记为n，第k个频率点记为fk，k＝1,2,...,n，(频率点总数量与gan hemt晶体管的工作频率范围密切相关，矢量网络分析仪测试s参数时会根据用户设定的频率范围自动生成频率点数，用户也可以根据实际需求设置频率点数，本发明的频率点数采用矢量网络分析仪自动生成的频率点数)将测试得到的冷夹断偏置状态下第k个频率点处gan hemt晶体管的s参数记为sk，sk为2
×
2的矩阵，sk中每个元素均为复数；
[0069]
对gan hemt晶体管的小信号等效电路模型进行简化得到第一个简化电路模型，如图4所示，第一个简化电路模型只包含寄生电容c1、c2、c3、c4、c5、c6和本征电容c7、c8、c9，此时，寄生电容c1的一端、寄生电容c2的一端、寄生电容c5的一端、寄生电容c6的一端、本征电容c7的一端和本征电容c9的一端连接且其连接端等效为gan hemt晶体管的源极，寄生电容c1的另一端、寄生电容c2的另一端、寄生电容c3的一端、寄生电容c4的一端、本征电容c7的另一端和本征电容c8的一端连接且其连接端等效为gan hemt晶体管的栅极，寄生电容c5的另一端、寄生电容c6的另一端、寄生电容c3的另一端、寄生电容c4的另一端、本征电容c8的另一端和本征电容c9的另一端连接且其连接端等效为gan hemt晶体管的漏极；将第一个简化电路模型建立在ads射频仿真软件中；
[0070]
将寄生电容c1、寄生电容c2和本征电容c7的电容值之和作为冷夹断偏置状态下gan hemt晶体管的栅-源分支总电容c
ga
，将寄生电容c3、寄生电容c4和本征电容c8的电容值之和作为冷夹断偏置状态下gan hemt晶体管的栅-漏分支总电容c
gdc
，将寄生电容c5、寄生电容c6和本征电容c9的电容值之和作为冷夹断偏置状态下ganhemt晶体管的漏-源分支总电容c
db
，得到式(1)至式(3)：
[0071]cga
＝c
pg
c
gsi
c
gs
ꢀꢀꢀ
(1)
[0072]cgdc
＝c
pgd
c
gd
c
gdi
ꢀꢀꢀ
(2)
[0073]cdb
＝c
pd
c
dsi
c
ds
ꢀꢀꢀ
(3)
[0074]
步骤s103：将测试得到的冷夹断偏置状态下每个频率点处gan hemt晶体管的s参数分别转化为y参数，得到冷夹断偏置状态下每个频率点处gan hemt晶体管的y参数，将冷夹断偏置状态下第k个频率点处gan hemt晶体管的y参数记为yk，yk为2
×
2的矩阵，yk中每个元素均为复数，将yk的第1行第1列元素记为y
11,k
，第1行第2列元素记为y
12,k
，第2行第1列元素记为y
21,k
，第2行第2列元素记为y
22,k
；此时，栅-源分支总电容c
ga
、栅-漏分支总电容c
gdc
、漏-源分支总电容c
db
与冷夹断偏置状态下gan hemt晶体管的y参数关系式为：
[0075]
im(y
11,k
) im(y
12,k
)＝c
ga,k
ωkꢀꢀꢀ
(4)
[0076]
im(y
12,k
)＝-c
gdc,k
ωkꢀꢀꢀ
(5)
[0077]
im(y
22,k
) im(y
12,k
)＝c
db,k
ωkꢀꢀꢀ
(6)
[0078]
其中，im()表示复数的虚部，ωk＝2*π*fk，ωk为第k个频率点fk对应的角频率；im(y
ij,k
)表示y
ij,k
的虚部值，i＝1,2，j＝1,2；c
ga,k
表示第k个频率点处c
ga
的取值，c
gdc,k
表示第k个频率点处c
gdc
的取值，c
db,k
表示第k个频率点处c
db
的取值；
[0079]
根据式(4)至式(6)拟合冷夹断偏置状态下第k个频率点处gan hemt晶体管的y参数中各元素虚部与角频率ωk的线性关系，即拟合im(y
11,k
) im(y
12,k
)与ωk、im(y
12,k
)与ωk、
hemt晶体管的栅极，寄生电感l1的另一端和寄生电阻r1的一端连接，寄生电阻r1的另一端和本征电容cg的一端连接，本征电容cg的另一端、本征电容cd的一端和本征电容cs的一端连接，本征电容cs的另一端和寄生电阻r3的一端连接，寄生电阻r3的另一端和寄生电感l3的一端连接，寄生电感l3的另一端等效为gan hemt晶体管的源极，本征电容cd的另一端和寄生电阻r2的一端连接，寄生电阻r2的另一端和寄生电感l2的一端连接，寄生电感l2的另一端等效为gan hemt晶体管的漏极；
[0086]
步骤s202：先从弱前向偏置状态下测试的gan hemt晶体管的s参数剥离步骤s1提取的一组最优寄生电容值(即c
pg
、c
gsi
、c
pgd
、c
gdi
、c
pd
、c
dsi
的最优电容值)，然后转换为z参数，得到弱前向偏置状态下每个频率点处剥离了一组最优寄生电容值的ganhemt晶体管的z参数，将弱前向偏置状态下第k个频率点处剥离了一组最优寄生电容值的gan hemt晶体管的z参数记为zk，zk为2
×
2的矩阵，zk中每个元素均为复数，re()表示复数的实部，im()表示复数的虚部，将zk的第1行第1列元素记为z
11,k
，第1行第2列元素记为z
12,k
，第2行第1列元素记为z
21,k
，第2行第2列元素记为z
22,k
；弱前向偏置状态下剥离了一组最优寄生电容值的z参数与寄生电感l1、l2、l3的电感值和寄生电阻r1、r2、r3的电阻值的关系表示为：
[0087]z11,k
＝r
g,k
r
s,k
jωk(l
g,k
l
s,k
)
ꢀꢀꢀ
(8)
[0088]z12,k
＝r
s,k
jωkl
s,k
ꢀꢀꢀ
(9)
[0089]z22,k
＝r
d,k
r
s,k
jωk(l
d,k
l
s,k
)
ꢀꢀꢀ
(10)
[0090]
其中，r
g,k
表示第k个频率点处rg的取值，r
s,k
表示第k个频率点处rs的取值，r
d,k
表示第k个频率点处rd的取值，l
g,k
表示第k个频率点处lg的取值，l
s,k
表示第k个频率点处ls的取值，l
d,k
表示第k个频率点处ld的取值；
[0091]
公式(8)、(9)、(10)中，在每个频率点处，有等式re(z
11,k
)＝r
g,k
r
s,k
、re(z
12,k
)＝r
s,k
和re(z
22,k
)＝r
d,k
r
s,k
成立，根据上述三个等式可求解出每个频率点处的r
g,k
，r
s,k
和r
d,k
值，将r
g,1
至r
g,n
的平均值作为rg，将r
s,1
至r
s,n
的平均值作为rs，将r
d,1
至r
d,n
的平均值作为rd，进一步降低测试不确定性造成的误差，寄生电阻的电阻值计算公式为：
[0092][0093][0094][0095]
其中，re(z
ij,k
)表示zk的第i行第j列元素的实部值，i＝1,2，j＝1,2；
[0096]
步骤s203:求解出寄生电阻后，在式(8)、(9)、(10)两端同时乘以角频率ωk，根据弱前向偏置状态下gan hemt晶体管剥离了最优寄生电容值的z参数虚部拟合ω
k im(z
ij
)与ω
k2
的线性关系，即拟合ωk*im(z
12,k
)与ω
k2
、ωk*im(z
11,k
)-ωk*im(z
12,k
)与ω
k2
、ωk*im(z
22,k
)-ωk*im(z
12,k
)与ω
k2
的线性关系，该拟合借助matlab软件实现，直接得到ωk*im(z
12,k
)与ω
k2
、ωk*im(z
11,k
)-ωk*im(z
12,k
)与ω
k2
、ωk*im(z
22,k
)-ωk*im(z
12,k
)与ω
k2
这三个拟合线性直线的斜率，分别做为lg、ls和ld的值，用公式表示为：
[0097]
ω
k im(z
12,k
)＝l
s,k
ω
k2
ꢀꢀꢀ
(14)
[0098]
ω
k im(z
11,k
)-ω
k im(z
12,k
)＝l
g,k
ω
k2
ꢀꢀꢀ
(15)
[0099]
ω
k im(z
22,k
)-ω
k im(z
12,k
)＝l
d,k
ω
k2
ꢀꢀꢀ
(16)
[0100]
其中，im(z
ij,k
)表示zk的第i行第j列元素的虚部值，i＝1,2，j＝1,2。
[0101]
步骤s3：提取本征本征电容的电容值、本征电阻的电阻值、本征电导的电导值，具体过程为：
[0102]
步骤s301：在热偏置状态下测试gan hemt晶体管的s参数，得到在热偏置状态下每个频率点处gan hemt晶体管的s参数，将热偏置状态下每个频率点处gan hemt晶体管的s参数剥离寄生电容、寄生电感和寄生电阻，得到剥离了寄生电容、寄生电感和寄生电阻的s参数，将剥离了寄生电容、寄生电感和寄生电阻的s参数转换为y参数，得到热偏置状态下剥离了寄生电容、寄生电感和寄生电阻的gan hemt晶体管每个频率点处的y参数，将热偏置状态下剥离了寄生电容、寄生电感和寄生电阻的ganhemt晶体管第k个频率点处的y参数记为为2
×
2的矩阵，中每个元素均为复数，将的第1行第1列元素记为第1行第2列元素记为第2行第1列元素记为第2行第2列元素记为此时本征电导g1、g2在第k个频率点的值g
gsf,k
、g
gdf,k
用表示为：
[0103][0104][0105]
针对每个频率点对应的本征电导g
gsf,k
、g
gdf,k
，取所有频率点下的该本征电导均值为该热偏置状态下本征电导的终值，由此得到热偏置状态下本征电导g
gsf
、g
gdf
的值；
[0106]
步骤s302：从中剥离本征电导g
gsf
、g
gdf
，得到剥离了本征电导g
gsf
、g
gdf
的y参数，记为为2
×
2的矩阵，的第1行第1列元素记为第1行第2列元素记为第2行第1列元素记为第2行第2列元素记为则剩余8个本征元件在第k个频率点的计算公式为：
[0107][0108][0109][0110]
[0111][0112][0113][0114][0115]
其中，re()表示复数的实部，im()表示复数的虚部，abs()为取绝对值，c
gs,k
表示第k个频率点处c
gs
的取值，c
gd,k
表示第k个频率点处c
gd
的取值，c
ds,k
表示第k个频率点处c
ds
的取值，g
m,k
表示第k个频率点处gm的取值，g
ds,k
表示第k个频率点处g
ds
的取值，r
i,k
表示第k个频率点处ri的取值，r
gd,k
表示第k个频率点处r
gd
的取值，τk表示第k个频率点处τ的取值；针对每个频率点处本征元件的参数值，取该本征元件所有频率点处取值的均值为该本征元件在热偏置状态下的参数终值，由此得到热偏置状态下每个本征元件的参数终值；
[0116]
步骤s4：将求解的一组最优寄生电容值、寄生电感、寄生电阻和本征元件的值均代入到ads射频仿真软件建立的小信号等效电路模型中，从而得到gan hemt晶体管完整的小信号模型。
[0117]
在热偏置状态下(选取gan hemt晶体管在栅压vgs＝-2.8v，漏压vds＝28v的偏置状态)，对本发明方法得到的gan hemt晶体管完整的小信号模型进行仿真，得到其仿真s参数，并将该仿真s参数与gan hemt晶体管的实测s参数进行对比，其中对比曲线如图6所示(实线为实测s参数数据，圆圈为仿真s参数数据)。图6中仿真s参数和实测s参数在史密斯原图中具有很好的一致性，验证了本发明提出的ganhemt晶体管小信号模型建模方法的可靠性。