CN111596544A - 一种基于模糊pi的频率跟踪搜索方法 - Google Patents
一种基于模糊pi的频率跟踪搜索方法 Download PDFInfo
- Publication number
- CN111596544A CN111596544A CN202010708046.4A CN202010708046A CN111596544A CN 111596544 A CN111596544 A CN 111596544A CN 202010708046 A CN202010708046 A CN 202010708046A CN 111596544 A CN111596544 A CN 111596544A
- Authority
- CN
- China
- Prior art keywords
- ultrasonic transducer
- frequency
- phase difference
- fuzzy
- change
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 30
- 238000004364 calculation method Methods 0.000 claims description 13
- 238000002604 ultrasonography Methods 0.000 claims description 7
- 150000001875 compounds Chemical class 0.000 claims description 6
- 230000005484 gravity Effects 0.000 claims description 3
- 230000006870 function Effects 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 3
- 238000010408 sweeping Methods 0.000 description 3
- 238000013461 design Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000007476 Maximum Likelihood Methods 0.000 description 1
- 230000032683 aging Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000002925 chemical effect Effects 0.000 description 1
- 238000004140 cleaning Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000012905 input function Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000003466 welding Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B11/00—Automatic controllers
- G05B11/01—Automatic controllers electric
- G05B11/36—Automatic controllers electric with provision for obtaining particular characteristics, e.g. proportional, integral, differential
- G05B11/42—Automatic controllers electric with provision for obtaining particular characteristics, e.g. proportional, integral, differential for obtaining a characteristic which is both proportional and time-dependent, e.g. P. I., P. I. D.
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B13/00—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
- G05B13/02—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
- G05B13/04—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators
- G05B13/042—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators in which a parameter or coefficient is automatically adjusted to optimise the performance
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N5/00—Computing arrangements using knowledge-based models
- G06N5/04—Inference or reasoning models
- G06N5/048—Fuzzy inferencing
Landscapes
- Engineering & Computer Science (AREA)
- Software Systems (AREA)
- Automation & Control Theory (AREA)
- General Physics & Mathematics (AREA)
- Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- Artificial Intelligence (AREA)
- Theoretical Computer Science (AREA)
- Medical Informatics (AREA)
- Health & Medical Sciences (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Fuzzy Systems (AREA)
- Computational Linguistics (AREA)
- Data Mining & Analysis (AREA)
- Computing Systems (AREA)
- General Engineering & Computer Science (AREA)
- Mathematical Physics (AREA)
- Apparatuses For Generation Of Mechanical Vibrations (AREA)
Abstract
本发明公开了一种基于模糊PI的频率跟踪搜索方法,包括:根据超声换能器两端电压电流的相位差计算当前时刻超声换能器的相位差变化率;根据当前时刻的所述相位差和相位差变化率模糊推理得到频率变化标识符;根据所述频率变化标识符判断变步长级别,并在所述变步长级别内指数化PI参数;根据所述PI参数确定超声换能器下一步步长的变化量,并计算超声换能器下一步的输出频率。本发明采用模糊PI全频域搜索和跟踪,使用者不需要设定初始频率,也不需要了解换能器的谐振频率,即能直接快速搜索到谐振频率。
Description
技术领域
本发明属于超声频率跟踪技术领域,特别是涉及一种基于模糊PI的频率跟踪搜索方法。
背景技术
超声波是指频率大于20kHz的声波,因为超出人的听觉上限,故名超声波。目前,由于超声波可以产生机械效应、空化效应、化学效应及热效应,被广泛应用于清洗、检测、加工、焊接等领域。
当超声换能器工作在谐振状态时,其输出功率最大,阻抗匹配技术和频率跟踪技术成为超声换能器工作于谐振状态的两个关键技术,超声技术正是利用这一原理产生超声波。由于超声换能器负载变化、换能器老化等诸多因素的影响,超声换能器的动态参数或静态参数会发生变化,超声换能器处于失谐状态。因此,为了使超声换能器长期工作于谐振状态,对超声换能器的频率跟踪成为超声研究领域一大热点问题。
胡武林、刘丽晨等人使用PI-DDS算法实现对频率的跟踪,这种跟踪方法速度快,性能好,但跟踪范围较小。李夏林、屈百达等人采用模糊控制算法实现对频率的跟踪,成贵等人采用极大似然估计方法对换能器的参数进行估计,进而实现对频率的跟踪这种跟踪算法由于采样的局限性进而导致参数不准确。夏旭峰、黄秋霖等人提出基于模糊PI的频率跟踪算法,但由于模糊语言清晰化的设计不足,导致了初始频率只能通过快速扫频的方式实现。彭呈祥等人提出的二分法、李长有等人提出的变步长扫频方法虽然解决了扫频问题,但是扫频速度过慢。
发明内容
本发明的目的在于提供一种基于模糊PI的频率跟踪搜索方法,克服现有技术在频率跟踪速度和频率跟踪范围内的不足,打破频率跟踪速度与频率跟踪范围之间的相互制约。
本发明的目的是通过以下技术方案来实现的:一种基于模糊PI的频率跟踪搜索方法,包括:
根据超声换能器两端电压电流的相位差计算当前时刻超声换能器的相位差变化率;
根据当前时刻的所述相位差和相位差变化率模糊推理得到频率变化标识符;
根据所述频率变化标识符判断变步长级别,并在所述变步长级别内指数化PI参数;
根据所述PI参数确定超声换能器下一步步长的变化量,并计算超声换能器下一步的输出频率。
优选的,所述超声换能器的相位差变化率的计算公式为:
优选的,所述频率变化标识符的获取方法为:根据预设的模糊规则,采用重心法推理得到。
优选的,所述预设的模糊规则为
优选的,所述指数化PI参数包括:计算PI调节的比例系数和PI调节的微分系数。
优选的,所述PI调节的比例系数的计算公式为:
所述PI调节的微分系数的计算公式为:
式中,Kp为PI调节的比例系数,Ki为PI调节的微分系数,flag为频率变化标识符,flag的变化范围为[0,7.5],k为变化方向标识符, Z为调节系数;当flag>4时,k为2;当flag<4时,k为1,当flag=4时,,停止频率搜索;flag在[0,2.5]U[5.5,7.5]内时,Z为5;flag在[2.5,5.5]内时,Z为4。
优选的,所述超声换能器下一步步长的变化量的计算公式为:
优选的,所述超声换能器下一步的输出频率的计算公式为:
本发明的有益效果是:
(1)本发明采用模糊PI全频域搜索和跟踪,使用者不需要设定初始频率,也不需要了解换能器的谐振频率,即能直接快速搜索到谐振频率;
(2)本发明通过相位差及相位差变化率识别换能器所处状态,进而调节步长,不会出现跟踪到并联谐振点的情况;
(3)本发明能够在18步内搜索到谐振频率附近,然后对谐振频率进行动态跟踪;从任意频率开始,其搜索效率总比普通PI变步长搜索高约40%;
(4)本发明方法的性能优于普通PI变步长跟踪,这是由于本发明的参数是处在动态调整中的,其兼具了调节步长和调节精度两方面的优点;而普通PI变步长跟踪算法在一定范围内是定参数的,因此跟踪步长变化率有限;
(5)本专利的指数化变化步长,可以使变化步长在1Hz至100kHz之间变化。
附图说明
图1为基于模糊PI的频率跟踪搜索方法的流程图;
图2为模糊PI控制的原理框图;
图3为超声换能器两端电压电流的相位差的模糊函数曲线图;
图4为超声换能器的相位差变化率的模糊函数曲线图。
具体实施方式
下面将结合实施例,对本发明的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域技术人员在没有付出创造性劳动的前提下所获得的所有其他实施例,都属于本发明保护的范围。
参阅图1-4,本发明提供一种基于模糊PI的频率跟踪搜索方法:
如图1所示,一种基于模糊PI的频率跟踪搜索方法,包括:
S1. 根据超声换能器两端电压电流的相位差计算当前时刻超声换能器的相位差变化率。
所述超声换能器的相位差变化率的计算公式为:
S2.根据当前时刻的所述相位差和相位差变化率模糊推理得到频率变化标识符。
如图2所示,本实施例中输入模糊化使用典型函数设计,根据超声换能器的相位变化特征对输入函数进行模糊化并在实际调试中对参数进行微调。超声换能器两端电压电流的相位差的变化范围为[-1.57,1.57],模糊论域为[NB NM NS ZO PS PB],其模糊函数曲线如图3所示。超声换能器的相位差变化率放大后,其变化范围为[-2,2],模糊论域为[NB NMNS ZO PS PB],其模糊函数曲线如图4所示。
所述频率变化标识符的获取方法为:根据预设的模糊规则,采用重心法推理得到。
所述预设的模糊规则为
模糊规则采用模糊条件推理if E is NM and EC is NS,Then U is NM的形式。为了避免搜索到并联谐振频率,制定模糊规则时,判断了并联谐振频率特征,即当模糊控制器追踪到并联谐振频率时,能够避开这一频率范围。
S3.根据所述频率变化标识符判断变步长级别,并在所述变步长级别内指数化PI参数。
所述指数化PI参数包括:计算PI调节的比例系数和PI调节的微分系数。
所述PI调节的比例系数的计算公式为:
所述PI调节的微分系数的计算公式为:
式中,Kp为PI调节的比例系数,Ki为PI调节的微分系数,flag为频率变化标识符,flag的变化范围为[0,7.5],k为变化方向标识符, Z为调节系数;当flag>4时,k为2;当flag<4时,k为1,当flag=4时,,停止频率搜索;flag在[0,2.5]U[5.5,7.5]内时,Z为5;flag在[2.5,5.5]内时,Z为4。
本实施中超声换能器下一步步长的变化量可实现更大化(df的范围可取为1~106)。当需要大步长搜索时,步长会取的足够大;当需要小步长跟踪时,步长可以取得足够小。因此,既可实现全频率快速搜索,又可以实现谐振频率邻域的跟踪。
S4.根据所述PI参数确定超声换能器下一步步长的变化量,并计算超声换能器下一步的输出频率。
所述超声换能器下一步步长的变化量的计算公式为:
所述超声换能器下一步的输出频率的计算公式为:
以上所述仅是本发明的优选实施方式,应当理解本发明并非局限于本文所披露的形式,不应看作是对其他实施例的排除,而可用于各种其他组合、修改和环境,并能够在本文所述构想范围内,通过上述教导或相关领域的技术或知识进行改动。而本领域人员所进行的改动和变化不脱离本发明的精神和范围,则都应在本发明所附权利要求的保护范围内。
Claims (8)
1.一种基于模糊PI的频率跟踪搜索方法,其特征在于,包括:
根据超声换能器两端电压电流的相位差计算当前时刻超声换能器的相位差变化率;
根据当前时刻的所述相位差和相位差变化率模糊推理得到频率变化标识符;
根据所述频率变化标识符判断变步长级别,并在所述变步长级别内指数化PI参数;
根据所述PI参数确定超声换能器下一步步长的变化量,并计算超声换能器下一步的输出频率。
3.根据权利要求1所述的一种基于模糊PI的频率跟踪搜索方法,其特征在于,所述频率变化标识符的获取方法为:根据预设的模糊规则,采用重心法推理得到。
5.根据权利要求1所述的一种基于模糊PI的频率跟踪搜索方法,其特征在于,所述指数化PI参数包括:计算PI调节的比例系数和PI调节的微分系数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010708046.4A CN111596544B (zh) | 2020-07-22 | 2020-07-22 | 一种基于模糊pd的频率跟踪搜索方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010708046.4A CN111596544B (zh) | 2020-07-22 | 2020-07-22 | 一种基于模糊pd的频率跟踪搜索方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111596544A true CN111596544A (zh) | 2020-08-28 |
CN111596544B CN111596544B (zh) | 2020-10-27 |
Family
ID=72192279
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010708046.4A Active CN111596544B (zh) | 2020-07-22 | 2020-07-22 | 一种基于模糊pd的频率跟踪搜索方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111596544B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113934137A (zh) * | 2021-10-20 | 2022-01-14 | 江南大学 | 一种超声波电源谐振频率跟踪方法及系统 |
CN116484155A (zh) * | 2022-10-20 | 2023-07-25 | 广东工业大学 | 一种压电换能器的全状态频率跟踪方法及系统 |
WO2023213049A1 (zh) * | 2022-05-05 | 2023-11-09 | 以诺康医疗科技(苏州)有限公司 | 超声波换能器频率跟踪方法、系统及超声设备 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6102865A (en) * | 1996-02-29 | 2000-08-15 | Acuson Corporation | Multiple ultrasound image registration system, method and transducer |
CN201304720Y (zh) * | 2008-12-05 | 2009-09-09 | 东莞市长江超声波机有限公司 | 一种超声波塑料焊接机中频率跟踪电路结构 |
US20130345566A1 (en) * | 2012-05-09 | 2013-12-26 | The Regents Of The University Of Michigan | Linear magnetic drive transducer for ultrasound imaging |
CN109290672A (zh) * | 2018-10-09 | 2019-02-01 | 湖南机电职业技术学院 | 一种基于超声波焊接的金属材料焊接方法及控制系统 |
CN109670139A (zh) * | 2019-01-29 | 2019-04-23 | 杭州国彪超声设备有限公司 | 一种超声波车削装置跟踪频率的确定方法及系统 |
CN110149056A (zh) * | 2019-05-27 | 2019-08-20 | 西安石油大学 | 基于模糊pi控制技术的超声波电源输出信号频率跟踪系统 |
CN110646518A (zh) * | 2019-09-26 | 2020-01-03 | 杭州电力设备制造有限公司 | 一种超声换能器的输出控制方法及相关设备 |
CN111420308A (zh) * | 2020-04-16 | 2020-07-17 | 重庆邮电大学 | 一种基于频率自动跟踪的超声理疗仪及频率自动跟踪方法 |
-
2020
- 2020-07-22 CN CN202010708046.4A patent/CN111596544B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6102865A (en) * | 1996-02-29 | 2000-08-15 | Acuson Corporation | Multiple ultrasound image registration system, method and transducer |
CN201304720Y (zh) * | 2008-12-05 | 2009-09-09 | 东莞市长江超声波机有限公司 | 一种超声波塑料焊接机中频率跟踪电路结构 |
US20130345566A1 (en) * | 2012-05-09 | 2013-12-26 | The Regents Of The University Of Michigan | Linear magnetic drive transducer for ultrasound imaging |
CN109290672A (zh) * | 2018-10-09 | 2019-02-01 | 湖南机电职业技术学院 | 一种基于超声波焊接的金属材料焊接方法及控制系统 |
CN109670139A (zh) * | 2019-01-29 | 2019-04-23 | 杭州国彪超声设备有限公司 | 一种超声波车削装置跟踪频率的确定方法及系统 |
CN110149056A (zh) * | 2019-05-27 | 2019-08-20 | 西安石油大学 | 基于模糊pi控制技术的超声波电源输出信号频率跟踪系统 |
CN110646518A (zh) * | 2019-09-26 | 2020-01-03 | 杭州电力设备制造有限公司 | 一种超声换能器的输出控制方法及相关设备 |
CN111420308A (zh) * | 2020-04-16 | 2020-07-17 | 重庆邮电大学 | 一种基于频率自动跟踪的超声理疗仪及频率自动跟踪方法 |
Non-Patent Citations (6)
Title |
---|
HONGJIE ZHANG: "A new automatic resonance frequency tracking method forpiezoelectric ultrasonic transducers used in thermosonic wirebonding", 《SENSORS AND ACTUATORS A: PHYSICAL》 * |
LANG-TE XU: "AN IMPROVED FREQUENCY TRACKING STRATEGY IN ULTRASONIC TRANSDUCER", 《2016 SYMPOSIUM ON PIEZOELECTRICITY, ACOUSTIC WAVES, AND DEVICE APPLICATIONS》 * |
夏旭峰: "基于模糊PI并联控制的超声传输系统频率跟踪", 《微特电机》 * |
廖晓辉: "基于模糊PI控制的超声波换能器谐振频率", 《河南科学》 * |
李夏林: "超声电源频率自动跟踪的模糊控制算法研究", 《应用声学》 * |
黄秋霖: "基于模糊PI的超声电源频率跟踪系统研究", 《现代计算机》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113934137A (zh) * | 2021-10-20 | 2022-01-14 | 江南大学 | 一种超声波电源谐振频率跟踪方法及系统 |
CN113934137B (zh) * | 2021-10-20 | 2022-07-15 | 江南大学 | 一种超声波电源谐振频率跟踪方法及系统 |
WO2023213049A1 (zh) * | 2022-05-05 | 2023-11-09 | 以诺康医疗科技(苏州)有限公司 | 超声波换能器频率跟踪方法、系统及超声设备 |
CN116484155A (zh) * | 2022-10-20 | 2023-07-25 | 广东工业大学 | 一种压电换能器的全状态频率跟踪方法及系统 |
CN116484155B (zh) * | 2022-10-20 | 2024-03-15 | 广东工业大学 | 一种压电换能器的全状态频率跟踪方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN111596544B (zh) | 2020-10-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111596544B (zh) | 一种基于模糊pd的频率跟踪搜索方法 | |
US10179346B2 (en) | Method of forming a transducer controller and circuit therefor | |
CN110572091A (zh) | 一种优化永磁同步电机无传感器控制方法 | |
Wang et al. | A high-tolerance matching method against load fluctuation for ultrasonic transducers | |
CN103500574B (zh) | 一种智能电子喇叭及其实现方法 | |
Shi et al. | Chaotic operation and chaos control of travelling wave ultrasonic motor | |
CN109546991B (zh) | 应用于宽频超声波电源的频率搜索匹配跟踪方法及系统 | |
CN112865341B (zh) | Lc串联拓扑无线充电系统的异物检测方法 | |
JP3949350B2 (ja) | 連系装置 | |
CN117240244A (zh) | 一种用于超声电源的t型动态匹配网络及其匹配方法 | |
CN103064447B (zh) | 基于窄域论的激光器pid温控参数近似三维整定方法 | |
CN110995110B (zh) | 一种单相永磁直线压缩机抗扰动控制系统及方法 | |
CN110446290B (zh) | 一种用于感应加热的准谐振控制方法 | |
CN115173834B (zh) | 基于陷波参数量化分析的多级串联陷波器优化设计方法 | |
CN103713688B (zh) | 一种自适应变步长mppt控制方法 | |
CN112671115A (zh) | 基于模糊自适应pi的wpt系统及频率跟踪方法 | |
CN116111861A (zh) | 一种数据中心单相Boost型PFC低频振荡快速抑制方法 | |
CN108155836A (zh) | 基于全局在线启发式动态规划永磁同步电机矢量控制方法 | |
KR102180166B1 (ko) | 트랜스듀서 컨트롤러 및 그를 위한 회로의 형성 방법 | |
CN117340421A (zh) | 一种用于超声焊接的谐振频率追踪方法 | |
CN106253744B (zh) | 超声波电机力矩滞回控制方法 | |
Soltani et al. | Adaptive fuzzy sliding-mode control of speed sensorless universal field oriented induction motor drive with on-line stator resistance tuning | |
CN105184077B (zh) | 过近距离下共振电能传输系统效率寻优粒子群指数方法 | |
CN116667275B (zh) | 一种智能矩阵式医疗清洗系统设备过流保护系统 | |
CN115826669B (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 |