CN107895079B - 全频带超宽带天线电路模型 - Google Patents

全频带超宽带天线电路模型 Download PDF

Info

Publication number
CN107895079B
CN107895079B CN201711119188.1A CN201711119188A CN107895079B CN 107895079 B CN107895079 B CN 107895079B CN 201711119188 A CN201711119188 A CN 201711119188A CN 107895079 B CN107895079 B CN 107895079B
Authority
CN
China
Prior art keywords
basic model
resistor
type
band
radio frequency
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201711119188.1A
Other languages
English (en)
Other versions
CN107895079A (zh
Inventor
唐路
王瀚升
何伟梁
张明辉
曹子建
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Southeast University
Original Assignee
Southeast University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Southeast University filed Critical Southeast University
Priority to CN201711119188.1A priority Critical patent/CN107895079B/zh
Publication of CN107895079A publication Critical patent/CN107895079A/zh
Application granted granted Critical
Publication of CN107895079B publication Critical patent/CN107895079B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/30Circuit design
    • G06F30/36Circuit design at the analogue level
    • G06F30/367Design verification, e.g. using simulation, simulation program with integrated circuit emphasis [SPICE], direct methods or relaxation methods

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Microelectronics & Electronic Packaging (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
  • Secondary Cells (AREA)
  • Measurement Of Resistance Or Impedance (AREA)

Abstract

本发明公开了一种全频带超宽带天线电路模型,由N个A类基本模型单元和M个B类基本模型单元串联而成。基本模型单元总数目K(K=N+M)与测试频段内天线射频阻抗参数实部的频域冲击数目相等。A类基本模型单元数目N与频段划分后各频段Z参数经矢量拟合所得多项式中的电容型多项式数目相等。B类基本模型单元数目M与频段划分后各频段Z参数(射频阻抗参数)实部经矢量拟合所得多项式中电感型多项式数目相等。本发明实现在全频带超宽带条件下,对天线输入阻抗特性的高精度仿真。还公开了基于矢量拟合方法的利用测量射频阻抗参数提取所述全频带超宽带天线模型元件参数值的方法。

Description

全频带超宽带天线电路模型
技术领域
本发明涉及一种天线等效模型及参数提取方法,特别涉及一种全频带超宽带天线等效电路模型及参数提取方法。
背景技术
天线是无线通信系统中射频信号接收与发送的重要元件,在各种不同领域具有广泛应用。在天线和电路系统的协同设计中,根据阻抗匹配原则,常将天线的输入阻抗调整到50欧,来匹配发射机和接收机的其他电路部分。对于窄带天线,这种方法具有很好的可行性,且能达到较高的精度要求。但在全频带超宽带天线应用,由于天线的输入阻抗无法在一段较宽频带上均匹配到50欧,故上述方法对于全频带天线与电路的协同设计不具有可行性。
现阶段,在全频带天线设计的过程中,常采用建立电路模型的方法研究天线的电磁学特性。然而,目前对全频带天线的电路模型建模研究还不是很完善,成为全频带天线和电路协同设计的一个主要难点。
目前国际上常用的天线电路模型,是类似于论文(Shin D R,Lee G,Park WS.Simplified Vector Potential and Circuit Equivalent Model for a Normal-ModeHelical Antenna.IEEE Antennas&Wireless Propagation Letters,2013,12(11):1037-1040.参考文献1)提出的电阻-电容-电感串联电路结构来表征天线的谐振特性,但该电路模型只适用于窄带天线电路建模,并且无法精确地表征天线的阻抗特性。
论文(Badawy M M,El-Azem Malhat H A,Zainud-Deen S H,et al.ASimpleEquivalent Circuit Model for Plasma Dipole Antenna.IEEE Transactions onPlasma Science,2015,43(12):4092-4098.参考文献2),提出了一种包含5个元件的电路模型,同参考文献1相比,增加考虑了辐射电路和辐射电容,相比于参考文献1,该论文的模型精度虽然有了改善,但仍是只适用于窄带天线的建模,无法模拟全频带超宽带天线的输入阻抗特性。
参考文献2利用迭代拟合的方法,利用天线的散射参数对天线电路模型进行参数提取。这种方法具有简单和适用面广等优势,尤其适用于大规模的传统方法求解困难或不可解的问题。然而,迭代拟合的方法仍面临不可避免的收敛问题。由于迭代初值的不同,迭代拟合方法的收敛值可能仅为局部最小值而非全局最小值,由此导致提取结果非最佳结果;或者,由于迭代初值不合理,导致迭代拟合会出现不收敛的情况。
为了避免上述两种方法的缺陷,可以采用矢量拟合方法进行天线电路模型建模。矢量方法核心在于先利用矢量拟合方法这一计算数学方法对所测天线射频阻抗参数或导纳参数进行处理,得到该射频参数对应的经拉普拉斯变换后的多项式。再根据多项式结果推导出拓扑模型,并完成参数提取。目前已有利用矢量拟合方法,进行片上电感等效模型构建和参数提取的文献(Liang Y,Wang Y,Li L.Rational modeling of on-chip inductorby vector fitting[J].Analog Integrated Circuits&Signal Processing,2010,65(2):253-258.参考文献3)。该论文首先利用矢量拟合方法对片上电感实测射频导纳参数值进行处理,提出了一个4元件的等效模型,并实现参数提取。
发明内容
鉴于全频带超宽带天线电路模型在天线电路的协同设计中的重要性,本发明提供一种全频带超宽带天线电路模型以及基于和矢量拟合算法的该电路模型的参数提取方法。本发明针对全频带超宽带天线电路模型构建,利用了矢量拟合方法,结合所测天线等效的射频阻抗参数,提出基于两类基本模型单元新结构的一种新型全频带超宽带天线电路模型,并创新地将矢量拟合的方法运用于天线建模领域,对提出的电路模型所包含各参数元件值的参数提取。本发明对宽带划分子频带,对每个子频带都建模,而且每个子带模型串联构成的总模型,还能很好拟合天线宽带特性。
本发明为解决上述技术问题采用以下技术方案:
本发明提供一种全频带超宽带天线电路模型,包括N个A类基本模型单元和M个B类基本模型单元,各基本模型单元间以串联形式级联,M≥0,N≥0,且M和N不同时为0;其中,A类基本模型单元由以串联形式级联的第一电阻和第一、第二并联结构组成,第一并联结构由第二电阻和第一电容并联组成,第二并联结构由第三电阻和第一电感的串联结构以及第四电阻、第二电容并联组成;B类基本模型单元由以串联形式级联的第五电阻和第三、第四并联结构组成,第三并联结构由第六电阻和第二电感并联组成,第四并联结构由第七电阻和第三电感的串联结构以及第八电阻、第三电容并联组成。
作为本发明的进一步技术方案,每个基本模型单元中的电路元件值不同。
作为本发明的进一步技术方案,基于矢量拟合方法的利用测量射频阻抗参数确定A类基本模型单元的数目N和B类基本模型单元的数目M:
1)对天线进行测试,扫描一组频率,得到不同频率点上天线输入端口射频散射参数;
2)将射频散射参数转换为天线输入端口射频阻抗参数,天线输入端口射频阻抗参数的实部在频域上有若干个极大值,按照极大值的分布将测试频段划分为若干个子频段,每个子频段有且仅有一个极大值;
3)在每个子频段上利用式(1)的电容型多项式和式(2)的电感型多项式,采用矢量拟合的方法对天线输入端口射频阻抗参数和复频率s=jω进行拟合,其中,A1、B1为复数,A* 1、B* 1分别A1、B1的共轭,A2、B2、D值均为实数;
Figure GDA0002957856610000031
Figure GDA0002957856610000032
对于每个子频段,首先采用电容型多项式对天线输入端口射频阻抗参数和复频率s=jω进行拟合,如果所有系数均收敛,则称该子频段收敛于电容型多项式;否则采用电感型多项式对天线输入端口射频阻抗参数和复频率s=jω进行拟合,此时所有系数均收敛,则称该子频段收敛于电感型多项式;
4)基本模型单元的总数目K与测试频段内天线射频阻抗参数实部的频域冲击数目相等,A类基本模型单元的数目N与收敛于电容型多项式的子频段数目相等,B类基本模型单元的数目M与收敛于电感型多项式的子频段数目相等。
作为本发明的进一步技术方案,A类基本模型单元中的电路元件值的确认方法为:
1)对收敛于电容型多项式的子频段,利用式(1)的电容型多项式对天线输入端口射频阻抗参数和复频率s=jω进行拟合,对应于第n个A类基本模型单元的子频段满足式(3):
Figure GDA0002957856610000033
其中,An1、Bn1为复数,A* n1、B* n1分别An1、Bn1的共轭,An2、Bn2、Dn值均为实数,An1、An2、Bn1、Bn2、Dn值均通过矢量拟合得到;
将式(3)改写为式(4):
Figure GDA0002957856610000034
2)根据A类基本模型单元的电路结构,得到第n个A类基本模型单元的射频阻抗参数函数:
Figure GDA0002957856610000035
其中,RAn1、RAn2、RAn3、RAn4、CAn1、CAn2、LAn1分别为第n个A类基本模型单元中第四电阻、第三电阻、第二电阻、第一电阻、第二电容、第一电容、第一电感的元件值;
3)根据式(4)和式(5)的系数对应关系,得到式(6)至式(12):
Figure GDA0002957856610000041
Figure GDA0002957856610000042
Figure GDA0002957856610000043
Figure GDA0002957856610000044
Figure GDA0002957856610000045
Figure GDA0002957856610000046
RAn4=Dn (12)
4)联立式(6)至式(12),计算得到RAn1、RAn2、RAn3、RAn4、CAn1、CAn2、LAn1的值。
作为本发明的进一步技术方案,B类基本模型单元中的电路元件值的确认方法为:
1)对收敛于电感型多项式的子频段,利用式(2)的电感型多项式对天线输入端口射频阻抗参数和复频率s=jω进行拟合,对应于第m个B类基本模型单元的子频段满足式(13):
Figure GDA0002957856610000047
其中,Am1、Bm1为复数,A* m1、B* m1分别Am1、Bm1的共轭,Am2、Bm2、Dm值均为实数,Am1、Am2、Bm1、Bm2、Dm值均通过矢量拟合得到;
将式(13)改写为式(14):
Figure GDA0002957856610000048
2)根据B类基本模型单元的电路结构,得到第m个B类基本模型单元的射频阻抗参数函数:
Figure GDA0002957856610000051
其中,RBm1、RBm2、RBm3、RBm4、CBm1、LBm1、LBm2分别为第m个B类基本模型单元中第八电阻、第七电阻、第六电阻、第五电阻、第三电容、第三电感、第二电感的元件值;
3)根据式(14)和式(15)的系数对应关系,得到式(16)至式(22):
Figure GDA0002957856610000052
Figure GDA0002957856610000053
Figure GDA0002957856610000054
Figure GDA0002957856610000055
RBm3=Bm2(20)
Figure GDA0002957856610000056
RBm4=Dm (22)
4)联立式(16)至式(22),计算得到RBm1、RBm2、RBm3、RBm4、CBm1、LBm1、LBm2的值。
本发明采用以上技术方案与现有技术相比,具有以下技术效果:本发明针对全频带超宽带天线电路模型构建,利用了矢量拟合方法,结合所测天线等效的射频阻抗参数,提出基于两类基本模型单元新结构的一种新型全频带超宽带天线电路模型,并创新地将矢量拟合的方法运用于天线建模领域,对提出的电路模型所包含各参数元件值的参数提取。本发明对宽带划分子频带,对每个子频带都建模,而且每个子带模型串联构成的总模型,还能很好拟合天线宽带特性。本发明所提出的一种适用于全频带超宽带天线电路模型适用于所有全频带超宽带天线,包括但不限于低频单陷波天线、高频单陷波天线、窄带单陷波天线、全频带单陷波天线。
附图说明
图1是基本模型单元数目为K的全频带超宽带天线电路模型的示意图。
图2是适用于全频带超宽带天线电路模型的A类基本模型单元示意图。
图3是适用于全频带超宽带天线电路模型的B类基本模型单元示意图。
图4是基本模型单元数目为5的全频带超宽带天线电路模型的示意图。
图5是超宽带天线Z参数(射频阻抗参数)测量值,其中,(a)是超宽带天线Z参数(射频阻抗参数)测量值的实部,(b)是超宽带天线Z参数(射频阻抗参数)测量值的虚部。
图6是超宽带天线S参数(射频散射参数)测量值与全频带超宽带新型天线电路模型仿真值的比较,其中,(a)是超宽带天线S参数(射频散射参数)测量值实部与全频带超宽带新型天线电路模型仿真值实部的比较,(b)是超宽带天线S参数(射频散射参数)测量值虚部与全频带超宽带新型天线电路模型仿真值虚部的比较。
其中:
1n00-第n个A类基本模型单元;
1n10-第n个A类基本模型单元的并联结构1;
1n20-第n个A类基本模型单元的并联结构2;
1n30-第n个A类基本模型单元的第一电阻(调节电阻)RAn4
1n11-第n个A类基本模型单元并联结构1的第四电阻RAn1
1n12-第n个A类基本模型单元并联结构1的第三电阻RAn2
1n13-第n个A类基本模型单元并联结构1的第二电容CAn1
1n14-第n个A类基本模型单元并联结构1的第一电感LAn1
1n21-第n个A类基本模型单元并联结构2的第二电阻RAn3
1n22-第n个A类基本模型单元并联结构2的第一电容CAn2
2m00-第m个B类基本模型单元;
2m10-第m个B类基本模型单元的并联结构1;
2m20-第m个B类基本模型单元的并联结构2;
2m30-第m个B类基本模型单元的第五电阻(调节电阻)RBm4
2m11-第m个B类基本模型单元并联结构单元1的第八电阻RBm1
2m12-第m个B类基本模型单元并联结构单元1的第七电阻RBm2
2m13-第m个B类基本模型单元并联结构单元1的第三电容CBm1
2m14-第m个B类基本模型单元并联结构单元1的第三电感LBm1
2m21-第m个B类基本模型单元并联结构单元2的第六电阻RBm3
2m22-第m个B类基本模型单元并联结构单元2的第二电感LBm2
Ln-第n个A类基本模型单元左侧结点;
Rn-第n个A类基本模型单元右侧结点;
Lm-第m个B类基本模型单元左侧结点;
Rm-第m个B类基本模型单元右侧结点;
51-全频带超宽带天线Z参数(射频阻抗参数)测试值;
61-全频带超宽带天线S参数(射频散射参数)测量值;
62-全频带超宽带天线电路模型的S参数(射频散射参数)仿真值。
具体实施方式
下面详细描述本发明的实施方式,所述实施方式的示例在附图中示出,其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施方式是示例性的,仅用于解释本发明,而不能解释为对本发明的限制。
本技术领域技术人员可以理解的是,除非另外定义,这里使用的所有术语(包括技术术语和科学术语)具有与本发明所属领域中的普通技术人员的一般理解相同的意义。还应该理解的是,诸如通用字典中定义的那些术语应该被理解为具有与现有技术的上下文中的意义一致的意义,并且除非像这里一样定义,不会用理想化或过于正式的含义来解释。
下面结合附图对本发明的技术方案做进一步的详细说明:
如图1所示,本发明的全频带超宽带天线电路模型,包括N个A类基本模型单元和M个B类基本模型单元,各基本模型单元间以串联形式级联,M≥0,N≥0,且M和N不同时为0。其中基本模型单元总数K(K=M+N)与所测频段内全频带天线Z参数实部的频域冲击个数相等。其中,所测频段指在对天线实物进行实际测试时测试所用的频段;频域冲击指对所测S参数转换后的相应Z参数的实部在频域进行表示后,所测频段内的所有极大值点。A类基本模型单元数目N与频段划分后各频段Z参数实部经矢量拟合处理所得多项式中的电容型多项式数目相等。B类基本模型单元数目M与频段划分后各频段Z参数(射频阻抗参数)实部经矢量拟合方法处理所得多项式中的电感型多项式数目相等。其中,频段划分指根据Z参数(射频阻抗参数)实部的频域冲击,对所测频率范围进行的划分,使每一个频段内有且仅有一个冲击。
电容型多项式指划分频段经矢量拟合方法处理后所得多项式,在拉普拉斯变换后呈现等式(1)型结构或其对应的变换形式,其对应的A类基本模型单元中的双元件并联结构单元为电容电阻并联结构。
Figure GDA0002957856610000071
电感型多项式指划分频段经矢量拟合方法后所得多项式,在拉普拉斯变换后呈现等式(2)型结构或其对应的变换形式,其对应的B类基本模型单元中的双元件并联结构单元为电容电阻并联结构。
Figure GDA0002957856610000081
A类基本模型单元如图2所示,包括并联结构1、并联结构2和调节电阻RAn4,其中,并联结构单元1包括电阻RAn1、RAn2、电容CAn1、电感LAn1,并联结构2包括电阻RAn3、电容CAn2,其中n取1到N的所有整数,N为A类基本模型单元数目。A类基本模型单元由并联结构单元1、并联结构单元2和调节电阻RAn4串联组成位于结构单元左侧结点Ln和右侧结点Rn之间,并与左侧结点Ln和右侧结点Rn相连接。其中n取1到N的所有整数,N为A类基本模型单元数目。并联结构1由电阻RAn1,电感LAn1的串联整体与电阻RAn2,电容CAn1的并联整体并联组成,并联结构1的左侧结点与结构单元左侧结点Ln连接,并联结构1的右侧结点与并联结构2的左侧结点连接。并联结构2由电阻RAn3和电容CAn2并联组成,并联结构2的左侧结点与并联结构1的右侧结点连接,并联结构2的右侧结点与调节电阻RAn4左侧结点连接。调节电阻RAn4由一电阻构成,调节电阻RAn4左侧结点与并联结构2的右侧结点连接,调节电阻RAn4右侧结点与结构单元右侧结点Rn连接。
B类基本模型单元如图3所示,包括并联结构1、并联结构2和调节电阻RBm4,其中,并联结构1包括电阻RBm1、RBm2、电容CBm1、电感LBm1,并联结构单元2包括电阻RBm3、电感LBm2。其中m取0到M的所有整数,M为B类基本模型单元数目。B类基本模型单元由并联结构1、并联结构2和调节电阻RBm4串联组成位于结构单元左侧结点Lm和右侧结点Rm之间,并与左侧结点Lm和右侧结点Rm相连接。其中m取0到M的所有整数,M为B类基本模型单元数目。并联结构1由电阻RBm1,电感LBm1的串联整体与电阻RBm2,电容CBm1的并联整体并联组成,并联结构1的左侧结点与结构单元左侧结点Lm连接,并联结构1的右侧结点与并联结构2的左侧结点连接。并联结构2由电阻RBm3和电感LBm2并联组成,并联结构2的左侧结点与并联结构1的右侧结点连接,并联结构2的右侧结点与调节电阻RBm4左侧结点连接。调节电阻RBm4由一电阻构成,调节电阻RBm4左侧结点与并联结构2的右侧结点连接,调节电阻RBm4右侧结点与结构单元右侧结点Rm连接。
本发明全频带超宽带电路模型由各基本模型单元构成,节点连接方式为:基本模型单元左侧结点和前一个基本模型单元右侧结点连接,基本模型单元右侧结点和后一个基本模型单元左侧结点连接。第一个基本模型单元左侧结点为输入端I,最后一个基本模型单元右侧结点为输出端O。
A类基本模型单元的数目N和B类基本模型单元的数目M的确定方法为:
对天线进行测试,扫描一组频率,从而得到不同频率点上天线输入端口S参数(射频散射参数)。对S参数(射频散射参数)转换为天线输入端口Z参数(射频阻抗参数)。天线输入端口Z参数(射频阻抗参数)的实部在频域上会出现若干个极大值,即频域冲击。按照极大值的分布将测试频段划分为若干个子频段,每个子频段有且仅有一个极大值。在每个子频段上利用电容型多项式((1)式)和电感型多项式((2)式)采用矢量拟合的方法对天线输入端口Z参数(射频阻抗参数)和复频率s=jω进行拟合。
Figure GDA0002957856610000091
Figure GDA0002957856610000092
其中,A1、B1为复数,A* 1、B* 1分别A1、B1的共轭,A2、B2、D值均为实数。
对于每个子频段,首先采用电容型多项式((1)式)对天线输入端口Z参数(射频阻抗参数)和复频率s=jω进行拟合,如果所有系数均收敛,则称该子频段收敛于电容型多项式。否则采用电感型多项式((2)式)对天线输入端口Z参数(射频阻抗参数)和复频率s=jω进行拟合,此时所有系数均收敛,则称该子频段收敛于电感型多项式。
基本模型单元总数目K(即上述M与N之和)与测试频段内天线Z参数(射频阻抗参数)实部的频域冲击数目相等,即与频段划分后子频段的数目相等。A类基本模型单元数目N与收敛于电容型多项式的子频段数目相等;B类基本模型单元数目M与收敛于电感型多项式的子频段数目相等。
第n个A类基本模型单元中的电路元件值RAn1、RAn2、RAn3、RAn4、CAn1、CAn2、LAn1的确定方法为:
对天线进行测试,扫描一组频率,从而得到不同频率点上天线输入端口S参数(射频散射参数)。对S参数(射频散射参数)进行处理得到天线输入端口Z参数(射频阻抗参数)。对测试频段进行划分,使划分后的每个子频段中有且仅有一个Z参数(射频阻抗参数)实部的频域冲击。对收敛于电容型多项式的子频段,利用电容型多项式((1)式)对天线输入端口Z参数(射频阻抗参数)和复频率s=jω进行拟合。对应于第n个A类基本模型单元的子频段满足(3)式。
Figure GDA0002957856610000093
其中,An1、Bn1为复数,A* n1、B* n1分别An1、Bn1的共轭,An2、Bn2、Dn值均为实数。An1、An2、Bn1、Bn2、Dn值均可通过矢量拟合得到。
对(3)式前两项进行通分运算,得到(4)式。
Figure GDA0002957856610000101
根据A类基本模型单元的电路结构,得到第n个A类基本模型单元的Z参数(射频阻抗参数)函数。
Figure GDA0002957856610000102
根据(4)式和(5)式的系数对应关系,得到以下方程组:
Figure GDA0002957856610000103
Figure GDA0002957856610000104
Figure GDA0002957856610000105
Figure GDA0002957856610000106
Figure GDA0002957856610000107
Figure GDA0002957856610000108
RAn4=Dn (12)
联立以上方程组,计算得到RAn1、RAn2、RAn3、RAn4、CAn1、CAn2、LAn1的值。对于每个A类基本模型单元所包括的电路元件值均可以基于对应的子频段的测试数据按照以上过程确定。
第m个B类基本模型单元中的电路元件值RBm1、RBm2、RBm3、RBm4、CBm1、LBm1、LBm2的确定方法为:
对天线进行测试,扫描一组频率,从而得到不同频率点上天线输入端口S参数(射频散射参数)。对S参数(射频散射参数)转换为天线输入端口Z参数(射频阻抗参数)。对测试频段进行划分,使划分后的每个子频段中有且仅有一个Z参数(射频阻抗参数)实部的频域冲击。对收敛于电感型多项式的子频段,利用电感型多项式((2)式)对天线输入端口Z参数(射频阻抗参数)和复频率s=jω进行拟合。对应于第m个B类基本模型单元的子频段满足(13)式。
Figure GDA0002957856610000111
其中,Am1、Bm1为复数,A* m1、B* m1分别Am1、Bm1的共轭,Am2、Bm2、Dm值均为实数。Am1、Am2、Bm1、Bm2、Dm值均可通过矢量拟合得到。
对(13)式前两项进行通分运算,得到(14)式。
Figure GDA0002957856610000112
根据B类基本模型单元的电路结构,得到第m个B类基本模型单元的Z参数(射频阻抗参数)函数。
Figure GDA0002957856610000113
根据(14)式和(15)式的系数对应关系,得到以下方程组。
Figure GDA0002957856610000114
Figure GDA0002957856610000115
Figure GDA0002957856610000116
Figure GDA0002957856610000117
RBm3=Bm2 (20)
Figure GDA0002957856610000118
RBm4=Dm (22)
联立以上方程组,计算得到RBm1、RBm2、RBm3、RBm4、CBm1、LBm1、LBm2的值。对于每个B类基本模型单元所包括的电路元件值均可以基于对应的子频段的测试数据按照以上过程确定。
下面通过具体实施例对本发明的技术方案做进一步阐述:
如图4所示,全频带超宽带天线电路模型的一实施例包括五个基本结构单元,其中,基本结构单元数与所测频段内全频带天线射频阻抗参数Z参数实部的频域冲击个数相等。五个基本结构单元以串联的形式进行级联。五个基本结构单元中包含四个A类基本模型单元(其分别为1100、1200、1300、1400)和一个B类基本模型单元(2100)。
A类基本模型单元1100内部包括并联结构1110、并联结构1120和调节电阻1130(RA14),并联结构1110包括电阻1111(RA11)、电阻1112(RA12)、电容1113(CA11)、电感1114(LA11),并联结构1120包括电阻1121(RA13)、电容1122(CA12)。
A类基本模型单元1200内部包括并联结构1210、并联结构1220和调节电阻1230(RA24),并联结构1210包括电阻1211(RA21)、电阻1212(RA22)、电容1213(CA21)、电感1214(LA21),并联结构1220包括电阻1221(RA23)、电容1222(CA22)。
A类基本模型单元1300内部包括并联结构1310、并联结构1320和调节电阻1330(RA34),并联结构1310包括电阻1311(RA31)、电阻1312(RA32)、电容1313(CA31)、电感1314(LA31),并联结构1320包括电阻1321(RA33)、电容1322(CA32)。
A类基本模型单元1400内部包括并联结构1410、并联结构1420和调节电阻1430(RA44),并联结构1410包括电阻1411(RA41)、电阻1412(RA42)、电容1413(CA41)、电感1414(LA41),并联结构1420包括电阻1421(RA43)、电容1422(CA42)。
B类基本模型单元2100内部包括并联结构2110、并联结构2120和调节电阻2130(RB14),并联结构2110包括电阻2111(RB11)、电阻2112(RB12)、电容2113(CB11)、电感2114(LB11),并联结构1120包括电阻1121(RB13)、电感2122(LB12)。
第n个A类基本模型单元1n00(n取1、2、3、4)的并联结构1n10、并联结构1n20和调节电阻1n30以串联形式级联。其中,并联结构1n10的左侧结点与结构单元左侧结点Ln连接,并联结构1n10的右侧结点与并联结构1n20的左侧结点连接;并联结构1n20的左侧结点与并联结构1n10的右侧结点连接,并联结构1n20的右侧结点与调节电阻1n30左侧结点连接;调节电阻1n30左侧结点与并联结构1n20的右侧结点连接,调节电阻1n30右侧结点与结构单元右侧结点Rn连接。
B类基本模型单元2100的并联结构2110、并联结构2120和调节电阻2130以串联形式级联。其中,并联结构2110的左侧结点与结构单元左侧结点L5连接,并联结构2110的右侧结点与并联结构2120的左侧结点连接;并联结构2120的左侧结点与并联结构2110的右侧结点连接,并联结构2120的右侧结点与调节电阻2130左侧结点连接;调节电阻1n30左侧结点与并联结构2120的右侧结点连接,调节电阻2130右侧结点与结构单元右侧结点R5连接。
本实施例的全频带超宽带天线电路模型由五个基本结构单元串联构成,节点连接方式为:
第1个A类基本模型单元的左侧结点L1和输入端I连接,右侧结点R1与第2个A类基本模型单元的左侧结点L2连接;
第2个A类基本模型单元的左侧结点L2和第1个A类基本模型单元的右侧结点R1连接,右侧结点R2与第3个A类基本模型单元的左侧结点L3连接;
第3个A类基本模型单元的左侧结点L2和第2个A类基本模型单元的右侧结点R2连接,右侧结点R3与第4个A类基本模型单元的左侧结点L4连接;
第4个A类基本模型单元的左侧结点L4和第3个A类基本模型单元的右侧结点R3连接,右侧结点R4与B类基本模型单元的左侧结点L5连接;
B类基本模型单元的左侧结点L5和第4个A类基本模型单元的右侧结点R4连接,右侧结点R5与输出端O连接。
本发明的一种适用于全频带超宽带天线电路模型建模具体包括以下步骤:
步骤一,对天线进行测试,扫描一组频率,得到不同频率点上两端口的Z参数(射频阻抗参数)。
步骤二,根据Z参数(射频阻抗参数)实部的频域冲击,对所测频率范围进行划分,使每一个频段内有且仅有一个冲击。
步骤三,每一个频段对应一个等效电路模型结构单元。根据已测得的不同频率点上两端口的Y参数(射频导纳参数),对每一个频段进行建模,运用特征函数法进行各结构单元的元件参数提取。
步骤四,在完成各结构单元的参数提取后,将各结构单元串联,并利用电路仿真软件进行优化。利用测试曲线与仿真结果之间曲线拟合的方法,迭代优化得到全频带超宽带天线电路模型中的元件值,最终可得该全频带超宽带天线实际电路模型。
具体地,本实施例的超宽带天线等效模型中的电路元件的参数提取方法是:
对超宽带天线进行测试,扫描一组频率(1.1GHz到11GHz),从而得到不同频率点上天线输入端口S参数(射频散射参数)。对S参数(射频散射参数)进行处理得到天线输入端口Z参数(射频阻抗参数)。如图5的(a)和(b)所示,根据处理得到的这一组频率中Z参数(射频阻抗参数)曲线。对扫描频段进行划分,其中1.1GHz到2.9GHz为第一频段,2.9GHz到4GHz为第二频段,4GHz到5.9GHz为第三频段,5.9GHz到7.4GHz为第四频段,7.4GHz到11GHz为第五频段,使划分后的每一频段中有且仅有一个Z参数(射频阻抗参数)实部的频域冲击。对于划分后的任意一个频段,均可以用所述A类或B类基本模型单元进行等效。
对划分后的第一个频段的Z参数(射频阻抗参数)曲线进行矢量拟合方法处理,具体处理方式如下:
首先,由于该频段有且仅有一个Z参数(射频阻抗参数)实部的频域冲击,故其Z参数(射频阻抗参数)可以写成如下形式:
Figure GDA0002957856610000141
其中,s=jω(ω为角频率,j为虚数单位),Ci和pi分别为常数系数值和极点值(i取1到L的所有整数,L为极点数目)。
由于最开始分母中各极点值未知,(23)为非线性方程,较难处理。于是为了简化问题,我们对各极点赋予合理的初值
Figure GDA0002957856610000142
使(23)成为线性方程。接下来,在(23)的基础上乘以一个未知的方程(24)得到等式(25),(24)与(23)有相同的极点初值:
Figure GDA0002957856610000143
Figure GDA0002957856610000144
其中(σZ)fit(s)≈σ(s)Z(s)。将(25)代入(26),得到(27):
Figure GDA0002957856610000145
对S个频域测量值sk(1≤i≤S,S>>L),我们有线性方程:
Akx=bk
其中:
Figure GDA0002957856610000146
bk=Z(sk)。
使用最小二乘法,可以求解上述方程:
Figure GDA0002957856610000151
Figure GDA0002957856610000152
Figure GDA0002957856610000153
其中,zi
Figure GDA0002957856610000154
是(σZ)fit(s)和σ(s)的零点。(30)式表明σ(s)的零点值和z(s)的极点值相等。于是,通过求解σ(s)的零点可以得到z(s)的下一轮迭代的极点初值。反复迭代,当|σ(s)-1|达到精度需要时,可视为σ(s)=1,
Figure GDA0002957856610000155
由此得到射频阻抗参数表达式。
经得矢量拟合方法处理后得到该频段射频阻抗参数(Z参数)的经拉普拉斯变换后的多项式函数。所得多项式函数如下:
Figure GDA0002957856610000156
其中,A1、A2、A3、B1、B* 1、B2、D值均可根据该频段Z参数(射频阻抗参数)计算得出,B1、B* 1共轭。
对(31)式前两项进行通分运算,得:
Figure GDA0002957856610000157
(31)为电容型多项式,则该频段可以用A类基本模型单元等效。根据所述A类基本模型单元结构,得其Z参数(射频阻抗参数)函数:
Figure GDA0002957856610000158
其中比较等式(32)和等式(33),得:
Figure GDA0002957856610000161
Figure GDA0002957856610000162
Figure GDA0002957856610000163
Figure GDA0002957856610000164
Figure GDA0002957856610000165
Figure GDA0002957856610000166
RA14=D
电阻RA11、RA12、RA13、RA14,电容CA11、CA12,电感LA11的值均可由上式联立计算得到。由此可得第1个A类基本模型单元的各元件参数值。
利用以上方法对上述第二、三、五频段进行矢量拟合方法处理,这三个频段均可以用A类基本模型单元等效,利用以上方法分别进行参数提取,可得第2个A类基本模型单元、第3个A类基本模型单元和第4个A类基本模型单元各元件参数值。
对于第四频段Z参数(射频阻抗参数)曲线进行矢量拟合方法处理,得到其Z参数(射频阻抗参数)的经拉普拉斯变换后的多项式函数。所得多项式函数如下:
Figure GDA0002957856610000167
其中,A1、A2、A3、B1、B* 1、B2、D值均可根据该频段射频阻抗参数(Z参数)计算得出,B1、B* 1共轭。
对(34)式前两项进行通分运算,得:
Figure GDA0002957856610000168
(34)为电感型多项式,则该频段可以用B类基本模型单元等效。根据所述B类基本模型单元结构,得其Z参数(射频阻抗参数)函数:
Figure GDA0002957856610000171
其中比较等式(35)和等式(36),得:
Figure GDA0002957856610000172
Figure GDA0002957856610000173
Figure GDA0002957856610000174
Figure GDA0002957856610000175
RB13=B2
Figure GDA0002957856610000176
RB14=D
电阻RB11、RB12、RB13、RB14,电容CB11,电感LB11、LB12值均可由上式联立计算得到。由此可得B类基本模型单元各元件参数值。
将五个基本模型单元串联,并结合等效模型的仿真数据与天线实测数据,进行比较,调整各调平电阻值,使仿真数据曲线与实测数据曲线相近,便得到优化后各调平电阻值。
以下是全频带超宽带天线S参数(射频散射参数)测量值和全频带超宽带天线电路模型的仿真值比较:
由图6中的(a)和(b)可以看出,当采用本发明的超宽带天线等效模型以及超宽带天线等效模型参数提取方法后可以同时实现S参数(射频散射参数)实部和虚部的同时高精度拟合。可以看出,利用本发明的本发明的超宽带天线等效模型以及本发明的超宽带天线等效模型参数提取方法可以实现同测试结果高精度的拟合。
以上所述,仅为本发明中的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉该技术的人在本发明所揭露的技术范围内,可理解想到的变换或替换,都应涵盖在本发明的包含范围之内,因此,本发明的保护范围应该以权利要求书的保护范围为准。

Claims (4)

1.全频带超宽带天线电路模型,其特征在于,包括N个A类基本模型单元和M个B类基本模型单元,各基本模型单元间以串联形式级联,M≥0,N≥0,且M和N不同时为0;其中,A类基本模型单元由以串联形式级联的第一电阻和第一、第二并联结构组成,第一并联结构由第二电阻和第一电容并联组成,第二并联结构由第三电阻和第一电感的串联结构以及第四电阻、第二电容并联组成;B类基本模型单元由以串联形式级联的第五电阻和第三、第四并联结构组成,第三并联结构由第六电阻和第二电感并联组成,第四并联结构由第七电阻和第三电感的串联结构以及第八电阻、第三电容并联组成;
基于矢量拟合方法的利用测量射频阻抗参数确定A类基本模型单元的数目N和B类基本模型单元的数目M:
1)对天线进行测试,扫描一组频率,得到不同频率点上天线输入端口射频散射参数;
2)将射频散射参数转换为天线输入端口射频阻抗参数,天线输入端口射频阻抗参数的实部在频域上有若干个极大值,按照极大值的分布将测试频段划分为若干个子频段,每个子频段有且仅有一个极大值;
3)在每个子频段上利用式(1)的电容型多项式和式(2)的电感型多项式,采用矢量拟合的方法对天线输入端口射频阻抗参数和复频率s=jω进行拟合,其中,A1、B1为复数,A* 1、B* 1分别A1、B1的共轭,A2、B2、D值均为实数;
Figure FDA0002957856600000011
Figure FDA0002957856600000012
对于每个子频段,首先采用电容型多项式对天线输入端口射频阻抗参数和复频率s=jω进行拟合,如果所有系数均收敛,则称该子频段收敛于电容型多项式;否则采用电感型多项式对天线输入端口射频阻抗参数和复频率s=jω进行拟合,此时所有系数均收敛,则称该子频段收敛于电感型多项式;
4)基本模型单元的总数目K与测试频段内天线射频阻抗参数实部的频域冲击数目相等,A类基本模型单元的数目N与收敛于电容型多项式的子频段数目相等,B类基本模型单元的数目M与收敛于电感型多项式的子频段数目相等。
2.根据权利要求1所述的全频带超宽带天线电路模型,其特征在于,每个基本模型单元中的电路元件值不同。
3.根据权利要求1所述的全频带超宽带天线电路模型,其特征在于,A类基本模型单元中的电路元件值的确认方法为:
1)对收敛于电容型多项式的子频段,利用式(1)的电容型多项式对天线输入端口射频阻抗参数和复频率s=jω进行拟合,对应于第n个A类基本模型单元的子频段满足式(3):
Figure FDA0002957856600000021
其中,An1、Bn1为复数,A* n1、B* n1分别An1、Bn1的共轭,An2、Bn2、Dn值均为实数,An1、An2、Bn1、Bn2、Dn值均通过矢量拟合得到;
将式(3)改写为式(4):
Figure FDA0002957856600000022
2)根据A类基本模型单元的电路结构,得到第n个A类基本模型单元的射频阻抗参数函数:
Figure FDA0002957856600000023
其中,RAn1、RAn2、RAn3、RAn4、CAn1、CAn2、LAn1分别为第n个A类基本模型单元中第四电阻、第三电阻、第二电阻、第一电阻、第二电容、第一电容、第一电感的元件值;
3)根据式(4)和式(5)的系数对应关系,得到式(6)至式(12):
Figure FDA0002957856600000024
Figure FDA0002957856600000025
Figure FDA0002957856600000026
Figure FDA0002957856600000027
Figure FDA0002957856600000028
Figure FDA0002957856600000029
RAn4=Dn (12)
4)联立式(6)至式(12),计算得到RAn1、RAn2、RAn3、RAn4、CAn1、CAn2、LAn1的值。
4.根据权利要求1所述的全频带超宽带天线电路模型,其特征在于,B类基本模型单元中的电路元件值的确认方法为:
1)对收敛于电感型多项式的子频段,利用式(2)的电感型多项式对天线输入端口射频阻抗参数和复频率s=jω进行拟合,对应于第m个B类基本模型单元的子频段满足式(13):
Figure FDA0002957856600000031
其中,Am1、Bm1为复数,A* m1、B* m1分别Am1、Bm1的共轭,Am2、Bm2、Dm值均为实数,Am1、Am2、Bm1、Bm2、Dm值均通过矢量拟合得到;
将式(13)改写为式(14):
Figure FDA0002957856600000032
2)根据B类基本模型单元的电路结构,得到第m个B类基本模型单元的射频阻抗参数函数:
Figure FDA0002957856600000033
其中,RBm1、RBm2、RBm3、RBm4、CBm1、LBm1、LBm2分别为第m个B类基本模型单元中第八电阻、第七电阻、第六电阻、第五电阻、第三电容、第三电感、第二电感的元件值;
3)根据式(14)和式(15)的系数对应关系,得到式(16)至式(22):
Figure FDA0002957856600000034
Figure FDA0002957856600000035
Figure FDA0002957856600000036
Figure FDA0002957856600000037
RBm3=Bm2 (20)
Figure FDA0002957856600000041
RBm4=Dm (22)
4)联立式(16)至式(22),计算得到RBm1、RBm2、RBm3、RBm4、CBm1、LBm1、LBm2的值。
CN201711119188.1A 2017-11-14 2017-11-14 全频带超宽带天线电路模型 Active CN107895079B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711119188.1A CN107895079B (zh) 2017-11-14 2017-11-14 全频带超宽带天线电路模型

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711119188.1A CN107895079B (zh) 2017-11-14 2017-11-14 全频带超宽带天线电路模型

Publications (2)

Publication Number Publication Date
CN107895079A CN107895079A (zh) 2018-04-10
CN107895079B true CN107895079B (zh) 2021-04-27

Family

ID=61804379

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711119188.1A Active CN107895079B (zh) 2017-11-14 2017-11-14 全频带超宽带天线电路模型

Country Status (1)

Country Link
CN (1) CN107895079B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110188381B (zh) * 2019-04-18 2023-04-18 中国北方车辆研究所 一种用于电磁干扰预测的仿真模型的构建方法及系统
CN117712670B (zh) * 2023-12-11 2024-07-16 中国人民解放军海军工程大学 一种磁电天线的串并联堆叠阵列的等效电路模型建立方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101237244A (zh) * 2008-03-11 2008-08-06 浙江大学 天线等效电路模型的建立方法
CN101310494A (zh) * 2005-11-14 2008-11-19 艾比奎蒂数字公司 用于am带内同频无线接收器的均衡器
CN104200041A (zh) * 2014-09-17 2014-12-10 电子科技大学 一种宽带特征模的快速提取方法
CN106777483A (zh) * 2016-11-18 2017-05-31 东南大学 用于集成电路的片上电感等效电路模型及参数提取方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20080201658A1 (en) * 2006-06-06 2008-08-21 Ivi Smart Technologies, Inc. Wireless Media Player Device and System, and Method for Operating the Same

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101310494A (zh) * 2005-11-14 2008-11-19 艾比奎蒂数字公司 用于am带内同频无线接收器的均衡器
CN101237244A (zh) * 2008-03-11 2008-08-06 浙江大学 天线等效电路模型的建立方法
CN104200041A (zh) * 2014-09-17 2014-12-10 电子科技大学 一种宽带特征模的快速提取方法
CN106777483A (zh) * 2016-11-18 2017-05-31 东南大学 用于集成电路的片上电感等效电路模型及参数提取方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
A low-phase-noise and low-power crystal oscillator for RF tuner;Tang Lu等;《Journal of Southeast University》;20121231;第28卷(第1期);第21-24页 *
Rational modeling of on-chip inductor by vector fitting;Yuanjun Liang等;《Analog Integr Circ Sig Process》;20101231;第253-258页 *
一种带陷频带可重构超宽带天线的设计;刘健;《延安大学学报(自然科学版)》;20111231;第3O卷(第4期);第79-83页 *

Also Published As

Publication number Publication date
CN107895079A (zh) 2018-04-10

Similar Documents

Publication Publication Date Title
CN107076789B (zh) 测量和确定噪声参数的系统和方法
CN107895079B (zh) 全频带超宽带天线电路模型
CN106777483B (zh) 用于集成电路的片上电感等效电路模型及参数提取方法
Song Non-Foster impedance matching and loading networks for electrically small antennas
CN106533382B (zh) 应用于超高频dc/dc功率变换器的t型阻抗匹配电路及该匹配电路参数设计方法
CN103562736A (zh) 用于感测并计算负载阻抗的简单且微创的方法及系统
CN108509671A (zh) 微波滤波器耦合参数提取方法
Atilla et al. Computer aided Darlington synthesis of an all purpose immittance function
CN108832907A (zh) 数传电台宽带阻抗匹配网络及其设计方法
JP2004235279A (ja) インダクタ素子のシミュレーション方法及びその等価回路
Xia et al. A non-foster matching circuit for an ultra-wideband electrically small monopole antenna
KR100876414B1 (ko) 멀티 레이어 칩 인덕터의 시뮬레이션 장치 및 방법
Goto et al. New on-chip de-embedding for accurate evaluation of symmetric devices
CN113675597A (zh) 一种改善匹配效果的宽带偶极天线及其附加匹配网络结构参数确定方法
Jacob et al. Noise in non-Foster antenna matching circuits
CN111382534B (zh) 一种基于混合时域有限元算法的低频电小系统设计方法
Mukhtar et al. Brune's algorithm for circuit synthesis
Almahroug et al. Antenna matching network using new class of non-Foster reactive elements
CN112394232B (zh) 一种压电晶体、电感器等效参数以及电容容值测量方法
Zhang et al. Rectenna Design with Space Mapping Optimization
Fernandez-Lopez et al. A methodology to design an EMC filter layout providing optimal response based on simulation and considering the inter-component couplings
Fan Multiband RF/Microwave Front-end Receiver Design for Multi-Standard Applications
Liu et al. Design of Broadband Matching Method for Square-law Detection Circuit of Microwave Imaging
Steenput A Spice circuit can be synthesised with a specified set of S-parameters
Mukhtar et al. Brune's synthesis of linear lossy distributed one-port and symmetric two-port microwave circuits

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant