CN106444390A - 一种基于fir滤波器和分数阶重复控制器的磁悬浮转子谐波电流抑制方法 - Google Patents
一种基于fir滤波器和分数阶重复控制器的磁悬浮转子谐波电流抑制方法 Download PDFInfo
- Publication number
- CN106444390A CN106444390A CN201611113486.5A CN201611113486A CN106444390A CN 106444390 A CN106444390 A CN 106444390A CN 201611113486 A CN201611113486 A CN 201611113486A CN 106444390 A CN106444390 A CN 106444390A
- Authority
- CN
- China
- Prior art keywords
- rotor
- harmonic
- current
- repetitive controller
- fir filter
- 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
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
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Computation (AREA)
- Medical Informatics (AREA)
- Software Systems (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Magnetic Bearings And Hydrostatic Bearings (AREA)
Abstract
本发明公开了一种基于FIR滤波器和分数阶重复控制器的磁悬浮转子谐波电流抑制方法,首先建立含质量不平衡和传感器谐波的磁悬浮转子动力学模型,其次FIR滤波器是具有线性相位特性的低通滤波器,分数阶重复控制器由分数延时滤波器近似得到,通过分数延时滤波器系数的在线更新,本发明能实现定转速下的谐波电流抑制,适用于存在质量不平衡和传感器谐波的磁悬浮转子谐波电流抑制。将FIR低通滤波器和分数阶重复控制器结合使用,应用FIR滤波器的线性相位特性和分数阶重复控制器系数在线更新实现谐波电流的抑制。
Description
技术领域
本发明涉及磁悬浮转子谐波电流抑制的技术领域,具体涉及一种基于FIR滤波器和分数阶重复控制器的磁悬浮转子谐波电流抑制方法,用于对磁悬浮飞轮转子系统中的谐波电流进行抑制,为磁悬浮飞轮在“超静”卫星平台上的应用提供技术支持。
背景技术
磁轴承转子系统具有无摩擦、长寿命和主动振动可控等特性,在磁悬浮飞轮、磁悬浮陀螺以及磁悬浮分子泵等方面有良好的应用。由于加工精度有限,磁悬浮转子不可避免地会存在质量不平衡,在高速转动过程中会产生与转速频率相同的同频控制电流;另一方面,由于传感器检测面、被检测面上电或者磁特性不一致,通过实验发现即使没有不平衡量,传感器检测信号中仍然含有同频及倍频噪声,这就是传感器谐波,传感器谐波会引发谐波控制电流。由磁悬浮转子质量不平衡和传感器谐波所产生的谐波控制电流,会使磁轴承产生谐波振动力,进而传递到基座影响超静卫星平台的姿态控制精度。
谐波振动抑制可以分为零电流、零位移和零振动三类,其中零电流可以用最少的计算量和功耗抑制大部分的振动,本发明将磁悬浮转子的谐波电流视为扰动进行谐波电流抑制,实现零电流。现有技术主要针对单一频率的干扰进行抑制,对于谐波扰动抑制研究相对较少,主要可以分为两类。第一类方法针对不同频率的振动并联多个滤波器,如并联多个陷波器或多个LMS(Least Mean Square,LMS)滤波器等。该方法不能针对所有振动同时抑制,计算量大,且需要考虑不同滤波器间的收敛速度问题,设计起来比较复杂。第二类方法无需并联多个滤波器便可实现对不同频率成分振动的同时抑制,如重复控制算法。现有的应用于磁悬浮转子控制系统的重复控制算法均没有频率自适应能力,只能实现某些固定转速下的谐波电流抑制,本发明提出一种频率自适应重复控制器以实现磁悬浮转子的谐波电流抑制。
发明内容
本发明的目的为:克服现有技术的不足,发明一种基于FIR滤波器和分数阶重复控制器的磁悬浮转子谐波电流抑制方法,通过在线更新分数阶重复控制器系数实现谐波电流抑制。
本发明采用的技术方案为:一种基于FIR滤波器和分数阶重复控制器的磁悬浮转子谐波电流抑制方法,包括以下步骤:
步骤(1)建立含质量不平衡和传感器谐波的磁悬浮转子动力学模型
磁悬浮飞轮通过主动磁轴承控制转子实现稳定悬浮,令O表示磁轴承定子的几何中心,G表示转子的几何中心,I表示转子的质心,以O为中心建立惯性坐标系OXYZ,(x,y)表示转子几何中心G在惯性坐标系下的坐标值,由于转子具有对称结构,故在X方向上对其谐波扰动来源以及控制算法进行分析与研究,
根据牛顿第二定律,磁悬浮转子在X方向的动力学方程如下:
其中,表示转子在X方向的加速度,m表示转子质量,fx表示磁轴承在X方向的轴承力,fu表示转子的不平衡力,可表示如下:
fu=meΩ2cos(Ωt+φ)
其中e表示转子几何中心与质心之间的偏差,Ω表示转子转速,φ表示转子不平衡质量的初始相位;
当转子在磁轴承中心位置小范围悬浮时,磁轴承的电磁力可近似表示为线性化方程:
fx≈Kxx+Kii
其中Kx和Ki分别表示磁轴承位移刚度和电流刚度,i表示磁轴承线圈控制电流;
由于机械加工精度和材料的不均匀因素的影响,磁悬浮转子的位移传感器检测面会出现圆度不理想、材质不均匀、剩磁特性不同的现象,位移传感器的输出将会出现同频和倍频的多谐波信号:
xs(t)=x(t)+xd(t)
其中x(t)表示转子几何中心真实的坐标值,xs(t)表示传感器的输出值,xd(t)为传感器输出值与真实值的误差,可表示如下:
其中l表示谐波次数,cl表示谐波系数,θl表示谐波初始相位;
将i、xd(t)、fu依次进行拉普拉斯变换可得I(s)、Xd(s)、Fu(s),则磁轴承电流I(s)的传递函数可表示如下:
其中,Gc(s)表示反馈控制器的传递函数,Gw(s)表示功放环节的传递函数,Gp(s)表示磁悬浮转子的传递函数,R(s)表示参考输入信号,Ks表示传感器增益;
综合以上分析可得,转子质量不平衡以及传感器误差会使磁轴承产生谐波控制电流,从而产生谐波振动力;
步骤(2)设计基于FIR滤波器和分数阶重复控制器的谐波电流抑制算法进行磁悬浮转子谐波电流抑制
以谐波电流为控制目标,将谐波电流i输入至重复控制器,重复控制器是基于内模原理实现输入信号中谐波量的消除,在实际磁轴承控制系统中,采样频率与谐波电流基频的比值一般不为整数,现有的针对磁轴承谐波电流抑制的重复控制器只能对其整数部分进行补偿,谐波电流抑制精度大大降低,针对这一不足,采用FIR滤波器和分数阶重复控制器相结合的方式实现谐波电流的抑制,FIR滤波器是具有线性相位特性的低通滤波器,分数阶重复控制器由分数延时滤波器近似得到,通过分数延时滤波器系数的在线更新,可以实现谐波电流抑制。
更进一步的,所述的步骤(2)谐波电流抑制算法为:
利用整数延时环节、FIR滤波器和分数阶延时环节相串联,Ts是采样周期;是采样周期的整数延时环节;F(s)是具有线性相位的FIR滤波器,滞后相位与频率ω成比例关系,可以获知并进行补偿,解决了低通滤波器相位滞后带来的影响;是采样周期的分数阶延时环节,由分数延时滤波器近似得到,采用以上三者相串联的形式,一方面可以消除低通滤波器相位滞后带来的误差,增大系统带宽;另一方面可以消除采样频率与谐波扰动信号基频的比值不为整数情况下带来的误差,实现磁轴承谐波电流抑制。
更进一步的,该谐波电流抑制算法以参考输入信号R(s)和等效谐波扰动信号D(s)作为输入,磁轴承线圈电流I(s)作为输出的灵敏度函数S2(s)可表示如下:
其中,表示除去重复控制器时以电流I(s)作为输出的灵敏度函数,C(s)表示相位补偿器。低通滤波器F(s)的截止频率ωc大于有效谐波扰动的最高频率ωmax,在ω∈(0,ωmax)范围内F(s)的幅值衰减很小,可近似认为|F(s)|=1。
本发明基本原理:对磁悬浮飞轮来讲,高频振动会降低卫星平台的指向精度和稳定度,必须加以抑制。其中,振动的主要来源是质量不平衡和传感器谐波。本发明针对谐波电流进行抑制,减小谐波振动。由于质量不平衡和传感器谐波的存在,导致控制电流中含有谐波,即谐波电流,从而使磁悬浮飞轮中含有谐波振动。通过建立含质量不平衡和传感器谐波的磁悬浮转子动力学模型,分析谐波电流,提出一种频率自适应重复控制器以实现磁悬浮转子高转速下的谐波电流抑制,重点从三个方面进行研究:与延时环节串联的具有线性相位特性的FIR低通滤波器设计,FIR滤波器滞后相位可以获知,以便于后续进行分析和补偿;分数阶延时滤波器设计,当转子转速改变时,可通过改变分数阶延时滤波器的系数实现小数部分补偿;相位补偿环节设计,采用重构谱和小增益理论进行稳定性分析,通过设计相位补偿环节以保证稳定性,最终实现磁悬浮转子谐波电流的抑制,
本发明与现有技术相比的优点在于:
(1)为了有效抑制磁悬浮转子系统中的谐波电流,本发明提出一种基于FIR滤波器和分数阶重复控制器的磁悬浮转子谐波电流抑制方法,能实现定转速下的谐波电流抑制,适用于存在质量不平衡和传感器谐波的磁悬浮转子谐波电流抑制。
(2)本发明将FIR低通滤波器和分数阶重复控制器结合使用,应用FIR滤波器的线性相位特性和分数阶重复控制器系数在线更新实现谐波电流的抑制。
附图说明
图1为本发明的流程图;
图2为磁悬浮飞轮结构示意图,其中,1为径轴一体位移传感器,2为径向磁轴承,3为轴向磁轴承,4为惯性轴,5为几何轴,6为磁悬浮转子;
图3为传感器谐波示意图;
图4为X通道磁轴承转子控制系统框图;
图5为改进插入式重复控制器系统框图;
图6为拉格朗日插值多项式幅频特性曲线,其中,图6(a)为n=2时拉格朗日插值多项式幅频特性曲线,图6(b)为n=3时拉格朗日插值多项式幅频特性曲线;
图7为简化后的改进重复控制器系统方框图。
具体实施方式
下面结合附图以及具体实施例进一步说明本发明。
如图1所示,一种基于FIR滤波器和分数阶重复控制器的磁悬浮转子谐波电流抑制方法的实施过程是:首先建立含质量不平衡和传感器谐波的磁悬浮转子动力学模型,然后设计一种基于FIR滤波器和分数阶重复控制器的磁悬浮转子谐波电流抑制方法。
(1)建立含质量不平衡和传感器谐波的磁悬浮转子动力学模型
磁悬浮转子系统结构示意图如图2所示,O表示磁轴承定子的几何中心,G表示转子的几何中心,I表示转子的质心。以O为中心建立惯性坐标系OXYZ,(x,y)表示转子几何中心G在惯性坐标系下的坐标值。由于转子具有对称结构,故在X方向上对其谐波扰动来源以及控制算法进行分析与研究。
根据牛顿第二定律,磁悬浮转子在X方向的动力学方程如下:
其中,表示转子在X方向的加速度,m表示转子质量,fx表示磁轴承在X方向的轴承力,fu表示转子的不平衡力,可表示如下:
fu=meΩ2cos(Ωt+φ)
其中e表示转子几何中心与质心之间的偏差,Ω表示转子转速,φ表示转子不平衡质量的初始相位。
当转子在磁轴承中心位置小范围悬浮时,磁轴承的电磁力可近似表示为线性化方程:
fx≈Kxx+Kii
其中Kx和Ki分别表示磁轴承位移刚度和电流刚度,i表示磁轴承线圈控制电流。
由于机械加工精度和材料的不均匀等因素的影响,磁悬浮转子的位移传感器检测面会出现圆度不理想、材质不均匀、剩磁特性不同的现象,如图3所示,位移传感器的输出将会出现同频和倍频的多谐波信号:
xs(t)=x(t)+xd(t)
其中x(t)表示转子几何中心真实的坐标值,xs(t)表示传感器的输出值,xd(t)为传感器输出值与真实值的误差,可表示如下:
其中l表示谐波次数,cl表示谐波系数,θl表示谐波初始相位。
磁轴承X方向平动控制系统如图4所示,将i、xd(t)、fu依次进行拉普拉斯变换可得I(s)、Xd(s)、Fu(s),则磁轴承电流I(s)的传递函数可表示如下:
其中,Gc(s)表示反馈控制器的传递函数,Gw(s)表示功放环节的传递函数,Gp(s)表示磁悬浮转子的传递函数,R(s)表示参考输入信号,Ks表示传感器增益;
结合以上分析可得,转子质量不平衡以及传感器误差会使磁轴承产生谐波控制电流,从而产生谐波振动力。
(2)设计基于FIR滤波器和分数阶重复控制器的磁悬浮转子谐波电流抑制方法
针对步骤(1)线圈电流中存在谐波电流这一问题,本发明采用一种基于FIR滤波器和分数阶重复控制器的磁悬浮转子谐波电流抑制方法。如图5所示,其中Ts是采样周期;是采样周期的整数延时环节;F(s)是具有线性相位的FIR滤波器,滞后相位与频率ω成比例关系,可以获知并进行补偿,解决了低通滤波器相位滞后带来的影响;是采样周期的分数阶延时环节,由分数延时滤波器近似得到。采用以上三者相串联的形式,一方面可以消除低通滤波器相位滞后带来的误差,增大系统带宽;另一方面可以消除采样频率与谐波扰动信号基频的比值不为整数情况下带来的误差,实现定转速下的磁轴承谐波电流抑制。
以参考输入信号R(s)和等效谐波扰动信号D(s)作为输入,磁轴承线圈电流I(s)作为输出的灵敏度函数S2(s)可表示如下:
其中,表示除去重复控制器时以电流I(s)作为输出的灵敏度函数,C(s)表示相位补偿器。低通滤波器F(s)的截止频率ωc大于有效谐波扰动的最高频率ωmax,在ω∈(0,ωmax)范围内F(s)的幅值衰减很小,可近似认为|F(s)|=1。
1.FIR低通滤波器分析
将FIR滤波器替换普通的低通滤波器,是因为FIR滤波器既能满足低通滤波器的特性,并且其相位与频率成线性关系,FIR滤波器的通用表达式如下:
其中M表示FIR滤波器的阶数,ai表示FIR滤波器的系数。当ai满足偶对称时,FIR滤波器的相位可表示如下:
其中,表示比例系数,Ts=0.0002s表示系统采样周期。
设计过程如下:
a)、本发明应用的磁悬浮转子系统,转子最高转速4800rpm,有效谐波扰动主要表现在同频、二倍频…七倍频,也就是有效谐波扰动的最高频率ωmax=3519rad/s。为了保证一定的稳定裕度和谐波抑制精度,FIR滤波器的截止频率ωc应该大于有效谐波扰动的最大频率,但是截止频率ωc过大会破坏系统稳定性,本发明中采用其中表示系统采样频率,则有ωc=6283rad/s。
b)、FIR滤波器的过渡带宽应该尽可能小,由于在FIR设计方法中,三角窗函数可以实现最小的过渡带宽,在本发明中,采用三角窗函数方法设计FIR滤波器,其过渡带宽BW可表示如下:
从上式可以看出,FIR滤波器阶次M越高,过渡带宽越小。但是,随着M的升高,FIR滤波器的系数增多,计算量也将急剧升高。本发明中过渡带宽要求不大于对应的选取FIR滤波器阶次M=9。
2.分数阶延时环节分析
在工程应用中,分数阶延时环节无法直接应用,需要找到一种替换形式。分数阶延时环节可用一种拉格朗日插值多项式来近似表示:
其中系数Dk可表示如下:
多项式与分数阶延时环节的差值Rn可表示如下:
其中ξ∈[Tk,Tk+1],Tk和Tk+1分别表示第k个和第k+1个采样时刻。从上式可以看出,随着拉格朗日插值多项式阶数n的增大,Rn逐渐减小,即拉格朗日插值多项式的近似程度逐渐升高,但是,随着n的增大,算法计算量大幅度增大。在实际工程中,应该综合考虑差值Rn和算法计算量两个因素,下面分别给出n=2和n=3两种情况下A从0到0.9变化时的拉格朗日插值多项式幅频特性曲线。
图6(a)和图6(b)两图分别给出在n=2和n=3两种情况下拉格朗日插值多项式的幅频特性曲线,拉格朗日插值多项式的截止频率均高于系统截止频率,并且在本系统截止频率范围(0,ωc)内,当n=2时,多项式最大幅值衰减仅为-0.565dB,与分数延时环节的近似程度极高,完全可以满足差值Rn尽量小的要求,并且计算量也相对较小,所以在本发明中选取n=2。
3.稳定性分析
根据对FIR滤波器和分数阶延时环节的分析可知,当ω∈(0,ωc)时,FIR滤波器的相位幅值|F(z)|≈1,分数阶延时环节通过拉格朗日插值多项式近似表示后,多项式的相位为-A·Tsω,幅值在对加入算法后的系统进行稳定性分析时,可作如下近似:
其中T表示磁悬浮转子的转动周期。
当ω∈(0,ωc)时,可将图5进行简化,如图7所示:其中相位补偿函数C(s)可以表示为:
其中Krc表示改进重复控制器增益,Kf(s)表示在低频段和中频段的相位补偿函数,表示高频段的相位补偿函数。
由图7可以得出加入改进重复控制器后系统的闭环特征方程,表示如下:
M(s)-N(s)e-Ts=0
其中,M(s)=1+Gc(s)Gw(s)Gp(s)Ks,N(s)=1+C(s)Gw(s)+Gc(s)Gw(s)Gp(s)Ks。
加入改进重复控制器后的系统重构谱函数R(ω)可表示如下:
重构谱函数可以作为一种判断系统稳定性的依据:根据最小增益理论可知,对于一个稳定系统,如果加入重复控制器后系统重构谱函数能满足下式所述的条件,则新的系统也是稳定的。
R(ω)<1,ω∈(0,ωc)
定义函数G(s):
其中G(s)|s=jω=L(ω)eiθ(ω),加入重复控制器后系统的稳定条件可等效为:
其中令λ(ω)=θ(ω)+θb(ω)+NhTsω,将上式通过欧拉公式展开如下:
|KrcL(ω)·Kb(ω)cosλ(ω)+jKrcL(ω)·Kb(ω)sinλ(ω)+1|<1
将上式等号两边的式子同时平方可得:
[KrcL(ω)·Kb(ω)]2<-2KrcL(ω)·Kb(ω)cosλ(ω)
由于同时满足Krc>0、L(ω)>0、Kb(ω)>0,所以上式可化简如下:
KrcL(ω)·Kb(ω)<-2cosλ(ω)
为了保证上式是有解的,必须保证cosλ(ω)<0,即
90°<λ(ω)<270°
综上所述,通过串联合适的相位补偿函数和增益系数,可以保证加入算法后系统的稳定性。
本发明说明书中未作详细描述的内容属于本领域专业技术人员公知的现有技术。
Claims (3)
1.一种基于FIR滤波器和分数阶重复控制器的磁悬浮转子谐波电流抑制方法,其特征在于:包括以下步骤:
步骤(1)建立含质量不平衡和传感器谐波的磁悬浮转子动力学模型
磁悬浮飞轮通过主动磁轴承控制转子实现稳定悬浮,令O表示磁轴承定子的几何中心,G表示转子的几何中心,I表示转子的质心,以O为中心建立惯性坐标系OXYZ,(x,y)表示转子几何中心G在惯性坐标系下的坐标值,由于转子具有对称结构,故在X方向上对其谐波扰动来源以及控制算法进行分析与研究,
根据牛顿第二定律,磁悬浮转子在X方向的动力学方程如下:
其中,表示转子在X方向的加速度,m表示转子质量,fx表示磁轴承在X方向的轴承力,fu表示转子的不平衡力,可表示如下:
fu=meΩ2cos(Ωt+φ)
其中e表示转子几何中心与质心之间的偏差,Ω表示转子转速,φ表示转子不平衡质量的初始相位;
当转子在磁轴承中心位置小范围悬浮时,磁轴承的电磁力可近似表示为线性化方程:
fx≈Kxx+Kii
其中Kx和Ki分别表示磁轴承位移刚度和电流刚度,i表示磁轴承线圈控制电流;
由于机械加工精度和材料的不均匀因素的影响,磁悬浮转子的位移传感器检测面会出现圆度不理想、材质不均匀、剩磁特性不同的现象,位移传感器的输出将会出现同频和倍频的多谐波信号:
xs(t)=x(t)+xd(t)
其中x(t)表示转子几何中心真实的坐标值,xs(t)表示传感器的输出值,xd(t)为传感器输出值与真实值的误差,可表示如下:
其中l表示谐波次数,cl表示谐波系数,θl表示谐波初始相位;
将i、xd(t)、fu依次进行拉普拉斯变换可得I(s)、Xd(s)、Fu(s),则磁轴承电流I(s)的传递函数可表示如下:
其中,Gc(s)表示反馈控制器的传递函数,Gw(s)表示功放环节的传递函数,Gp(s)表示磁悬浮转子的传递函数,R(s)表示参考输入信号,Ks表示传感器增益;
综合以上分析可得,转子质量不平衡以及传感器误差会使磁轴承产生谐波控制电流,从而产生谐波振动力;
步骤(2)设计基于FIR滤波器和分数阶重复控制器的谐波电流抑制算法进行磁悬浮转子谐波电流抑制
以谐波电流为控制目标,将谐波电流i输入至重复控制器,重复控制器是基于内模原理实现输入信号中谐波量的消除,在实际磁轴承控制系统中,采样频率与谐波电流基频的比值一般不为整数,现有的针对磁轴承谐波电流抑制的重复控制器只能对其整数部分进行补偿,谐波电流抑制精度大大降低,针对这一不足,采用FIR滤波器和分数阶重复控制器相结合的方式实现谐波电流的抑制,FIR滤波器是具有线性相位特性的低通滤波器,分数阶重复控制器由分数延时滤波器近似得到,通过分数延时滤波器系数的在线更新,可以实现谐波电流抑制。
2.根据权利要求1所述的一种基于FIR滤波器和分数阶重复控制器的磁悬浮转子谐波电流抑制方法,其特征在于:所述的步骤(2)谐波电流抑制算法为:
利用整数延时环节、FIR滤波器和分数阶延时环节相串联,Ts是采样周期;是采样周期的整数延时环节;F(s)是具有线性相位的FIR滤波器,滞后相位与频率ω成比例关系,可以获知并进行补偿,解决了低通滤波器相位滞后带来的影响;是采样周期的分数阶延时环节,由分数延时滤波器近似得到,采用以上三者相串联的形式,一方面可以消除低通滤波器相位滞后带来的误差,增大系统带宽;另一方面可以消除采样频率与谐波扰动信号基频的比值不为整数情况下带来的误差,实现磁轴承谐波电流抑制。
3.根据权利要求1所述的一种基于FIR滤波器和分数阶重复控制器的磁悬浮转子谐波电流抑制方法,其特征在于:该谐波电流抑制算法以参考输入信号R(s)和等效谐波扰动信号D(s)作为输入,磁轴承线圈电流I(s)作为输出的灵敏度函数S2(s)可表示如下:
其中,表示除去重复控制器时以电流I(s)作为输出的灵敏度函数,C(s)表示相位补偿器,低通滤波器F(s)的截止频率ωc大于有效谐波扰动的最高频率ωmax,在ω∈(0,ωmax)范围内F(s)的幅值衰减很小,可近似认为|F(s)|=1。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611113486.5A CN106444390B (zh) | 2016-12-06 | 2016-12-06 | 一种基于fir滤波器和分数阶重复控制器的磁悬浮转子谐波电流抑制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611113486.5A CN106444390B (zh) | 2016-12-06 | 2016-12-06 | 一种基于fir滤波器和分数阶重复控制器的磁悬浮转子谐波电流抑制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106444390A true CN106444390A (zh) | 2017-02-22 |
CN106444390B CN106444390B (zh) | 2019-03-29 |
Family
ID=58217520
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611113486.5A Active CN106444390B (zh) | 2016-12-06 | 2016-12-06 | 一种基于fir滤波器和分数阶重复控制器的磁悬浮转子谐波电流抑制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106444390B (zh) |
Cited By (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106873655A (zh) * | 2017-04-13 | 2017-06-20 | 北京航空航天大学 | 一种基于有限维重复控制的磁轴承系统多谐波振动抑制方法 |
CN106886152A (zh) * | 2017-04-21 | 2017-06-23 | 北京航空航天大学 | 一种基于二阶奇次重复控制器的磁悬浮转子奇次谐波电流抑制方法 |
CN106896726A (zh) * | 2017-04-21 | 2017-06-27 | 北京航空航天大学 | 一种基于并联式二阶重复控制dpsorc的数字控制器的磁悬浮转子谐波电流抑制方法 |
CN107844052A (zh) * | 2017-10-14 | 2018-03-27 | 河南工业大学 | 基于并联式有限维重复控制的磁轴承系统谐波电流抑制方法 |
CN107870568A (zh) * | 2017-12-21 | 2018-04-03 | 北京航空航天大学 | 一种基于二阶双模重复控制的磁悬浮转子谐波电流抑制方法 |
CN107966906A (zh) * | 2017-11-13 | 2018-04-27 | 南京航空航天大学 | 基于采样控制分离原理的分数阶延时实现方法 |
CN108227499A (zh) * | 2018-01-16 | 2018-06-29 | 北京航空航天大学 | 一种谐波振荡器并联分数阶重复控制器的磁悬浮转子谐波电流抑制方法 |
CN108490777A (zh) * | 2018-03-19 | 2018-09-04 | 北京航空航天大学 | 一种基于改进奇次重复控制的磁悬浮转子谐波振动力抑制方法 |
CN108535541A (zh) * | 2018-03-19 | 2018-09-14 | 湖北工业大学 | 一种带有数字工频滤波的相位比较装置及方法 |
CN109164700A (zh) * | 2018-08-07 | 2019-01-08 | 中国科学院光电技术研究所 | 一种抑制全频扰动的控制方法 |
CN111473049A (zh) * | 2020-04-17 | 2020-07-31 | 河海大学 | 一种实心定子磁悬浮励磁电流的控制方法 |
CN111708278A (zh) * | 2020-06-24 | 2020-09-25 | 北京航空航天大学 | 一种基于最小均方差前馈补偿算法的同频振动力抑制方法 |
CN111752153A (zh) * | 2020-06-24 | 2020-10-09 | 北京航空航天大学 | 一种基于1.5阶混合重复控制器的谐波电流抑制方法 |
CN112467744A (zh) * | 2020-12-11 | 2021-03-09 | 东北电力大学 | 面向配网频率偏移的apf抗频扰谐波指令电流预测方法 |
CN112945165A (zh) * | 2021-03-23 | 2021-06-11 | 华驰动能(北京)科技有限公司 | 一种基于谐波小波的磁悬浮转子动不平衡位移检测方法 |
CN114442485A (zh) * | 2022-01-12 | 2022-05-06 | 武汉理工大学 | 一种通用的分数阶多周期特征谐波重复控制器设计方法及控制器 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20010030474A1 (en) * | 2000-04-14 | 2001-10-18 | Tokin Corporation | Multi-functional vibration actuator capable of suppressing an unstable operation around a resonance frequency |
CN101638061A (zh) * | 2009-08-26 | 2010-02-03 | 上海磁浮交通发展有限公司 | 可回馈制动能量的磁浮列车牵引供电系统抑制网侧谐波的方法 |
JP2013115970A (ja) * | 2011-11-30 | 2013-06-10 | Panasonic Corp | 整流子モータおよびこれを備えた電気機器 |
CN104503237A (zh) * | 2014-12-15 | 2015-04-08 | 北京航空航天大学 | 一种基于傅里叶变换的磁悬浮转子谐波振动控制方法 |
CN104503238A (zh) * | 2014-12-15 | 2015-04-08 | 北京航空航天大学 | 一种基于自适应重复控制器的磁悬浮转子系统电流谐波抑制方法 |
CN104660139A (zh) * | 2014-11-12 | 2015-05-27 | 江苏大学 | 基于矩阵变换器的交流磁轴承电主轴运行控制系统及方法 |
CN105007060A (zh) * | 2015-08-07 | 2015-10-28 | 广西星宇智能电气有限公司 | 一种清除有源滤波器中非选定谐波的方法 |
CN105159342A (zh) * | 2015-09-06 | 2015-12-16 | 北京航空航天大学 | 一种基于并联相移滤波器的磁悬浮转子谐波电流抑制方法 |
-
2016
- 2016-12-06 CN CN201611113486.5A patent/CN106444390B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20010030474A1 (en) * | 2000-04-14 | 2001-10-18 | Tokin Corporation | Multi-functional vibration actuator capable of suppressing an unstable operation around a resonance frequency |
CN101638061A (zh) * | 2009-08-26 | 2010-02-03 | 上海磁浮交通发展有限公司 | 可回馈制动能量的磁浮列车牵引供电系统抑制网侧谐波的方法 |
JP2013115970A (ja) * | 2011-11-30 | 2013-06-10 | Panasonic Corp | 整流子モータおよびこれを備えた電気機器 |
CN104660139A (zh) * | 2014-11-12 | 2015-05-27 | 江苏大学 | 基于矩阵变换器的交流磁轴承电主轴运行控制系统及方法 |
CN104503237A (zh) * | 2014-12-15 | 2015-04-08 | 北京航空航天大学 | 一种基于傅里叶变换的磁悬浮转子谐波振动控制方法 |
CN104503238A (zh) * | 2014-12-15 | 2015-04-08 | 北京航空航天大学 | 一种基于自适应重复控制器的磁悬浮转子系统电流谐波抑制方法 |
CN105007060A (zh) * | 2015-08-07 | 2015-10-28 | 广西星宇智能电气有限公司 | 一种清除有源滤波器中非选定谐波的方法 |
CN105159342A (zh) * | 2015-09-06 | 2015-12-16 | 北京航空航天大学 | 一种基于并联相移滤波器的磁悬浮转子谐波电流抑制方法 |
Non-Patent Citations (3)
Title |
---|
JIQIANG TANG 等: "Suppression of vibration caused by residual unbalance of rotor for magnetically suspended flywheel", 《JOURNAL OF VIBRATION AND CONTROL》 * |
单鸿涛 等: "基于FIR 滤波器的逆变器重复控制技术研究", 《中国电工技术学会电力电子学会第十一届学术年会》 * |
吕世轩 等: "基于电流观测与重复控制的逆变器多环控制策略研究", 《电气传动》 * |
Cited By (24)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106873655A (zh) * | 2017-04-13 | 2017-06-20 | 北京航空航天大学 | 一种基于有限维重复控制的磁轴承系统多谐波振动抑制方法 |
CN106886152A (zh) * | 2017-04-21 | 2017-06-23 | 北京航空航天大学 | 一种基于二阶奇次重复控制器的磁悬浮转子奇次谐波电流抑制方法 |
CN106896726A (zh) * | 2017-04-21 | 2017-06-27 | 北京航空航天大学 | 一种基于并联式二阶重复控制dpsorc的数字控制器的磁悬浮转子谐波电流抑制方法 |
CN106896726B (zh) * | 2017-04-21 | 2020-06-19 | 北京航空航天大学 | 一种磁悬浮转子谐波电流抑制方法 |
CN107844052A (zh) * | 2017-10-14 | 2018-03-27 | 河南工业大学 | 基于并联式有限维重复控制的磁轴承系统谐波电流抑制方法 |
CN107844052B (zh) * | 2017-10-14 | 2020-11-10 | 河南工业大学 | 基于并联式有限维重复控制的磁轴承系统谐波电流抑制方法 |
CN107966906B (zh) * | 2017-11-13 | 2020-04-21 | 南京航空航天大学 | 基于采样控制分离原理的分数阶延时实现方法 |
CN107966906A (zh) * | 2017-11-13 | 2018-04-27 | 南京航空航天大学 | 基于采样控制分离原理的分数阶延时实现方法 |
CN107870568A (zh) * | 2017-12-21 | 2018-04-03 | 北京航空航天大学 | 一种基于二阶双模重复控制的磁悬浮转子谐波电流抑制方法 |
CN108227499A (zh) * | 2018-01-16 | 2018-06-29 | 北京航空航天大学 | 一种谐波振荡器并联分数阶重复控制器的磁悬浮转子谐波电流抑制方法 |
CN108490777B (zh) * | 2018-03-19 | 2020-06-19 | 北京航空航天大学 | 一种基于改进奇次重复控制的磁悬浮转子谐波振动力抑制方法 |
CN108535541A (zh) * | 2018-03-19 | 2018-09-14 | 湖北工业大学 | 一种带有数字工频滤波的相位比较装置及方法 |
CN108490777A (zh) * | 2018-03-19 | 2018-09-04 | 北京航空航天大学 | 一种基于改进奇次重复控制的磁悬浮转子谐波振动力抑制方法 |
CN109164700B (zh) * | 2018-08-07 | 2021-06-18 | 中国科学院光电技术研究所 | 一种抑制全频扰动的控制方法 |
CN109164700A (zh) * | 2018-08-07 | 2019-01-08 | 中国科学院光电技术研究所 | 一种抑制全频扰动的控制方法 |
CN111473049A (zh) * | 2020-04-17 | 2020-07-31 | 河海大学 | 一种实心定子磁悬浮励磁电流的控制方法 |
CN111473049B (zh) * | 2020-04-17 | 2021-08-20 | 河海大学 | 一种实心定子磁悬浮励磁电流的控制方法 |
CN111708278B (zh) * | 2020-06-24 | 2021-07-13 | 北京航空航天大学 | 一种基于最小均方差前馈补偿算法的同频振动力抑制方法 |
CN111752153A (zh) * | 2020-06-24 | 2020-10-09 | 北京航空航天大学 | 一种基于1.5阶混合重复控制器的谐波电流抑制方法 |
CN111708278A (zh) * | 2020-06-24 | 2020-09-25 | 北京航空航天大学 | 一种基于最小均方差前馈补偿算法的同频振动力抑制方法 |
CN112467744A (zh) * | 2020-12-11 | 2021-03-09 | 东北电力大学 | 面向配网频率偏移的apf抗频扰谐波指令电流预测方法 |
CN112945165A (zh) * | 2021-03-23 | 2021-06-11 | 华驰动能(北京)科技有限公司 | 一种基于谐波小波的磁悬浮转子动不平衡位移检测方法 |
CN112945165B (zh) * | 2021-03-23 | 2022-09-23 | 华驰动能(北京)科技有限公司 | 一种基于谐波小波的磁悬浮转子动不平衡位移检测方法 |
CN114442485A (zh) * | 2022-01-12 | 2022-05-06 | 武汉理工大学 | 一种通用的分数阶多周期特征谐波重复控制器设计方法及控制器 |
Also Published As
Publication number | Publication date |
---|---|
CN106444390B (zh) | 2019-03-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106444390A (zh) | 一种基于fir滤波器和分数阶重复控制器的磁悬浮转子谐波电流抑制方法 | |
CN104503238B (zh) | 一种基于自适应重复控制器的磁悬浮转子系统电流谐波抑制方法 | |
CN106773673A (zh) | 一种基于频率自适应的分数补偿重复控制器的磁悬浮转子谐波电流抑制方法 | |
CN106647843B (zh) | 一种基于复合分数重复控制器的磁悬浮转子谐波电流抑制方法 | |
CN108227499B (zh) | 一种谐波振荡器并联分数阶重复控制器的磁悬浮转子谐波电流抑制方法 | |
CN106886152B (zh) | 一种基于二阶奇次重复控制器的磁悬浮转子奇次谐波电流抑制方法 | |
CN104503237B (zh) | 一种基于傅里叶变换的磁悬浮转子谐波振动控制方法 | |
CN105159342B (zh) | 一种基于并联相移滤波器的磁悬浮转子谐波电流抑制方法 | |
CN106330042A (zh) | 一种永磁同步电机谐波电流抑制方法及装置 | |
CN108897218B (zh) | 一种基于混合奇次重复控制器的磁悬浮转子奇次谐波振动抑制方法 | |
CN108490777A (zh) | 一种基于改进奇次重复控制的磁悬浮转子谐波振动力抑制方法 | |
CN112432634A (zh) | 一种基于多同步旋转坐标变换的谐波振动力抑制方法 | |
CN107276120A (zh) | 一种双馈风机同步机系统以及抑制同步发电机振荡的方法 | |
CN107591809A (zh) | 多周期并联重复控制谐波和间谐波指定次补偿方法 | |
CN106873655A (zh) | 一种基于有限维重复控制的磁轴承系统多谐波振动抑制方法 | |
CN101795105B (zh) | 无轴承永磁同步电机悬浮转子等效扰动电流补偿控制装置 | |
CN107844052A (zh) | 基于并联式有限维重复控制的磁轴承系统谐波电流抑制方法 | |
CN106289208B (zh) | 一种基于非线性自适应算法的磁轴承系统惯性轴辨识方法 | |
Cui et al. | Harmonic vibration control of MSCMG based on multisynchronous rotating frame transformation | |
CN105790269B (zh) | 通过双馈风机定子侧变流器抑制次同步谐振的方法及装置 | |
CN111752153A (zh) | 一种基于1.5阶混合重复控制器的谐波电流抑制方法 | |
CN113741181A (zh) | 一种转速自适应的磁悬浮转子系统奇次谐波电流抑制方法 | |
CN109062274A (zh) | 一种基于复变量有限维重复控制的磁轴承振动力矩抑制方法 | |
Cui et al. | Harmonic current suppression of magnetically suspended rotor system via odd-harmonic fractional RC | |
Peng et al. | Parallel mode notch filters for vibration control of magnetically suspended flywheel in the full speed range |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |