CN107203495B - 一种s参数圆形插值方法 - Google Patents
一种s参数圆形插值方法 Download PDFInfo
- Publication number
- CN107203495B CN107203495B CN201710477870.1A CN201710477870A CN107203495B CN 107203495 B CN107203495 B CN 107203495B CN 201710477870 A CN201710477870 A CN 201710477870A CN 107203495 B CN107203495 B CN 107203495B
- Authority
- CN
- China
- Prior art keywords
- interpolation
- interpolation point
- parameter
- parameters
- circular
- 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
- 238000000034 method Methods 0.000 title claims abstract description 23
- 238000005259 measurement Methods 0.000 claims description 15
- 238000004364 calculation method Methods 0.000 claims description 4
- 238000010586 diagram Methods 0.000 claims 1
- 230000009466 transformation Effects 0.000 claims 1
- 238000006243 chemical reaction Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/17—Function evaluation by approximation methods, e.g. inter- or extrapolation, smoothing, least mean square method
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Computational Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Mathematical Physics (AREA)
- Theoretical Computer Science (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Complex Calculations (AREA)
Abstract
本发明提出了一种S参数的圆形插值方法,利用被测件S参数在Smith圆图上连续且呈现圆形的特点,将被测件S参数转化到Smith圆图的UV坐标系下;利用UV坐标下的数据进行圆形参数的求解;然后利用插值点频率对应的相位数据,作为圆形插值的依据;将插值点临近2点的数据,作为求解范围的限定;进而在Smith圆图上进行插值求解。本发明利用S参数在Smith圆图上的特性进行插值,能够得到更加精确的插值结果。
Description
技术领域
本发明涉及测试技术领域,特别涉及一种S参数圆形插值方法。
背景技术
在进行网络参数测试、仿真时,通常使用S参数进行计算。若使用的数据频率与测量到的数据频率不一致,则需要进行插值运算。虽然现有的插值方法种类很多,如线性插值、样条插值、拉格朗日插值等,但是这些都是对直角坐标系下的数据进行插值。如何更准确的得到插值数据,还需要根据被测件的S参数性质进行处理。
S参数有多种格式:实部+虚部、线性幅度+相位、对数幅度+相位、Smith圆图。在需要进行插值的情况下,最简单的是将S参数的实部、虚部数据进行分别插值。但是被测件的S参数在实部、虚部格式下没有呈现特定关系,无法表现被测件的特性。更精确的方法是将S参数转化到线性幅度、线性相位格式下,然后再分别对线性幅度、相位进行插值。这种情况下,可以将被测件的相位线性关系表现出来。但是由于被测件S参数在线性幅度下没有特定的线性关系,所以通过此种方式会与实际值有偏差。
由于对S参数的实部+虚部、线性幅度+相位插值方法,不能完整的表现被测件的特性,所以插值结果会有偏差。
真实的被测件可以分解成电感、电容、电阻模型的组合,它们在Smith圆图上是连续的,通常呈现圆形趋势,所以在Smith圆图上对被测件采用圆形插值,将会得到与实际值更接近的插值数据。
如何利用S参数在Smith圆图上的特性进行插值,是目前亟待解决的问题。
发明内容
本发明提出一种S参数圆形插值方法,解决了现有S参数插值方法插值结果有偏差的问题。
本发明的技术方案是这样实现的:
一种S参数的圆形插值方法,利用被测件S参数在Smith圆图上连续且呈现圆形的特点,将被测件S参数转化到Smith圆图的UV坐标系下;利用UV坐标下的数据进行圆形参数的求解;然后利用插值点频率对应的相位数据,作为圆形插值的依据;将插值点临近2点的数据,作为求解范围的限定;进而在Smith圆图上进行插值求解。
可选地,上述S参数的圆形插值方法,具体包括以下步骤:
第二步:所有插值点循环,若所有插值点都运算完成,则转到第八步;否则,转到第三步;
第四步:将插值点临近的3个测量值的实部虚部值,转化为Smith圆图UV坐标下的值;
第五步:利用下述公式(6),得到圆形的a、b、c参数;
对于直角坐标系而言,圆形的公式为:
(u-a)2+(v-b)2=c (6)
由第四步中的3点U、V值确定此圆的a、b、c参数;
第六步:在插值点临近2点范围内,求解满足公式(6)且相位为插值点相位的解;
第七步:插值点i增加1,转到第二步;
第八步:插值结束。
可选地,所述第四步中,S参数实部虚部值到Smith圆图UV坐标下的转化过程具体如下:
S参数与反射系数的计算公式如下:
其中,Γ为反射系数;Zn为S参数复数;
Zn=r+jx (3)
其中,Zn为S参数复数;r为S参数实部数据;x为S参数虚部数据;
Smith圆图是由公式(2)所描述的r和x在Γ复平面上的轨迹,将Γ分离为实部U和虚部V,得到:
其中:
变换后的实部U和虚部V构成了新的直角坐标系。
本发明的有益效果是:
本发明利用S参数在Smith圆图上的特性进行插值,能够得到更加精确的插值结果。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为Smith圆图下的UV坐标系示意图;
图2为本发明的Smith圆图下的圆形插值流程图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明提出了一种S参数的圆形插值方法,利用被测件S参数在Smith圆图上连续且呈现圆形的特点,将被测件S参数转化到Smith圆图的UV坐标系下;利用UV坐标下的数据进行圆形参数的求解;然后利用插值点频率对应的相位数据,作为圆形插值的依据;将插值点临近2点的数据,作为求解范围的限定,上述插值点临近2点是指插值点左侧最近一点和右侧最近一点;进而在Smith圆图上进行插值求解。利用本发明的S参数圆形插值方法能够显著提高对被测件S参数的插值效果。
下面结合说明书附图对本发明的S参数的圆形插值方法进行详细说明。
如图1所示,本发明的S参数的圆形插值方法包括以下步骤:
第二步:所有插值点循环,若所有插值点都运算完成,则转到第八步;否则,转到第三步。
第四步:将插值点临近的3个测量值的实部虚部值,转化为Smith圆图UV坐标下的值。上述插值点临近的3个测量值可以是选取插值点左侧最近的两个测量值和插值点右侧最近的一个测量值,也可以是选取插值点左侧最近的一个测量值和插值点右侧最近的两个测量值。
S参数实部+虚部到Smith圆图UV坐标下的转化过程具体如下:
S参数与反射系数的计算公式如下:
其中,Γ为反射系数;Zn为S参数复数;
Zn=r+jx (3)
其中,Zn为S参数复数;r为S参数实部数据;x为S参数虚部数据;
Smith圆图是由公式(2)所描述的r和x在Γ复平面上的轨迹,将Γ分离为实部U和虚部V,便得到
其中:
图2所示即为Smith圆图,变换后的实部U和虚部V构成了新的直角坐标系。
第五步:利用下述公式(6),得到圆形的a、b、c参数;
对于直角坐标系而言,圆形的公式为:
(u-a)2+(v-b)2=c (6)
由3点即可确定此圆的a、b、c。
第六步:在插值点临近2点范围内,求解满足公式(6)且相位为插值点相位的解。上述插值点临近2点是指插值点左侧最近一个测量值和插值点右侧最近一个测量值。
第七步:插值点i增加1,转到第二步。
第八步:插值结束。
至此,完成了对于S参数在Smith圆图上的圆形插值。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (2)
1.一种S参数的圆形插值方法,其特征在于,利用被测件S参数在Smith圆图上连续且呈现圆形的特点,将被测件S参数转化到Smith圆图的UV坐标系下;利用UV坐标下的数据进行圆形参数的求解;然后利用插值点频率对应的相位数据,作为圆形插值的依据;将插值点临近2点的数据,作为求解范围的限定;进而在Smith圆图上进行插值求解;
具体包括以下步骤:
第二步:所有插值点循环,若所有插值点都运算完成,则转到第八步;否则,转到第三步;
第三步:利用当前插值点i的频率进行插值,得到插值点的相位;
第四步:将插值点临近的3个测量值的实部虚部值,转化为Smith圆图UV坐标下的值,上述插值点临近的3个测量值是选取插值点左侧最近的两个测量值和插值点右侧最近的一个测量值,或者是选取插值点左侧最近的一个测量值和插值点右侧最近的两个测量值;
第五步:利用下述公式(6),得到圆形的a、b、c参数;
对于直角坐标系而言,圆形的公式为:
(u-a)2+(v-b)2=c (6)
由第四步中的3点U、V值确定此圆的a、b、c参数;
第六步:在插值点临近2点范围内,求解满足公式(6)且相位为插值点相位的解,上述插值点临近2点是指插值点左侧最近一个测量值和插值点右侧最近一个测量值;
第七步:插值点i增加1,转到第二步;
第八步:插值结束。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710477870.1A CN107203495B (zh) | 2017-06-09 | 2017-06-09 | 一种s参数圆形插值方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710477870.1A CN107203495B (zh) | 2017-06-09 | 2017-06-09 | 一种s参数圆形插值方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107203495A CN107203495A (zh) | 2017-09-26 |
CN107203495B true CN107203495B (zh) | 2020-12-22 |
Family
ID=59907941
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710477870.1A Active CN107203495B (zh) | 2017-06-09 | 2017-06-09 | 一种s参数圆形插值方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107203495B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108845182B (zh) * | 2018-05-29 | 2019-11-22 | 中国电子科技集团公司第四十一研究所 | 一种应用于夹具参数提取的相位计算方法 |
CN112084455B (zh) * | 2020-08-04 | 2024-03-01 | 中电科思仪科技股份有限公司 | 一种s参数的最小二乘插值方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106053945A (zh) * | 2016-07-05 | 2016-10-26 | 中国电子科技集团公司第四十研究所 | 一种短时s函数史密斯圆图分析装置及方法 |
CN106372306A (zh) * | 2016-08-31 | 2017-02-01 | 浪潮电子信息产业股份有限公司 | 确定负载端阻抗在史密斯圆图中的位置的方法及装置 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8319504B2 (en) * | 2009-05-29 | 2012-11-27 | Freescale Semiconductor, Inc. | Tuner characterization methods and apparatus |
US8904330B2 (en) * | 2013-03-15 | 2014-12-02 | Synopsys, Inc. | Phase determination for extraction from scattering parameters |
-
2017
- 2017-06-09 CN CN201710477870.1A patent/CN107203495B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106053945A (zh) * | 2016-07-05 | 2016-10-26 | 中国电子科技集团公司第四十研究所 | 一种短时s函数史密斯圆图分析装置及方法 |
CN106372306A (zh) * | 2016-08-31 | 2017-02-01 | 浪潮电子信息产业股份有限公司 | 确定负载端阻抗在史密斯圆图中的位置的方法及装置 |
Non-Patent Citations (4)
Title |
---|
列表点双圆弧样条曲线拟合计算法;钱一晨,陈虎娣;《机械研究与应用》;20121231(第119期);第19-26页 * |
基于负载牵引法的微波大功率自动调配器的控制软件的研究及其应用;安明;《中国优秀硕士学位论文全文数据库信息科技辑》;20110815;I135-184 * |
应用史密斯圆图提取慢波微带线特征阻抗方法;董峦;《电子测试》;20110531(第5期);第94-96页 * |
短波天线宽带匹配网络研究;张奇;《中国优秀硕士学位论文全文数据库信息科技辑》;20170315;I136-114 * |
Also Published As
Publication number | Publication date |
---|---|
CN107203495A (zh) | 2017-09-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110333389B (zh) | 基于插值dft的正弦信号频率估计方法 | |
US9627886B2 (en) | State estimation for power system using hybrid measurements | |
CN106777836B (zh) | 一种变压器直流偏磁仿真模拟方法及装置 | |
CN107085140B (zh) | 基于改进的SmartDFT算法的非平衡***频率估计方法 | |
CN108037361A (zh) | 一种基于滑动窗dft的高精度谐波参数估计方法 | |
CN103399213B (zh) | 一种基于基波频率修正的高压设备介质损耗因数计算方法 | |
CN107203495B (zh) | 一种s参数圆形插值方法 | |
CN110572779B (zh) | 一种基于ToF模型的基站的***误差估计方法 | |
CN104102836A (zh) | 一种电力***快速抗差状态估计方法 | |
CN114460527B (zh) | 关联度延拓Hilbert移相电子式互感器校验仪溯源方法和*** | |
CN103475369A (zh) | 基于信号源误差一次性校准识别的高精度adc测试方法 | |
CN114280526B (zh) | 一种电子式互感器校验仪数字微差溯源***及方法 | |
CN110188480B (zh) | 一种直流偏磁条件下铁磁材料的磁滞特性模拟分析***和方法 | |
CN103575979B (zh) | 一种数字化测量交流电频率的方法 | |
CN104407197A (zh) | 一种基于三角函数迭代的信号相量测量的方法 | |
CN108111448B (zh) | 预失真查找表的生成方法、装置与预失真校准设备 | |
Dutra et al. | High-precision frequency estimation of real sinusoids with reduced computational complexity using a model-based matched-spectrum approach | |
CN109638811B (zh) | 基于模型等值的配电网电压功率灵敏度鲁棒估计方法 | |
CN116680882A (zh) | 基于分潮重组与验潮约束的海洋深度基准构建方法及装置 | |
CN106405287A (zh) | 适合日常电能质量监测***应用的主导谐波源追踪方法 | |
CN106324315B (zh) | 一种空心罗氏线圈电流互感器的谐波特性补偿方法 | |
CN112084455B (zh) | 一种s参数的最小二乘插值方法 | |
CN107800434A (zh) | 一种基于参数提取的快速adc测试方法 | |
CN109327026B (zh) | 一种低压配电网区间状态估计方法 | |
CN103713686B (zh) | 一种小信号功率校正电路 |
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 |