CN112054778A - 一种截止频率自调整低通滤波器及设计方法 - Google Patents

一种截止频率自调整低通滤波器及设计方法 Download PDF

Info

Publication number
CN112054778A
CN112054778A CN202010897502.4A CN202010897502A CN112054778A CN 112054778 A CN112054778 A CN 112054778A CN 202010897502 A CN202010897502 A CN 202010897502A CN 112054778 A CN112054778 A CN 112054778A
Authority
CN
China
Prior art keywords
frequency
cut
pass filter
signal
integrator
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
CN202010897502.4A
Other languages
English (en)
Other versions
CN112054778B (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.)
Rocket Force University of Engineering of PLA
Original Assignee
Rocket Force University of Engineering of PLA
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 Rocket Force University of Engineering of PLA filed Critical Rocket Force University of Engineering of PLA
Priority to CN202010897502.4A priority Critical patent/CN112054778B/zh
Publication of CN112054778A publication Critical patent/CN112054778A/zh
Application granted granted Critical
Publication of CN112054778B publication Critical patent/CN112054778B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H03ELECTRONIC CIRCUITRY
    • H03HIMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
    • H03H7/00Multiple-port networks comprising only passive electrical elements as network components
    • H03H7/01Frequency selective two-port networks
    • HELECTRICITY
    • H03ELECTRONIC CIRCUITRY
    • H03HIMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
    • H03H17/00Networks using digital techniques
    • H03H17/02Frequency selective networks

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Mathematical Physics (AREA)
  • Control Of Electric Motors In General (AREA)

Abstract

本发明公开一种截止频率自调整低通滤波器的设计方法,该方法包括以下步骤:步骤S1:从频域的角度分析量化噪声;步骤S2:根据步骤S1分析的结果设计出截止频率自调整的低通滤波器。本发明设计的截止频率自调整低通滤波器,可以根据伺服电机转速实时地调整截止频率,有效地抑制量化噪声对闭环系统的影响,位移信号和速度信号估计精度更高,控制信号振动减弱,系统能耗降低。

Description

一种截止频率自调整低通滤波器及设计方法
技术领域
本发明属于伺服电机控制技术领域,具体涉及一种截止频率自调整低通滤波器及设计方法。
背景技术
对于伺服电机控制系统,高精度的位置传感器,是实现高精度位置控制的基础,对于旋转位置的检测与控制,常使用光电编码器。光电编码器根据能否检测绝对位置,可分为绝对式和增量式,无论哪种形式,测量的基础都是编码器码盘上的刻线,刻线的间隔造成了位置测量中的分辨率限制,这也从原理上决定了光电编码器的输出信号对于真实的位置信息具有“量化效应”,检测得到的信号与真实位置之间具有“量化误差”,也被称为量化噪声。
将带有量化噪声的编码器信号用于反馈控制,会降低位移和速度的控制精度,甚至造成极限环振荡。特别是在伺服电机转速较低的情况下,编码器测角的量化噪声对控制系统系统的影响非常明显,影响位移信号和速度信号的估计精度,同时会导致控制信号的振动,增加系统的能耗。
目前的研究主要是将量化噪声近似为均匀分布的白噪声,建立量化噪声的统计模型,通过Kalman滤波来抑制量化噪声,这种近似将量化噪声视为随机信号,虽然便于处理但没有反映出量化噪声的本质。同时有研究认为,量化效应是一种确定性和无记忆的非线性环节,并采用基于观测器理论的状态估计器,用于从量化的信号中提取有用的信息,虽然便于实现但是适应性较差,当伺服电机转速变化较大时,闭环系统性能下降严重。
发明内容
本发明针对上述存在的技术问题,提出了一种截止频率自调整低通滤波器,可以根据伺服电机转速实时地调整截止频率,有效地抑制量化噪声对闭环系统的影响,位移信号和速度信号估计精度更高,控制信号振动减弱,系统能耗降低。
本发明采用的技术方案:
一种截止频率自调整低通滤波器的设计方法,
该方法包括以下步骤:
步骤S1:从频域的角度分析量化噪声;
(1)编码器量化噪声是关于位移输入的奇函数,随位移输入的变化呈锯齿状周期变化,周期为Δq,Δq表示量化阶,其为编码器能够检测到的最小角度位置变化量;
(2)则令x表示位移输入,f(x)表示关于位移输入的量化误差,利用傅里叶级数对f(x)展开:
Figure BDA0002658932790000021
其中,
Figure BDA0002658932790000022
Figure BDA0002658932790000023
Figure BDA0002658932790000024
计算可知:
Figure BDA0002658932790000031
则公式1可写为
Figure BDA0002658932790000032
以时刻t作为初始时刻,伺服电机转速为ω,则有x=ωt,公式3可写为
Figure BDA0002658932790000033
根据公式4可知定义fq=ω/Δq,则量化噪声f(x)是由频率为nfq幅值为
Figure BDA0002658932790000034
的谐波分量构成;
步骤S2:根据步骤S1分析的结果设计出截止频率自调整的低通滤波器。
优选的,上述步骤S2中,截止频率自调整的低通滤波器设计过程如下:
(1)伺服电机转速为ω时,对应的角频率为
Figure BDA0002658932790000035
可推得:
Figure BDA0002658932790000036
根据公式5,将低通滤波器的截止频率设置在fe和fq之间,可以顺利地滤除掉量化噪声,保留有用的电机转速信号;
(2)上述步骤(1)给出了滤除量化噪声的截止频率的设置范围,在步骤(1)的基础上,截止频率自调整低通滤波器可表达
Figure BDA0002658932790000037
式中,ξ=0.707,ωn为其3dB截止频率,则
Figure BDA0002658932790000041
设xe表示有量化噪声的实际测量输出,xf表示滤波器输出,那么式(6)可写为
Figure BDA0002658932790000042
其中,xf表示估计的位移信号,
Figure BDA0002658932790000043
表示估计的速度信号。
一种截止频率自调整低通滤波器,包括第一放大器、第一积分器、第二放大器和第二积分器;所述第一放大器的输出端和第二放大器的输出端均与第一积分器的输入端相连,所述第一积分器的输出端分别与第二放大器和第二积分器的输入端相连,所述第二积分器的输出端的一路与第一放大器的输入端相连;
使用时,含有量化噪声的实际位移信号xe与反馈的滤波器输出位移信号xf做差,得到的误差信号通过第一放大器进行放大,然后将放大后的误差信号与反馈的通过第二放大器放大的滤波器的估计速度信号做差,得到误差信号通过第一积分器积分,获得滤波器估计的速度信号
Figure BDA0002658932790000044
,然后通过第二积分器积分,得到滤波器输出位移信号xf
与现有技术相比,本发明的有益效果是:本发明设计的截止频率自调整低通滤波器可以有效地抑制编码器量化噪声对伺服电机位置信号和速度信息估计精度的影响,有效地减小估计的位置信号和速度信号以及解算的控制信号的振动,减少闭环系统的能耗。
附图说明
为了更清楚地说明本发明实施例的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,应当理解,以下附图仅示出了本发明的某些实施例,因此不应被看作是对范围的限定,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他相关的附图。
图1为编码器量化输出和量化噪声图;(a)编码器量化输出;(b)编码器量化噪声;
图2为截止频率自调整低通滤波器结构图;
图3为本发明设计的截止频率自校正低通滤波器的仿真结果;
图4为固定截止频率低通滤波器的仿真结果。
其中,1-第一放大器;2-第一积分器;3-第二放大器;4-第二积分器。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。通常在此处附图中描述和示出的本发明实施例的组件可以以各种不同的配置来布置和设计。
因此,以下对在附图中提供的本发明的实施例的详细描述并非旨在限制要求保护的本发明的范围,而是仅仅表示本发明的选定实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明具体提供了一种截止频率自调整低通滤波器的设计方法,该方法包括以下步骤:
步骤S1:从频域的角度分析量化噪声;
(1)编码器量化输出和量化噪声分别如图(a)和(b)所示,编码器量化噪声是关于位移输入的奇函数,随位移输入的变化呈锯齿状周期变化,周期为Δq,Δq表示量化阶,其为编码器能够检测到的最小角度位置变化量;
(2)则令x表示位移输入,f(x)表示关于位移输入的量化误差,考虑函数f(x)的周期函数特性,利用傅里叶级数对f(x)展开:
Figure BDA0002658932790000061
其中,
Figure BDA0002658932790000062
Figure BDA0002658932790000063
Figure BDA0002658932790000064
计算可知:
Figure BDA0002658932790000065
则公式1可写为
Figure BDA0002658932790000066
以时刻t作为初始时刻,伺服电机转速为ω,则有x=ωt,公式3可写为
Figure BDA0002658932790000067
根据公式4可知定义fq=ω/Δq,则量化噪声f(x)是由频率为nfq幅值为
Figure BDA0002658932790000071
的谐波分量构成;
步骤S2:根据步骤S1分析的结果设计出截止频率自调整的低通滤波器;
(1)伺服电机转速为ω时,对应的角频率为
Figure BDA0002658932790000072
可推得:
Figure BDA0002658932790000073
根据公式5,将低通滤波器的截止频率设置在fe和fq之间,可以顺利地滤除掉量化噪声,保留有用的电机转速信号;当伺服电机转速较高时,截止频率设计地适当大,当伺服电机转速较低时,截止频率设计地适当小;相比于截止频率固定的低通滤波器,截止频率自调整低通滤波器在伺服电机高速时,可以提高位移信号和速度信号的精度并减小相位滞后,在伺服电机低速时,有效滤除掉量化噪声,进而提高位移信号和速度信号的估计精度,降低控制输入的振动;
(2)上述步骤(1)给出了滤除量化噪声的截止频率的设置范围,在步骤(1)的基础上,截止频率自调整低通滤波器可表达
Figure BDA0002658932790000074
式中,ξ=0.707,ωn为其3dB截止频率,则
Figure BDA0002658932790000075
设xe表示有量化噪声的实际测量输出,xf表示滤波器输出,那么式(6)可写为
Figure BDA0002658932790000076
其中,xf表示估计的位移信号,
Figure BDA0002658932790000077
表示估计的速度信号。
本发明还提出了
一种截止频率自调整低通滤波器,包括第一放大器1、第一积分器2、第二放大器3和第二积分器4;所述第一放大器1的输出端和第二放大器3的输出端均与第一积分器2的输入端相连,所述第一积分器1的输出端分别与第二放大器3和第二积分器4的输入端相连,所述第二积分器的输出端的一路与第一放大器的输入端相连;
使用时,含有量化噪声的实际位移信号xe与反馈的滤波器输出位移信号xf做差,得到的误差信号通过第一放大器1进行放大,然后将放大后的误差信号与反馈的通过第二放大器3放大的滤波器的估计速度信号做差,得到误差信号通过第一积分器2积分,获得滤波器估计的速度信号
Figure BDA0002658932790000081
然后通过第二积分器4积分,得到滤波器输出位移信号xf
仿真试验
通过数值仿真验证本发明提供的截止频率自调整低通滤波器的可实施性。考虑伺服电机
Figure BDA0002658932790000082
其中,m表示负载质量,x表示电机转轴的位移,
Figure BDA0002658932790000083
表示受到的干扰,K表示控制增益系数,y(t)表示传感器输出,f(x)表示光电编码器的量化噪声。
设b=K/m,xr(t)表示期望信号,为了实现伺服电机位移x对xr(t)的跟踪,可设计PID控制器
Figure BDA0002658932790000084
其中,ef(t)=xr(t)-xf(t),
Figure BDA0002658932790000085
xf(t)和
Figure BDA0002658932790000086
表示设计的截止频率自校正的低通滤波器的输出且满足
Figure BDA0002658932790000087
数值仿真中,式(9)表示的伺服电机的参数b=1,干扰d设计为
Figure BDA0002658932790000091
量化阶Δq=2π/2500。式(10)设计的PID控制器参数kp=100,kd=14.14,ki=150,α=0.1.为了防止滤波器截止频率过大过小影响系统的稳定性,对截止频率设置饱和界
Figure BDA0002658932790000092
为了直观体现设计的低通滤波器的优势,与利用固定截止频率的低通滤波器对比,截止频率固定为100。当伺服转速较高时,
Figure BDA0002658932790000093
设计的低通滤波器退化为截止频率固定的低通滤波器,因此,当伺服电机转速较低时,设计的低通滤波器的优势体现得更明显。设伺服电机的期望转速为1r/min,此时
Figure BDA0002658932790000094
对应的仿真结果如图3-4所示。
根据图3-4的仿真结果可知,相比于固定截止频率低通滤波器,利用设计的低通滤波器的控制器在实现精确位移跟踪和调速的同时,可有效地滤除量化噪声,速度观测更加准确,控制输入波动幅值更小,闭环系统更加节能。
以上所述,仅是本发明的较佳实施例,并非对本发明作任何限制。凡是根据发明技术实质对以上实施例所作的任何简单修改、变更以及等效变化,均仍属于本发明技术方案的保护范围内。

Claims (3)

1.一种截止频率自调整低通滤波器的设计方法,其特征在于,该方法包括以下步骤:
步骤S1:从频域的角度分析量化噪声;
(1)编码器量化噪声是关于位移输入的奇函数,随位移输入的变化呈锯齿状周期变化,周期为Δq,Δq表示量化阶,其为编码器能够检测到的最小角度位置变化量;
(2)则令x表示位移输入,f(x)表示关于位移输入的量化误差,利用傅里叶级数对f(x)展开:
Figure FDA0002658932780000011
其中,
Figure FDA0002658932780000012
Figure FDA0002658932780000013
Figure FDA0002658932780000014
计算可知:
Figure FDA0002658932780000015
则公式1可写为
Figure FDA0002658932780000021
以时刻t作为初始时刻,伺服电机转速为ω,则有x=ωt,公式3可写为
Figure FDA0002658932780000022
根据公式4可知定义fq=ω/Δq,则量化噪声f(x)是由频率为nfq幅值为
Figure FDA0002658932780000023
的谐波分量构成;
步骤S2:根据步骤S1分析的结果设计出截止频率自调整的低通滤波器。
2.根据权利要求1所述的一种截止频率自调整低通滤波器的设计方法,其特征在于,上述步骤S2中,截止频率自调整的低通滤波器设计过程如下:
(1)伺服电机转速为ω时,对应的角频率为
Figure FDA0002658932780000024
可推得:
Figure FDA0002658932780000025
根据公式5,将低通滤波器的截止频率设置在fe和fq之间,可以顺利地滤除掉量化噪声,保留有用的电机转速信号;
(2)上述步骤(1)给出了滤除量化噪声的截止频率的设置范围,在步骤(1)的基础上,截止频率自调整低通滤波器可表达
Figure FDA0002658932780000026
式中,ξ=0.707,ωn为其3dB截止频率,则
Figure FDA0002658932780000027
设xe表示有量化噪声的实际测量输出,xf表示滤波器输出,那么式(6)可写为
Figure FDA0002658932780000031
其中,xf表示估计的位移信号,
Figure FDA0002658932780000032
表示估计的速度信号。
3.一种截止频率自调整低通滤波器,其特征在于,包括第一放大器、第一积分器、第二放大器和第二积分器;所述第一放大器的输出端和第二放大器的输出端均与第一积分器的输入端相连,所述第一积分器的输出端分别与第二放大器和第二积分器的输入端相连,所述第二积分器的输出端的一路与第一放大器的输入端相连;
使用时,含有量化噪声的实际位移信号xe与反馈的滤波器输出位移信号xf做差,得到的误差信号通过第一放大器进行放大,然后将放大后的误差信号与反馈的通过第二放大器放大的滤波器的估计速度信号做差,得到误差信号通过第一积分器积分,获得滤波器估计的速度信号
Figure FDA0002658932780000033
然后通过第二积分器积分,得到滤波器输出位移信号xf
CN202010897502.4A 2020-08-31 2020-08-31 一种截止频率自调整低通滤波器及设计方法 Active CN112054778B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010897502.4A CN112054778B (zh) 2020-08-31 2020-08-31 一种截止频率自调整低通滤波器及设计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010897502.4A CN112054778B (zh) 2020-08-31 2020-08-31 一种截止频率自调整低通滤波器及设计方法

Publications (2)

Publication Number Publication Date
CN112054778A true CN112054778A (zh) 2020-12-08
CN112054778B CN112054778B (zh) 2024-03-19

Family

ID=73606638

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010897502.4A Active CN112054778B (zh) 2020-08-31 2020-08-31 一种截止频率自调整低通滤波器及设计方法

Country Status (1)

Country Link
CN (1) CN112054778B (zh)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS62258221A (ja) * 1986-05-02 1987-11-10 Mitsubishi Heavy Ind Ltd 磁気軸受の制御方式
US20020150003A1 (en) * 2001-04-12 2002-10-17 Chi-Mou Chao Multimode filter
CN109714015A (zh) * 2018-12-28 2019-05-03 电子科技大学 一种基于磁介复合材料的叠层低通滤波器
WO2019129838A1 (de) * 2017-12-29 2019-07-04 Avl List Gmbh Verfahren zum filtern eines periodischen, verrauschten messsignals mit einer grundfrequenz und harmonischen schwingungsanteilen
CN111049501A (zh) * 2020-01-14 2020-04-21 西安清泰科新能源技术有限责任公司 一种电压信号采样滤波电路

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS62258221A (ja) * 1986-05-02 1987-11-10 Mitsubishi Heavy Ind Ltd 磁気軸受の制御方式
US20020150003A1 (en) * 2001-04-12 2002-10-17 Chi-Mou Chao Multimode filter
WO2019129838A1 (de) * 2017-12-29 2019-07-04 Avl List Gmbh Verfahren zum filtern eines periodischen, verrauschten messsignals mit einer grundfrequenz und harmonischen schwingungsanteilen
CN109714015A (zh) * 2018-12-28 2019-05-03 电子科技大学 一种基于磁介复合材料的叠层低通滤波器
CN111049501A (zh) * 2020-01-14 2020-04-21 西安清泰科新能源技术有限责任公司 一种电压信号采样滤波电路

Also Published As

Publication number Publication date
CN112054778B (zh) 2024-03-19

Similar Documents

Publication Publication Date Title
CN107482977B (zh) 一种永磁同步电机转子位置和转速检测方法
Su et al. A simple improved velocity estimation for low-speed regions based on position measurements only
JP6197923B1 (ja) 制御システム
CN112532133A (zh) 一种适用于永磁同步电机滤波补偿滑模自抗扰控制方法
CN113791240B (zh) 基于高阶滑膜跟踪微分器的加速度估计方法、系统、设备和介质
Ohishi et al. Robust high speed tracking servo system for optical disk system
CN112054778A (zh) 一种截止频率自调整低通滤波器及设计方法
Nguyen et al. Auto-calibration and noise reduction for the sinusoidal signals of magnetic encoders
CN210380707U (zh) 一种基于参数已知的电机控制系统
CN116520684A (zh) 基于自抗扰控制和Youla参数化的控制器优化方法
CN116594302A (zh) 一种基于干扰观测的压电定位平台降噪控制方法
CN110879527A (zh) 一种基于改进型自抗扰的位置角控制器
Deur et al. A comparative study of servosystems with acceleration feedback
US6810366B2 (en) Method and apparatus for filtering a signal using a deadband
Yan et al. An improved position detection method for permanent magnet linear motor using linear hall sensors
CN115313931A (zh) 基于aekf的永磁同步电动机无传感器矢量控制方法
Ruan et al. Modeling research of MEMS gyro drift Based on Kalman filter
CN110022137B (zh) 一种简易互补融合滤波及微分估计方法
Amalkar et al. Design and implementation of sensorless speed control for induction motor drive using an optimized extended kalman filter
CN116301081B (zh) 一种惯性测试设备的速率控制方法、装置、设备及介质
CN110750047A (zh) 一种自适应互耦pid协同控制理论新方法
CN115102448B (zh) 一种永磁同步电机位置传感器谐波干扰的提取和补偿方法
US20240151279A1 (en) Method for standstill control of a multi-body system
Jarzebowicz et al. Angle tracking observer for filtering rotor position estimates in sensorless electric drives
Yu et al. Discrete-Time Sliding Mode Filter with Parameter adaptation

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