CN108107392B - 多线trl校准方法及终端设备 - Google Patents
多线trl校准方法及终端设备 Download PDFInfo
- Publication number
- CN108107392B CN108107392B CN201711159640.7A CN201711159640A CN108107392B CN 108107392 B CN108107392 B CN 108107392B CN 201711159640 A CN201711159640 A CN 201711159640A CN 108107392 B CN108107392 B CN 108107392B
- Authority
- CN
- China
- Prior art keywords
- calibration
- transmission line
- line
- constant
- error
- 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
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R35/00—Testing or calibrating of apparatus covered by the other groups of this subclass
- G01R35/005—Calibrating; Standards or reference devices, e.g. voltage or resistance standards, "golden" references
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
- Measurement Of Resistance Or Impedance (AREA)
Abstract
本发明提供了多线TRL校准方法及终端设备,该方法包括:对TRL校准过程中的误差进行分析,建立用于求解传播常数和校准常数的误差分析模型;利用多根、冗余的传输线作为标准覆盖每一个频点,根据有效相移规则选取公共线,并将公共线与其它每个传输线组成线对,每组线对之间形成独立测量,并根据所述误差分析模型得到多组传播常数和校准常数的观测值;通过预处理方法对传输线的测量结果进行处理,并根据处理结果更新公共传输线。上述方法及终端设备,能够提高在片S参数测试准确度。
Description
技术领域
本发明属于晶原级半导体器件微波特性测量技术领域,尤其涉及多线TRL校准方法及终端设备。
背景技术
微电子行业中配备的大量“在片S参数测试系统”在使用前,需要使用在片校准件进行矢量校准,校准件的类型包括SOLT(Short-Open-Load-Thru)、TRL(Thru-Reflect-Line)、LRRM(Line-Reflect-Reflect-Match)等。影响在片矢网校准准确度的原因主要有两种:一是系统参考阻抗引入的系统误差,二是探针与被测件接触的重复性误差。
SOLT校准参考阻抗为负载(Load),并设计直通线的特征阻抗与之相等,但事实是SOLT校准覆盖频段较宽,直通线的特征阻抗具有一定的频响,加之短路、开路校准件的定义方式不够完善,导致校准结果精度不高,一直停留在工业应用方面。
TRL校准中使用了易于加工制作的传输线标准,并且校准件的定义中采用长度,其精度得到一定程度提高。TRL校准的参考阻抗为传输线的特征阻抗,在设计传输线尺寸时以仿真结果50欧姆作为标准,但这忽略了仿真模型和频响带来的影响,其校准准确程度受制于传输线特征阻抗与50欧姆接近程度。TRL校准无法得到传输线特征阻抗,也就无法将测量的被测件S参数变换到某一阻抗下。
发明内容
有鉴于此,本发明实施例提供了多线TRL校准方法及终端设备,以解决现有技术中多线TRL校准S参数测量精度较低的问题。
本发明实施例的第一方面提供了一种多线TRL校准方法,包括:
对TRL校准过程中的误差进行分析,建立用于求解传播常数和校准常数的误差分析模型;
利用多根、冗余的传输线作为标准覆盖每一个频点,根据有效相移规则选取公共线,并将公共线与其它每个传输线组成线对,每组线对之间形成独立测量,并根据所述误差分析模型得到多组传播常数和校准常数的观测值;
通过预处理方法对传输线的测量结果进行处理,并根据处理结果更新公共传输线。
可选的,所述对TRL校准过程中的误差进行分析,建立用于求解传播常数和校准常数的误差分析模型的过程为:
矢网测量的第i个校准件的级联传输矩阵Mi为
在理想情况下,第i条传输线标准的传输矩阵Ti为
式中,γ是传播常数,li为第i个传输线标准的长度;
考虑到探针与校准件接触重复性等随机误差,对Ti修正为:
其中,δ1i为端口1不理想引起的随机误差,δ2i为端口2不理想引起的随机误差,且δ1i,δ2i中的元素值远小于1;
给定任意两个传输线标准的测量结果,根据级联传输矩阵的计算公式可得:
MijX=XTij (4)
其中,
Mij≡Mj(Mi)-1 (5)
Tij≡Tj(Ti)-1 (6)
若随机误差δ1i、δ2i不存在,则Tij简化为Lij:
由于随机误差的存在,Tij不再是对角矩阵,Mij的特征值和特征向量也无法直接求解传播常数和校准常数,假定Vij,Λij分别为Tij的特征向量和特征值,则
TijVij=VijΛij (8)
MijUij=UijΛij (9)
其中,Uij为Mij的特征向量;
Uij=XVij (10)
由公式(9)可知,Mij和Tij具有相同的特征值,并且它们特征向量关系如公式(10);
Tij与Lij的一阶线性误差方程为:
Tij≈Lij+εij (11)
根据公式(4),理想Tij的对角元素是Mij,Tij的特征值,Mij特征值为:
其中,δ=0表示δ1i=δ2i=δ1j=δ2j=0;经计算得出:
由公式(17)可得
分配完成后,定义λij为
结合公式(13)~(16)可得到:
其中,λij为Mij特征值;
分析随机误差对传输线传播常数γ的影响,定量给出观测值与待估量传播常数γ和随机误差的关系,并求解测量误差Δγij的协方差矩阵为:
其中,i,m,n分别表示校准件的序号;
其中,
接着求解出Δαij和Δβij的协方差为:
在求取B1和C1/A1时,根据Mij求出四组特征向量,再与估算的B1和C1/A1进行计算,求取差值最小的解为最优解;在得到校准常数X中B1和C1/A1的基础上,根据测量校准件中的一对短路器求解比例系数R1和A1的量值;
对于校准常数Y,求解过程与校准常数X一致。
可选的,所述根据有效相移规则选取公共线,并将公共线与其它每个传输线组成线对,每组线对之间形成独立测量,并根据所述误差分析模型得到多组传播常数和校准常数的观测值为:
确定传输线精确的有效介电常数εeff,计算各个频点下的传播常数γest,定义有效相移θeff为:
在各个频点测量所有传输线标准,根据θeff最大原则确定公共传输线,同时把得到的公共线为最长线和最短线变换为中间长度线;
将所述最终公共传输线与其它传输线标准组成多个线对,并对各个线对进行多次独立测量,根据所述误差分析模型得到多组传播常数和校准常数的观测值。
可选的,在所述通过预处理方法对传输线的测量结果进行处理之前,还包括:
根据线性噪声测量理论,采用最小二乘法评估,对得到的多组传播常数和校准常数的观测值进行处理,得到传播常数和校准常数的最佳估计值。
可选的,所述通过预处理方法对传输线的测量结果进行处理,并根据处理结果更新公共传输线为:
系统校准完成后,对于每个频点,测量任一条传输线标准的S参数;
若测量得到的S参数为奇异解,则重新确定公共传输线。
本发明实施例的第二方面提供了一种多线TRL校准终端设备,包括存储器、处理器,所述存储器中存储有可在所述处理器上运行的计算机程序,所述处理器执行所述计算机程序时实现如下步骤:
对TRL校准过程中的误差进行分析,建立用于求解传播常数和校准常数的误差分析模型;
利用多根、冗余的传输线作为标准覆盖每一个频点,根据有效相移规则选取公共线,并将公共线与其它每个传输线组成线对,每组线对之间形成独立测量,并根据所述误差分析模型得到多组传播常数和校准常数的观测值;
通过预处理方法对传输线的测量结果进行处理,并根据处理结果更新公共传输线。
可选的,所述对TRL校准过程中的误差进行分析,建立用于求解传播常数和校准常数的误差分析模型的过程为:
矢网测量的第i个校准件的级联传输矩阵Mi为
在理想情况下,第i条传输线标准的传输矩阵Ti为
式中,γ是传播常数,li为第i个传输线标准的长度;
考虑到探针与校准件接触重复性等随机误差,对Ti修正为:
其中,δ1i为端口1不理想引起的随机误差,δ2i为端口2不理想引起的随机误差,且δ1i,δ2i中的元素值远小于1;
给定任意两个传输线标准的测量结果,根据级联传输矩阵的计算公式可得:
MijX=XTij (4)
其中,
Mij≡Mj(Mi)-1 (5)
Tij≡Tj(Ti)-1 (6)
若随机误差δ1i、δ2i不存在,则Tij简化为Lij:
由于随机误差的存在,Tij不再是对角矩阵,Mij的特征值和特征向量也无法直接求解传播常数和校准常数,假定Vij,Λij分别为Tij的特征向量和特征值,则
TijVij=VijΛij (8)
MijUij=UijΛij (9)
其中,Uij为Mij的特征向量;
Uij=XVij (10)
由公式(9)可知,Mij和Tij具有相同的特征值,并且它们特征向量关系如公式(10);
Tij与Lij的一阶线性误差方程为:
Tij≈Lij+εij (11)
根据公式(4),理想Tij的对角元素是Mij,Tij的特征值,Mij特征值为:
其中,δ=0表示δ1i=δ2i=δ1j=δ2j=0;经计算得出:
由公式(17)可得
分配完成后,定义λij为
结合公式(13)~(16)可得到:
其中,λij为Mij特征值;
分析随机误差对传输线传播常数γ的影响,定量给出观测值与待估量传播常数γ和随机误差的关系,并求解测量误差Δγij的协方差矩阵为:
其中,i,m,n分别表示校准件的序号;
其中,
接着求解出Δαij和Δβij的协方差为:
在求取B1和C1/A1时,根据Mij求出四组特征向量,再与估算的B1和C1/A1进行计算,求取差值最小的解为最优解;在得到校准常数X中B1和C1/A1的基础上,根据测量校准件中的一对短路器求解比例系数R1和A1的量值;
对于校准常数Y,求解过程与校准常数X一致。
可选的,所述根据有效相移规则选取公共线,并将公共线与其它每个传输线组成线对,每组线对之间形成独立测量,并根据所述误差分析模型得到多组传播常数和校准常数的观测值为:
确定传输线精确的有效介电常数εeff,计算各个频点下的传播常数γest,定义有效相移θeff为:
在各个频点测量所有传输线标准,根据θeff最大原则确定公共传输线,同时把得到的公共线为最长线和最短线变换为中间长度线;
将所述最终公共传输线与其它传输线标准组成多个线对,并对各个线对进行多次独立测量,根据所述误差分析模型得到多组传播常数和校准常数的观测值。
可选的,所述处理器执行所述计算机程序时还实现如下步骤:
根据线性噪声测量理论,采用最小二乘法评估,对得到的多组传播常数和校准常数的观测值进行处理,得到传播常数和校准常数的最佳估计值。
可选的,所述通过预处理方法对传输线的测量结果进行处理,并根据处理结果更新公共传输线为:
系统校准完成后,对于每个频点,测量任一条传输线标准的S参数;
若测量得到的S参数为奇异解,则重新确定公共传输线。
本发明实施例的第三方面提供了一种计算机可读存储介质,所述计算机可读存储介质存储有计算机程序,所述计算机程序被处理器执行时实现如上述任务执行控制方法的步骤。
本发明实施例与现有技术相比存在的有益效果是:本发明实施例,通过对TRL校准过程中的误差进行分析,建立用于求解传播常数和校准常数的误差分析模型,优化了传播常数分配原则,并对公共线的选取进行优化,减少对校准件有效介电常数准确度的敏感度,解决S参数测量结果奇异性的问题,从而提高在片S参数精确测量,能够达到较好的指标,满足市场上商用的在片S参数校准工作。
附图说明
为了更清楚地说明本发明实施例中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1是本发明实施例提供的多线TRL校准方法的流程图;
图2是本发明实施例提供的特征值优化分配流程图;
图3是本发明实施例提供的参数奇异点解决方案的流程图;
图4是本发明实施例提供的传播常数优化前后的测试结果对比;
图5是本发明实施例提供的S参数优化前后的测试结果对比;
图6是本发明实施例提供的对奇异点优化前后的对比;
图7是本发明实施例提供的多线TRL校准程序的运行环境示意图;
图8是本发明实施例提供的多线TRL校准程序的程序模块图。
具体实施方式
以下描述中,为了说明而不是为了限定,提出了诸如特定系统结构、技术之类的具体细节,以便透彻理解本发明实施例。然而,本领域的技术人员应当清楚,在没有这些具体细节的其它实施例中也可以实现本发明。在其它情况中,省略对众所周知的系统、装置、电路以及方法的详细说明,以免不必要的细节妨碍本发明的描述。
为了说明本发明所述的技术方案,下面通过具体实施例来进行说明。
实施例一
图1示出了本发明实施例一提供的多线TRL校准方法的实现流程,详述如下:
步骤S101,对TRL校准过程中的误差进行分析,建立用于求解传播常数和校准常数的误差分析模型。
本步骤中,所述对TRL校准过程中的误差进行分析,建立用于求解传播常数和校准常数的误差分析模型的过程为:
矢网测量的第i个校准件的级联传输矩阵Mi为
其中,Ti为校准件i的实际传输矩阵,X、Y为待求的误差网络传输矩阵,即校准常数; 表示将信号传输方向与Y的信号传输方向反向。例如,若Y表示从端口2到端口1的传输矩阵,则表示从端口1到端口2的传输矩阵。
在理想情况下(如果传输线标准是理想的,并且它与探针的连接也是理想的情况下),第i条传输线标准的传输矩阵Ti为
式中,γ是传播常数,li为第i个传输线标准的长度;
考虑到探针与校准件接触重复性等随机误差,对Ti修正为:
其中,δ1i为端口1不理想引起的随机误差,δ2i为端口2不理想引起的随机误差,且δ1i,δ2i中的元素值远小于1;
给定任意两个传输线标准的测量结果,根据公式(1)可得:
MijX=XTij (4)
其中,
Mij≡Mj(Mi)-1 (5)
Tij≡Tj(Ti)-1 (6)
若随机误差δ1i、δ2i不存在,则Tij简化为Lij:
由于Lij为对角矩阵,由公式(4)可知,此时传播常数和校准常数的求解便转化为矩阵的特征值和特征向量问题,即Mij的特征值对应Tij的对角值 由此可推导出传播常数;Mij的特征向量为X的列向量,由此可推导出校准常数。以上为传统的TRL校准中的求解方法。
由于随机误差的存在,使传播常数和校准常数的求解复杂化。Tij不再是对角矩阵,Mij的特征值和特征向量也无法直接求解传播常数和校准常数。在实际计算过程中,计算Mij的特征值和特征向量是比较容易的,它们与Tij的特征值和特征向量具有一定的关系。假定Vij,Λij分别为Tij的特征向量和特征值,则
TijVij=VijΛij (8)
MijUij=UijΛij (9)
其中,Uij为Mij的特征向量,
Uij=XVij (10)
由公式(9)可知,Mij和Tij具有相同的特征值,并且它们特征向量关系如公式(10);因此,可以通过随机误差对Tij的影响来间接分析对Mij的影响,即对Mij的特征值和特征向量的影响,而Mij的特征值和特征向量分别对应着传播常数和校准常数。
Tij与Lij的一阶线性误差方程为:
Tij≈Lij+εij (11)
根据公式(4),理想Tij的对角元素是Mij,Tij的特征值,Mij特征值为:
其中,δ=0表示δ1i=δ2i=δ1j=δ2j=0;经计算得出:
由公式(17)可得
事实上,考虑到随机误差δ1i、δ2i,Tij不再是对角阵,(13)、(14)与(15)、(16)近似相等。因此特征值对分配,即如何判断 与或相等是传播常数估计的关键,特别是当衰减或相位差引起的测量结果变化相比测量噪声很小时。本实施例中给出的解决方案如图2所示。
分配完成后,定义λij为
结合公式(13)~(16)可得到:
其中,λij为Mij特征值;
分析随机误差对传输线传播常数γ的影响,定量给出观测值与待估量传播常数γ和随机误差的关系,并求解测量误差Δγij的协方差矩阵为:
其中,i,m,n分别表示校准件的序号;
推导出传输线的传播常数后,再按照传统的方法提取出线电容,并计算得到传输线的特征阻抗,从而可将S参数变换到任一阻抗下。
其中,
接着求解出Δαij和Δβij的协方差为:
在求取B1和C1/A1时,同样要面对特征向量的分配,本实施例中采用的方案为:根据Mij求出四组特征向量,再与估算的B1和C1/A1进行计算,求取差值最小的解为最优解。在得到校准常数X中B1和C1/A1的基础上,根据测量校准件中的一对短路器求解比例系数R1和A1的量值;
对于校准常数Y,求解过程与校准常数X一致。
步骤S102,利用多根、冗余的传输线作为标准覆盖每一个频点,根据有效相移规则选取公共线,并将公共线与其它每个传输线组成线对,每组线对之间形成独立测量,并根据所述误差分析模型得到多组传播常数和校准常数的观测值。
本步骤的具体过程为:
确定传输线精确的有效介电常数εeff,包含实部虚部,然后计算各个频点下的传播常数γest,对于低损耗的介质而言,传播常数实部非常小。由公式(29)、(30)可知,对Δαij、Δβij的影响很大,而实际需要的模值尽可能大,这样αij、βij的误差可以尽可能地小,从而使得B1,C1/A1更准确。定义有效相移θeff为:
其中,有效相移θeff越接近90°,越大,测量误差越小。在各个频点测量所有传输线标准,根据θeff最大原则,确定一条传输线为公共传输线标准。在后续计算过程中,以此公共传输线标准与其它传输线标准组成N-1个线对,对其测量结果进行计算。公共线的选取保证N-1个线对对应待估量的N-1次独立测量。结合前面所述传播常数和校准常数的求解方案,可分别得到N-1个观测值。
但在这一步中,NIST(National Institute of Standard Technology)认为,当大于1,是噪声影响,从而把该项设为0。而实际上对于选定公共线而言,有可能接近1,此时加上噪声后其值大于1,把该项设为0存在较大的误判风险,导致该频点在所有传输线为公共线时计算的有效相移都很小,增大测量误差,从而导致后续计算的传播常数、特征阻抗和S参数出现不连续点,极端情况下甚至出现错误。
本实施例中,在各个频点测量所有传输线标准,根据θeff最大原则确定公共传输线,同时把最初得到的最长公共线和最短公共线之间的中间长度的传输线,作为最终公共传输线。将所述最终公共传输线与其它传输线标准组成多个线对,并对各个线对进行多次独立测量,根据所述误差分析模型得到多组传播常数和校准常数的观测值,从而可有效避免传播常数和阻抗不连续导致的测量系统误差。
步骤S103,通过预处理方法对传输线的测量结果进行处理,并根据处理结果更新公共传输线。
本步骤中,所述通过预处理方法对传输线的测量结果进行处理,并根据处理结果更新公共传输线为:
系统校准完成后,对于每个频点,测量任一条传输线标准的S参数;
若测量得到的S参数为奇异解,则重新确定公共传输线。
具体的,系统校准完成后,自动测量计算某一条传输线标准,得到的S参数有可能是奇异解,即其S参数大于1(无源器件四个S参数均小于1),此时采用入图3所述的判断的方式进行处理。对每个频点,测量任一条传输线标准的S参数,若测量得到的Sij参数大于1,则重新确定公共传输线,直至在该频点下测量得到的任一条传输线标准的S参数小于等于1;若测量得到的Sij参数小于等于1,则对下一频点进行上述过程。
可选的,在步骤S103之前,该多线TRL校准方法还可以包括:根据线性噪声测量理论,采用最小二乘法评估,对得到的多组传播常数和校准常数的观测值进行处理,得到传播常数和校准常数的最佳估计值。而步骤S103中,则通过预处理方法对得到传播常数和校准常数的最佳估计值进行处理,并根据处理结果更新公共传输线。
公式(23)、(24)和(25)构成了使用统计方法处理多个传输线标准测量结果的理论基础。设想要求解某一参数x,对anx做了多次测量,得到观测值bn,an为事先已知的常数,则每次测量结果偏离anx约定真值的量值为en:
bn=anx+en (33)
其中,en是一个随机变量,其期望为零,整个测量过程没有系统误差存在。
设定b和a分别为bn和an的列向量,根据Guass-Markou定律,在每次独立测量并不等权的前提下,x的最佳无偏估计为:
x=(σx)2aHV-1b (34)
其中,“aH”表示a的共轭转置变换,en协方差矩阵V定义为:
其中,σx是x的标准偏差,定义为
这种估计是最优的,因为它使得x标准偏差最小。
以传播常数γ为例,可以得到:
x=γ (37)
bij=ln(λij) (38)
aij=li-lj (39)
eij=(li-lj)Δγij (40)
最后,对上述多线TRL校准方法进行了实验,实验结果如图4至图6所示。
上述多线TRL校准方法,优化了传播常数分配原则,对公共线的选取进行优化,减少对校准件有效介电常数准确度的敏感度,解决S参数测量结果奇异性的问题,从而提高在片S参数精确测量,能够达到较好的指标,满足市场上商用的在片S参数校准工作。
应理解,上述实施例中各步骤的序号的大小并不意味着执行顺序的先后,各过程的执行顺序应以其功能和内在逻辑确定,而不应对本发明实施例的实施过程构成任何限定。
实施例二
对应于上文实施例所述的多线TRL校准方法,图7示出了本发明实施例提供的多线TRL校准程序的运行环境示意图。为了便于说明,仅示出了与本实施例相关的部分。
在本实施例中,所述的多线TRL校准程序200安装并运行于终端设备20中。该终端设备20可包括,但不仅限于,存储器201和处理器202。图7仅示出了具有组件201-202的终端设备20,但是应理解的是,并不要求实施所有示出的组件,可以替代的实施更多或者更少的组件。
所述存储器201在一些实施例中可以是所述终端设备20的内部存储单元,例如该终端设备20的硬盘或内存。所述存储器201在另一些实施例中也可以是所述终端设备20的外部存储设备,例如所述终端设备20上配备的插接式硬盘,智能存储卡(Smart MediaCard,SMC),安全数字(Secure Digital,SD)卡,闪存卡(Flash Card)等。进一步地,所述存储器201还可以既包括所述终端设备20的内部存储单元也包括外部存储设备。所述存储器201用于存储安装于所述终端设备20的应用软件及各类数据,例如所述多线TRL校准程序200的程序代码等。所述存储器201还可以用于暂时地存储已经输出或者将要输出的数据。
所述处理器202在一些实施例中可以是一中央处理器(Central ProcessingUnit,CPU),微处理器或其他数据处理芯片,用于运行所述存储器201中存储的程序代码或处理数据,例如执行所述多线TRL校准程序200等。
该终端设备20还可包括显示器,所述显示器在一些实施例中可以是LED显示器、液晶显示器、触控式液晶显示器以及OLED(Organic Light-Emitting Diode,有机发光二极管)触摸器等。
请参阅图8,是本发明实施例提供的多线TRL校准程序200的程序模块图。在本实施例中,所述的多线TRL校准程序200可以被分割成一个或多个模块,所述一个或者多个模块被存储于所述存储器201中,并由一个或多个处理器(本实施例为所述处理器202)所执行,以完成本发明。例如,在图8中,所述的多线TRL校准程序200可以被分割成模型建立模块301、公共传输线选取测量模块302和奇异点处理模块303。本发明所称的模块是指能够完成特定功能的一系列计算机程序指令段,比程序更适合于描述所述多线TRL校准程序200在所述终端设备20中的执行过程。以下描述将具体介绍所述模块301-303的功能。
其中,模型建立模块301,用于对TRL校准过程中的误差进行分析,建立用于求解传播常数和校准常数的误差分析模型。
公共传输线选取测量模块302,用于利用多根、冗余的传输线作为标准覆盖每一个频点,根据有效相移规则选取公共线,并将公共线与其它每个传输线组成线对,每组线对之间形成独立测量,并根据所述误差分析模型得到多组传播常数和校准常数的观测值。
奇异点处理模块303,用于通过预处理方法对传输线的测量结果进行处理,并根据处理结果更新公共传输线。
作为一种可实施方式,模型建立模块301对TRL校准过程中的误差进行分析,建立用于求解传播常数和校准常数的误差分析模型的过程为:
矢网测量的第i个校准件的级联传输矩阵Mi为
在理想情况下,第i条传输线标准的传输矩阵Ti为
式中,γ是传播常数,li为第i个传输线标准的长度;
考虑到探针与校准件接触重复性等随机误差,对Ti修正为:
其中,δ1i为端口1不理想引起的随机误差,δ2i为端口2不理想引起的随机误差,且δ1i,δ2i中的元素值远小于1;
给定任意两个传输线标准的测量结果,根据级联传输矩阵的计算公式可得:
MijX=XTij (4)
其中,
Mij≡Mj(Mi)-1 (5)
Tij≡Tj(Ti)-1 (6)
若随机误差δ1i、δ2i不存在,则Tij简化为Lij:
由于随机误差的存在,Tij不再是对角矩阵,Mij的特征值和特征向量也无法直接求解传播常数和校准常数,假定Vij,Λij分别为Tij的特征向量和特征值,则
TijVij=VijΛij (8)
MijUij=UijΛij (9)
其中,Uij为Mij的特征向量。
Uij=XVij (10)
由公式(9)可知,Mij和Tij具有相同的特征值,并且它们特征向量关系如公式(10);
Tij与Lij的一阶线性误差方程为:
Tij≈Lij+εij (11)
根据公式(4),理想Tij的对角元素是Mij,Tij的特征值,Mij特征值为:
其中,δ=0表示δ1i=δ2i=δ1j=δ2j=0;经计算得出:
由公式(17)可得
分配完成后,定义λij为
结合公式(13)~(16)可得到:
其中,λij为Mij特征值;
分析随机误差对传输线传播常数γ的影响,定量给出观测值与待估量传播常数γ和随机误差的关系,并求解测量误差Δγij的协方差矩阵为:
其中,i,m,n分别表示校准件的序号;
其中,
接着求解出Δαij和Δβij的协方差为:
在求取B1和C1/A1时,根据Mij求出四组特征向量,再与估算的B1和C1/A1进行计算,求取差值最小的解为最优解;在得到校准常数X中B1和C1/A1的基础上,根据测量校准件中的一对短路器求解比例系数R1和A1的量值;
对于校准常数Y,求解过程与校准常数X一致。
作为另一种可实施方式,公共传输线选取测量模块302具体用于:
确定传输线精确的有效介电常数εeff,计算各个频点下的传播常数γest,定义有效相移θeff为:
在各个频点测量所有传输线标准,根据θeff最大原则确定公共传输线,同时把得到的公共线为最长线和最短线变换为中间长度线。
将所述最终公共传输线与其它传输线标准组成多个线对,并对各个线对进行多次独立测量,根据所述误差分析模型得到多组传播常数和校准常数的观测值。
可选的,多线TRL校准程序200还可以被分割出优化模块。所述优化模块,用于根据线性噪声测量理论,采用最小二乘法评估,对得到的多组传播常数和校准常数的观测值进行处理,得到传播常数和校准常数的最佳估计值。所述优化模块得出的最佳估计值传输给奇异点处理模块303。奇异点处理模块303通过预处理方法对优化模块发来的传输线的测量结果进行处理,并根据处理结果更新公共传输线。
进一步的,奇异点处理模块303具体用于:
系统校准完成后,对于每个频点,测量任一条传输线标准的S参数;
若测量得到的S参数为奇异解,则重新确定公共传输线。
所属领域的技术人员可以清楚地了解到,为了描述的方便和简洁,仅以上述各功能单元、模块的划分进行举例说明,实际应用中,可以根据需要而将上述功能分配由不同的功能单元、模块完成,即将所述装置的内部结构划分成不同的功能单元或模块,以完成以上描述的全部或者部分功能。实施例中的各功能单元、模块可以集成在一个处理单元中,也可以是各个单元单独物理存在,也可以两个或两个以上单元集成在一个单元中,上述集成的单元既可以采用硬件的形式实现,也可以采用软件功能单元的形式实现。另外,各功能单元、模块的具体名称也只是为了便于相互区分,并不用于限制本申请的保护范围。上述系统中单元、模块的具体工作过程,可以参考前述方法实施例中的对应过程,在此不再赘述。
在上述实施例中,对各个实施例的描述都各有侧重,某个实施例中没有详述或记载的部分,可以参见其它实施例的相关描述。
本领域普通技术人员可以意识到,结合本文中所公开的实施例描述的各示例的单元及算法步骤,能够以电子硬件、或者计算机软件和电子硬件的结合来实现。这些功能究竟以硬件还是软件方式来执行,取决于技术方案的特定应用和设计约束条件。专业技术人员可以对每个特定的应用来使用不同方法来实现所描述的功能,但是这种实现不应认为超出本发明的范围。
在本发明所提供的实施例中,应该理解到,所揭露的装置/终端设备和方法,可以通过其它的方式实现。例如,以上所描述的装置/终端设备实施例仅仅是示意性的,例如,所述模块或单元的划分,仅仅为一种逻辑功能划分,实际实现时可以有另外的划分方式,例如多个单元或组件可以结合或者可以集成到另一个系统,或一些特征可以忽略,或不执行。另一点,所显示或讨论的相互之间的耦合或直接耦合或通讯连接可以是通过一些接口,装置或单元的间接耦合或通讯连接,可以是电性,机械或其它的形式。
所述作为分离部件说明的单元可以是或者也可以不是物理上分开的,作为单元显示的部件可以是或者也可以不是物理单元,即可以位于一个地方,或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本实施例方案的目的。
另外,在本发明各个实施例中的各功能单元可以集成在一个处理单元中,也可以是各个单元单独物理存在,也可以两个或两个以上单元集成在一个单元中。上述集成的单元既可以采用硬件的形式实现,也可以采用软件功能单元的形式实现。
所述集成的模块/单元如果以软件功能单元的形式实现并作为独立的产品销售或使用时,可以存储在一个计算机可读取存储介质中。基于这样的理解,本发明实现上述实施例方法中的全部或部分流程,也可以通过计算机程序来指令相关的硬件来完成,所述的计算机程序可存储于一计算机可读存储介质中,该计算机程序在被处理器执行时,可实现上述各个方法实施例的步骤。其中,所述计算机程序包括计算机程序代码,所述计算机程序代码可以为源代码形式、对象代码形式、可执行文件或某些中间形式等。所述计算机可读介质可以包括:能够携带所述计算机程序代码的任何实体或装置、记录介质、U盘、移动硬盘、磁碟、光盘、计算机存储器、只读存储器(ROM,Read-Only Memory)、随机存取存储器(RAM,Random Access Memory)、电载波信号、电信信号以及软件分发介质等。需要说明的是,所述计算机可读介质包含的内容可以根据司法管辖区内立法和专利实践的要求进行适当的增减,例如在某些司法管辖区,根据立法和专利实践,计算机可读介质不包括电载波信号和电信信号。
以上所述实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围,均应包含在本发明的保护范围之内。
Claims (8)
1.一种多线TRL校准方法,其特征在于,包括:
对TRL校准过程中的误差进行分析,建立用于求解传播常数和校准常数的误差分析模型;
利用多根、冗余的传输线作为标准覆盖每一个频点,根据有效相移规则选取公共传输线,并将公共传输线与其它每个传输线组成线对,每组线对之间形成独立测量,并根据所述误差分析模型得到多组传播常数和校准常数的观测值;
通过预处理方法对传输线的测量结果进行处理,并根据处理结果更新公共传输线;
其中,所述通过预处理方法对传输线的测量结果进行处理,并根据处理结果更新公共传输线包括:对每个频点,测量任一条传输线的S参数,若测量得到的S参数大于1,则重新确定公共传输线,直至在该频点下测量得到的任一条传输线的S参数小于等于1;若测量得到的S参数小于等于1,则对下一频点进行上述过程。
2.如权利要求1所述的多线TRL校准方法,其特征在于,所述对TRL校准过程中的误差进行分析,建立用于求解传播常数和校准常数的误差分析模型的过程为:
矢网测量的第i个校准件的级联传输矩阵Mi为
在理想情况下,第i条传输线的传输矩阵Ti为
式中,γ是传播常数,li为第i个传输线的长度;
考虑到探针与校准件接触重复性等随机误差,对Ti修正为:
其中,δ1i为端口1不理想引起的随机误差,δ2i为端口2不理想引起的随机误差,且δ1i,δ2i中的元素值远小于1;
给定任意两个传输线的测量结果,根据公式(1)可得:
MijX=XTij(4)
其中,
Mij≡Mj(Mi)-1(5)
Tij≡Tj(Ti)-1(6)
若随机误差δ1i、δ2i不存在,则Tij简化为Lij:
由于随机误差的存在,Tij不再是对角矩阵,Mij的特征值和特征向量也无法直接求解传播常数和校准常数,假定Vij,Λij分别为Tij的特征向量和特征值,则
TijVij=VijΛij(8)
MijUij=UijΛij(9)
其中,Uij为Mij的特征向量,
Uij=XVij(10)
由公式(9)可知,Mij和Tij具有相同的特征值,并且它们特征向量关系如公式(10);
Tij与Lij的一阶线性误差方程为:
Tij≈Lij+εij(11)
根据公式(4),理想Tij的对角元素是Mij,Tij的特征值,Mij特征值为:
其中,δ=0表示δ1i=δ2i=δ1j=δ2j=0;经计算得出:
由公式(17)可得
分配完成后,定义λij为
结合公式(13)~(16)可得到:
其中,λij为Mij特征值;
分析随机误差对传输线传播常数γ的影响,定量给出观测值与待估量传播常数γ和随机误差的关系,并求解测量误差Δγij的协方差矩阵为:
其中,i,m,n分别表示校准件的序号;
其中,
接着求解出Δαij和Δβij的协方差为:
在求取B1和C1/A1时,根据Mij求出四组特征向量,再与估算的B1和C1/A1进行计算,求取差值最小的解为最优解;在得到校准常数X中B1和C1/A1的基础上,根据测量校准件中的一对短路器求解比例系数R1和A1的量值;
对于校准常数Y,求解过程与校准常数X一致。
3.如权利要求1所述的多线TRL校准方法,其特征在于,所述根据有效相移规则选取公共传输线,并将公共传输线与其它每个传输线组成线对,每组线对之间形成独立测量,并根据所述误差分析模型得到多组传播常数和校准常数的观测值为:
确定传输线精确的有效介电常数εeff,计算各个频点下的传播常数γest,定义有效相移θeff为:
在各个频点测量所有传输线,根据θeff最大原则确定公共传输线,同时把得到的公共传输线为最长线和最短线变换为中间长度线;
将所述公共传输线与其它传输线组成多个线对,并对各个线对进行多次独立测量,根据所述误差分析模型得到多组传播常数和校准常数的观测值。
4.如权利要求1所述的多线TRL校准方法,其特征在于,在所述通过预处理方法对传输线的测量结果进行处理之前,还包括:
根据线性噪声测量理论,采用最小二乘法评估,对得到的多组传播常数和校准常数的观测值进行处理,得到传播常数和校准常数的最佳估计值。
5.一种多线TRL校准终端设备,其特征在于,包括存储器、处理器,所述存储器中存储有可在所述处理器上运行的计算机程序,所述处理器执行所述计算机程序时实现如下步骤:
对TRL校准过程中的误差进行分析,建立用于求解传播常数和校准常数的误差分析模型;
利用多根、冗余的传输线作为标准覆盖每一个频点,根据有效相移规则选取公共传输线,并将公共传输线与其它每个传输线组成线对,每组线对之间形成独立测量,并根据所述误差分析模型得到多组传播常数和校准常数的观测值;
通过预处理方法对传输线的测量结果进行处理,并根据处理结果更新公共传输线;
其中,所述通过预处理方法对传输线的测量结果进行处理,并根据处理结果更新公共传输线包括:对每个频点,测量任一条传输线的S参数,若测量得到的S参数大于1,则重新确定公共传输线,直至在该频点下测量得到的任一条传输线的S参数小于等于1;若测量得到的S参数小于等于1,则对下一频点进行上述过程。
6.如权利要求5所述的多线TRL校准终端设备,其特征在于,所述对TRL校准过程中的误差进行分析,建立用于求解传播常数和校准常数的误差分析模型的过程为:
矢网测量的第i个校准件的级联传输矩阵Mi为
在理想情况下,第i条传输线的传输矩阵Ti为
式中,γ是传播常数,li为第i个传输线的长度;
考虑到探针与校准件接触重复性等随机误差,对Ti修正为:
其中,δ1i为端口1不理想引起的随机误差,δ2i为端口2不理想引起的随机误差,且δ1i,δ2i中的元素值远小于1;
给定任意两个传输线的测量结果,根据级联传输矩阵的计算公式可得:
MijX=XTij(4)
其中,
Mij≡Mj(Mi)-1(5)
Tij≡Tj(Ti)-1(6)
若随机误差δ1i、δ2i不存在,则Tij简化为Lij:
由于随机误差的存在,Tij不再是对角矩阵,Mij的特征值和特征向量也无法直接求解传播常数和校准常数,假定Vij,Λij分别为Tij的特征向量和特征值,则
TijVij=VijΛij(8)
MijUij=UijΛij(9)
其中,Uij为Mij的特征向量;
Uij=XVij(10)
由公式(9)可知,Mij和Tij具有相同的特征值,并且它们特征向量关系如公式(10);
Tij与Lij的一阶线性误差方程为:
Tij≈Lij+εij(11)
根据公式(4),理想Tij的对角元素是Mij,Tij的特征值,Mij特征值为:
其中,δ=0表示δ1i=δ2i=δ1j=δ2j=0;经计算得出:
由公式(17)可得
分配完成后,定义λij为
结合公式(13)~(16)可得到:
其中,λij为Mij特征值;
分析随机误差对传输线传播常数γ的影响,定量给出观测值与待估量传播常数γ和随机误差的关系,并求解测量误差Δγij的协方差矩阵为:
其中,i,m,n分别表示校准件的序号;
其中,
接着求解出Δαij和Δβij的协方差为:
在求取B1和C1/A1时,根据Mij求出四组特征向量,再与估算的B1和C1/A1进行计算,求取差值最小的解为最优解;在得到校准常数X中B1和C1/A1的基础上,根据测量校准件中的一对短路器求解比例系数R1和A1的量值;
对于校准常数Y,求解过程与校准常数X一致。
7.如权利要求5所述的多线TRL校准终端设备,其特征在于,所述根据有效相移规则选取公共传输线,并将公共传输线与其它每个传输线组成线对,每组线对之间形成独立测量,并根据所述误差分析模型得到多组传播常数和校准常数的观测值为:
确定传输线精确的有效介电常数εeff,计算各个频点下的传播常数γest,定义有效相移θeff为:
在各个频点测量所有传输线,根据θeff最大原则确定公共传输线,同时把得到的公共传输线为最长线和最短线变换为中间长度线;
将所述公共传输线与其它传输线组成多个线对,并对各个线对进行多次独立测量,根据所述误差分析模型得到多组传播常数和校准常数的观测值。
8.一种计算机可读存储介质,所述计算机可读存储介质存储有计算机程序,其特征在于,所述计算机程序被处理器执行时实现如权利要求1至4任一项所述方法的步骤。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711159640.7A CN108107392B (zh) | 2017-11-20 | 2017-11-20 | 多线trl校准方法及终端设备 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711159640.7A CN108107392B (zh) | 2017-11-20 | 2017-11-20 | 多线trl校准方法及终端设备 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108107392A CN108107392A (zh) | 2018-06-01 |
CN108107392B true CN108107392B (zh) | 2020-03-24 |
Family
ID=62206502
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711159640.7A Active CN108107392B (zh) | 2017-11-20 | 2017-11-20 | 多线trl校准方法及终端设备 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108107392B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109270479B (zh) * | 2018-09-21 | 2021-03-19 | 临海市云谱光电有限公司 | 一种基于张量分解的多线trl校准方法 |
CN112098793B (zh) * | 2020-08-14 | 2023-02-28 | 中国电子科技集团公司第十三研究所 | 单端口在片校准件模型确定的方法及终端设备 |
CN112098794B (zh) * | 2020-08-14 | 2023-02-28 | 中国电子科技集团公司第十三研究所 | 在片校准件模型中参数确定的方法及终端设备 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH10197577A (ja) * | 1997-01-10 | 1998-07-31 | Kyocera Corp | 高周波測定の校正標準器および校正法ならびに高周波用伝送線路の伝送損失の測定方法 |
US7030625B1 (en) * | 2005-01-18 | 2006-04-18 | Agilent Technologies, Inc. | Method and apparatus for performing a minimum connection multiport through-reflect-line calibration and measurement |
CN101339213A (zh) * | 2008-08-07 | 2009-01-07 | 中国电子科技集团公司第十研究所 | 毫米波频段准确测量波导宽边缝隙导纳的方法 |
CN104502878A (zh) * | 2014-12-26 | 2015-04-08 | 中国电子科技集团公司第十三研究所 | 微波GaAs衬底在片S参数微带线TRL校准件 |
CN204666731U (zh) * | 2015-06-12 | 2015-09-23 | 中国电子科技集团公司第十三研究所 | 太赫兹在片散射参数测量校准件 |
CN105891759A (zh) * | 2016-05-10 | 2016-08-24 | 北京无线电计量测试研究所 | 一种传输线标准器件的确定方法、设备以及校准设备 |
CN106383327A (zh) * | 2016-08-26 | 2017-02-08 | 工业和信息化部电子工业标准化研究院 | 一种微波器件标准样片的校准方法 |
CN106405462A (zh) * | 2016-08-30 | 2017-02-15 | 中国电子科技集团公司第十三研究所 | 在片散射参数的溯源及不确定度评估方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7068046B2 (en) * | 2004-11-18 | 2006-06-27 | Anritsu Company | Calibration techniques for simplified high-frequency multiport differential measurements |
-
2017
- 2017-11-20 CN CN201711159640.7A patent/CN108107392B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH10197577A (ja) * | 1997-01-10 | 1998-07-31 | Kyocera Corp | 高周波測定の校正標準器および校正法ならびに高周波用伝送線路の伝送損失の測定方法 |
US7030625B1 (en) * | 2005-01-18 | 2006-04-18 | Agilent Technologies, Inc. | Method and apparatus for performing a minimum connection multiport through-reflect-line calibration and measurement |
CN101339213A (zh) * | 2008-08-07 | 2009-01-07 | 中国电子科技集团公司第十研究所 | 毫米波频段准确测量波导宽边缝隙导纳的方法 |
CN104502878A (zh) * | 2014-12-26 | 2015-04-08 | 中国电子科技集团公司第十三研究所 | 微波GaAs衬底在片S参数微带线TRL校准件 |
CN204666731U (zh) * | 2015-06-12 | 2015-09-23 | 中国电子科技集团公司第十三研究所 | 太赫兹在片散射参数测量校准件 |
CN105891759A (zh) * | 2016-05-10 | 2016-08-24 | 北京无线电计量测试研究所 | 一种传输线标准器件的确定方法、设备以及校准设备 |
CN106383327A (zh) * | 2016-08-26 | 2017-02-08 | 工业和信息化部电子工业标准化研究院 | 一种微波器件标准样片的校准方法 |
CN106405462A (zh) * | 2016-08-30 | 2017-02-15 | 中国电子科技集团公司第十三研究所 | 在片散射参数的溯源及不确定度评估方法 |
Non-Patent Citations (6)
Title |
---|
A Multiline Method of Network Analyzer Calibration;Roger B.Marks et.al;《IEEE TRANSACTIONS ON MICROWAVE THEORY AND TECHNIQUES》;19910731;第39卷(第7期);第1205-1215页 * |
An improved Multiline TRL method;J. E. ZUfiiga-Juarez et.al;《67th ARFTG Conference》;20061231;第139-142页 * |
Multi-Frequency Approach to the Coaxial Multiline Through-Reflect-Line Calibration;Arkadiusz Lewandowski et.al;《18-th INTERNATIONAL CONFERENCE ON MICROWAVES, RADAR AND WIRELESS COMMUNICATIONS》;20101231;正文第1-4页 * |
关于S参数;汪进进;《中国集成电路》;20151231(第9期);第61-65页 * |
基于Multi_TRL算法的传输线特征阻抗定标;王一帮 等;《计量学报》;20170331;第38卷(第2期);第225-229页 * |
怎样设计和验证TRL 校准件及具体过程;刘迪;《电子产品世界》;20081231(第3期);第123-126页 * |
Also Published As
Publication number | Publication date |
---|---|
CN108107392A (zh) | 2018-06-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110286345B (zh) | 一种矢量网络分析仪在片s参数的校准方法、系统及设备 | |
CN109444721B (zh) | 检测s参数的方法及终端设备 | |
US11340286B2 (en) | On-wafer S-parameter calibration method | |
CN108107392B (zh) | 多线trl校准方法及终端设备 | |
CN111142057B (zh) | 太赫兹频段在片s参数的校准方法及终端设备 | |
US6823276B2 (en) | System and method for determining measurement errors of a testing device | |
CN112098791B (zh) | 在片校准件模型及在片校准件模型中参数确定的方法 | |
CN108646208A (zh) | 一种多端口夹具自动去嵌入方法 | |
US20080195344A1 (en) | Method for determining measurement errors in scattering parameter measurements | |
Peng | Estimating and testing quantile-based process capability indices for processes with skewed distributions | |
CN107991537A (zh) | 基于共面波导传输线的介电常数提取方法及终端设备 | |
CN108563895A (zh) | 一种考虑相关性的区间模型修正方法 | |
CN112526470A (zh) | 标定雷达参数的方法和装置、电子设备、存储介质 | |
CN115684780A (zh) | 电磁信号测量方法、装置、计算机设备和存储介质 | |
CN113849958A (zh) | 在片s参数测量系统串扰误差修正方法及电子设备 | |
CN111579869A (zh) | 互易二端口网络s参数测量方法、装置及终端设备 | |
CN110955862B (zh) | 一种用于设备模型趋势相似度的评估方法及装置 | |
CN115728445A (zh) | 传感器校准方法、装置、设备及计算机可读存储介质 | |
CN109444786B (zh) | 提高在片负载牵引测量准确度的方法、系统及终端设备 | |
CN111025214B (zh) | 获取功率校准模型的方法及终端设备 | |
CN111027190A (zh) | 一种用于设备模型数值相似度的评估方法及装置 | |
JP7153309B2 (ja) | ベクトルネットワークアナライザを用いた反射係数の測定方法 | |
CN108362650B (zh) | 液体色度检测方法及装置 | |
CN116306890A (zh) | 一种反演蜂窝等效介电常数的神经网络的训练方法 | |
CN109270479B (zh) | 一种基于张量分解的多线trl校准方法 |
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 |