CN115455745B - 自适应频点采样的扫频方法、系统及相关设备 - Google Patents

自适应频点采样的扫频方法、系统及相关设备 Download PDF

Info

Publication number
CN115455745B
CN115455745B CN202211401494.5A CN202211401494A CN115455745B CN 115455745 B CN115455745 B CN 115455745B CN 202211401494 A CN202211401494 A CN 202211401494A CN 115455745 B CN115455745 B CN 115455745B
Authority
CN
China
Prior art keywords
frequency
sampling
point
frequency point
points
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
CN202211401494.5A
Other languages
English (en)
Other versions
CN115455745A (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.)
Shenzhen Volans Technology Co Ltd
Original Assignee
Shenzhen Volans Technology Co Ltd
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 Shenzhen Volans Technology Co Ltd filed Critical Shenzhen Volans Technology Co Ltd
Priority to CN202211401494.5A priority Critical patent/CN115455745B/zh
Publication of CN115455745A publication Critical patent/CN115455745A/zh
Application granted granted Critical
Publication of CN115455745B publication Critical patent/CN115455745B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • 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
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02D—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
    • Y02D30/00—Reducing energy consumption in communication networks
    • Y02D30/70—Reducing energy consumption in communication networks in wireless communication networks

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Data Mining & Analysis (AREA)
  • General Engineering & Computer Science (AREA)
  • Software Systems (AREA)
  • Computational Mathematics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Databases & Information Systems (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Algebra (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • Feedback Control In General (AREA)

Abstract

本发明属于无线通信技术领域，尤其涉及一种自适应频点采样的扫频方法、系统及相关设备，所述扫频方法运用于射频器件的仿真。本发明采用了一种不依赖于具体数值方法自适应频点采样方法，能够对频率响应快速变化处进行精确仿真，而在其他不感兴趣的区域进行粗略仿真，从而提高射频器件设计过程的精度与效率；同时，针对频率响应中Y参数与S参数的特点使用不同的评价指标，从而减少扫频所需要的采样点，提高扫频速度。

Description

自适应频点采样的扫频方法、系统及相关设备
技术领域
本发明属于无线通信技术领域，尤其涉及一种自适应频点采样的扫频方法、系统及相关设备。
背景技术
随着无线通讯技术的发展，各类电子消费产品对射频器件的需求越来越大，相对应的，在射频器件设计过程中，设计人员对仿真软件扫频效率与准确度的要求也越来越高。通常情况下，密集的扫频点可以提高扫频曲线的精确度，但同样也会大大增加计算所需的时间。因此需要加速算法在保证计算精度的前提下尽可能的减小计算所需时间。
主流的仿真软件，例如COMSOL，通常会采用渐进波形估计方法（AsymptoticWaveform Evaluation，简称AWE）与自适应兰索斯-帕德扫频方法（ALPS）等方法来提高扫频的计算效率。在有限单元法中，以上方法的基本原理是通过计算几个典型频点下的频率响应，然后通过泰勒级数、Pade有理逼近等数学方法对整体矩阵进行近似展开，从而对相邻频点的频率响应以一种很小的计算开销来估计。问题在于，AWE等方法所基于的泰勒级数展开往往只能在采样频点周围具有较高的精度，而在距离采样频点较远处以及频率响应快速变化处存在频率相应估计精度不高的问题，而频率相应快速变化的区域往往是采样过程中比较关心的区域。
在实际仿真中，往往只对频带部分区域，如通带内的频率响应的精度具有较高要求，而对其他区域的频率响应的精度则没有太高的要求。传统的估计方法则会不加区分的将全部频带划分成若干个较窄的频段，然后对每个频带分别进行扫频分析，因此仍然存在效率较低的问题，因此，传统的估计方法一方面存在频率响应快速变化处可能精度不高的问题，另一方面也存在对于采样中不关心的区域消耗了大量计算资源的问题。并且，传统的估计方法往往依赖于具体的数值模拟算法，如渐进波形估计方法需要将有限元方法的整体矩阵进行分解，从而得到电场值随频率变化的解析表达式，这样的方法对于实验测量、机器学习等不存在显式控制方程的环境不适用。
发明内容
本发明实施例提供一种自适应频点采样的扫频方法、系统及相关设备，旨在解决传统的扫频方法依赖具体数值、不适用于实验测量、机器学习等场景的问题。
第一方面，本发明实施例提供一种自适应频点采样的扫频方法，运用于射频器件的仿真，所述扫频方法包括以下步骤：
S1、确定频点采样数n、频点最终采样数和采样范围，并使用预设仿真方法在所述采样范围中进行采样，并计算n处采样点的频率响应，其中，采样点处的所述频率响应为；
S2、计算所述频率响应的一阶导数和二阶导数；
S3、根据所述一阶导数、所述二阶导数对n处所述采样点进行频率响应曲线估计，得到条估计曲线，其中，采样点对应的估计曲线为，适用范围为；
S4、依次获取相邻的所述估计曲线与之间的最大误差频点，直至得到个所述最大误差频点，每一所述最大误差频点具有误差值；
S5、根据预设评价方法对每一所述最大误差频点的所述误差值进行计算，得到每一所述最大误差频点的评价指标；
S6、根据所述评价指标进行排序，并选择前m个所述最大误差频点作为已采样点，其中，m小于n，m、n均为大于0的整数，且：
若所述已采样点的累计数量未达到所述频点最终采样数，则返回步骤S2进行迭代执行；
若所述已采样点的累计数量达到所述频点最终采样数，则停止迭代；
S7、输出所有所述已采样点作为最终采样点。
更进一步地，步骤S1中，所述预设仿真方法在所述采样范围中进行采样的方式为均匀采样。
更进一步地，步骤S2中，计算所述频率响应的一阶导数和二阶导数所使用的方法为中心差分法。
更进一步地，步骤S3中，所述估计曲线满足以下关系式（1）：
（1）。
更进一步地，步骤S4中，所述最大误差频点xp满足以下关系式（2）：
（2）；
所述误差值满足以下关系式（3）：
（3）。
更进一步地，步骤S5中，当所述频率响应的类型为Y参数时，所述预设评价方法为：
S51、在所述采样点中，选择Y参数取值最大的频点最大值和频点最小值，并计算出对应的最大值估计曲线和最小值估计曲线；
S52、根据所述最大值估计曲线和所述最小值估计曲线得到谐振点和反谐振点，并将所述谐振点和所述反谐振点加入迭代使用的所述采样点中；
S53、定义所述最大误差频点的所述评价指标为，所述评价指标满足以下关系式（4）：
（4）。
更进一步地，步骤S5中，当所述频率响应的类型为S参数时，所述预设评价方法为：
S51、定义所述频率响应的S参数包括第一S参数和第二S参数；
S52、分别计算所述第一S参数和所述第二S参数的估计曲线的评价指标和，其中，所述第一S参数的评价指标满足以下关系式（5）：
（5）；
所述第二S参数的评价指标F21(xp)满足以下关系式（6）：
（6）；
S53、定义所述最大误差频点的所述评价指标为，所述评价指标满足以下关系式（7）：
（7）。
第二方面，本发明实施例还提供一种自适应频点采样的扫频系统，运用于射频器件的仿真，包括：
初始化模块，用于确定频点采样数n、频点最终采样数和采样范围，并使用预设仿真方法在所述采样范围中进行采样，并计算n处采样点的频率响应，其中，采样点处的所述频率响应为；
求导模块，用于计算所述频率响应的一阶导数和二阶导数；
频率响应估计模块，用于根据所述一阶导数、所述二阶导数对n处所述采样点进行频率响应曲线估计，得到条估计曲线，其中，采样点对应的估计曲线为，适用范围为；
最大误差计算模块，用于依次获取相邻的所述估计曲线与之间的最大误差频点，直至得到个所述最大误差频点，每一所述最大误差频点具有误差值；
评价模块，用于根据预设评价方法对每一所述最大误差频点的所述误差值进行计算，得到每一所述最大误差频点的评价指标；
迭代模块，用于根据所述评价指标进行排序，并选择前m个所述最大误差频点作为已采样点，其中，m小于n，m、n均为大于0的整数，且：
若所述已采样点的累计数量未达到所述频点最终采样数，则返回所述频率响应估计模块进行迭代执行；
若所述已采样点的累计数量达到所述频点最终采样数，则停止迭代；
输出模块，用于输出所有所述已采样点作为最终采样点。
第三方面，本发明实施例还提供一种计算机设备，包括：存储器、处理器及存储在所述存储器上并可在所述处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现如上述实施例中任意一项所述的自适应频点采样的方法中的步骤。
第四方面，本发明实施例还提供一种计算机可读存储介质，所述计算机可读存储介质上存储有计算机程序，所述计算机程序被处理器执行时实现如上述实施例中任意一项所述的自适应频点采样的方法中的步骤。
本发明所达到的有益效果，由于采用了一种不依赖于具体数值方法自适应频点采样方法，能够对频率响应快速变化处进行精确仿真，而在其他不感兴趣的区域进行粗略仿真，从而提高射频器件设计过程的精度与效率；同时，针对频率响应中Y参数与S参数的特点使用不同的评价指标，从而减少扫频所需要的采样点，提高扫频速度。
附图说明
图1是本发明实施例提供的自适应频点采样的扫频方法的步骤流程框图；
图2是本发明实施例提供的自适应频点采样的扫频方法在S参数中采样的示意图；
图3是本发明实施例提供的自适应频点采样的扫频方法在Y参数中采样的示意图；
图4是本发明实施例提供的自适应频点采样的扫频系统的结构示意图；
图5是本发明实施例提供的计算机设备的结构示意图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。
请参照图1，图1是本发明实施例提供的自适应频点采样的扫频方法的步骤流程框图，所述扫频方法运用于射频器件的仿真，所述扫频方法包括以下步骤：
S1、确定频点采样数n、频点最终采样数和采样范围，并使用预设仿真方法在所述采样范围中进行采样，并计算n处采样点的频率响应，其中，采样点处的所述频率响应为。
更进一步地，步骤S1中，所述预设仿真方法在所述采样范围中进行采样的方式为均匀采样。
具体的，所述预设仿真方法可以是任意一种数值仿真计算方法，例如Matlab等。
S2、计算所述频率响应的一阶导数和二阶导数。
更进一步地，步骤S2中，计算所述频率响应的一阶导数和二阶导数所使用的方法为中心差分法。
S3、根据所述一阶导数、所述二阶导数对n处所述采样点进行频率响应曲线估计，得到条估计曲线，其中，采样点对应的估计曲线为，适用范围为。
具体的，根据所述一阶导数、所述二阶导数对n处所述采样点进行频率响应曲线估计时，按照每两个相邻的所述采样点进行，并且可以进一步通过泰勒公式进行估计，从而得到所述估计曲线。
更进一步地，步骤S3中，所述估计曲线满足以下关系式（1）：
（1）。
S4、依次获取相邻的所述估计曲线与之间的最大误差频点，直至得到个所述最大误差频点，每一所述最大误差频点具有误差值。
具体的，对于任意范围中的一点，其都可以通过其两侧采样点的估计曲线与进行预测，即与。
同时，为了便于后续评价指标及迭代的计算，本发明实施例通过数据验证，将所述最大误差频点取值为其计算使用的所述估计曲线对应的所述采样点的中点。
更进一步地，步骤S4中，所述最大误差频点xp满足以下关系式（2）：
   （2）；
所述误差值满足以下关系式（3）：
（3）。
S5、根据预设评价方法对每一所述最大误差频点的所述误差值进行计算，得到每一所述最大误差频点的评价指标。
具体的，所述频率响应的参数类型可以为多种，包括S参数和Y参数，不同参数虽然本质上都等效为二端口网络之后的二端口网络参数，但得到数值的过程不同，在本发明实施例中，根据数值类型的不同，所述评价指标的计算方法也不同。
更进一步地，步骤S5中，当所述频率响应的类型为Y参数时，所述预设评价方法为：
S51、在所述采样点中，选择Y参数取值最大的频点最大值和频点最小值，并计算出对应的最大值估计曲线和最小值估计曲线；
S52、根据所述最大值估计曲线和所述最小值估计曲线得到谐振点和反谐振点，并将所述谐振点和所述反谐振点加入迭代使用的所述采样点中；
S53、定义所述最大误差频点的所述评价指标为，所述评价指标满足以下关系式（4）：
 （4）。
所述谐振点和所述反谐振点加入迭代使用的所述采样点后，采样点所体现的采集方式会逐渐向所述谐振点和所述反谐振点的位置靠近，从而实现目标兴趣频带的扫频。
对应的，本发明还提供另一种实施例，步骤S5中，当所述频率响应的类型为S参数时，所述预设评价方法为：
S51、定义所述频率响应的S参数包括第一S参数和第二S参数；
S52、分别计算所述第一S参数和所述第二S参数的估计曲线的评价指标和，其中，所述第一S参数的评价指标满足以下关系式（5）：
（5）；
所述第二S参数的评价指标F21(xp)满足以下关系式（6）：
（6）；
S53、定义所述最大误差频点的所述评价指标为，所述评价指标满足以下关系式（7）：
 （7）。
S6、根据所述评价指标进行排序，并选择前m个所述最大误差频点作为已采样点，其中，m小于n，m、n均为大于0的整数，且：
若所述已采样点的累计数量未达到所述频点最终采样数，则返回步骤S2进行迭代执行；
若所述已采样点的累计数量达到所述频点最终采样数，则停止迭代。
具体的，本发明实施例设定所述频点最终采样数，并将步骤S2至步骤S6作为一次完整的迭代过程，其中，在每次迭代后，在原本划分的所述采样点的基础上再进行频率估计，会得到比迭代前一次使用的所述采样点的数量更多的用于计算的所述采样点，通过多次迭代，并将每次迭代中满足所述评价指标的所述最大误差频点作为选用的采样点，从而在较少频点情况下对目标频带进行高效地采样。
S7、输出所有所述已采样点作为最终采样点。
示例性的，请参照图2和图3，其分别是本发明实施例提供的自适应频点采样的扫频方法在S参数和Y参数中采样的示意图，其中，对于S参数，在通带内采样点的连线与真实曲线具有较高的一致性；对于Y参数，采样点的连线与真实曲线整体具有较高的一致性；图2、图3中下方的散点为每次采样时选择的点，其采样批次呈层级变化，可以看出，第一批次采样点与第二批次采样点均体现为均匀采样，但第三批次及之后的采样点则集中在了较为集中的区域，在图2及图3中，第三批次及之后的采样点均向S参数的带内区域以及Y参数的谐振点、反谐振点处靠近。
本发明所达到的有益效果，由于采用了一种不依赖于具体数值方法自适应频点采样方法，能够对频率响应快速变化处进行精确仿真，而在其他不感兴趣的区域进行粗略仿真，从而提高射频器件设计过程的精度与效率；同时，针对频率响应中Y参数与S参数的特点使用不同的评价指标，从而减少扫频所需要的采样点，提高扫频速度。
本发明实施例还提供一种自适应频点采样的扫频系统，请参照图4，图4是本发明实施例提供的自适应频点采样的扫频系统的结构示意图，所述自适应频点采样的扫频系统200包括：
初始化模块201，用于确定频点采样数n、频点最终采样数和采样范围，并使用预设仿真方法在所述采样范围中进行采样，并计算n处采样点的频率响应，其中，采样点处的所述频率响应为；
求导模块202，用于计算所述频率响应的一阶导数和二阶导数；
频率响应估计模块203，用于根据所述一阶导数、所述二阶导数对n处所述采样点进行频率响应曲线估计，得到条估计曲线，其中，采样点对应的估计曲线为，适用范围为；
最大误差计算模块204，用于依次获取相邻的所述估计曲线与之间的最大误差频点，直至得到个所述最大误差频点，每一所述最大误差频点具有误差值；
评价模块205，用于根据预设评价方法对每一所述最大误差频点的所述误差值进行计算，得到每一所述最大误差频点的评价指标；
迭代模块206，用于根据所述评价指标进行排序，并选择前m个所述最大误差频点作为已采样点，其中，m小于n，m、n均为大于0的整数，且：
若所述已采样点的累计数量未达到所述频点最终采样数，则返回所述频率响应估计模块进行迭代执行；
若所述已采样点的累计数量达到所述频点最终采样数，则停止迭代；
输出模块207，用于输出所有所述已采样点作为最终采样点。
所述自适应频点采样的扫频系统200能够实现如上述实施例中的自适应频点采样的扫频方法中的步骤，且能实现同样的技术效果，参上述实施例中的描述，此处不再赘述。
本发明实施例还提供一种计算机设备，请参照图5，图5是本发明实施例提供的计算机设备的结构示意图，所述计算机设备300包括：存储器302、处理器301及存储在所述存储器302上并可在所述处理器301上运行的计算机程序。
所述处理器301调用所述存储器302存储的计算机程序，执行本发明实施例提供的自适应频点采样的扫频方法中的步骤，请结合图1，具体包括：
S1、确定频点采样数n、频点最终采样数和采样范围，并使用预设仿真方法在所述采样范围中进行采样，并计算n处采样点的频率响应，其中，采样点处的所述频率响应为；
S2、计算所述频率响应的一阶导数和二阶导数；
S3、根据所述一阶导数、所述二阶导数对n处所述采样点进行频率响应曲线估计，得到条估计曲线，其中，采样点对应的估计曲线为，适用范围为；
S4、依次获取相邻的所述估计曲线与之间的最大误差频点，直至得到个所述最大误差频点，每一所述最大误差频点具有误差值；
S5、根据预设评价方法对每一所述最大误差频点的所述误差值进行计算，得到每一所述最大误差频点的评价指标；
S6、根据所述评价指标进行排序，并选择前m个所述最大误差频点作为已采样点，其中，m小于n，m、n均为大于0的整数，且：
若所述已采样点的累计数量未达到所述频点最终采样数，则返回步骤S2进行迭代执行；
若所述已采样点的累计数量达到所述频点最终采样数，则停止迭代；
S7、输出所有所述已采样点作为最终采样点。
更进一步地，步骤S1中，所述预设仿真方法在所述采样范围中进行采样的方式为均匀采样。
更进一步地，步骤S2中，计算所述频率响应的一阶导数和二阶导数所使用的方法为中心差分法。
更进一步地，步骤S3中，所述估计曲线满足以下关系式（1）：
（1）。
更进一步地，步骤S4中，所述最大误差频点满足以下关系式（2）：
  （2）；
所述误差值满足以下关系式（3）：
  （3）。
更进一步地，步骤S5中，当所述频率响应的类型为Y参数时，所述预设评价方法为：
S51、在所述采样点中，选择Y参数取值最大的频点最大值和频点最小值，并计算出对应的最大值估计曲线和最小值估计曲线；
S52、根据所述最大值估计曲线和所述最小值估计曲线得到谐振点和反谐振点，并将所述谐振点和所述反谐振点加入迭代使用的所述采样点中；
S53、定义所述最大误差频点的所述评价指标为，所述评价指标满足以下关系式（4）：
  （4）。
更进一步地，步骤S5中，当所述频率响应的类型为S参数时，所述预设评价方法为：
S51、定义所述频率响应的S参数包括第一S参数和第二S参数；
S52、分别计算所述第一S参数和所述第二S参数的估计曲线的评价指标和，其中，所述第一S参数的评价指标满足以下关系式（5）：
（5）；
所述第二S参数的评价指标满足以下关系式（6）：
（6）；
S53、定义所述最大误差频点的所述评价指标为，所述评价指标满足以下关系式（7）：
 （7）。
本发明实施例提供的计算机设备300能够实现如上述实施例中的自适应频点采样的扫频方法中的步骤，且能实现同样的技术效果，参上述实施例中的描述，此处不再赘述。
本发明实施例还提供一种计算机可读存储介质，所述计算机可读存储介质上存储有计算机程序，该计算机程序被处理器执行时实现本发明实施例提供的自适应频点采样的扫频方法中的各个过程及步骤，且能实现相同的技术效果，为避免重复，这里不再赘述。
本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程，是可以通过计算机程序来指令相关的硬件来完成，所述的程序可存储于一计算机可读取存储介质中，该程序在执行时，可包括如上述各方法的实施例的流程。其中，所述的存储介质可为磁碟、光盘、只读存储记忆体（Read-Only Memory，ROM）或随机存取存储器（Random AccessMemory，简称RAM）等。
需要说明的是，在本文中，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、物品或者装置不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、物品或者装置所固有的要素。在没有更多限制的情况下，由语句“包括一个……”限定的要素，并不排除在包括该要素的过程、方法、物品或者装置中还存在另外的相同要素。
通过以上的实施方式的描述，本领域的技术人员可以清楚地了解到上述实施例方法可借助软件加必需的通用硬件平台的方式来实现，当然也可以通过硬件，但很多情况下前者是更佳的实施方式。基于这样的理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质（如ROM/RAM、磁碟、光盘）中，包括若干指令用以使得一台终端（可以是手机，计算机，服务器，空调器，或者网络设备等）执行本发明各个实施例所述的方法。
上面结合附图对本发明的实施例进行了描述，所揭露的仅为本发明较佳实施例而已，但是本发明并不局限于上述的具体实施方式，上述的具体实施方式仅仅是示意性的，而不是限制性的，本领域的普通技术人员在本发明的启示下，在不脱离本发明宗旨和权利要求所保护的范围情况下，还可做出很多形式用等同变化，均属于本发明的保护之内。

Claims (10)

1.一种自适应频点采样的扫频方法，运用于射频器件的仿真，其特征在于，所述扫频方法包括以下步骤：
S1、确定频点采样数n、频点最终采样数n`和采样范围[fmin，fmax]，并使用预设仿真方法在所述采样范围中进行采样，并计算n处采样点的频率响应，其中，采样点xn处的所述频率响应为f(x)；
S2、计算所述频率响应f(x)的一阶导数和二阶导数；
S3、根据所述一阶导数、所述二阶导数对n处所述采样点进行频率响应曲线估计，得到n-1条估计曲线，其中，采样点xn对应的估计曲线为gxn(x)，适用范围为[xn-1, xn+1]；
S4、依次获取相邻的所述估计曲线gxn(x)与gxn+1(x)之间的最大误差频点xp，直至得到n-1个所述最大误差频点，每一所述最大误差频点具有误差值；
S5、根据预设评价方法对每一所述最大误差频点的所述误差值进行计算，得到每一所述最大误差频点的评价指标；
S6、根据所述评价指标进行排序，并选择前m个所述最大误差频点作为已采样点，其中，m小于n，m、n均为大于0的整数，且：
若所述已采样点的累计数量未达到所述频点最终采样数n`，则返回步骤S2进行迭代执行；
若所述已采样点的累计数量达到所述频点最终采样数n`，则停止迭代；
S7、输出所有所述已采样点作为最终采样点。
2.如权利要求1所述的扫频方法，其特征在于，步骤S1中，所述预设仿真方法在所述采样范围中进行采样的方式为均匀采样。
3.如权利要求1所述的扫频方法，其特征在于，步骤S2中，计算所述频率响应f(x)的一阶导数和二阶导数所使用的方法为中心差分法。
4.如权利要求1所述的扫频方法，其特征在于，步骤S3中，所述估计曲线gxn(x)满足以下关系式（1）：
gxn(x)=f(xn)+f '(xn)*(x-xn)+1/2*f ''(xn)*(x-xn)2（1）。
5.如权利要求4所述的扫频方法，其特征在于，步骤S4中，所述最大误差频点xp满足以下关系式（2）：
xp=(xn+xn+1)/2   （2）；
定义所述误差值为Axp，所述误差值Axp满足以下关系式（3）：
Axp=Abs[gxn(xp)-gxn+1(xp)]   （3）。
6.如权利要求5所述的扫频方法，其特征在于，步骤S5中，当所述频率响应的类型为Y参数时，所述预设评价方法为：
S51、在所述采样点中，选择Y参数取值最大的频点最大值xmax和频点最小值xmin，并计算出对应的最大值估计曲线gxmax(x)和最小值估计曲线gxmin(x)；
S52、根据所述最大值估计曲线gxmax(x)和所述最小值估计曲线gxmin(x)得到谐振点xr和反谐振点xa，并将所述谐振点xr和所述反谐振点xa加入迭代使用的所述采样点中；
S53、定义所述最大误差频点的所述评价指标为F(xp)，所述评价指标F(xp)满足以下关系式（4）：
F(xp)=Abs[gxn(xp)-gxn+1(xp)]   （4）。
7.如权利要求5所述的扫频方法，其特征在于，步骤S5中，当所述频率响应的类型为S参数时，所述预设评价方法为：
S51、定义所述频率响应的S参数包括第一S参数和第二S参数；
S52、分别计算所述第一S参数和所述第二S参数的估计曲线的评价指标F11(xp)和F21(xp)，其中，所述第一S参数的评价指标F11(xp)满足以下关系式（5）：
F11(xp)=Abs[gxn(xp)-gxn+1(xp)]/Abs[(gxn(xp)+gxn+1(xp))]2（5）；
所述第二S参数的评价指标F21(xp)满足以下关系式（6）：
F21(xp)=Abs[gxn(xp)-gxn+1(xp)]*Abs[(gxn(xp)+gxn+1(xp))]2（6）；
S53、定义所述最大误差频点的所述评价指标为F(xp)，所述评价指标F(xp)满足以下关系式（7）：
F(xp)=F11(xp)+F21(xp)   （7）。
8.一种自适应频点采样的扫频系统，运用于射频器件的仿真，其特征在于，包括：
初始化模块，用于确定频点采样数n、频点最终采样数n`和采样范围[fmin，fmax]，并使用预设仿真方法在所述采样范围中进行采样，并计算n处采样点的频率响应，其中，采样点xn处的所述频率响应为f(x)；
求导模块，用于计算所述频率响应f(x)的一阶导数和二阶导数；
频率响应估计模块，用于根据所述一阶导数、所述二阶导数对n处所述采样点进行频率响应曲线估计，得到n-1条估计曲线，其中，采样点xn对应的估计曲线为gxn(x)，适用范围为[xn-1, xn+1]；
最大误差计算模块，用于依次获取相邻的所述估计曲线gxn(x)与gxn+1(x)之间的最大误差频点xp，直至得到n-1个所述最大误差频点，每一所述最大误差频点具有误差值；
评价模块，用于根据预设评价方法对每一所述最大误差频点的所述误差值进行计算，得到每一所述最大误差频点的评价指标；
迭代模块，用于根据所述评价指标进行排序，并选择前m个所述最大误差频点作为已采样点，其中，m小于n，m、n均为大于0的整数，且：
若所述已采样点的累计数量未达到所述频点最终采样数n`，则返回所述频率响应估计模块进行迭代执行；
若所述已采样点的累计数量达到所述频点最终采样数n`，则停止迭代；
输出模块，用于输出所有所述已采样点作为最终采样点。
9.一种计算机设备，其特征在于，包括：存储器、处理器及存储在所述存储器上并可在所述处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现如权利要求1至7中任意一项所述的自适应频点采样的扫频方法中的步骤。
10.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质上存储有计算机程序，所述计算机程序被处理器执行时实现如权利要求1至7中任意一项所述的自适应频点采样的扫频方法中的步骤。
CN202211401494.5A 2022-11-09 2022-11-09 自适应频点采样的扫频方法、系统及相关设备 Active CN115455745B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211401494.5A CN115455745B (zh) 2022-11-09 2022-11-09 自适应频点采样的扫频方法、系统及相关设备

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211401494.5A CN115455745B (zh) 2022-11-09 2022-11-09 自适应频点采样的扫频方法、系统及相关设备

Publications (2)

Publication Number Publication Date
CN115455745A CN115455745A (zh) 2022-12-09
CN115455745B true CN115455745B (zh) 2023-04-07

Family

ID=84310907

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211401494.5A Active CN115455745B (zh) 2022-11-09 2022-11-09 自适应频点采样的扫频方法、系统及相关设备

Country Status (1)

Country Link
CN (1) CN115455745B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116136938B (zh) * 2023-04-18 2023-07-04 深圳飞骧科技股份有限公司 声表面波器件仿真参数的快速拟合方法、系统及相关设备

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2014059890A1 (zh) * 2012-10-15 2014-04-24 苏州上声电子有限公司 在线迭代的声重放系统频响均衡方法和装置
WO2017084350A1 (zh) * 2015-11-19 2017-05-26 中兴通讯股份有限公司 一种确定第一反射点距离射频单元的距离的方法和装置

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106199185B (zh) * 2016-08-18 2019-04-16 苏州上声电子股份有限公司 一种基于连续对数扫频的线性脉冲响应测量方法及系统
CN107577872B (zh) * 2017-09-05 2021-02-19 湖南时变通讯科技有限公司 一种时域频率不变波束形成器设计方法及装置
CN112019317B (zh) * 2019-05-28 2022-10-28 中兴通讯股份有限公司 频率校准方法及装置、存储介质、电子装置
CN111585658B (zh) * 2020-05-25 2022-09-06 电子科技大学 一种带宽交替采样系统频率交叠带相位补偿方法
CN114252767A (zh) * 2021-12-10 2022-03-29 歌尔股份有限公司 扫频特性曲线的生成方法、设备及存储介质
CN114996969A (zh) * 2022-07-01 2022-09-02 无锡飞谱电子信息技术有限公司 电磁仿真并行计算方法、装置、计算机设备及存储介质

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2014059890A1 (zh) * 2012-10-15 2014-04-24 苏州上声电子有限公司 在线迭代的声重放系统频响均衡方法和装置
WO2017084350A1 (zh) * 2015-11-19 2017-05-26 中兴通讯股份有限公司 一种确定第一反射点距离射频单元的距离的方法和装置

Also Published As

Publication number Publication date
CN115455745A (zh) 2022-12-09

Similar Documents

Publication Publication Date Title
Sun et al. A novel APSO-aided maximum likelihood identification method for Hammerstein systems
CN115455745B (zh) 自适应频点采样的扫频方法、系统及相关设备
CN112597610B (zh) 机械臂结构轻量化设计的优化方法、装置及设备
CN112231985B (zh) 一种射频滤波器建模方法
CN110032706A (zh) 一种低阶时滞系统的两阶段参数估计方法及系统
CN115481549A (zh) 圆筒直线电机多目标优化方法、设备及存储介质
CN113221298B (zh) 一种机电暂态过程的仿真方法及系统
WO2009049703A1 (en) Method and device for generating a model of a multiparameter system
CN112685841B (zh) 一种带有连接关系的结构的有限元建模、修正方法及系统
CN115542173B (zh) 电池电量估计方法及装置
CN115563840B (zh) 减小单元矩阵的级联误差的仿真方法、系统及相关设备
WO2024099061A1 (zh) 自适应频点采样的扫频方法、系统及相关设备
CN111242379A (zh) 一种基于随机傅里叶特征的核递归最大相关熵时间序列在线预测方法
CN111274665A (zh) 一种模拟裂纹尖端应力场的数值方法和系统
CN111210877A (zh) 一种推断物性参数的方法及装置
CN112147685A (zh) 一种基于波动方程的正演模拟方法及装置
CN115828819B (zh) 一种高精度传输线高频耦合电流计算方法及系统
CN116592816B (zh) 基于旋转机构的数据获取方法、装置、设备及存储介质
CN107247828A (zh) 一种基于逆kriging函数的结构有限元模型修正方法
CN113449382B (zh) 一种整车排放结果的仿真方法及系统
CN116796681A (zh) 一种基于fdtd的低频拟合算法及系统
CN116702571A (zh) 基于多重光滑度量因子的数值模拟方法及装置
CN116455427A (zh) 电力线噪声的时频融合分析方法、装置、设备及存储介质
JP2023100053A (ja) 計算装置、計算プログラム及び計算方法
CN116738820A (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