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

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

Info

Publication number
CN115455745A
CN115455745A CN202211401494.5A CN202211401494A CN115455745A CN 115455745 A CN115455745 A CN 115455745A CN 202211401494 A CN202211401494 A CN 202211401494A CN 115455745 A CN115455745 A CN 115455745A
Authority
CN
China
Prior art keywords
frequency
sampling
points
frequency point
point
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
CN202211401494.5A
Other languages
English (en)
Other versions
CN115455745B (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
Priority to PCT/CN2023/126000 priority patent/WO2024099061A1/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)
  • Algebra (AREA)
  • Geometry (AREA)
  • Computational Mathematics (AREA)
  • Evolutionary Computation (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Computer Hardware Design (AREA)
  • Pure & Applied Mathematics (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • Feedback Control In General (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

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

Description

自适应频点采样的扫频方法、系统及相关设备
技术领域
本发明属于无线通信技术领域，尤其涉及一种自适应频点采样的扫频方法、系统及相关设备。
背景技术
随着无线通讯技术的发展，各类电子消费产品对射频器件的需求越来越大，相对应的，在射频器件设计过程中，设计人员对仿真软件扫频效率与准确度的要求也越来越高。通常情况下，密集的扫频点可以提高扫频曲线的精确度，但同样也会大大增加计算所需的时间。因此需要加速算法在保证计算精度的前提下尽可能的减小计算所需时间。
主流的仿真软件，例如COMSOL，通常会采用渐进波形估计方法（AsymptoticWaveform Evaluation，简称AWE）与自适应兰索斯-帕德扫频方法（ALPS）等方法来提高扫频的计算效率。在有限单元法中，以上方法的基本原理是通过计算几个典型频点下的频率响应，然后通过泰勒级数、Pade有理逼近等数学方法对整体矩阵进行近似展开，从而对相邻频点的频率响应以一种很小的计算开销来估计。问题在于，AWE等方法所基于的泰勒级数展开往往只能在采样频点周围具有较高的精度，而在距离采样频点较远处以及频率响应快速变化处存在频率相应估计精度不高的问题，而频率相应快速变化的区域往往是采样过程中比较关心的区域。
在实际仿真中，往往只对频带部分区域，如通带内的频率响应的精度具有较高要求，而对其他区域的频率响应的精度则没有太高的要求。传统的估计方法则会不加区分的将全部频带划分成若干个较窄的频段，然后对每个频带分别进行扫频分析，因此仍然存在效率较低的问题，因此，传统的估计方法一方面存在频率响应快速变化处可能精度不高的问题，另一方面也存在对于采样中不关心的区域消耗了大量计算资源的问题。并且，传统的估计方法往往依赖于具体的数值模拟算法，如渐进波形估计方法需要将有限元方法的整体矩阵进行分解，从而得到电场值随频率变化的解析表达式，这样的方法对于实验测量、机器学习等不存在显式控制方程的环境不适用。
发明内容
本发明实施例提供一种自适应频点采样的扫频方法、系统及相关设备，旨在解决传统的扫频方法依赖具体数值、不适用于实验测量、机器学习等场景的问题。
第一方面，本发明实施例提供一种自适应频点采样的扫频方法，运用于射频器件的仿真，所述扫频方法包括以下步骤：
S1、确定频点采样数n、频点最终采样数
Figure 293681DEST_PATH_IMAGE001
和采样范围
Figure 895564DEST_PATH_IMAGE002
，并使用预设仿 真方法在所述采样范围中进行采样，并计算n处采样点的频率响应，其中，采样点
Figure 214550DEST_PATH_IMAGE003
处的所 述频率响应为
Figure 319909DEST_PATH_IMAGE004
；
S2、计算所述频率响应
Figure 66148DEST_PATH_IMAGE004
的一阶导数和二阶导数；
S3、根据所述一阶导数、所述二阶导数对n处所述采样点进行频率响应曲线估计， 得到
Figure 624168DEST_PATH_IMAGE005
条估计曲线，其中，采样点
Figure 746845DEST_PATH_IMAGE006
对应的估计曲线为
Figure 706711DEST_PATH_IMAGE007
，适用范围为
Figure 889431DEST_PATH_IMAGE008
；
S4、依次获取相邻的所述估计曲线
Figure 200326DEST_PATH_IMAGE009
与
Figure 861115DEST_PATH_IMAGE010
之间的最大误差频点
Figure 675487DEST_PATH_IMAGE011
，直至得 到
Figure 29108DEST_PATH_IMAGE012
个所述最大误差频点，每一所述最大误差频点具有误差值；
S5、根据预设评价方法对每一所述最大误差频点的所述误差值进行计算，得到每一所述最大误差频点的评价指标；
S6、根据所述评价指标进行排序，并选择前m个所述最大误差频点作为已采样点，其中，m小于n，m、n均为大于0的整数，且：
若所述已采样点的累计数量未达到所述频点最终采样数
Figure 561720DEST_PATH_IMAGE013
，则返回步骤S2进行迭 代执行；
若所述已采样点的累计数量达到所述频点最终采样数
Figure 26200DEST_PATH_IMAGE013
，则停止迭代；
S7、输出所有所述已采样点作为最终采样点。
更进一步地，步骤S1中，所述预设仿真方法在所述采样范围中进行采样的方式为均匀采样。
更进一步地，步骤S2中，计算所述频率响应
Figure 695079DEST_PATH_IMAGE014
的一阶导数和二阶导数所使用的方 法为中心差分法。
更进一步地，步骤S3中，所述估计曲线
Figure 485180DEST_PATH_IMAGE015
满足以下关系式（1）：
Figure 505089DEST_PATH_IMAGE016
（1）。
更进一步地，步骤S4中，所述最大误差频点xp满足以下关系式（2）：
Figure 507680DEST_PATH_IMAGE017
（2）；
所述误差值满足以下关系式（3）：
Figure 296644DEST_PATH_IMAGE018
（3）。
更进一步地，步骤S5中，当所述频率响应的类型为Y参数时，所述预设评价方法为：
S51、在所述采样点中，选择Y参数取值最大的频点最大值
Figure 6716DEST_PATH_IMAGE019
和频点最小值
Figure 248342DEST_PATH_IMAGE020
， 并计算出对应的最大值估计曲线
Figure 54624DEST_PATH_IMAGE021
和最小值估计曲线
Figure 963674DEST_PATH_IMAGE022
；
S52、根据所述最大值估计曲线
Figure 829999DEST_PATH_IMAGE021
和所述最小值估计曲线
Figure 824499DEST_PATH_IMAGE022
得到谐振点
Figure 168893DEST_PATH_IMAGE023
和反谐振点
Figure 932450DEST_PATH_IMAGE024
，并将所述谐振点
Figure 704097DEST_PATH_IMAGE023
和所述反谐振点
Figure 185894DEST_PATH_IMAGE024
加入迭代使用的所述采样点中；
S53、定义所述最大误差频点的所述评价指标为
Figure 333978DEST_PATH_IMAGE025
，所述评价指标
Figure 952041DEST_PATH_IMAGE025
满足以下 关系式（4）：
Figure 425748DEST_PATH_IMAGE026
（4）。
更进一步地，步骤S5中，当所述频率响应的类型为S参数时，所述预设评价方法为：
S51、定义所述频率响应的S参数包括第一S参数和第二S参数；
S52、分别计算所述第一S参数和所述第二S参数的估计曲线的评价指标
Figure 129262DEST_PATH_IMAGE027
和
Figure 815458DEST_PATH_IMAGE028
，其中，所述第一S参数的评价指标
Figure 553607DEST_PATH_IMAGE027
满足以下关系式（5）：
Figure 667056DEST_PATH_IMAGE029
（5）；
所述第二S参数的评价指标F21(xp)满足以下关系式（6）：
Figure 857866DEST_PATH_IMAGE030
（6）；
S53、定义所述最大误差频点的所述评价指标为
Figure 347754DEST_PATH_IMAGE031
，所述评价指标
Figure 940409DEST_PATH_IMAGE031
满足以 下关系式（7）：
Figure 755918DEST_PATH_IMAGE032
（7）。
第二方面，本发明实施例还提供一种自适应频点采样的扫频系统，运用于射频器件的仿真，包括：
初始化模块，用于确定频点采样数n、频点最终采样数
Figure 434024DEST_PATH_IMAGE033
和采样范围
Figure 462023DEST_PATH_IMAGE034
，并 使用预设仿真方法在所述采样范围中进行采样，并计算n处采样点的频率响应，其中，采样 点
Figure 909185DEST_PATH_IMAGE035
处的所述频率响应为
Figure 630016DEST_PATH_IMAGE036
；
求导模块，用于计算所述频率响应
Figure 795418DEST_PATH_IMAGE036
的一阶导数和二阶导数；
频率响应估计模块，用于根据所述一阶导数、所述二阶导数对n处所述采样点进行 频率响应曲线估计，得到
Figure 627108DEST_PATH_IMAGE037
条估计曲线，其中，采样点
Figure 928776DEST_PATH_IMAGE038
对应的估计曲线为
Figure 820509DEST_PATH_IMAGE039
，适用范 围为
Figure 741716DEST_PATH_IMAGE040
；
最大误差计算模块，用于依次获取相邻的所述估计曲线
Figure 111518DEST_PATH_IMAGE041
与
Figure 533272DEST_PATH_IMAGE042
之间的最大 误差频点
Figure 595906DEST_PATH_IMAGE043
，直至得到
Figure 470321DEST_PATH_IMAGE044
个所述最大误差频点，每一所述最大误差频点具有误差值；
评价模块，用于根据预设评价方法对每一所述最大误差频点的所述误差值进行计算，得到每一所述最大误差频点的评价指标；
迭代模块，用于根据所述评价指标进行排序，并选择前m个所述最大误差频点作为已采样点，其中，m小于n，m、n均为大于0的整数，且：
若所述已采样点的累计数量未达到所述频点最终采样数
Figure 909392DEST_PATH_IMAGE045
，则返回所述频率响应 估计模块进行迭代执行；
若所述已采样点的累计数量达到所述频点最终采样数
Figure 185653DEST_PATH_IMAGE045
，则停止迭代；
输出模块，用于输出所有所述已采样点作为最终采样点。
第三方面，本发明实施例还提供一种计算机设备，包括：存储器、处理器及存储在所述存储器上并可在所述处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现如上述实施例中任意一项所述的自适应频点采样的方法中的步骤。
第四方面，本发明实施例还提供一种计算机可读存储介质，所述计算机可读存储介质上存储有计算机程序，所述计算机程序被处理器执行时实现如上述实施例中任意一项所述的自适应频点采样的方法中的步骤。
本发明所达到的有益效果，由于采用了一种不依赖于具体数值方法自适应频点采样方法，能够对频率响应快速变化处进行精确仿真，而在其他不感兴趣的区域进行粗略仿真，从而提高射频器件设计过程的精度与效率；同时，针对频率响应中Y参数与S参数的特点使用不同的评价指标，从而减少扫频所需要的采样点，提高扫频速度。
附图说明
图1是本发明实施例提供的自适应频点采样的扫频方法的步骤流程框图；
图2是本发明实施例提供的自适应频点采样的扫频方法在S参数中采样的示意图；
图3是本发明实施例提供的自适应频点采样的扫频方法在Y参数中采样的示意图；
图4是本发明实施例提供的自适应频点采样的扫频系统的结构示意图；
图5是本发明实施例提供的计算机设备的结构示意图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。
请参照图1，图1是本发明实施例提供的自适应频点采样的扫频方法的步骤流程框图，所述扫频方法运用于射频器件的仿真，所述扫频方法包括以下步骤：
S1、确定频点采样数n、频点最终采样数
Figure 684767DEST_PATH_IMAGE046
和采样范围
Figure 46479DEST_PATH_IMAGE047
，并使用预设仿真 方法在所述采样范围中进行采样，并计算n处采样点的频率响应，其中，采样点
Figure 758083DEST_PATH_IMAGE048
处的所述 频率响应为
Figure 154429DEST_PATH_IMAGE049
。
更进一步地，步骤S1中，所述预设仿真方法在所述采样范围中进行采样的方式为均匀采样。
具体的，所述预设仿真方法可以是任意一种数值仿真计算方法，例如Matlab等。
S2、计算所述频率响应
Figure 558865DEST_PATH_IMAGE049
的一阶导数和二阶导数。
更进一步地，步骤S2中，计算所述频率响应
Figure 407873DEST_PATH_IMAGE049
的一阶导数和二阶导数所使用的方 法为中心差分法。
S3、根据所述一阶导数、所述二阶导数对n处所述采样点进行频率响应曲线估计， 得到
Figure 923168DEST_PATH_IMAGE050
条估计曲线，其中，采样点
Figure 439600DEST_PATH_IMAGE051
对应的估计曲线为
Figure 14938DEST_PATH_IMAGE052
，适用范围为
Figure 351241DEST_PATH_IMAGE053
。
具体的，根据所述一阶导数、所述二阶导数对n处所述采样点进行频率响应曲线估计时，按照每两个相邻的所述采样点进行，并且可以进一步通过泰勒公式进行估计，从而得到所述估计曲线。
更进一步地，步骤S3中，所述估计曲线
Figure 670227DEST_PATH_IMAGE054
满足以下关系式（1）：
Figure 775586DEST_PATH_IMAGE055
（1）。
S4、依次获取相邻的所述估计曲线
Figure 521825DEST_PATH_IMAGE056
与
Figure 79845DEST_PATH_IMAGE057
之间的最大误差频点
Figure 468102DEST_PATH_IMAGE058
，直至得 到
Figure 427967DEST_PATH_IMAGE059
个所述最大误差频点，每一所述最大误差频点具有误差值。
具体的，对于任意范围
Figure 345108DEST_PATH_IMAGE060
中的一点
Figure 656003DEST_PATH_IMAGE061
，其都可以通过其两侧采样点的估计 曲线
Figure 316792DEST_PATH_IMAGE062
与
Figure 131164DEST_PATH_IMAGE063
进行预测，即
Figure 484785DEST_PATH_IMAGE064
与
Figure 285906DEST_PATH_IMAGE065
。
同时，为了便于后续评价指标及迭代的计算，本发明实施例通过数据验证，将所述 最大误差频点
Figure 484807DEST_PATH_IMAGE066
取值为其计算使用的所述估计曲线对应的所述采样点的中点。
更进一步地，步骤S4中，所述最大误差频点xp满足以下关系式（2）：
Figure 419264DEST_PATH_IMAGE067
（2）；
所述误差值满足以下关系式（3）：
Figure 943787DEST_PATH_IMAGE068
（3）。
S5、根据预设评价方法对每一所述最大误差频点的所述误差值进行计算，得到每一所述最大误差频点的评价指标。
具体的，所述频率响应的参数类型可以为多种，包括S参数和Y参数，不同参数虽然本质上都等效为二端口网络之后的二端口网络参数，但得到数值的过程不同，在本发明实施例中，根据数值类型的不同，所述评价指标的计算方法也不同。
更进一步地，步骤S5中，当所述频率响应的类型为Y参数时，所述预设评价方法为：
S51、在所述采样点中，选择Y参数取值最大的频点最大值
Figure 698116DEST_PATH_IMAGE069
和频点最小值
Figure 966286DEST_PATH_IMAGE070
， 并计算出对应的最大值估计曲线
Figure 755251DEST_PATH_IMAGE071
和最小值估计曲线
Figure 185095DEST_PATH_IMAGE072
；
S52、根据所述最大值估计曲线
Figure 692300DEST_PATH_IMAGE071
和所述最小值估计曲线
Figure 764161DEST_PATH_IMAGE072
得到谐振点
Figure 407632DEST_PATH_IMAGE073
和反谐振点
Figure 273957DEST_PATH_IMAGE074
，并将所述谐振点
Figure 2878DEST_PATH_IMAGE073
和所述反谐振点
Figure 612851DEST_PATH_IMAGE074
加入迭代使用的所述采样点中；
S53、定义所述最大误差频点的所述评价指标为
Figure 376408DEST_PATH_IMAGE075
，所述评价指标
Figure 148055DEST_PATH_IMAGE075
满足以 下关系式（4）：
Figure 895431DEST_PATH_IMAGE076
（4）。
所述谐振点
Figure 43516DEST_PATH_IMAGE077
和所述反谐振点
Figure 661579DEST_PATH_IMAGE078
加入迭代使用的所述采样点后，采样点所体现的 采集方式会逐渐向所述谐振点
Figure 604127DEST_PATH_IMAGE077
和所述反谐振点
Figure 307641DEST_PATH_IMAGE078
的位置靠近，从而实现目标兴趣频带的 扫频。
对应的，本发明还提供另一种实施例，步骤S5中，当所述频率响应的类型为S参数时，所述预设评价方法为：
S51、定义所述频率响应的S参数包括第一S参数和第二S参数；
S52、分别计算所述第一S参数和所述第二S参数的估计曲线的评价指标
Figure 259416DEST_PATH_IMAGE079
和
Figure 731986DEST_PATH_IMAGE080
，其中，所述第一S参数的评价指标
Figure 111015DEST_PATH_IMAGE079
满足以下关系式（5）：
Figure 301825DEST_PATH_IMAGE081
（5）；
所述第二S参数的评价指标F21(xp)满足以下关系式（6）：
Figure 57291DEST_PATH_IMAGE082
（6）；
S53、定义所述最大误差频点的所述评价指标为
Figure 649946DEST_PATH_IMAGE031
，所述评价指标
Figure 199876DEST_PATH_IMAGE031
满足以 下关系式（7）：
Figure 615333DEST_PATH_IMAGE083
（7）。
S6、根据所述评价指标进行排序，并选择前m个所述最大误差频点作为已采样点，其中，m小于n，m、n均为大于0的整数，且：
若所述已采样点的累计数量未达到所述频点最终采样数
Figure 908911DEST_PATH_IMAGE084
，则返回步骤S2进行迭 代执行；
若所述已采样点的累计数量达到所述频点最终采样数
Figure 356073DEST_PATH_IMAGE084
，则停止迭代。
具体的，本发明实施例设定所述频点最终采样数
Figure 76904DEST_PATH_IMAGE084
，并将步骤S2至步骤S6作为一 次完整的迭代过程，其中，在每次迭代后，在原本划分的所述采样点的基础上再进行频率估 计，会得到比迭代前一次使用的所述采样点的数量更多的用于计算的所述采样点，通过多 次迭代，并将每次迭代中满足所述评价指标的所述最大误差频点作为选用的采样点，从而 在较少频点情况下对目标频带进行高效地采样。
S7、输出所有所述已采样点作为最终采样点。
示例性的，请参照图2和图3，其分别是本发明实施例提供的自适应频点采样的扫频方法在S参数和Y参数中采样的示意图，其中，对于S参数，在通带内采样点的连线与真实曲线具有较高的一致性；对于Y参数，采样点的连线与真实曲线整体具有较高的一致性；图2、图3中下方的散点为每次采样时选择的点，其采样批次呈层级变化，可以看出，第一批次采样点与第二批次采样点均体现为均匀采样，但第三批次及之后的采样点则集中在了较为集中的区域，在图2及图3中，第三批次及之后的采样点均向S参数的带内区域以及Y参数的谐振点、反谐振点处靠近。
本发明所达到的有益效果，由于采用了一种不依赖于具体数值方法自适应频点采样方法，能够对频率响应快速变化处进行精确仿真，而在其他不感兴趣的区域进行粗略仿真，从而提高射频器件设计过程的精度与效率；同时，针对频率响应中Y参数与S参数的特点使用不同的评价指标，从而减少扫频所需要的采样点，提高扫频速度。
本发明实施例还提供一种自适应频点采样的扫频系统，请参照图4，图4是本发明实施例提供的自适应频点采样的扫频系统的结构示意图，所述自适应频点采样的扫频系统200包括：
初始化模块201，用于确定频点采样数n、频点最终采样数
Figure DEST_PATH_IMAGE085
和采样范围
Figure 773465DEST_PATH_IMAGE086
，并使用预设仿真方法在所述采样范围中进行采样，并计算n处采样点的频率响应，其中，采 样点
Figure 605155DEST_PATH_IMAGE087
处的所述频率响应为
Figure 906823DEST_PATH_IMAGE088
；
求导模块202，用于计算所述频率响应
Figure 798555DEST_PATH_IMAGE088
的一阶导数和二阶导数；
频率响应估计模块203，用于根据所述一阶导数、所述二阶导数对n处所述采样点 进行频率响应曲线估计，得到
Figure 185674DEST_PATH_IMAGE089
条估计曲线，其中，采样点
Figure 821055DEST_PATH_IMAGE090
对应的估计曲线为
Figure 242809DEST_PATH_IMAGE091
，适 用范围为
Figure 305443DEST_PATH_IMAGE092
；
最大误差计算模块204，用于依次获取相邻的所述估计曲线
Figure 179858DEST_PATH_IMAGE093
与
Figure 618930DEST_PATH_IMAGE094
之间的 最大误差频点
Figure 895190DEST_PATH_IMAGE095
，直至得到
Figure 128726DEST_PATH_IMAGE096
个所述最大误差频点，每一所述最大误差频点具有误差值；
评价模块205，用于根据预设评价方法对每一所述最大误差频点的所述误差值进行计算，得到每一所述最大误差频点的评价指标；
迭代模块206，用于根据所述评价指标进行排序，并选择前m个所述最大误差频点作为已采样点，其中，m小于n，m、n均为大于0的整数，且：
若所述已采样点的累计数量未达到所述频点最终采样数
Figure 490437DEST_PATH_IMAGE097
，则返回所述频率响应 估计模块进行迭代执行；
若所述已采样点的累计数量达到所述频点最终采样数
Figure 467620DEST_PATH_IMAGE097
，则停止迭代；
输出模块207，用于输出所有所述已采样点作为最终采样点。
所述自适应频点采样的扫频系统200能够实现如上述实施例中的自适应频点采样的扫频方法中的步骤，且能实现同样的技术效果，参上述实施例中的描述，此处不再赘述。
本发明实施例还提供一种计算机设备，请参照图5，图5是本发明实施例提供的计算机设备的结构示意图，所述计算机设备300包括：存储器302、处理器301及存储在所述存储器302上并可在所述处理器301上运行的计算机程序。
所述处理器301调用所述存储器302存储的计算机程序，执行本发明实施例提供的自适应频点采样的扫频方法中的步骤，请结合图1，具体包括：
S1、确定频点采样数n、频点最终采样数
Figure 863966DEST_PATH_IMAGE097
和采样范围
Figure 268403DEST_PATH_IMAGE098
，并使用预设仿真 方法在所述采样范围中进行采样，并计算n处采样点的频率响应，其中，采样点
Figure 851831DEST_PATH_IMAGE099
处的所述 频率响应为
Figure 632705DEST_PATH_IMAGE100
；
S2、计算所述频率响应
Figure 149137DEST_PATH_IMAGE100
的一阶导数和二阶导数；
S3、根据所述一阶导数、所述二阶导数对n处所述采样点进行频率响应曲线估计， 得到
Figure 724475DEST_PATH_IMAGE101
条估计曲线，其中，采样点
Figure 795199DEST_PATH_IMAGE102
对应的估计曲线为
Figure 114185DEST_PATH_IMAGE103
，适用范围为
Figure 219544DEST_PATH_IMAGE104
；
S4、依次获取相邻的所述估计曲线
Figure 965783DEST_PATH_IMAGE105
与
Figure 46173DEST_PATH_IMAGE106
之间的最大误差频点
Figure 903271DEST_PATH_IMAGE107
，直至得 到
Figure 128716DEST_PATH_IMAGE108
个所述最大误差频点，每一所述最大误差频点具有误差值；
S5、根据预设评价方法对每一所述最大误差频点的所述误差值进行计算，得到每一所述最大误差频点的评价指标；
S6、根据所述评价指标进行排序，并选择前m个所述最大误差频点作为已采样点，其中，m小于n，m、n均为大于0的整数，且：
若所述已采样点的累计数量未达到所述频点最终采样数
Figure 45856DEST_PATH_IMAGE109
，则返回步骤S2进行迭 代执行；
若所述已采样点的累计数量达到所述频点最终采样数
Figure 91173DEST_PATH_IMAGE109
，则停止迭代；
S7、输出所有所述已采样点作为最终采样点。
更进一步地，步骤S1中，所述预设仿真方法在所述采样范围中进行采样的方式为均匀采样。
更进一步地，步骤S2中，计算所述频率响应
Figure 751961DEST_PATH_IMAGE110
的一阶导数和二阶导数所使用的方 法为中心差分法。
更进一步地，步骤S3中，所述估计曲线
Figure 831912DEST_PATH_IMAGE111
满足以下关系式（1）：
Figure 919954DEST_PATH_IMAGE112
（1）。
更进一步地，步骤S4中，所述最大误差频点
Figure 452567DEST_PATH_IMAGE113
满足以下关系式（2）：
Figure 917046DEST_PATH_IMAGE114
（2）；
所述误差值满足以下关系式（3）：
Figure 117083DEST_PATH_IMAGE115
（3）。
更进一步地，步骤S5中，当所述频率响应的类型为Y参数时，所述预设评价方法为：
S51、在所述采样点中，选择Y参数取值最大的频点最大值
Figure 641605DEST_PATH_IMAGE116
和频点最小值
Figure 661514DEST_PATH_IMAGE117
， 并计算出对应的最大值估计曲线
Figure 929684DEST_PATH_IMAGE118
和最小值估计曲线
Figure 718649DEST_PATH_IMAGE119
；
S52、根据所述最大值估计曲线
Figure 148493DEST_PATH_IMAGE118
和所述最小值估计曲线
Figure 655698DEST_PATH_IMAGE119
得到谐振点
Figure 727559DEST_PATH_IMAGE120
和反谐振点
Figure 371030DEST_PATH_IMAGE121
，并将所述谐振点
Figure 237355DEST_PATH_IMAGE120
和所述反谐振点
Figure 966276DEST_PATH_IMAGE121
加入迭代使用的所述采样点中；
S53、定义所述最大误差频点的所述评价指标为
Figure 576249DEST_PATH_IMAGE122
，所述评价指标
Figure 74227DEST_PATH_IMAGE122
满足以 下关系式（4）：
Figure 377032DEST_PATH_IMAGE123
（4）。
更进一步地，步骤S5中，当所述频率响应的类型为S参数时，所述预设评价方法为：
S51、定义所述频率响应的S参数包括第一S参数和第二S参数；
S52、分别计算所述第一S参数和所述第二S参数的估计曲线的评价指标
Figure 593250DEST_PATH_IMAGE124
和
Figure 741334DEST_PATH_IMAGE125
，其中，所述第一S参数的评价指标
Figure 359398DEST_PATH_IMAGE126
满足以下关系式（5）：
Figure 570455DEST_PATH_IMAGE127
（5）；
所述第二S参数的评价指标
Figure 539548DEST_PATH_IMAGE125
满足以下关系式（6）：
Figure 491323DEST_PATH_IMAGE128
（6）；
S53、定义所述最大误差频点的所述评价指标为
Figure 229472DEST_PATH_IMAGE129
，所述评价指标
Figure 874080DEST_PATH_IMAGE129
满足以下 关系式（7）：
Figure 64890DEST_PATH_IMAGE130
（7）。
本发明实施例提供的计算机设备300能够实现如上述实施例中的自适应频点采样的扫频方法中的步骤，且能实现同样的技术效果，参上述实施例中的描述，此处不再赘述。
本发明实施例还提供一种计算机可读存储介质，所述计算机可读存储介质上存储有计算机程序，该计算机程序被处理器执行时实现本发明实施例提供的自适应频点采样的扫频方法中的各个过程及步骤，且能实现相同的技术效果，为避免重复，这里不再赘述。
本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程，是可以通过计算机程序来指令相关的硬件来完成，所述的程序可存储于一计算机可读取存储介质中，该程序在执行时，可包括如上述各方法的实施例的流程。其中，所述的存储介质可为磁碟、光盘、只读存储记忆体（Read-Only Memory，ROM）或随机存取存储器（Random AccessMemory，简称RAM）等。
需要说明的是，在本文中，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、物品或者装置不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、物品或者装置所固有的要素。在没有更多限制的情况下，由语句“包括一个……”限定的要素，并不排除在包括该要素的过程、方法、物品或者装置中还存在另外的相同要素。
通过以上的实施方式的描述，本领域的技术人员可以清楚地了解到上述实施例方法可借助软件加必需的通用硬件平台的方式来实现，当然也可以通过硬件，但很多情况下前者是更佳的实施方式。基于这样的理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质（如ROM/RAM、磁碟、光盘）中，包括若干指令用以使得一台终端（可以是手机，计算机，服务器，空调器，或者网络设备等）执行本发明各个实施例所述的方法。
上面结合附图对本发明的实施例进行了描述，所揭露的仅为本发明较佳实施例而已，但是本发明并不局限于上述的具体实施方式，上述的具体实施方式仅仅是示意性的，而不是限制性的，本领域的普通技术人员在本发明的启示下，在不脱离本发明宗旨和权利要求所保护的范围情况下，还可做出很多形式用等同变化，均属于本发明的保护之内。

Claims (10)

1.一种自适应频点采样的扫频方法，运用于射频器件的仿真，其特征在于，所述扫频方法包括以下步骤：
S1、确定频点采样数n、频点最终采样数
Figure 217846DEST_PATH_IMAGE001
和采样范围
Figure 680793DEST_PATH_IMAGE002
，并使用预设仿真方法 在所述采样范围中进行采样，并计算n处采样点的频率响应，其中，采样点
Figure 726110DEST_PATH_IMAGE003
处的所述频率 响应为
Figure 652478DEST_PATH_IMAGE004
；
S2、计算所述频率响应
Figure 732429DEST_PATH_IMAGE004
的一阶导数和二阶导数；
S3、根据所述一阶导数、所述二阶导数对n处所述采样点进行频率响应曲线估计，得到
Figure 820471DEST_PATH_IMAGE005
条估计曲线，其中，采样点
Figure 353083DEST_PATH_IMAGE006
对应的估计曲线为
Figure 817563DEST_PATH_IMAGE007
，适用范围为
Figure 752021DEST_PATH_IMAGE008
；
S4、依次获取相邻的所述估计曲线
Figure 10964DEST_PATH_IMAGE009
与
Figure 296451DEST_PATH_IMAGE010
之间的最大误差频点
Figure 299042DEST_PATH_IMAGE011
，直至得到
Figure 88007DEST_PATH_IMAGE012
个所述最大误差频点，每一所述最大误差频点具有误差值；
S5、根据预设评价方法对每一所述最大误差频点的所述误差值进行计算，得到每一所述最大误差频点的评价指标；
S6、根据所述评价指标进行排序，并选择前m个所述最大误差频点作为已采样点，其中，m小于n，m、n均为大于0的整数，且：
若所述已采样点的累计数量未达到所述频点最终采样数
Figure 783430DEST_PATH_IMAGE013
，则返回步骤S2进行迭代执 行；
若所述已采样点的累计数量达到所述频点最终采样数
Figure 290635DEST_PATH_IMAGE013
，则停止迭代；
S7、输出所有所述已采样点作为最终采样点。
2.如权利要求1所述的扫频方法，其特征在于，步骤S1中，所述预设仿真方法在所述采样范围中进行采样的方式为均匀采样。
3.如权利要求1所述的扫频方法，其特征在于，步骤S2中，计算所述频率响应
Figure 831338DEST_PATH_IMAGE014
的一阶 导数和二阶导数所使用的方法为中心差分法。
4.如权利要求1所述的扫频方法，其特征在于，步骤S3中，所述估计曲线
Figure 740388DEST_PATH_IMAGE015
满足以下 关系式（1）：
Figure 606713DEST_PATH_IMAGE016
（1）。
5.如权利要求4所述的扫频方法，其特征在于，步骤S4中，所述最大误差频点
Figure 601214DEST_PATH_IMAGE017
满足以 下关系式（2）：
Figure 945607DEST_PATH_IMAGE018
（2）；
所述误差值满足以下关系式（3）：
Figure 709164DEST_PATH_IMAGE019
（3）。
6.如权利要求5所述的扫频方法，其特征在于，步骤S5中，当所述频率响应的类型为Y参数时，所述预设评价方法为：
S51、在所述采样点中，选择Y参数取值最大的频点最大值
Figure 11969DEST_PATH_IMAGE020
和频点最小值
Figure 228187DEST_PATH_IMAGE021
，并计 算出对应的最大值估计曲线
Figure 376272DEST_PATH_IMAGE022
和最小值估计曲线
Figure 994335DEST_PATH_IMAGE023
；
S52、根据所述最大值估计曲线
Figure 202462DEST_PATH_IMAGE024
和所述最小值估计曲线
Figure 905976DEST_PATH_IMAGE025
得到谐振点
Figure 592172DEST_PATH_IMAGE026
和反谐振点
Figure 598830DEST_PATH_IMAGE027
，并将所述谐振点
Figure 977859DEST_PATH_IMAGE026
和所述反谐振点
Figure 168669DEST_PATH_IMAGE027
加入迭代使用的所述采样点中；
S53、定义所述最大误差频点的所述评价指标为
Figure 658556DEST_PATH_IMAGE028
，所述评价指标
Figure 251211DEST_PATH_IMAGE029
满足以下关 系式（4）：
Figure 801141DEST_PATH_IMAGE030
（4）。
7.如权利要求5所述的扫频方法，其特征在于，步骤S5中，当所述频率响应的类型为S参数时，所述预设评价方法为：
S51、定义所述频率响应的S参数包括第一S参数和第二S参数；
S52、分别计算所述第一S参数和所述第二S参数的估计曲线的评价指标
Figure 479247DEST_PATH_IMAGE031
和
Figure 772825DEST_PATH_IMAGE032
，其中，所述第一S参数的评价指标
Figure 219987DEST_PATH_IMAGE033
满足以下关系式（5）：
Figure 675239DEST_PATH_IMAGE034
（5）；
所述第二S参数的评价指标F21(xp)满足以下关系式（6）：
Figure 106221DEST_PATH_IMAGE035
（6）；
S53、定义所述最大误差频点的所述评价指标为
Figure 937911DEST_PATH_IMAGE036
，所述评价指标
Figure 505158DEST_PATH_IMAGE036
满足以下关 系式（7）：
Figure 131312DEST_PATH_IMAGE037
（7）。
8.一种自适应频点采样的扫频系统，运用于射频器件的仿真，其特征在于，包括：
初始化模块，用于确定频点采样数n、频点最终采样数
Figure 784010DEST_PATH_IMAGE038
和采样范围
Figure 419390DEST_PATH_IMAGE039
，并使 用预设仿真方法在所述采样范围中进行采样，并计算n处采样点的频率响应，其中，采样点
Figure 575565DEST_PATH_IMAGE040
处的所述频率响应为
Figure 903778DEST_PATH_IMAGE041
；
求导模块，用于计算所述频率响应
Figure 778194DEST_PATH_IMAGE041
的一阶导数和二阶导数；
频率响应估计模块，用于根据所述一阶导数、所述二阶导数对n处所述采样点进行频率 响应曲线估计，得到
Figure 482844DEST_PATH_IMAGE042
条估计曲线，其中，采样点
Figure 24684DEST_PATH_IMAGE043
对应的估计曲线为
Figure 258219DEST_PATH_IMAGE044
，适用范围 为
Figure 619931DEST_PATH_IMAGE045
；
最大误差计算模块，用于依次获取相邻的所述估计曲线
Figure 597114DEST_PATH_IMAGE046
与
Figure 727881DEST_PATH_IMAGE047
之间的最大误 差频点
Figure 132317DEST_PATH_IMAGE048
，直至得到
Figure 981325DEST_PATH_IMAGE049
个所述最大误差频点，每一所述最大误差频点具有误差值；
评价模块，用于根据预设评价方法对每一所述最大误差频点的所述误差值进行计算，得到每一所述最大误差频点的评价指标；
迭代模块，用于根据所述评价指标进行排序，并选择前m个所述最大误差频点作为已采样点，其中，m小于n，m、n均为大于0的整数，且：
若所述已采样点的累计数量未达到所述频点最终采样数
Figure 765129DEST_PATH_IMAGE050
，则返回所述频率响应估计 模块进行迭代执行；
若所述已采样点的累计数量达到所述频点最终采样数
Figure 15981DEST_PATH_IMAGE051
，则停止迭代；
输出模块，用于输出所有所述已采样点作为最终采样点。
9.一种计算机设备，其特征在于，包括：存储器、处理器及存储在所述存储器上并可在所述处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现如权利要求1至7中任意一项所述的自适应频点采样的方法中的步骤。
10.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质上存储有计算机程序，所述计算机程序被处理器执行时实现如权利要求1至7中任意一项所述的自适应频点采样的方法中的步骤。
CN202211401494.5A 2022-11-09 2022-11-09 自适应频点采样的扫频方法、系统及相关设备 Active CN115455745B (zh)

Priority Applications (2)

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

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 true CN115455745A (zh) 2022-12-09
CN115455745B 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 (2)

Country Link
CN (1) CN115455745B (zh)
WO (1) WO2024099061A1 (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116136938A (zh) * 2023-04-18 2023-05-19 深圳飞骧科技股份有限公司 声表面波器件仿真参数的快速拟合方法、系统及相关设备
WO2024099061A1 (zh) * 2022-11-09 2024-05-16 深圳飞骧科技股份有限公司 自适应频点采样的扫频方法、系统及相关设备

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2014059890A1 (zh) * 2012-10-15 2014-04-24 苏州上声电子有限公司 在线迭代的声重放系统频响均衡方法和装置
CN106199185A (zh) * 2016-08-18 2016-12-07 苏州上声电子有限公司 一种基于连续对数扫频的线性脉冲响应测量方法及系统
WO2017084350A1 (zh) * 2015-11-19 2017-05-26 中兴通讯股份有限公司 一种确定第一反射点距离射频单元的距离的方法和装置
CN107577872A (zh) * 2017-09-05 2018-01-12 湖南时变通讯科技有限公司 一种频率不变波束形成器设计方法及装置
CN111585658A (zh) * 2020-05-25 2020-08-25 电子科技大学 一种带宽交替采样系统频率交叠带相位补偿方法
WO2020238780A1 (zh) * 2019-05-28 2020-12-03 中兴通讯股份有限公司 频率校准方法及装置、存储介质、电子装置
CN114252767A (zh) * 2021-12-10 2022-03-29 歌尔股份有限公司 扫频特性曲线的生成方法、设备及存储介质
CN114996969A (zh) * 2022-07-01 2022-09-02 无锡飞谱电子信息技术有限公司 电磁仿真并行计算方法、装置、计算机设备及存储介质

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2013010266A1 (en) * 2011-07-21 2013-01-24 Enbala Power Networks Inc. Method and system for providing local primary frequency response
CN111896802B (zh) * 2020-06-30 2023-03-31 成都四威功率电子科技有限公司 一种频率自适应的采样方法
CN114139372A (zh) * 2021-11-29 2022-03-04 杭州电子科技大学 一种使用高效自适应频率扫描的宽带电磁仿真方法
CN114236383A (zh) * 2021-12-10 2022-03-25 歌尔股份有限公司 扫频特性曲线的生成方法、设备及存储介质
CN115455745B (zh) * 2022-11-09 2023-04-07 深圳飞骧科技股份有限公司 自适应频点采样的扫频方法、系统及相关设备

Patent Citations (8)

* 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 中兴通讯股份有限公司 一种确定第一反射点距离射频单元的距离的方法和装置
CN106199185A (zh) * 2016-08-18 2016-12-07 苏州上声电子有限公司 一种基于连续对数扫频的线性脉冲响应测量方法及系统
CN107577872A (zh) * 2017-09-05 2018-01-12 湖南时变通讯科技有限公司 一种频率不变波束形成器设计方法及装置
WO2020238780A1 (zh) * 2019-05-28 2020-12-03 中兴通讯股份有限公司 频率校准方法及装置、存储介质、电子装置
CN111585658A (zh) * 2020-05-25 2020-08-25 电子科技大学 一种带宽交替采样系统频率交叠带相位补偿方法
CN114252767A (zh) * 2021-12-10 2022-03-29 歌尔股份有限公司 扫频特性曲线的生成方法、设备及存储介质
CN114996969A (zh) * 2022-07-01 2022-09-02 无锡飞谱电子信息技术有限公司 电磁仿真并行计算方法、装置、计算机设备及存储介质

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
何十全等: "基于稳健自适应频率采样的宽频电磁数值仿真", 《电波科学学报》 *
韩丽娟等: "基于PSpice的变压器绕组变形参数仿真与分析", 《宁夏工程技术》 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2024099061A1 (zh) * 2022-11-09 2024-05-16 深圳飞骧科技股份有限公司 自适应频点采样的扫频方法、系统及相关设备
CN116136938A (zh) * 2023-04-18 2023-05-19 深圳飞骧科技股份有限公司 声表面波器件仿真参数的快速拟合方法、系统及相关设备
CN116136938B (zh) * 2023-04-18 2023-07-04 深圳飞骧科技股份有限公司 声表面波器件仿真参数的快速拟合方法、系统及相关设备

Also Published As

Publication number Publication date
CN115455745B (zh) 2023-04-07
WO2024099061A1 (zh) 2024-05-16

Similar Documents

Publication Publication Date Title
CN115455745B (zh) 自适应频点采样的扫频方法、系统及相关设备
Sun et al. A novel APSO-aided maximum likelihood identification method for Hammerstein systems
CN115577604B (zh) 快速构建整体矩阵的仿真方法、系统及相关设备
CN112597610B (zh) 机械臂结构轻量化设计的优化方法、装置及设备
CN116136938B (zh) 声表面波器件仿真参数的快速拟合方法、系统及相关设备
CN112231985B (zh) 一种射频滤波器建模方法
US20100280801A1 (en) Method and device for generating a model of a multiparameter system
CN107562714B (zh) 一种语句相似度计算方法及装置
CN111460734A (zh) 一种结合先进自适应采样和人工神经网络的微波器件自动建模方法
CN115542173B (zh) 电池电量估计方法及装置
CN115563840B (zh) 减小单元矩阵的级联误差的仿真方法、系统及相关设备
CN111967184A (zh) 一种基于序列模型的多目标天线设计方法
CN116776695A (zh) 基于超高阶有限元技术的一维电磁计算方法、系统及设备
CN114722341A (zh) 一种用于集成电路的汉克尔变换滤波器的改进方法及装置
CN111210877A (zh) 一种推断物性参数的方法及装置
CN115828819B (zh) 一种高精度传输线高频耦合电流计算方法及系统
CN111274665A (zh) 一种模拟裂纹尖端应力场的数值方法和系统
CN113449382B (zh) 一种整车排放结果的仿真方法及系统
CN107247828A (zh) 一种基于逆kriging函数的结构有限元模型修正方法
CN116055340B (zh) 分布式网络未知参数估计方法、装置及电子设备
CN111062159B (zh) 一种无人机虚拟试验方法
CN116796681A (zh) 一种基于fdtd的低频拟合算法及系统
CN117610386A (zh) 基于有限体积法的高频电磁场仿真方法、系统及计算机
CN117829048A (zh) 基于神经网络的数据插值方法、系统、设备及介质
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