CN103616848B - 改善速度前馈的自适应滤波器的滤波方法和系统 - Google Patents

改善速度前馈的自适应滤波器的滤波方法和系统 Download PDF

Info

Publication number
CN103616848B
CN103616848B CN201310648244.6A CN201310648244A CN103616848B CN 103616848 B CN103616848 B CN 103616848B CN 201310648244 A CN201310648244 A CN 201310648244A CN 103616848 B CN103616848 B CN 103616848B
Authority
CN
China
Prior art keywords
sef
filter
speed
adapting filter
velocity
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.)
Expired - Fee Related
Application number
CN201310648244.6A
Other languages
English (en)
Other versions
CN103616848A (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.)
GOOGOL TECHNOLOGY (SHENZHEN) Ltd
Zhejiang Normal University CJNU
Original Assignee
GOOGOL TECHNOLOGY (SHENZHEN) Ltd
Zhejiang Normal University CJNU
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 GOOGOL TECHNOLOGY (SHENZHEN) Ltd, Zhejiang Normal University CJNU filed Critical GOOGOL TECHNOLOGY (SHENZHEN) Ltd
Priority to CN201310648244.6A priority Critical patent/CN103616848B/zh
Publication of CN103616848A publication Critical patent/CN103616848A/zh
Application granted granted Critical
Publication of CN103616848B publication Critical patent/CN103616848B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Feedback Control In General (AREA)

Abstract

本发明涉及一种改善速度前馈的自适应滤波器的滤波方法和系统。所述方法包括以下步骤:获取速度变化量;根据所述速度变化量计算自适应滤波器的滤波系数,并根据所述自适应滤波器的滤波系数对速度进行滤波处理。上述改善速度前馈的自适应滤波器的滤波方法和系统,根据速度变化量获取滤波系数,并根据滤波系数对速度进行滤波处理,改善速度的跳变,平滑速度,从而使得速度前馈的输出更加平滑,减少机床的振动。

Description

改善速度前馈的自适应滤波器的滤波方法和系统
技术领域
本发明涉及自动控制领域,特别是涉及一种改善速度前馈的自适应滤波器的滤波方法和系统。
背景技术
在工业生产中,通过多轴连续轨迹数控设备进行加工得到广泛的应用。该多轴连续轨迹数控设备由于运动轴的加减速或轨迹矢量方向变化导致单轴速度跳变,从而使速度前馈的输出跳变,容易引起机床的振动。
发明内容
基于此,有必要针对速度跳变引起机床振动的问题,提供一种改善速度前馈的自适应滤波器的滤波方法,能使速度前馈输出平滑,减少机床的振动。
此外,还有必要提供一种改善速度前馈的自适应滤波器的滤波方法,能使速度前馈输出平滑,减少机床的振动。
一种改善速度前馈的自适应滤波器的滤波方法,包括以下步骤:
获取速度变化量;
根据所述速度变化量计算自适应滤波器的滤波系数,并根据所述自适应滤波器的滤波系数对速度进行滤波处理。
在其中一个实施例中,所述根据所述速度变化量计算自适应滤波器的滤波系数,并根据所述自适应滤波器的滤波系数对速度进行滤波处理的步骤包括:
当所述速度变化量为零或预设加速度,则自适应滤波器的滤波系数为1,对速度不进行滤波;
当所述速度变化量不为零或预设加速度,则自适应滤波器的滤波系数小于1,对速度进行滤波处理。
在其中一个实施例中,所述根据所述速度变化率计算自适应滤波器的滤波系数,并根据所述自适应滤波器的滤波系数对速度进行滤波处理的计算公式为:
dn=|xn-yn-1|
f = 1 1 + k * | 4 * d n ( d n - a ) a 2 |
yn=f*xn+(1-f)*yn-1
其中,a为预设加速度,xn为n时刻的原始速度,yn为n时刻经过滤波后的速度,k为可调系数,dn为n时刻的速度变化量,f为自适应滤波器的滤波系数。
一种改善速度前馈的自适应滤波器的滤波系统,包括:
获取模块,用于获取速度变化量;
处理模块,用于根据所述速度变化量计算自适应滤波器的滤波系数,并根据所述自适应滤波器的滤波系数对速度进行滤波处理。
在其中一个实施例中,所述处理模块根据所述速度变化量计算自适应滤波器的滤波系数,并根据所述自适应滤波器的滤波系数对速度进行滤波处理,包括:
当所述速度变化量为零或预设加速度,则自适应滤波器的滤波系数为1,对速度不进行滤波;
当所述速度变化量不为零或预设加速度,则自适应滤波器的滤波系数小于1,对速度进行滤波处理。
在其中一个实施例中,所述根据所述速度变化率计算自适应滤波器的滤波系数,并根据所述自适应滤波器的滤波系数对速度进行滤波处理的计算公式为:
dn=|xn-yn-1|
f = 1 1 + k * | 4 * d n ( d n - a ) a 2 |
yn=f*xn+(1-f)*yn-1
其中,a为预设加速度,xn为n时刻的原始速度,yn为n时刻经过滤波后的速度,k为可调系数,dn为n时刻的速度变化量,f为自适应滤波器的滤波系数。
上述改善速度前馈的自适应滤波器的滤波方法和系统,根据速度变化量获取滤波系数,并根据滤波系数对速度进行滤波处理,改善速度的跳变,平滑速度,从而使得速度前馈的输出更加平滑,减少机床的振动。
附图说明
图1为一个实施例中改善速度前馈的自适应滤波器的滤波方法的流程图;
图2为当k为0.5,a为1时,随dn变化的f值示意图;
图3为一个实施例中改善速度前馈的自适应滤波器的滤波系统的结构框图;
图4为未经过自适应滤波器滤波的X轴时间-速度的示意图;
图5为经过自适应滤波器滤波的X轴时间-速度的示意图;
图6为未经过自适应滤波器滤波的X轴时间-加速度的示意图;
图7为经过自适应滤波器滤波的X轴时间-加速度的示意图;
图8为未经过自适应滤波器滤波的Y轴时间-速度的示意图;
图9为经过自适应滤波器滤波的Y轴时间-速度的示意图;
图10为未经过自适应滤波器滤波的Y轴时间-加速度的示意图;
图11为经过自适应滤波器滤波的Y轴时间-加速度的示意图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
图1为一个实施例中改善速度前馈的自适应滤波器的滤波方法的流程图。该改善速度前馈的自适应滤波器的滤波方法,包括以下步骤:
步骤102,获取速度变化量。
步骤104,根据该速度变化量计算自适应滤波器的滤波系数,并根据该自适应滤波器的滤波系数对速度进行滤波处理。
本实施例中,该根据该速度变化量计算自适应滤波器的滤波系数,并根据该自适应滤波器的滤波系数对速度进行滤波处理的步骤包括:
(1)当该速度变化量为零或预设加速度,则自适应滤波器的滤波系数为1,对速度不进行滤波。
具体的,速度变化量为零,即匀速运动。预设加速度可在实际加工中由用户根据设定,主要根据加工机床的性能和具体加工材料来确定。
(2)当该速度变化量不为零或预设加速度,则自适应滤波器的滤波系数小于1,对速度进行滤波处理。
本实施例中,该根据该速度变化率计算自适应滤波器的滤波系数,并根据该自适应滤波器的滤波系数对速度进行滤波处理的计算公式为:
dn=|xn-yn-1|
f = 1 1 + k * | 4 * d n ( d n - a ) a 2 |
yn=f*xn+(1-f)*yn-1
其中,a为预设加速度,xn为n时刻的原始速度,yn为n时刻经过滤波后的速度,k为可调系数,dn为n时刻的速度变化量,f为自适应滤波器的滤波系数,yn-1为(n-1)时刻经过滤波后的速度。
图2为当k为0.5,a为1时,随dn变化的f值。由图2可知,当速度变化量为0或a时,滤波系数为1,即滤波输出为原始输入;当速度变化量小于a时,滤波系数当速度变化量大于a时,随着速度跳变量的增大,滤波系数以二阶函数形式不断减小,即滤波强度不断变大,使速度更加平滑,即速度变化量偏离0和a越远,滤波效果越明显。
上述改善速度前馈的自适应滤波器的滤波方法,根据速度变化量获取滤波系数,并根据滤波系数对速度进行滤波处理,改善速度的跳变,平滑速度,从而使得速度前馈的输出更加平滑,减少机床的振动。此外,通过预设加速度,并设计对应的滤波系数算法,更易调整滤波系数,控制速度平滑输出。
图3为一个实施例中改善速度前馈的自适应滤波器的滤波系统的结构框图。该改善速度前馈的自适应滤波器的滤波系统,包括获取模块320和处理模块340。其中:
获取模块320用于获取速度变化量。
处理模块340用于根据该速度变化量计算自适应滤波器的滤波系数,并根据该自适应滤波器的滤波系数对速度进行滤波处理。
该处理模块340根据该速度变化量计算自适应滤波器的滤波系数,并根据该自适应滤波器的滤波系数对速度进行滤波处理,包括:
(1)当该速度变化量为零或预设加速度,则自适应滤波器的滤波系数为1,对速度不进行滤波。
(2)当该速度变化量不为零或预设加速度,则自适应滤波器的滤波系数小于1,对速度进行滤波处理。
本实施例中,该根据该速度变化率计算自适应滤波器的滤波系数,并根据该自适应滤波器的滤波系数对速度进行滤波处理的计算公式为:
dn=|xn-yn-1|
f = 1 1 + k * | 4 * d n ( d n - a ) a 2 |
yn=f*xn+(1-f)*yn-1
其中,a为预设加速度,xn为n时刻的原始速度,yn为n时刻经过滤波后的速度,k为可调系数,dn为n时刻的速度变化量,f为自适应滤波器的滤波系数,yn-1为(n-1)时刻经过滤波后的速度。
上述改善速度前馈的自适应滤波器的滤波方法和系统应用于加工半径为10毫米圆形,取a=500mm/s2(毫米/秒2),k为0.01为例对速度和加速度值进行修正,得到对应的修正结果。图4为未经过自适应滤波器滤波的X轴时间-速度的示意图;图5为经过自适应滤波器滤波的X轴时间-速度的示意图;图6为未经过自适应滤波器滤波的X轴时间-加速度的示意图;图7为经过自适应滤波器滤波的X轴时间-加速度的示意图;图8为未经过自适应滤波器滤波的Y轴时间-速度的示意图;图9为经过自适应滤波器滤波的Y轴时间-速度的示意图;图10为未经过自适应滤波器滤波的Y轴时间-加速度的示意图;图11为经过自适应滤波器滤波的Y轴时间-加速度的示意图。图4、图5、图8和图9的横坐标为插补时间,单位为插补周期250微秒,纵坐标为速度值,单位为毫米/插补周期;图6、图7、图10和图11的横坐标为插补时间,单位为250微秒,纵坐标为加速度值,单位为毫米/插补周期2。从图6和图7可知,单轴速度的跳变即加速度明显减小。
本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程,是可以通过计算机程序来指令相关的硬件来完成,所述的程序可存储于一计算机可读取存储介质中,该程序在执行时,可包括如上述各方法的实施例的流程。其中,所述的存储介质可为磁碟、光盘、只读存储记忆体(Read-Only Memory,ROM)或随机存储记忆体(Random Access Memory,RAM)等。
以上所述实施例仅表达了本发明的几种实施方式,其描述较为具体和详细,但并不能因此而理解为对本发明专利范围的限制。应当指出的是,对于本领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干变形和改进,这些都属于本发明的保护范围。因此,本发明专利的保护范围应以所附权利要求为准。

Claims (4)

1.一种改善速度前馈的自适应滤波器的滤波方法,包括以下步骤:
获取速度变化量;
根据所述速度变化量计算自适应滤波器的滤波系数,并根据所述自适应滤波器的滤波系数对速度进行滤波处理;
所述滤波系数根据预设加速度设计对应的滤波系数算法计算得出;
所述根据所述速度变化量计算自适应滤波器的滤波系数,并根据所述自适应滤波器的滤波系数对速度进行滤波处理的计算公式为:
dn=|xn-yn-1|
f = 1 1 + k * | 4 * d n ( d n - a ) a 2 |
yn=f*xn+(1-f)*yn-1
其中,a为预设加速度,xn为n时刻的原始速度,yn为n时刻经过滤波后的速度,k为可调系数,dn为n时刻的速度变化量,f为自适应滤波器的滤波系数。
2.根据权利要求1所述的改善速度前馈的自适应滤波器的滤波方法,其特征在于,所述根据所述速度变化量计算自适应滤波器的滤波系数,并根据所述自适应滤波器的滤波系数对速度进行滤波处理的步骤包括:
当所述速度变化量为零或预设加速度,则自适应滤波器的滤波系数为1,对速度不进行滤波;
当所述速度变化量不为零或预设加速度,则自适应滤波器的滤波系数小于1,对速度进行滤波处理。
3.一种改善速度前馈的自适应滤波器的滤波系统,其特征在于,包括:
获取模块,用于获取速度变化量;
处理模块,用于根据所述速度变化量计算自适应滤波器的滤波系数,并根据所述自适应滤波器的滤波系数对速度进行滤波处理;
所述滤波系数根据预设加速度设计对应的滤波系数算法计算得出;
所述根据所述速度变化量计算自适应滤波器的滤波系数,并根据所述自适应滤波器的滤波系数对速度进行滤波处理的计算公式为:
dn=|xn-yn-1|
f = 1 1 + k * | 4 * d n ( d n - a ) a 2 |
yn=f*xn+(1-f)*yn-1
其中,a为预设加速度,xn为n时刻的原始速度,yn为n时刻经过滤波后的速度,k为可调系数,dn为n时刻的速度变化量,f为自适应滤波器的滤波系数。
4.根据权利要求3所述的改善速度前馈的自适应滤波器的滤波系统,其特征在于,所述处理模块根据所述速度变化量计算自适应滤波器的滤波系数,并根据所述自适应滤波器的滤波系数对速度进行滤波处理,包括:
当所述速度变化量为零或预设加速度,则自适应滤波器的滤波系数为1,对速度不进行滤波;
当所述速度变化量不为零或预设加速度,则自适应滤波器的滤波系数小于1,对速度进行滤波处理。
CN201310648244.6A 2013-12-04 2013-12-04 改善速度前馈的自适应滤波器的滤波方法和系统 Expired - Fee Related CN103616848B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310648244.6A CN103616848B (zh) 2013-12-04 2013-12-04 改善速度前馈的自适应滤波器的滤波方法和系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310648244.6A CN103616848B (zh) 2013-12-04 2013-12-04 改善速度前馈的自适应滤波器的滤波方法和系统

Publications (2)

Publication Number Publication Date
CN103616848A CN103616848A (zh) 2014-03-05
CN103616848B true CN103616848B (zh) 2016-08-17

Family

ID=50167552

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310648244.6A Expired - Fee Related CN103616848B (zh) 2013-12-04 2013-12-04 改善速度前馈的自适应滤波器的滤波方法和系统

Country Status (1)

Country Link
CN (1) CN103616848B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106597402A (zh) * 2016-11-23 2017-04-26 天津津航计算技术研究所 一种基于统计规律的自适应动目标显示方法

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP7035733B2 (ja) * 2018-03-30 2022-03-15 ブラザー工業株式会社 工作機械及び平滑化方法
CN109254563B (zh) * 2018-10-22 2021-04-06 大族激光科技产业集团股份有限公司 一种数控速度滤波方法及其滤波系统
CN109240087B (zh) * 2018-10-23 2022-03-01 固高科技股份有限公司 实时改变指令规划频率抑制振动的方法和系统

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0376250A1 (fr) * 1988-12-30 1990-07-04 Alcatel Telspace Dispositif d'égalisation auto-adaptative pour installation de démodulation différentiellement cohérente
CN101246353A (zh) * 2008-03-17 2008-08-20 湖南宇晶机器实业有限公司 Ic材料线切割机床多电机系统速度同步自适应逆控制方法
CN101702091A (zh) * 2009-10-19 2010-05-05 大连海事大学 一种基于自适应滤波器的电液伺服系统随机振动控制方法
CN102175889A (zh) * 2011-01-24 2011-09-07 长春工业大学 一种伺服转台角加速度自适应测量方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2006059273A1 (en) * 2004-12-03 2006-06-08 Arima Devices Corporation Self-adaptive filter circuit
CN101470434B (zh) * 2007-12-28 2010-08-25 中国科学院沈阳计算技术研究所有限公司 基于滤波技术的速度前瞻控制方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0376250A1 (fr) * 1988-12-30 1990-07-04 Alcatel Telspace Dispositif d'égalisation auto-adaptative pour installation de démodulation différentiellement cohérente
CN101246353A (zh) * 2008-03-17 2008-08-20 湖南宇晶机器实业有限公司 Ic材料线切割机床多电机系统速度同步自适应逆控制方法
CN101702091A (zh) * 2009-10-19 2010-05-05 大连海事大学 一种基于自适应滤波器的电液伺服系统随机振动控制方法
CN102175889A (zh) * 2011-01-24 2011-09-07 长春工业大学 一种伺服转台角加速度自适应测量方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
基于滤波技术的数控系统加减速研究;于东 等;《中国机械工程》;20080415;第19卷(第7期);第804-807页 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106597402A (zh) * 2016-11-23 2017-04-26 天津津航计算技术研究所 一种基于统计规律的自适应动目标显示方法
CN106597402B (zh) * 2016-11-23 2019-03-19 天津津航计算技术研究所 一种基于统计规律的自适应动目标显示方法

Also Published As

Publication number Publication date
CN103616848A (zh) 2014-03-05

Similar Documents

Publication Publication Date Title
CN103616848B (zh) 改善速度前馈的自适应滤波器的滤波方法和系统
Liu et al. Adaptive interpolation scheme for NURBS curves with the integration of machining dynamics
Xi et al. Improving CNC contouring accuracy by robust digital integral sliding mode control
CN111267857B (zh) 一种自动驾驶车辆的变道轨迹生成方法及装置
US6823235B2 (en) Controller for machining gears
CN107340732B (zh) 伺服控制装置以及伺服控制方法
CN111015738A (zh) 一种工业机器人振动抑制方法
CN105811817A (zh) 电动机控制装置
CN109643102B (zh) 指令值生成装置
CN104204977A (zh) 轨迹控制装置
CN104977895A (zh) 伺服控制装置
CN102707671A (zh) 应用于工具机的加工路径最佳化方法
CN114779721B (zh) 基于位置点进行速度滤波的方法、设备和存储介质
CN102360198A (zh) 数控系统中操作设备的速度规划方法、装置及数控机床
JP2935713B2 (ja) 数値制御装置
CN105425730A (zh) 一种NURBS曲线Taylor迭代的插补算法
CN103777573B (zh) 减小轮廓加工误差的方法和系统
CN107544247A (zh) 一种抑制机械运动结构振动的方法及系统
CN103454979A (zh) 一种封装成PLCOpen指令的变速曲线圆弧快速插补方法
Chang et al. Design and implementation of look-ahead linear jerk filter for a computerized numerical controlled machine
Tsai et al. Integration of input shaping technique with interpolation for vibration suppression of servo-feed drive system
CN108762064B (zh) 一种伺服驱动器的速度平滑方法
JP2010182157A (ja) 補正装置及び数値制御装置
CN104106016B (zh) 指令生成装置以及指令生成方法
CN111813053A (zh) 数控系统中针对速度滤波导致转角处速度超限进行控制的方法、装置及其存储介质

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C41 Transfer of patent application or patent right or utility model
CB03 Change of inventor or designer information

Inventor after: Wu Hong

Inventor after: Zhu Xinzhong

Inventor after: Lv Shu

Inventor after: Wang Rui

Inventor after: Zhang Tinghua

Inventor after: Xu Huiying

Inventor after: Zhao Jianmin

Inventor before: Wu Hong

Inventor before: Lv Shu

Inventor before: Wang Rui

Inventor before: Zhang Tinghua

COR Change of bibliographic data
TA01 Transfer of patent application right

Effective date of registration: 20160516

Address after: Room two, building W211 block 518057 West Guangdong province Shenzhen city Nanshan District high tech Zone ier

Applicant after: GOOGOL TECHNOLOGY (SHENZHEN) Ltd.

Applicant after: ZHEJIANG NORMAL University

Address before: Room two, building W211 block 518057 West Guangdong province Shenzhen city Nanshan District high tech Zone ier

Applicant before: Googol Technology (Shenzhen) Ltd.

C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20160817

Termination date: 20211204

CF01 Termination of patent right due to non-payment of annual fee