CN107291998B - 一种基于滑模思想的新型低通滤波器设计方法 - Google Patents

一种基于滑模思想的新型低通滤波器设计方法 Download PDF

Info

Publication number
CN107291998B
CN107291998B CN201710402551.4A CN201710402551A CN107291998B CN 107291998 B CN107291998 B CN 107291998B CN 201710402551 A CN201710402551 A CN 201710402551A CN 107291998 B CN107291998 B CN 107291998B
Authority
CN
China
Prior art keywords
sliding mode
filter
low
pass filter
signal
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
CN201710402551.4A
Other languages
English (en)
Other versions
CN107291998A (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.)
Yantai Nanshan University
Naval Aeronautical Engineering Institute of PLA
Original Assignee
Yantai Nanshan University
Naval Aeronautical Engineering Institute 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 Yantai Nanshan University, Naval Aeronautical Engineering Institute of PLA filed Critical Yantai Nanshan University
Priority to CN201710402551.4A priority Critical patent/CN107291998B/zh
Publication of CN107291998A publication Critical patent/CN107291998A/zh
Application granted granted Critical
Publication of CN107291998B publication Critical patent/CN107291998B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/30Circuit design
    • G06F30/36Circuit design at the analogue level
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00Details relating to CAD techniques
    • G06F2111/10Numerical modelling

Abstract

本发明公开了一种基于滑模思想的新型低通滤波器设计方法。首先按照一阶惯性环节传递函数,构造虚拟滤波被控对象;针对上述虚拟滤波被控对象,设计简单滑模控制器;对滑模控制微分方程组进行离散化处理,得到滤波器的离散化表达形式;输入信号,通过选取合适的滤波参数,并观察输出曲线,从而确定最终的低通滤波器参数,使得整个低通滤波器具有令人满意的低通特性。本发明的有益效果是基于非线性滑模的非线性数字滤波器,其能吸收滑模控制抗干扰强、响应快的优点,使得滤波器失真小,而且具有较强的抑制噪声能量。

Description

一种基于滑模思想的新型低通滤波器设计方法
技术领域
本发明属于滤波器设计与制造技术领域,涉及一种基于滑模思想的新型低通滤波器设计方法。
背景技术
滤波器是一种使有用频率信号通过而同时抑制无用频率的电子装置,其在工程上有着非常广泛的应用,主要应用于各类信号处理、数据传输与干扰抑制之中。而低通滤波器则广泛应用于抑制高频噪声。传统的有源滤波器基于线性传递函数理论,采用电阻电容运放等元器件组成滤波器电路;而随着计算机技术的发展,目前采用离散方法可以设计各种不同类型的数字滤波器,具有稳定性高,精度高一届设计灵活的特点,广泛应用于语音、数据传输等各个方面。
发明内容
本发明的目的在于提供一种基于滑模思想的新型低通滤波器设计方法。
本发明所采用的技术方案是按照以下步骤进行:
步骤一:按照一阶惯性环节传递函数,构造虚拟滤波被控对象;
选取基本滤波器时间常数T,构造基本的一阶惯性环节虚拟滤波被控对象,并用传递函数描述如下:
Figure BDA0001310066100000011
其中y为滤波器的输出,u为虚拟控制量,s为传递函数中的微分算子;
步骤二:针对上述虚拟滤波被控对象,设计简单滑模控制器;
将上述传递函数,定义误差变量为e=r-y,其中r为滤波器输入信号,定义柔化函数型非线性滑模面为
Figure BDA0001310066100000012
其中
Figure BDA0001310066100000021
c1、c2与ε为待设计的正常数;
构造虚拟控制量u如下:
Figure BDA0001310066100000022
步骤三:对滑模控制微分方程组进行离散化处理,得到滤波器的离散化表达形式;
采用欧拉法进行离散化,设定采样时间为Δt,则有
y(n)-y(n-1)=Δt/T(u(n-1)-y(n-1))
而其中
Figure BDA0001310066100000023
e(n-1)=r(n-1)-y(n-1),
Figure BDA0001310066100000024
Figure BDA0001310066100000025
Figure BDA0001310066100000026
以上数据中y(n)表示当前时刻的数据,y(n-1)表示相对当前时刻上一个采样时刻数据,而y(n-2)表示相当当前时刻上两个采样时刻的数据,最终所得到的y(n)为滤波器的输出。
进一步,设输入信号为低频有用信号、高频噪声信号和近似脉冲信号的混合,如下所示:
u(t)=sin(5t)+k1sin(100000t)+k2Δ
其中t为时间,sin(5t)为低频有用信号,k1sin(100000t)为高频噪声信号,Δ为近似脉冲信号,在t=t1时刻加入,左右时间长为Δt1左右;通过选取合适的滤波参数,并观察输出曲线,从而确定最终的低通滤波器参数,使得整个低通滤波器具有令人满意的低通特性。
本发明的有益效果是基于非线性滑模的非线性数字滤波器,其能吸收滑模控制抗干扰强、响应快的优点,使得滤波器失真小,而且具有较强的抑制噪声能量。
附图说明
图1是本发明提供的一种基于传递函数与数据处理的新型低通滤波器设计流程;
图2是本发明实施例所设定的混有噪声与脉冲信号的输入信号;
图3是本发明实施例提供的滤波器的输出低频信号;
图4是本发明实施例提供的滤波器的输出低频信号与原有用信号对比图;
图5是本发明实施例提供的有用信号与滤波信号相比的误差曲线。
具体实施方式
下面结合具体实施方式对本发明进行详细说明。
本发明按照滑模控制的思想,首先选取基本滤波时间常数,构建虚拟滤波被控对象,然后选取一类非线性柔化滑模面,构造简单积分滑模虚拟控制律。然后对整个闭环系统的输入输出关系微分方程组采用欧拉法进行离散化,得到新型低通滤波器的差分方程表达形式,从而实现对输入信号的低通滤波,抑制高频噪声,同时还可以剔除或减弱奇异异常数据的冲击影响。本发明如图1所示,步骤如下:
步骤一:按照一阶惯性环节传递函数,构造虚拟滤波被控对象;
选取基本滤波器时间常数T,构造基本的一阶惯性环节虚拟滤波被控对象,并用传递函数描述如下:
Figure BDA0001310066100000031
其中y为滤波器的输出,u为虚拟控制量,s为传递函数中的微分算子。
步骤二:针对上述虚拟滤波被控对象,设计简单滑模控制器;
将上述传递函数,定义误差变量为e=r-y,其中r为滤波器输入信号,定义柔化函数型非线性滑模面为
Figure BDA0001310066100000032
其中
Figure BDA0001310066100000041
c1、c2与ε为待设计的正常数。
构造虚拟控制量u如下:
Figure BDA0001310066100000042
步骤三:对滑模控制微分方程组进行离散化处理,得到滤波器的离散化表达形式;
采用欧拉法进行离散化,设定采样时间为Δt,则有
y(n)-y(n-1)=Δt/T(u(n-1)-y(n-1))
而其中
Figure BDA0001310066100000043
e(n-1)=r(n-1)-y(n-1),
Figure BDA0001310066100000044
Figure BDA0001310066100000045
Figure BDA0001310066100000046
以上数据中y(n)表示当前时刻的数据,y(n-1)表示相对当前时刻上一个采样时刻数据,而y(n-2)表示相当当前时刻上两个采样时刻的数据,其它变量含义与此相同。最终所得到的y(n)为滤波器的输出。
步骤四:输入信号的模拟;
为了说明本方法,在此假设输入信号为低频有用信号、高频噪声信号和近似脉冲信号的混合,如下所示:
u(t)=sin(5t)+k1sin(100000t)+k2Δ
其中t为时间,sin(5t)为低频有用信号,k1sin(100000t)为高频噪声信号,Δ为近似脉冲信号,在t=t1时刻加入,左右时间长为Δt1左右。上述输入信号将用于验证滤波器的低通滤波功能,滤波器的设计目标为抑制高频信号,消除脉冲信号的影响,同时保留低频有用信号。
步骤五:将步骤一至步骤三所得的滤波器,输入入步骤四所建立的输入信号,通过选取合适的滤波参数,并观察输出曲线,从而确定最终的低通滤波器参数,使得整个低通滤波器具有令人满意的低通特性。
案例实施与计算机仿真模拟结果分析
选取滤波器的基本时间常数为T=0.1,按照上述设计步骤进行仿真,选取输入信号为k1=2,k2=100,t1=5,Δt1=0.005,初始状态如下:y(0)=0,u(0)=0,选取滤波参数为:c1=10,c2=0.5,k1=0.2,ε=0.2,将步骤三所得到的滤波器输入步骤四的输入信号进行仿真,得到仿真结果图2至图5所示。
通过以上仿真结果与曲线图2至图5可以看出,图2中输入信号由于混有高频噪声信号而将有用信号完全淹没,而且幅值上,在5s时的近似脉冲信号也属于异常数据,滤波器应当将其剔除。图3为通过本设计的滤波器后保留下来的低频信号,高频信号完全得到了抑制,尤其是近似脉冲信号也得到了衰减,从峰值100减弱到4以下,而采用同样时间常数的传统滤波器时,该峰值为6左右。图4为有用信号和滤波器获得的低频信号对比图,可以看出低频有用信号得到了比较好的保留,波形有一定的相位差,幅值有一定的误差。图5给出了幅值误差曲线。由图4可以看出,滤波效果是令人满意的,波形基本没有太大的失真,而幅值消弱也较小。当然还可以通过加快滤波器响应来减少失真,但同时又会增大滤波器的通频带,使得噪声抑制功能减弱,在实际系统设计是可以根据工程设计而权衡选取参数。
从以上案例仿真结果可以看出,本发明提供的基于滑模思想的新型低通滤波器设计方法在原理上是完全正确与可行的,同时由于滑模控制的引入,使得滤波和传统线性滤波方法有了本质的不同,使得新的低通滤波器设计可调节参数更多,设计方法更为灵活,从而能够设计出更为丰富多样的低通滤波器,满足各种不同实际系统的滤波要求。同时滑模控制本身具有较好的快速性与抗干扰特性,也使得滤波器具有更好的性能。因此本发明具有很好的理论价值与实用价值,同时也丰富了滤波器设计方法。
以上所述仅是对本发明的较佳实施方式而已,并非对本发明作任何形式上的限制,凡是依据本发明的技术实质对以上实施方式所做的任何简单修改,等同变化与修饰,均属于本发明技术方案的范围内。

Claims (2)

1.一种基于滑模思想的新型低通滤波器设计方法,其特征在于按照以下步骤进行:
步骤一:按照一阶惯性环节传递函数,构造虚拟滤波被控对象;
选取基本滤波器时间常数T,构造基本的一阶惯性环节虚拟滤波被控对象,并用传递函数描述如下:
Figure FDA0002539381010000011
其中y为滤波器的输出,u为虚拟控制量,s为传递函数中的微分算子;
步骤二:针对上述虚拟滤波被控对象,设计简单滑模控制器;
将上述传递函数,定义误差变量为e=r-y,其中r为滤波器输入信号,定义柔化函数型非线性滑模面为
Figure FDA0002539381010000012
其中
Figure FDA0002539381010000013
c1、c2与ε为待设计的正常数;
构造虚拟控制量u如下:
Figure FDA0002539381010000014
步骤三:对滑模控制微分方程组进行离散化处理,得到滤波器的离散化表达形式;
采用欧拉法进行离散化,设定采样时间为Δt,则有
y(n)-y(n-1)=Δt/T(u(n-1)-y(n-1))
而其中
Figure FDA0002539381010000015
e(n-1)=r(n-1)-y(n-1),
Figure FDA0002539381010000016
Figure FDA0002539381010000021
Figure FDA0002539381010000022
y(n)表示当前时刻的数据,y(n-1)表示相对当前时刻上一个采样时刻数据,而y(n-2)表示相对当前时刻上两个采样时刻的数据,最终所得到的y(n)为滤波器的输出。
2.按照权利要求1所述一种基于滑模思想的新型低通滤波器设计方法,其特征在于:设输入信号为低频有用信号、高频噪声信号和近似脉冲信号的混合,如下所示:
u(t)=sin(5t)+k1 sin(100000t)+k2Δ
其中t为时间,sin(5t)为低频有用信号,k1 sin(100000t)为高频噪声信号,Δ为近似脉冲信号,在t=t1时刻加入,左右时间长为Δt1;通过选取滤波参数,并观察输出曲线,从而确定最终的低通滤波器参数,使得整个低通滤波器具有低通特性。
CN201710402551.4A 2017-06-01 2017-06-01 一种基于滑模思想的新型低通滤波器设计方法 Active CN107291998B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710402551.4A CN107291998B (zh) 2017-06-01 2017-06-01 一种基于滑模思想的新型低通滤波器设计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710402551.4A CN107291998B (zh) 2017-06-01 2017-06-01 一种基于滑模思想的新型低通滤波器设计方法

Publications (2)

Publication Number Publication Date
CN107291998A CN107291998A (zh) 2017-10-24
CN107291998B true CN107291998B (zh) 2020-08-07

Family

ID=60095436

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710402551.4A Active CN107291998B (zh) 2017-06-01 2017-06-01 一种基于滑模思想的新型低通滤波器设计方法

Country Status (1)

Country Link
CN (1) CN107291998B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108446491A (zh) * 2018-03-19 2018-08-24 成都欧督系统科技有限公司 一种无源低通滤波器的设计方法
CN108469572A (zh) * 2018-03-19 2018-08-31 成都欧督系统科技有限公司 一种低通自动滤波系统
DE102018222320A1 (de) * 2018-12-19 2020-06-25 Robert Bosch Gmbh Objekterkennungsvorrichtung für Fahrzeuge und Verfahren zur Erkennung eines Objektes für Fahrzeuge
CN111650345B (zh) * 2020-07-14 2021-02-19 中科三清科技有限公司 大气环境污染检测数据的处理方法、装置、设备及介质

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2002045831A1 (de) * 2000-12-07 2002-06-13 Aaflowsystems Gmbh & Co. Kg Verfahren zur herstellung eines membranmoduls
CN103034249A (zh) * 2012-12-27 2013-04-10 中国电子科技集团公司第二十六研究所 一种基于模糊决策的双轴光电转台复合控制系统和控制方法
CN103647490A (zh) * 2013-09-27 2014-03-19 天津大学 一种永磁电机的滑模控制策略
CN104468427A (zh) * 2014-10-25 2015-03-25 哈尔滨工业大学 基于fri时频域综合分析的信号高效采样及信号重构方法
CN106763983A (zh) * 2016-12-30 2017-05-31 天津市职业大学 一种新型无刷直流电机的电动执行机构控制装置和方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2002045831A1 (de) * 2000-12-07 2002-06-13 Aaflowsystems Gmbh & Co. Kg Verfahren zur herstellung eines membranmoduls
CN103034249A (zh) * 2012-12-27 2013-04-10 中国电子科技集团公司第二十六研究所 一种基于模糊决策的双轴光电转台复合控制系统和控制方法
CN103647490A (zh) * 2013-09-27 2014-03-19 天津大学 一种永磁电机的滑模控制策略
CN104468427A (zh) * 2014-10-25 2015-03-25 哈尔滨工业大学 基于fri时频域综合分析的信号高效采样及信号重构方法
CN106763983A (zh) * 2016-12-30 2017-05-31 天津市职业大学 一种新型无刷直流电机的电动执行机构控制装置和方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Cem Özgür 等.Slip cast forming of multilayer ceramic filter by fine particles migration.《Ceramics International》.2008, *
基于低通滤波器的液压挖掘机工作装置轨迹跟踪滑模控制;李发喜 等;《辽宁石油化工大学学报》;20160630;第50-53页 *

Also Published As

Publication number Publication date
CN107291998A (zh) 2017-10-24

Similar Documents

Publication Publication Date Title
CN107291998B (zh) 一种基于滑模思想的新型低通滤波器设计方法
US6181754B1 (en) System and method for modeling mixed signal RF circuits in a digital signal environment
CN108958041B (zh) 一种基于双曲正割吸引律的离散双周期重复控制方法
CN104796111A (zh) 一种用于动态迟滞系统建模与补偿的非线性自适应滤波器
CN106531145B (zh) 基于m估计器的递归有源噪声控制方法
CN105871356A (zh) 基于最大混合互相关熵准则的自适应滤波方法
CN109445275B (zh) 一种比例滞后控制方法、装置及设备
JP4292517B2 (ja) 衝撃応答スペクトルが現れる時間の差が短い衝撃波形の合成方法、及び衝撃応答内訳を求めるデジタルフィルタとそのデジタル逆フィルタ
CN113595527B (zh) 一种滤波参数确定方法、滤波方法及相关装置
CN107302349B (zh) 一种基于传递函数与数据处理的新型低通滤波器设计方法
CN109119061A (zh) 一种基于梯度矩阵的有源噪声控制方法
CN102724152B (zh) 基于Laguerre结构的多项式自适应有源噪声对消方法
Piskorowski Phase-compensated time-varying Butterworth filters
CN110989353B (zh) 一种周期扰动观测器的设计方法
Aarnoudse et al. Nonlinear iterative learning control: a frequency-domain approach for fast convergence and high accuracy
CN112198795B (zh) 机电伺服控制方法、系统、终端设备及存储介质
Karpaty Modeling Amplifiers as Analog Filters Increases SPICE Simulation Speed
CN103095251A (zh) 数字滤波高速化的改进
Pratheep et al. Model order reduction of interval systems by pole clustering technique using GA
Zhu et al. An input shaping method based on system output
Kanjanasurat et al. Gain equalizer by four parameters logistic bézier bernstein polynomial
CN111641400B (zh) 一种cic滤波器组的传递函数等效方法
CN114398755B (zh) 一种弹性滤波器设计方法
Zeinaly et al. Robust stability for bilateral teleoperation systems under time-varying delay by communication disturbance observer
Jalali et al. Line-of-sight stabilization by robust L1 controller based on Linear Matrix Inequality (LMI) approach

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