CN107203495A - 一种s参数圆形插值方法 - Google Patents

一种s参数圆形插值方法 Download PDF

Info

Publication number
CN107203495A
CN107203495A CN201710477870.1A CN201710477870A CN107203495A CN 107203495 A CN107203495 A CN 107203495A CN 201710477870 A CN201710477870 A CN 201710477870A CN 107203495 A CN107203495 A CN 107203495A
Authority
CN
China
Prior art keywords
mrow
parameter
interpolation
msup
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.)
Granted
Application number
CN201710477870.1A
Other languages
English (en)
Other versions
CN107203495B (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.)
CETC 41 Institute
Original Assignee
CETC 41 Institute
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 CETC 41 Institute filed Critical CETC 41 Institute
Priority to CN201710477870.1A priority Critical patent/CN107203495B/zh
Publication of CN107203495A publication Critical patent/CN107203495A/zh
Application granted granted Critical
Publication of CN107203495B publication Critical patent/CN107203495B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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参数性质进行处理。
S参数有多种格式：实部+虚部、线性幅度+相位、对数幅度+相位、Smith圆图。在需要进行插值的情况下，最简单的是将S参数的实部、虚部数据进行分别插值。但是被测件的S参数在实部、虚部格式下没有呈现特定关系，无法表现被测件的特性。更精确的方法是将S参数转化到线性幅度、线性相位格式下，然后再分别对线性幅度、相位进行插值。这种情况下，可以将被测件的相位线性关系表现出来。但是由于被测件S参数在线性幅度下没有特定的线性关系，所以通过此种方式会与实际值有偏差。
由于对S参数的实部+虚部、线性幅度+相位插值方法，不能完整的表现被测件的特性，所以插值结果会有偏差。
真实的被测件可以分解成电感、电容、电阻模型的组合，它们在Smith圆图上是连续的，通常呈现圆形趋势，所以在Smith圆图上对被测件采用圆形插值，将会得到与实际值更接近的插值数据。
如何利用S参数在Smith圆图上的特性进行插值，是目前亟待解决的问题。
发明内容
本发明提出一种S参数圆形插值方法，解决了现有S参数插值方法插值结果有偏差的问题。
本发明的技术方案是这样实现的：
一种S参数的圆形插值方法，利用被测件S参数在Smith圆图上连续且呈现圆形的特点，将被测件S参数转化到Smith圆图的UV坐标系下；利用UV坐标下的数据进行圆形参数的求解；然后利用插值点频率对应的相位数据，作为圆形插值的依据；将插值点临近2点的数据，作为求解范围的限定；进而在Smith圆图上进行插值求解。
可选地，上述S参数的圆形插值方法，具体包括以下步骤：
第一步：将所有测量值的实部虚部，通过公式(1)计算得到相位
其中，为S参数相位值；r为S参数实部数据；x为S参数虚部数据；
第二步：所有插值点循环，若所有插值点都运算完成，则转到第八步；否则，转到第三步；
第三步：利用当前插值点i的频率进行插值，得到插值点的相位
第四步：将插值点临近的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参数的圆形插值方法包括以下步骤：
第一步：将所有测量值的实部虚部，通过公式(1)计算得到相位
其中，为S参数相位值；r为S参数实部数据；x为S参数虚部数据。
第二步：所有插值点循环，若所有插值点都运算完成，则转到第八步；否则，转到第三步。
第三步：利用当前插值点i的频率进行插值，得到插值点的相位
第四步：将插值点临近的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 (3)

1.一种S参数的圆形插值方法，其特征在于，利用被测件S参数在Smith圆图上连续且呈现圆形的特点，将被测件S参数转化到Smith圆图的UV坐标系下；利用UV坐标下的数据进行圆形参数的求解；然后利用插值点频率对应的相位数据，作为圆形插值的依据；将插值点临近2点的数据，作为求解范围的限定；进而在Smith圆图上进行插值求解。
2.如权利要求1所述的一种S参数的圆形插值方法，其特征在于，具体包括以下步骤：
第一步：将所有测量值的实部虚部，通过公式(1)计算得到相位
其中，为S参数相位值；r为S参数实部数据；x为S参数虚部数据；
第二步：所有插值点循环，若所有插值点都运算完成，则转到第八步；否则，转到第三步；
第三步：利用当前插值点i的频率进行插值，得到插值点的相位
第四步：将插值点临近的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，转到第二步；
第八步：插值结束。
3.如权利要求2所述的一种S参数的圆形插值方法，其特征在于，所述第四步中，S参数实部虚部值到Smith圆图UV坐标下的转化过程具体如下：
S参数与反射系数的计算公式如下：
<mrow> <mi>&amp;Gamma;</mi> <mo>=</mo> <mfrac> <mrow> <msub> <mi>Z</mi> <mi>n</mi> </msub> <mo>-</mo> <mn>1</mn> </mrow> <mrow> <msub> <mi>Z</mi> <mi>n</mi> </msub> <mo>+</mo> <mn>1</mn> </mrow> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </mrow>
其中，Γ为反射系数；Zn为S参数复数；
Zn＝r+jx (3)
其中，Zn为S参数复数；r为S参数实部数据；x为S参数虚部数据；
Smith圆图是由公式(2)所描述的r和x在Γ复平面上的轨迹，将Γ分离为实部U和虚部V，得到：
<mrow> <mi>&amp;Gamma;</mi> <mo>=</mo> <mfrac> <mrow> <mi>r</mi> <mo>+</mo> <mi>j</mi> <mi>x</mi> <mo>-</mo> <mn>1</mn> </mrow> <mrow> <mi>r</mi> <mo>+</mo> <mi>j</mi> <mi>x</mi> <mo>+</mo> <mn>1</mn> </mrow> </mfrac> <mo>=</mo> <mi>U</mi> <mo>+</mo> <mi>j</mi> <mi>V</mi> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>4</mn> <mo>)</mo> </mrow> </mrow>
其中：
<mrow> <mtable> <mtr> <mtd> <mrow> <mi>U</mi> <mo>=</mo> <mfrac> <mrow> <msup> <mi>r</mi> <mn>2</mn> </msup> <mo>-</mo> <mn>1</mn> <mo>+</mo> <msup> <mi>x</mi> <mn>2</mn> </msup> </mrow> <mrow> <msup> <mrow> <mo>(</mo> <mi>r</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msup> <mi>x</mi> <mn>2</mn> </msup> </mrow> </mfrac> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>V</mi> <mo>=</mo> <mfrac> <mrow> <mn>2</mn> <mi>x</mi> </mrow> <mrow> <msup> <mrow> <mo>(</mo> <mi>r</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msup> <mi>x</mi> <mn>2</mn> </msup> </mrow> </mfrac> </mrow> </mtd> </mtr> </mtable> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>5</mn> <mo>)</mo> </mrow> </mrow> 1
变换后的实部U和虚部V构成了新的直角坐标系。
CN201710477870.1A 2017-06-09 2017-06-09 一种s参数圆形插值方法 Active CN107203495B (zh)

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 true CN107203495A (zh) 2017-09-26
CN107203495B 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)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108845182A (zh) * 2018-05-29 2018-11-20 中国电子科技集团公司第四十研究所 一种应用于夹具参数提取的相位计算方法
CN112084455A (zh) * 2020-08-04 2020-12-15 中电科仪器仪表有限公司 一种s参数的最小二乘插值方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100301875A1 (en) * 2009-05-29 2010-12-02 Freescale Semiconductor, Inc. Tuner characterization methods and apparatus
US20140282335A1 (en) * 2013-03-15 2014-09-18 Synopsys, Inc. Phase determination for extraction from scattering parameters
CN106053945A (zh) * 2016-07-05 2016-10-26 中国电子科技集团公司第四十研究所 一种短时s函数史密斯圆图分析装置及方法
CN106372306A (zh) * 2016-08-31 2017-02-01 浪潮电子信息产业股份有限公司 确定负载端阻抗在史密斯圆图中的位置的方法及装置

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100301875A1 (en) * 2009-05-29 2010-12-02 Freescale Semiconductor, Inc. Tuner characterization methods and apparatus
US20140282335A1 (en) * 2013-03-15 2014-09-18 Synopsys, Inc. Phase determination for extraction from scattering parameters
CN106053945A (zh) * 2016-07-05 2016-10-26 中国电子科技集团公司第四十研究所 一种短时s函数史密斯圆图分析装置及方法
CN106372306A (zh) * 2016-08-31 2017-02-01 浪潮电子信息产业股份有限公司 确定负载端阻抗在史密斯圆图中的位置的方法及装置

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
安明: "基于负载牵引法的微波大功率自动调配器的控制软件的研究及其应用", 《中国优秀硕士学位论文全文数据库信息科技辑》 *
张奇: "短波天线宽带匹配网络研究", 《中国优秀硕士学位论文全文数据库信息科技辑》 *
波塞: "《天线系统及其在电子战系统中的应用》", 31 January 2014 *
董峦: "应用史密斯圆图提取慢波微带线特征阻抗方法", 《电子测试》 *
钱一晨，陈虎娣: "列表点双圆弧样条曲线拟合计算法", 《机械研究与应用》 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108845182A (zh) * 2018-05-29 2018-11-20 中国电子科技集团公司第四十研究所 一种应用于夹具参数提取的相位计算方法
CN108845182B (zh) * 2018-05-29 2019-11-22 中国电子科技集团公司第四十一研究所 一种应用于夹具参数提取的相位计算方法
CN112084455A (zh) * 2020-08-04 2020-12-15 中电科仪器仪表有限公司 一种s参数的最小二乘插值方法
CN112084455B (zh) * 2020-08-04 2024-03-01 中电科思仪科技股份有限公司 一种s参数的最小二乘插值方法

Also Published As

Publication number Publication date
CN107203495B (zh) 2020-12-22

Similar Documents

Publication Publication Date Title
CN103986357B (zh) 抑制高压直流传输用模块化多电平换流器中的环流的方法
CN102508116B (zh) 基于电流平衡原理的智能变电站同杆双回线路选相方法
CN102865862B (zh) 全站仪免置平自由设站的测量方法
CN105223466B (zh) 一种利用模极大值比的特高压直流输电线路单端测距方法
CN107276660A (zh) 非正交多址空地协同通信系统中资源分配方法及装置
CN107203495A (zh) 一种s参数圆形插值方法
CN202281802U (zh) 一种三相三线电能计量表
CN105184027B (zh) 一种基于交互式多模型算法的电力负荷建模方法
CN103490772B (zh) 一种基于无功补偿的滑动加权单相软锁相方法
CN104237740A (zh) 架空线-电缆混合线路分布式故障测距方法
CN106451522B (zh) 一种高压直流回路谐波阻抗解析等效模型及频率分析方法
CN107681682A (zh) 一种基于ward等值的交直流系统等值方法
CN106505557A (zh) 一种遥测错误辨识方法及装置
CN105159063A (zh) 一种分数相位超前补偿重复控制器及控制方法
CN105652147A (zh) 一种单/双回输电线路全线高精度故障测距方法
CN110715657A (zh) 一种航空施药区域全覆盖路径规划方法
CN109507630A (zh) 接线判断方法及系统
CN109239533A (zh) 一种基于人工神经网络的特高压直流输电线路的故障定位方法
CN105929305A (zh) 一种非全程混压双回线故障区段识别与精确测距方法
CN102930140A (zh) 一种同塔多回复杂输电线路传输方程的求解方法
CN103401447B (zh) 一种级联型pwm整流器的调制方法
CN105301440A (zh) 一种线缆混合输电线路双端行波测距的方法和系统
CN103310124B (zh) 一种最优潮流模型的建模方法及装置
CN203299274U (zh) 一种高压开关柜相位测量电路
CN104266646B (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