基于直接泄漏滤波x最小均方算法的海水泵振动控制应用研究-刘锦春.pdf
《基于直接泄漏滤波x最小均方算法的海水泵振动控制应用研究-刘锦春.pdf》由会员分享,可在线阅读,更多相关《基于直接泄漏滤波x最小均方算法的海水泵振动控制应用研究-刘锦春.pdf(8页珍藏版)》请在淘文阁 - 分享文档赚钱的网站上搜索。
1、第48卷第3期2叭6年5月四川大学学报(工程科学版)JOURNAL OF SICHUAN UNIVERSITY(ENGINEERINC SCIENCE EDITl0N)V0148 No3Mav 2016文章编号:1009-3087(2016)03旬148旬8 DOI:1015961jjsuese201603020基于直接泄漏滤波X最小均方算法的海水泵振动控制应用研究刘锦春,朱石坚,何其伟,丁少春(海军工程大学动力工程学院,湖北武汉430033)摘要:针对参考信号频率与振源频率偏差造成往复海水泵自适应主动控制性能下降的问题,提出直接泄漏滤波x最小均方的自适应振动控制算法。通过一种几何关系来进行参
2、考信号的相位更新,并基于泄露滤波x最小均方误差方法,设计了直接进行频率函数控制更新和控制权系数更新的振动控制算法。通过泄露因子限制控制输出能量。为验证该方法的有效性,进行了仿真。仿真结果表明,该方法能够通过自调节参考信号频率,跟踪振源信号变化,取得良好的振动抑制效果。进行了水泵试验台架的振动控制试验,检验该自适应控制算法性能。试验结果表明,该方法能够对频率缓慢变化的水泵振动进行有效的自适应振动控制。关键词:振动;振动控制;往复泵;自适应算法中图分类号:032 文献标志码:AR姻咖ch AppH姐廿伽of A嘣Ve Water Pump Vibmtion Control B嬲ed仰Direct
3、Leaky Fntered-X Le嬲t Mean-square AlgorithmUU J诹hun,zHU Sh诤胍。HE Qiet,DING Shchun(Couege of P0wer Eng,Naval univof Eng,Wuh帅430033,China)Abstmct:In order to tackle the pe怕nflaIlce deFadation of active recipmcating pump vibr“on controL system caused by frequency errorbetween reference sigIlal and primar
4、y vibration source,a direct appmach_based adaptive vibration contml a190rithm waLs pmposed埘th adaptive semgenerating reference s培nalThe phase 0f refbrence signal was updated based on geometry relationship,besides tIle fequencyfunction recursion and tlle contml coefficients recursions were desigIled
5、directly based on tlle leaky 6ltered-x least meansquare algorithmControl efrort w踮limited by the leaky factorsSimulation aIld experiment were ca币ed out to validate tle efrectiveness of the pmposed methodE丘ective suppression of、ribmtion was achieved by adaptively tracking the vibration source f托quenc
6、y in simlllations,andtlle e矗ective pedomance of pmposed algorithm was coIlfiHned in experiments of slowly varying f诧quency water reciprocating pump、,ibKnionKey wortls:vibration;vibration control;recipmcating pump;adaptive a190rithms舰船机械设备是舰船线谱噪声的主要来源,严重影响其声隐身性能。隔离舰船线谱振动噪声传播的方法主要有两种,即被动控制方法和主动控制方法悼J。
7、被动控制方法主要是通过橡胶隔振器等被动器件来实现线谱振动的隔离衰减。目前,国内外船舶上已经广泛采用包括空气弹簧、钢丝绳弹簧等被动隔振元器件。虽然这些被动控制方法能够在一定程度上有效抑制高频振动噪声的传播,但是,在低频振动的隔离上难以取得有效的控制效果旧。J。振动主动控制方法作为一种有效抑制低频振动的手段,也得到广泛研究。自适应控制方法由于其计算简单方便,稳定性好,故实际应用广泛一6l。针对船舶电压等因素变化使旋转机械设备工况频率与实际振源线谱信号频率存在波动和偏差从而引起的振动主动控制性能下降的问题,传统方法通常是取靠近振源的信号作为理想的线谱振源信号,但由于受到环境限制,往往很难获得精确的振
8、源信收稿日期:2015一04一Ol基金项目:国家自然科学基金资助项目(51179197);高等学校全国优秀博士学位论文作者专项资金资助项目(201057)作者简介:刘锦春(1987一),男,博士生研究方向:舰船振动噪声控制E-mail:www333coca163com万方数据第3期 刘锦春,等:基于直接泄漏滤波x最小均方算法的海水泵振动控制应用研究 149号信息,故难以准确估计振源信号频率。将工况频率作为参考信号的频率,则必然与实际振源信号频率存在偏差,振动抑制效果就会降低、恶化6。1引。反馈式自适控制算法是解决这类问题的一种有效方法jn一7:文献6,11,13一15等自适应反馈振动控制方法,
9、都需要先拟合生成振源信号,间接通过频率估计方法估计频率和生成参考信号,再按照已知频率的自适应控制方法进行控制。因此,这些方法的性能差异主要取决于所采用的频率估计方法。文献6通过迭代最小无损响应谱方差的频率估计方法估计拟合的振源信号频率;文献1l,1415采用自适应陷波器方法估计频率,自适应陷波器频率收敛速度决定了其控制性能的好坏;文献13采用子空间辨识方法,辨识拟合振源信号频率,能够有效地跟踪识别未知波动的线谱信号。相比文献6,11,13一15,文献12,1718等自适应反馈控制方法无需拟合需要抑制的振源信号,可同时直接进行频率和信号的生成:文献1217通过余弦函数信号生成器来生成参考信号,并
10、结合传统算法解决了这类参考信号频率存在偏差引起的性能下降问题:其巾,文献12等提出利用2阶自回归模型的余弦函数信号生成器对参考信号进行预处理,能够在一定程度上鳃决频率不一致的问题,但理论和仿真分析表明该方法在线谱信号较低信噪比时系统不稳定。文献17对此进行了改进和发展,通过并行方式生成参考信号。文献18依据另一几何原理,通过正弦函数进行了参考信号生成,并形成了一种控制方法,但其并未指出由两种参考信号生成结构的差异引起的算法整体性能的差异;而且文献18中若仅按照文献12方法进行更新,则也存在与文献12相同的稳定性不佳的问题。作者在文献18算法基础上,指出了两种参考信号生成模式的差异与相互联系,并
11、主要针对文献18的3个部分进行改进。首先,针对文献18中频率更新方程计算较为复杂的缺点,通过近似化简处理,大大减小了计算量(这种改进同样适合文献17)。其次,针对文献18直接型带通滤波器带通参数调整较为复杂,采用格型带通器进行带通滤波,可分别控制调整带通频率和工作频宽,且对非目标线谱。具有更强的衰减。最后,针对文献18无法在控制过程中控制输出能量的饱和,通过采用泄漏滤波x最小均方算法(1eaky 6lteredx least meansquared)进行饱和抑制。针对时变频率振动信号进行J,自适应主动控制仿真,仿真结果表明,本文方法可以在存在频率偏差的情况下,收敛到振动信号频率,并跟踪其频率变
12、化,抑制时变振动。进行了海水泵振动主动控制试验,实验结果表明,本文方法对频率未精确已知且缓慢变化的振动设备,可以取得良好的振动抑制效果。1 参考信号生成方法分析文献17中提出一种基于余弦函数进行参考信号更新的方法,如图1所示。其估计频率。(n)与频率相关参数的关系为c。(n)=一2cos。(凡)。而文献18中提出了另一种方法,利用正弦函数进行参考信号生成,如图2所示。相应的频率与频率相关参数的关系为c,。(n)=一2sin。(n)。y+I,轨迹Ly+l,轨迹,图l 并行方式结构图171Fig1 ParaUel fbm17t lm图2交叉方式结构图Fig2 Cross updatedfb珊18分
13、析两种更新模式在计算过程中,c“(凡)、c。i(n)均为有界函数,且需满足稳定条件ki(n)I2,I c。(n)【2,。(凡)(o,1T2)。分析两者几何原理可知,两者存在三角函数映射关系。(n)一1T2一。(n)。即图l的低频段与图2的万方数据150 四川大学学报(工程科学版) 第48卷高频段更新性能类似。而且对同一值c“(n)=c。i(凡),两者频率值之和等于2。但对同一低频信号更新即。(n)-+0时,在频率相关参数更新方程的更新过程中, lim c,i(n)l一2, 而t(”,lim c。i(n)J=o2。故基于正弦函数的更新方mIn)哪 式具有更大的数值稳定裕度。为便于理解图2的基本原
14、理,现做如下阐述。设一个幅值为l的向量信号i为:以,i=石。(忍)+魄(忍) (1)式中,戈。(n)和戈。i(凡)分别为余弦和正弦信号。设与频率相关的系数为: 气如)=-等为了实现余弦、正弦信号的更新,则可按照图2中几何关系进行参考信号的更新:o=一c“(n)_1,ie2+-2。 (3),i=洳l以工o (4)式中,1J。l为竹,I0的模。由式(3)、(4)则有向量信号i的幅值为1,相位为妒i(n)。,i通过利用旧信号孔i、以-2i来对向量信号竹,i进行更新。同时,相应的余弦和正弦信号也实现了更新。如图l所示,从0、一c“(n)饮-1ie2和-2,。构成的几何关系可知,一c“(乃)-11ie2
15、和,i之间的相位差为一9i(n)+12,一c。(几)以_1,ie2和也。之间的相位差为9i(n一1)+1T2。故从图2可以推导得到如下关系:2sini(n)=sin妒;(n)+17。o sin9i(n一1)(5)当c“(凡)收敛后,则有lim妒。(乃)一i(n)。可认为向量信号镌i是频率为i(凡)、幅值为1的信号。即戈。(n)=cos di(n) (6)并6(n)=sin“。(凡) (7)式中,ai(n)为参考信号相位。式(3)和(4)相应的参考信号实数更新方程可以表示为:戈。i,o(n)=c,i(n)菇6。(n一1)+z。(n一2)(8)戈6。o(n)=一c,;(n)戈。i(n一1)+菇6。
16、(n一2)(9)菇。(凡)=戈o(n)戈孙(n)+戎,o(n)(10)菇6。(n)=z咄o(凡)菇:加(凡)+z;加(n)(11)式中,c叫(凡)=一2sini(n),戈。(n一1)、戈6。(凡一1)分别为n一1时刻余弦和正弦参考信号,z。棚(n)、戈“。(n)表示未进行归一化的参考信号。2 基于直接泄漏滤波X最小均方原理(1昀姆劬:eredu假)的自适应控制算法为了在自适应控制过程中,防止控制系统输出饱和,将残余振动误差信号和控制器输出的总能量作为性能误差函数口J:,(n)=e2(n)+p,2(n) (12)式中:届为一个泄漏系数,用于控制控制器的总能量输出;e(n)为振动主动控制中施加次级
17、作用力之后的残余振动信号,定义为:8(凡)=d(n)一yd(n) (13)式中:d(n)为需要抑制的振动信号,通常由多个谐波信号组成;次级控制力为y。(凡)为:yd(n)=s木1(n)x(n) (14)式中:s为次级通道模型;奉表示卷积运算;控制权系数,为I,(n)=五lo,o,q。1T;x为参考信号向量,由各谐波余弦、正弦函数组成x(n)=戈。(凡),戈。,(n),戈(n),髫k(n)T,g为需要抑制的谐波信号数目。21控制权系数更新按照梯度下降算法,控制权系数控制更新方程为2|:吣+1)=w(挖)一躲(15)联立方程(12)一(14),并代入式(15),可得:,(n+1)=(1一卢p),(
18、n)+肛e(n)X(n)(16)式中:肛为控制权系数,(n)的更新步长系数渺=肛1,p。,肛:,肛:,心,肛,7;滤波参考信号向量为:x(凡)=戈。如),菇s,(n),戈7(n),戈7(n)T。其中:其中:j为次级通道s的估计模型,其系数为;。珊一:;肘为模型阶数。22频率更新方法振源信号频率与参考信号频率通常存在偏差。自适应主动控制方法需要能够自动调节参考信号的频率,跟踪线谱振源信号频率。结合泄漏滤波x最小均方算法和第1节参考信号生成方法,则c“(n)mm一一nnq虬X菇mm“=删万方数据第3期 刘锦春,等:基于直接泄漏滤波x最小均方算法的海水泵振动控制应用研究 151史新万程口J为: 气如
19、+1)弘“旷K器一如)。咖)器(”)式中,肌为更新常数。结合式(13),ae(见)ac。(乃)化简为:器“州鲰豢+K如小木器(18)式中,先将油(n)、(n)进行近似处理,认为其在卷积运算中为定值;再将s近似为其估计模型i;缸q(凡)a。“(n)、以6。(n)aq,。(n)可以通过第1节中式(8)(11)计算并化简得:批。(n) z:,(n)z6(凡一1)+石。(聘)戈6(凡)戈。(凡一1)(19)批屯(n) z:。(n)zq(n一1)+戈。(n)z6。(n)。(n一1)加“(n) fi(20)式中,i=。加2(n)+菇味。2(n)。将式(19)、(20)代人(17),则可得:c叫(n+1)=
20、c“(凡)+p。,e(n)o(n)i:l=戈;。(n)戈屯(n1)+z。(n)戈a。(n)戈。(n一1)i一。(凡)i木 戈。2(n)戈。(n一1)+菇。(凡)戈6。(凡)戈6。(n一1)fi (21)观察式(21)可知,次级通道估计模型i阶数较高时,卷积运算过程的计算量较为复杂。基于参考信号的三角函数特性,对上述计算方法进行如下简化。首先,由于z。(凡)和戈。(n)分别为余弦和正弦信号,将式(6)、(7)代人(19)、(20),可得:桊s(她k如m。(22)昙芋生专兰:一c。(仅i(n)戈。(n)。 (23)石丽2一c08 L仅i Ln J)戈dt Ln)7。Lzj式中,di(n)=ai(凡
21、)一ai(n一1)。从式(22)和(23)可知:梯度分量均含有相位差信号仅i(n),并且在式(21)中均需要进行卷积运算;虽然卷积运算过程中cos(ai(n)有变化,但由于参考信号频率与需要抑制的线谱振动信号频率偏差通常较小。假设cos(d。(n)在卷积运算过程中为定值,令cos(仪i(n)一cosi(0),其中,。(O)可设置为参考信号的频率初值。同时,类似地将幅值ti假设为定值,令z,一1。故式(21)可以化简为:c;i(n+1)=c,i(n)+肛。,cos(oi(O)e(凡)危i,o(n)戈。7(n)一i,。(n)戈。,7(n)(24)式中,滤波器参考信号z。7(n)、戈。7(n)可以直
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 基于 直接 泄漏 滤波 最小 算法 海水 振动 控制 应用 研究 刘锦春
限制150内