CN101047368A - 一种章动频率自动跟踪的高通数字滤波方法 - Google Patents
一种章动频率自动跟踪的高通数字滤波方法 Download PDFInfo
- Publication number
- CN101047368A CN101047368A CN 200710064318 CN200710064318A CN101047368A CN 101047368 A CN101047368 A CN 101047368A CN 200710064318 CN200710064318 CN 200710064318 CN 200710064318 A CN200710064318 A CN 200710064318A CN 101047368 A CN101047368 A CN 101047368A
- Authority
- CN
- China
- Prior art keywords
- nutation frequency
- frequency
- rotor
- nutation
- motion tracking
- 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
Images
Abstract
一种章动频率自动跟踪的高通数字滤波方法,该方法主要包括Kφφ的测试、转子转速检测、计算章动频率、计算数字滤波器系数、数字滤波五个步骤;在不加交叉解耦控制时,测试Kφφ;转子的转速Ω由霍尔传感器检测,经调理电路调理后由微处理器获得;将测得的转速Ω代入章动频率计算公式在线计算章动频率fn;以fn为截止频率设计数字滤波器;最后对输入的转子位移信号进行在线滤波计算。该发明主要应用于磁悬浮控制力矩陀螺或磁悬浮飞轮系统中磁轴承的交叉解耦控制。应用该发明可以实时跟踪磁轴承的章动频率,并保持对章动频率的良好的滤波性能,有利于实现交叉耦合控制中对章动频率的有效抑制。
Description
技术领域
本发明涉及一种频率自跟踪数字滤波方法,特别是一种章动频率自动跟踪的高通数字滤波方法,可用于磁悬浮控制力矩陀螺或磁悬浮飞轮系统中磁轴承的交叉解耦控制。
背景技术
磁悬浮交叉解耦控制技术是一种解决磁轴承的稳定悬浮问题的优选方案。在这种方案中首先要把进动信号与章动信号分离出来,然后分别对进动和章动进行控制。进动信号的频率随转子转速的升高降低,章动信号的频率随转子转速的升高而升高,两者随转子转速的变化曲线没有交叉点。可以采用截止频率固定不变的交叉高通滤波器或多段交叉高通滤波器将章动信号分离出来。当采用截止频率固定不变的交叉高通滤波器时,随转速升高,章动频率升高,这种交叉高通滤波器的相位滞后加大,从而导致转速不高时磁轴承章动失稳。目前多采用多段交叉高通滤波求取频率较高的章动信号。这种方法随转子转速的升高不断提高交叉高通滤波器的截止频率,以克服截止频率固定不变的交叉高通滤波器在转速升高时因相位滞后过多而导致磁转子章动失稳的问题。采用分段高通滤波器方式主要有两个缺点:一是调试过程比较繁琐,需要转子反复升速以确定多段高通滤波器之间的切换频率;二是在相邻两段高通滤波器进行切换时,由于其截止频率是突然改变的,滤波后的章动信号的相位和幅值在切换前后也是突然改变的,这会给磁悬浮轴承系统带来冲击,增加了导致磁悬浮轴承失稳的不安全因素。
发明内容
本发明的技术解决问题是:提出一种磁悬浮轴承交叉解耦控制中自动跟踪章动频率,并实时改变截止频率的交叉高通滤波器的设计方法,以克服多段交叉高通滤波器对磁悬浮轴承系统的冲击,实现章动频率的平滑抑制。
本发明的技术解决方案:本发明主要有五个步骤组成,即K的测试、转子转速的检测、章动频率的跟踪计算、数字滤波器设计、数字滤波。
(1)K的测试。在不加交叉解耦控制时,将转子转速升到Ω1,测得此时的章动频率fn1,按公式
K=2πfn1[2πJdfn1-JzΩ1] (1)
求得K。其中:K是转子对广义坐标的广义角刚度,Jz、Jd分别是转子极转动惯量和赤道转动惯量。
(2)进行转子转速Ω的检测。转子的转速通过霍尔传感器检测,经调理电路进行滤波、整形、隔离后将其变为频率与转速成正比的矩形波,由微处理器得到转子的转速Ω。
(3)计算章动频率。章动频率fn的计算公式为:
(4)计算数字滤波器系数。以章动频率fn为截止频率先设计模拟高通滤波器,然后采用离散化方法将模拟高通滤波器离散化得到相应的数字滤波器系数。
(5)对输入的转子位移进行在线数字滤波计算。
先对输入的转子位移信号求差值,再对差值进行数字滤波,数字滤波器的计算公式为:
c_out(k)=b0×c_in(k)+Λ+bn×c_in(k-n)-[a1×c_out(k-1)+am×c_out(k-m)] (3)其中b0,Λ,bn,a1,Λam是数字滤波器的系数,c_in(k)是转子位移信号的第k次采样差值,c_out(k)是第k次滤波输出值。
本发明的原理是:在不加交叉解耦控制时,将转子转速升到Ω1,测得此时的章动频率fn1,按公式(1)求得K;转子的转速由霍尔传感器检测,经调理电路调理后送到微处理器得到转子的转速Ω;根据公式(2)在线计算章动频率fn,这样就实现了章动频率的自动跟踪;以fn为截止频率设计模拟高通滤波器,并采用合适的离散化方法将模拟高通滤波器离散化就得到相应的数字滤波器,然后对输入的位移信号进行数字滤波运算。根据滤波器在截止频率处的幅值和相位不随截止频率改变而改变的特点,可以保证该数字滤波器在章动频率fn处的幅值和相位不随fn的改变而改变。
本发明的方案与现有方案相比,主要优点在于:(1)高通截止频率可以实时跟踪章动频率的变化而变化,并能保持良好的滤波性能;(2)可以保证在章动频率变化时,高通滤波器在章动频率处的幅值和相位不变;(3)高通滤波器截止频率的改变是连续的,不会给系统带来切换冲击;(4)调试过程比较简单。
附图说明
图1为本发明的结构框图;
图2为本发明的流程图;
图3为磁轴承与转子示意图;
图4为本发明的交叉高通滤波器在交叉解耦控制中的应用。
具体实施方式
本发明主要包括五个步骤:转子转速检测、K的测试、章动频率计算、数字滤波器系数计算、数字滤波。如图1所示,在硬件上主要包括霍尔传感器4、调理电路5、DSP1共三个部分,转速检测6、章动频率计算2、数字滤器系数计算7、交叉高通滤波3是DSP1的内部算法。转子的转速由霍尔传感器4检测,经调理电路5调理后送到DSP1,由DSP1计算得到转子的转速Ω;将测得的转速Ω代入章动频率计算公式(1)在线计算章动频率fn,以实现章动频率的自动跟踪;以fn为截止频率计算得到数字滤波器的系数,最后对输入的转子位移信号进行在线滤波计算。
(1)K在实际系统中取决于系统参数和转速,且难以得出其表达式,但对于一个给定的系统,在控制参数不变时,其变化不大,可以用一个恒定的值来表示。对于同一个磁轴承控制系统在不加交叉解耦控制时,将转子转速升到Ω1,用示波器的FFT功能测得此时的章动频率fn1,代入公式(1)中可以解得K。
(2)转子的转速通过霍尔传感器4检测,经调理电路5进行滤波、整形、隔离将其变为频率与转速成正比的DSP(也可以采用微控制器ARM)可接收的矩形波,由DSP捕获,通过测周期或测频率的方法得到转子的转速Ω。
(3)按公式(2)在线计算章动频率fn。
(4)以fn为截止频率设计模拟高通滤波器,并采用合适的离散化方法(如离散化方法是向后差分法,或向前差分法、或双线性变换法,或预修正双线性变换法等)将模拟高通滤波器离散化得到相应的数字滤波器。模拟高通滤波器是一阶的,也可以是二阶的。本实施方案中,使用的是一阶模拟高通滤波器,其传递函数为
其中T是时间常数。取截止频率f0=fn,则
在f0变化时,可以保证在fn处的幅值
(5)对输入的转子位移信号进行高通数字滤波。如图3所示,转子8由两端径向磁轴承悬浮,一端为A端磁轴承9,另一端为B端磁轴承10。图4中的ax、ay、bx、by分别是A端磁轴承9在x轴、y轴及B端磁轴承10在x轴、y轴上的位移,KxL、KxH、KyL、KyH分别是x轴上位移的低通滤波系数、高通滤波系数及y轴上位移的低通滤波系数、高通滤波系数。对于x轴上的位移信号有:cx=ax-bx,其滤波输出为 其中:cx(k)、cx(k-1)是cx在当前采样周期和上一采样周期的采样值,cx_out(k)、cx_out(k-1)是当前采样周期和上一采样周期的滤波输出,并且T=2πfn。对于y轴上的位移信号有:cy=ay-by,其滤波输出为
其中:cy(k)、cy(k-1)是cy在当前采样周期和上一采样周期的采样值,cy_out(k)、cy_out(k-1)是当前采样周期和上一采样周期的滤波输出。
本发明说明书中未作详细描述的内容属于本领域专业技术人员公知的现有技术。
Claims (6)
1、一种章动频率自动跟踪的高通数字滤波方法,其特征在于包括以下步骤:
(1)进行K的测试,在不加交叉解耦控制时,将转子转速升到Ω1,测得此时的章动频率fn1,按公式:
K=2πfn1[2πJdfn1-JzΩ1]
求得K,其中:K是转子对广义坐标的广义角刚度,Jz、Jd分别是转子极转动惯量和赤道转动惯量;
(2)进行转子转速Ω的检测;
(3)根据K和转子转速Ω计算章动频率,章动频率fn的计算公式为:
(4)计算数字滤波器系数,以章动频率fn为截止频率先计算模拟高通滤波器的时间常数,然后采用离散化方法计算得到相应的数字滤波器系数;
(5)对输入的转子位移进行在线数字滤波计算。
2、根据权利要求1所述的章动频率自动跟踪的高通数字滤波方法,其特征在于:所述步骤(1)中Ω1的范围是376.8rad/s-3768rad/s。
3、根据权利要求1所述的章动频率自动跟踪的高通数字滤波方法,其特征在于:所述步骤(2)中的转子转速Ω的检测是:通过霍尔传感器检测,经调理电路进行滤波、整形、隔离后变为矩形波,由微处理器计算得到。
4、根据权利要求3所述的章动频率自动跟踪的高通数字滤波方法,其特征在于:所述的微处理器是DSP,或微控制器ARM。
5、根据权利要求1所述的章动频率自动跟踪的高通数字滤波方法,其特征在于:所述步骤(4)中的模拟高通滤波器为一阶,或二阶。
6、根据权利要求1所述的一种章动频率自动跟踪的高通数字滤波方法,其特征在于:所述步骤(4)中的离散化方法是向后差分法,或向前差分法、或双线性变换法,或预修正双线性变换法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2007100643186A CN100463367C (zh) | 2007-03-12 | 2007-03-12 | 一种章动频率自动跟踪的高通数字滤波方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2007100643186A CN100463367C (zh) | 2007-03-12 | 2007-03-12 | 一种章动频率自动跟踪的高通数字滤波方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101047368A true CN101047368A (zh) | 2007-10-03 |
CN100463367C CN100463367C (zh) | 2009-02-18 |
Family
ID=38771658
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNB2007100643186A Expired - Fee Related CN100463367C (zh) | 2007-03-12 | 2007-03-12 | 一种章动频率自动跟踪的高通数字滤波方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN100463367C (zh) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102424117A (zh) * | 2011-11-06 | 2012-04-25 | 北京航空航天大学 | 一种补偿磁悬浮控制力矩陀螺磁轴承相位滞后的方法 |
CN103472730A (zh) * | 2013-09-23 | 2013-12-25 | 北京航空航天大学 | 一种基于谐波减速器扭转刚度迟滞模型的双框架控制力矩陀螺高精度框架速率伺服系统 |
CN103543636A (zh) * | 2013-09-24 | 2014-01-29 | 北京航天控制仪器研究所 | 一种动力调谐陀螺仪解耦伺服控制回路 |
CN104568165A (zh) * | 2014-12-10 | 2015-04-29 | 中国航空工业集团公司北京长城计量测试技术研究所 | 通过自适应信号调理测量周期变化辐射温度的方法 |
CN104634363A (zh) * | 2014-11-27 | 2015-05-20 | 上海新跃仪表厂 | 陀螺章动频率的测试系统及其测试方法 |
CN106289776A (zh) * | 2016-07-29 | 2017-01-04 | 南京航空航天大学 | 一种柔性磁悬浮轴承转子刚度阻尼辨识方法 |
CN107992063A (zh) * | 2017-12-29 | 2018-05-04 | 哈尔滨工业大学 | 基于变参数章动阻尼的变速倾侧动量轮进动控制方法 |
CN110798101A (zh) * | 2019-11-12 | 2020-02-14 | 威科达(东莞)智能控制有限公司 | 一种永磁同步电机转矩控制下的抖动抑制方法 |
CN111256989A (zh) * | 2020-02-22 | 2020-06-09 | 南京凯奥思数据技术有限公司 | 一种变转速旋转设备振动特征值的采集方法 |
RU2766130C1 (ru) * | 2020-08-20 | 2022-02-08 | Федеральное государственное бюджетное учреждение науки Институт проблем управления им. В.А. Трапезникова Российской академии наук | Способ диагностики подшипников роторного оборудования на основе оценки микровариаций вращения вала |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
SE517525C2 (sv) * | 1999-09-07 | 2002-06-18 | Ericsson Telefon Ab L M | Förfarande och anordning för konstruktion av digitala filter |
US6976045B2 (en) * | 2001-08-08 | 2005-12-13 | Tektronix, Inc. | Variable sample rate recursive digital filter |
-
2007
- 2007-03-12 CN CNB2007100643186A patent/CN100463367C/zh not_active Expired - Fee Related
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102424117A (zh) * | 2011-11-06 | 2012-04-25 | 北京航空航天大学 | 一种补偿磁悬浮控制力矩陀螺磁轴承相位滞后的方法 |
CN103472730A (zh) * | 2013-09-23 | 2013-12-25 | 北京航空航天大学 | 一种基于谐波减速器扭转刚度迟滞模型的双框架控制力矩陀螺高精度框架速率伺服系统 |
CN103543636A (zh) * | 2013-09-24 | 2014-01-29 | 北京航天控制仪器研究所 | 一种动力调谐陀螺仪解耦伺服控制回路 |
CN103543636B (zh) * | 2013-09-24 | 2016-02-10 | 北京航天控制仪器研究所 | 一种动力调谐陀螺仪解耦伺服控制回路 |
CN104634363A (zh) * | 2014-11-27 | 2015-05-20 | 上海新跃仪表厂 | 陀螺章动频率的测试系统及其测试方法 |
CN104568165B (zh) * | 2014-12-10 | 2019-04-05 | 中国航空工业集团公司北京长城计量测试技术研究所 | 通过自适应信号调理测量周期变化辐射温度的方法 |
CN104568165A (zh) * | 2014-12-10 | 2015-04-29 | 中国航空工业集团公司北京长城计量测试技术研究所 | 通过自适应信号调理测量周期变化辐射温度的方法 |
CN106289776A (zh) * | 2016-07-29 | 2017-01-04 | 南京航空航天大学 | 一种柔性磁悬浮轴承转子刚度阻尼辨识方法 |
CN106289776B (zh) * | 2016-07-29 | 2020-07-07 | 南京航空航天大学 | 一种柔性磁悬浮轴承转子刚度阻尼辨识方法 |
CN107992063A (zh) * | 2017-12-29 | 2018-05-04 | 哈尔滨工业大学 | 基于变参数章动阻尼的变速倾侧动量轮进动控制方法 |
CN107992063B (zh) * | 2017-12-29 | 2020-08-04 | 哈尔滨工业大学 | 基于变参数章动阻尼的变速倾侧动量轮进动控制方法 |
CN110798101A (zh) * | 2019-11-12 | 2020-02-14 | 威科达(东莞)智能控制有限公司 | 一种永磁同步电机转矩控制下的抖动抑制方法 |
CN111256989A (zh) * | 2020-02-22 | 2020-06-09 | 南京凯奥思数据技术有限公司 | 一种变转速旋转设备振动特征值的采集方法 |
RU2766130C1 (ru) * | 2020-08-20 | 2022-02-08 | Федеральное государственное бюджетное учреждение науки Институт проблем управления им. В.А. Трапезникова Российской академии наук | Способ диагностики подшипников роторного оборудования на основе оценки микровариаций вращения вала |
Also Published As
Publication number | Publication date |
---|---|
CN100463367C (zh) | 2009-02-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101047368A (zh) | 一种章动频率自动跟踪的高通数字滤波方法 | |
CN102779238B (zh) | 一种基于自适应卡尔曼滤波的无刷直流电机系统辨识方法 | |
CN103269198B (zh) | 一种基于编码器自动调零的永磁同步电机控制方法及系统 | |
CN104601076B (zh) | 用于电动汽车电动机的无传感器滑模观测器设计方法 | |
CN201277910Y (zh) | 一种航空发动机试车参数测量系统 | |
CN101938246A (zh) | 无速度传感器电机转速的模糊融合辨识方法 | |
CN101917151B (zh) | 一种开关磁阻电机无位置传感器控制方法 | |
CN105634364B (zh) | 一种三相四开关变频调速系统中电容电压漂移的抑制方法 | |
CN104852662B (zh) | 永磁同步电机静态电感参数的测量方法及系统 | |
CN107769636B (zh) | 一种无位置传感器永磁同步电机转子位置检测方法 | |
CN110336507A (zh) | 共振频率的检测、机械共振的抑制方法、装置及介质 | |
CN103401243A (zh) | 一种指定次谐波补偿apf及其谐波检测和控制方法 | |
JP2021508825A (ja) | トルク発生装置の内部の実効トルクを推定する方法 | |
CN101972170A (zh) | 最小二乘支持向量机自适应滤波器及其滤波方法 | |
CN115479595A (zh) | 一种半球谐振陀螺分时复用参数辨识系统及方法 | |
CN104483502A (zh) | 一种基于单片机的电机实时转速宽范围精确测速方法 | |
CN102117374B (zh) | 一种基于扰动信号进行发电机组轴系扭振模态阻尼的计算方法 | |
CN103281031A (zh) | 基于可测电气量的鼠笼式异步电机等值电路参数辨识方法 | |
CN105759207B (zh) | 一种基于m-ekf算法的海流发电机绝缘系统故障检测方法 | |
CN105305898B (zh) | 一种开关磁阻电机无位置传感器的转子位置估计方法 | |
CN104579092B (zh) | 电机的控制方法、控制系统及电机电感的计算方法、装置 | |
CN103529697B (zh) | 一种确定轴向磁轴承pid参数鲁棒稳定域的方法 | |
CN107733313A (zh) | 旋转机械谐振自适应抑制方法 | |
CN112350635A (zh) | 高压大功率变频器无速度矢量控制磁链观测器及观测方法 | |
CN112968644B (zh) | 永磁同步电机参数在线识别方法、系统、终端及介质 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20090218 Termination date: 20180312 |