专利名称：正弦振动台失真度控制方法
技术领域：
本发明涉及一种失真度控制方法，具体地说是一种针对正弦振动台激励信号进行
闭环反馈控制的方法，应用于测量与控制领域。
背景技术：
由于振动台自身的结构设计原理，其支撑弹簧，气隙间磁场强度等非线性因素、振 动台的安装环境及负载等因素都会对振动台所产生的激励信号产生影响。振动台开环工作 时，尤其是在低频及超低频大振幅激励下，其加速度波形失真度将达到10 % ，甚至更高，很 难满足技术指标要求。因此标准振动台能否提供一个谐波失真最小的正弦振动激励信号， 直接关系到振动信号的测量与校准精度，这就使振动台失真度控制方法的研究显得尤为重 要。 现有振动台失真度控制方法主要有PID负反馈控制方法和重复控制方法。其中， PID负反馈控制方法是按偏差的比例、积分和微分进行控制的调节方法，该方法的实现需要 分析振动台波形失真产生机理，测量相关参数，建立振动台数学模型。当环境发生改变时， 振动台的实际模型会发生变化，控制性能会不同程度地下降，尤其对于建模误差较大时，会 出现很大偏差。重复控制方法基于内模原理，加到被控对象的输入信号除了新的偏差外，还 叠加了一个过去的偏差信号，该方法要通过多个周期的运行之后，才能提高系统的跟踪精 度，这对于低频或超低频信号，其反馈控制时间较长，影响测量速度。
发明内容
本发明的目的为解决现有技术振动台失真度控制方法建模难、误差较大、所需时 间较长的问题，本发明提供了一种无需对振动台建模、且可实现快速控制的正弦振动台失 真度闭环控制方法。
本发明的技术解决方案一种基于谐波陷波原理的正弦振动台失真度控制方法，
谐波陷波法控制方法是通过分析振动台输出信号的谐波成份，将其主要谐波成份负反馈回
振动台的驱动信号源，进行谐波补偿与抑制，达到降低振动台输出波形失真度的目的。
所述正弦振动台失真度控制方法包括如下步骤 步骤1由D/A产生正弦信号，经功率放大器驱动振动台振动； 步骤2计算振动台振动基频相位延迟； 步骤3对振动台的振动信号进行谐波分析； 步骤4采用谐波陷波算法，并通过谐波陷波补偿公式计算信号反馈量； 步骤5将反馈量与振动台输入信号相加，产生新的输入信号经功率放大器驱动振
动台振动； 步骤6重复步骤3至步骤5，直至将振动台的失真度控制在设定范围内。 所述振动信号为位移信号时，所述谐波陷波补偿公式为，
A (0 = - ， si一, +1 -(眠_ &))。 、
所述振动信号为速度信号时，所述谐波陷波补偿公式为，
A (0 = _ ， sin(枷,+ p。v4 - - & ) + (" 1) x 90。)。
所述振动信号为加速度信号时，所述谐波陷波补偿公式为，
A (0 = -，si一, + ^ —(眠—。+ (A: — 1) x 180。)。 所述步骤2进一步包括同步采集振动激光多普勒信号与振动台输入信号，对激光 干涉多普勒信号进行相位展开，获得振动信号，将测得的振动信号与振动台输入信号进行 相位比较，从而计算出振动台振动基频相位延迟。 所述步骤4中谐波陷波算法计算信号反馈量的过程为首先判断振动信号主要谐 波分量的谐波次数，其次通过查找表获得该谐波频率的振动台相位延迟，然后，采用谐波陷 波补偿公式计算反馈信号的相位值，从而确定反馈信号量。 本发明的原理将功率放大器、振动台、激光干涉仪与解调作为一个整体系统H来 对于该系统H，假设输入x(t) =Asin("t)时，其位移解调输出信号如式(1)所示。
分析。
相位。
= sin(W) + Jo2 sin(2W + p。2) +…
式中，A。" A。2……A。n为谐波分析得到各阶次信号幅值，^。2……P。"为各阶次信号 考虑到振动台输出与输入之间存在延迟信号，则可表示为如式(2)所示。
W)=」0l si咖/ - 。 +」。2 sin(2欲+ po2 - 2^ )
+ +」0" sin(""/ + p0 — ) + C70 (2)
9d为振动基频相位延迟。
为了消除y(t)的k次谐波，则在新的输入信号中加入反馈信号 ^^sinO^ + ^)，此时输出为如式(3)所示。
/ (0 =《sin一 - & ) +」o2 sin(2W + po2 _ 2^ ) +… +」OT sin("6^ + 一 ) 一sin(^^ +— ) _ <2 sin(2b, + 乂2 一 ) _… (3) _《sin(由/ + - O + C,
由式(3)可知，若要消除谐波，必须使基频的谐波响应与谐波的基频响应相等。因
此，可以得到反馈信号相位如式(4)所示。
%=%广购-。
故用位移信号反馈的k次谐波反馈量计算公式如式(5)所示。
<formula>formula see original document page 5</formula> 已知位移信号，可以通过一次求导的方式求得速度，二次求导的方式求得加速度。
则其速度与加速度信号可表示为
凡(O =」。vl sin(纽)+ 4v2 sin(2纽+ pov2)
+ + 0—, + ^,) + cv (6)
W)=《i sin(") + 42 sin(2w + %。2)
由于位移信号相位P。"与速度信号相位P,及加速度信号相位P。a"之间存在如下关
系
所示。
L = L - (m -1) x 90。 (8) P00 =P0"—("一l)xl80。 (9)
因此，用速度、加速度信号反馈的k次谐波反馈量计算公式分别如式ao)、式(ii)
A (0 = - , sin— + ^ - (d) + (A — 1) x 90。) (10)
A (0 = _ ， sin一 +1 —(眠—& ) + (A: -1) x 180°) (11) [oo4o] 式中c"(;为一常数。 本发明的有益效果本发明基于谐波陷波原理的正弦振动台失真度控制方法通过 分析振动台振动过程中的基频相位延迟和振动加速度谐波，利用谐波陷波算法计算反馈控 制数字信号，将该反馈信号与振动台的输入数字信号相加作为振动台新的输入信号，从而 不需要分析振动台的数学模型及各种影响参量，即可以快速、有效地补偿并抑制加速度波 形谐波失真度。该方法实现简单，可移植性强，是非常实用的振动台失真度控制方法，尤其
在低频或超低频测量与控制应用效果明显。
图1是本发明本发明正弦振动台失真度控制方法第一实施方式的原理框图；
图2是本发明本发明正弦振动台失真度控制方法第一实施方式中激光干涉仪的
5双路正交多普勒时域信号； 图3是本发明本发明正弦振动台失真度控制方法第三实施方式有无反馈时振动 台台面加速度时域信号的比较示意图， 图3a是图3中不加反馈时振动台台面加速度时域信号；
图3b是图3中加反馈后的振动台台面加速度时域信号； 其中，1-D/A电路、2-功率放大器、3-振动台、4-激光干涉仪、5-数据采集器、6-控 制算法模块、7_激光多普勒信号相位展开、8_谐波计算、9_反馈量计算、10-加法器。
具体实施例方式下面通过具体实施方式
对本发明作进一步的说明 请参阅图1，其是本发明正弦振动台失真度控制方法第一实施方式的原理框图。本 实施方式中，本发明正弦振动台失真度控制方法所对应的控制系统由D/A电路1、功率放大 器2、振动台3、激光干涉仪4、数据采集器5、控制算法模块6与加法器10组成。其中，所述 D/A电路1 一端与加法器10，另一端分别与功率放大器2和数据采集器5相接，所述功率放 大器2、振动台3、激光干涉仪4以及数据采集器5依次相接，所述数据采集器5经控制算法 模块6后与加法器10相连。而所述控制算法模块6包括顺次相连的激光多普勒信号相位 展开7、谐波计算8与谐波陷波法反馈量计算9。其中，所述谐波计算8采用正弦波三参量 最小二乘拟合算法拟合出谐波信息，所述谐波陷波法控制是通过分析振动台输出信号的谐 波成份，将其主要谐波成份负反馈回振动台的驱动信号源，进行谐波补偿与抑制，达到降低 振动台输出波形失真度的目的。 本实施方式中，所述功率放大器2、振动台3、激光干涉仪4与解调作为一个整体系
统H来分析，其输入为x(t) 二Asin("t)，其位移解调输出信号如式(1)所示。
= 4 sin(的+ 4>2 sin(2w + + ... + On— + L) + C。 (1) 式中，A", A。2……A。n为谐波分析得到各阶次信号幅值，^。2P。"为各阶次信号 相位。 考虑到振动台3输出与输入之间存在延迟信号，则可表示为如式(2)所示。
= sin(W D + 4>2 sin(2欣+ %2 - ) + " + j0 sin("^rf + 一 "<9rf ) + C。 (2) e d为振动基频相位延迟。 为了消除y(t)的k次谐波，本实施方式中，在新的输入信号中加入反馈信号
= —J^inO^ + ^)，此时振动台3的输出为如式(3)所示。<formula>formula see original document page 7</formula> ... (3) 由式(3)可知，若要消除谐波，必须使基频的谐波响应与谐波的基频响应相等，因
此％ +《d =P。A —，从而可以得到反馈信号相位如式(4)所示。
<formula>formula see original document page 7</formula> (4) 故用位移信号反馈的k次谐波反馈量计算公式如式(5)所示。
<formula>formula see original document page 7</formula> (5)、 式中Cs为一常数。然后将式(5)所得到的正弦谐波反馈量重新经D/A电路后负
反馈回到振动台的驱动信号源，以消除谐波，从而降低振动台输出波形失真度。 其具体实现步骤如下 步骤1 :将初始正弦信号源数字序列由D/A电路1转换为模拟信号，再经功率放大器2放大后驱动振动台3振动，再经激光干涉仪4干涉产生双路正交多普勒信号，其中，所述双路正交多普勒时域信号由图2所示； 步骤2 :信号的相位展开，同步采集激光干涉仪4的双路正交多普勒信号与D/A电路1输出信号，对激光干涉多普勒信号进行相位展开7，获得振动位移信号；将该振动信号与信号源信号进行相位比较，计算振动台基频相位延迟e d ; 步骤3 :谐波计算，用正弦波三参数最小二乘拟合算法对振动台的振动位移信号进行拟合，先拟合出基波成份的幅值相位信息，再将原始数据去除基波成份，拟合出2次谐波的幅值相位信息，依此类推，拟合出n次谐波信息，从而获得振动位移信号的数学表达式如式(1)所示; 步骤4:计算谐波反馈量，找出谐波分量中幅值最大的两个谐波，其中最大谐波分量为k次，第二大谐波分量为j次，则根据式(5)可计算出为消除k次和j次谐波所需要的反馈信号量，两次谐波的反馈信号量之和为总反馈量； 步骤5，将总反馈量叠加到信号源，产生新的信号源驱动振动台振动； 步骤6，重复步骤3 步骤5，直至将振动台的失真度控制在设定范围内。 本实施方式中，通过对振动台的振动位移信号进行拟合，得到其谐波信息，在计算
得到谐波反馈量，并将所计算得到的谐波反馈量经加法器重新输入，从而抑制或消除谐波
的影响，降低失真度。所述正弦振动台失真度控制方法不需要分析振动台的数学模型及各
种影响参量，即可以快速、有效地补偿并抑制加速度波形谐波失真度，具有实现简单，可移
植性强的特点，具有较大实际应用价值。 本发明正弦振动台失真度控制方法第二实施方式与第一实施方式相似，只是其所采取的谐波反馈量是速度信号，由于位移信号已知，因此通过一次求导就可以得到速度。
此时，与第一实施方式的位移信号相似，本实施方式中，速度输出可以表示为
凡(O = sin(W) +」。v2 sin(2W + pov2)
+... + sin(MC^ + p續)+ Gv (6) 且由于位移信号相位^。"与速度信号相位P,之间的关系可以式(8)表示，
L 二 L _ (" -1) x 90。 (8) 因此速度信号反馈的k次谐波反馈量具有式(10)的形式，
0V (0 = — , sin(b, + %v「 — &) + — x 90。)(10)[oo78] 式中c;为一常数。因此本实施方式中，通过对振动位移信号进行求导，计算其速
度信号，用最小二乘法计算速度信号基波与谐波信息，获得其数学表达式如式(6)所示，再根据式(10)计算出反馈信号量进行反馈，使振动台的失真度控制在设定范围内。
本发明正弦振动台失真度控制方法第三实施方式与第一实施方式相似，只是其所采取的谐波反馈量是加速度信号，由于位移信号已知，因此通过二次求导就可以得到加速度。 此时，与第一实施方式的位移信号相似，本实施方式中，所求的加速度输出可以表
示为
凡(O = 4。! sin(W) +」。a2 sin(2奴+ poa2)
+... + ^謹sin— + p麵)+ C。 (7) 且由于位移信号相位P。"与加速度信号相位^^之间的关系可以式(9)表示，
P鍾^^训-0-l)xl8(T (9) 因此，用加速度信号反馈的k次谐波反馈量计算公式，如式(11)所示
A (0 = - ，sin一 + ^ — — 。 + (" 1) x i,)(1 i) 式中c;为一常数。 本实施方式中，通过对振动位移信号进行二次求导，计算其加速度信号，用最小二乘法计算加速度信号基波与谐波信息，获得其数学表达式如式(7)所示，再根据式(11)计算出反馈信号量进行反馈，使振动台的失真度控制在设定范围内。 请同时参阅图3，其是本实施方式有无反馈时振动台台面加速度时域信号的比较示意图，其中，图3a是图3中不加反馈时振动台台面加速度时域信号；图3b是图3中加反馈后的振动台台面加速度时域信号。 综上所述本发明基于谐波陷波原理的正弦振动台失真度控制方法，通过对振动位移信号进行求导，计算其速度信号或加速度信号，用最小二乘法计算速度信号或加速度信号基波与谐波信息，再计算出反馈信号量并通过加法器作为新的输入进行反馈，使振动台的失真度控制在设定范围内。本发明基于谐波陷波原理的正弦振动台失真度控制方法实现
8C 简单，可移植性强，可以快速、有效地补偿并抑制加速度波形谐波失真度，具有较大的实际应用价值，特别是在低频或超低频测量与控制应用效果更为明显。
权利要求
一种正弦振动台失真度控制方法，其特征在于，包括如下步骤步骤1由D/A产生正弦信号，经功率放大器驱动振动台振动；步骤2计算振动台振动基频相位延迟；步骤3对振动台的振动信号进行谐波分析；步骤4采用谐波陷波算法，并通过谐波陷波补偿公式计算信号反馈量；步骤5将反馈量与振动台输入信号相加，产生新的输入信号经功率放大器驱动振动台振动；步骤6重复步骤3至步骤5，直至将振动台的失真度控制在设定范围内。
2. 根据权利要求1所述的正弦振动台失真度控制方法，其特征在于所述振动信号为 位移信号时，所述谐波陷波补偿公式为，<formula>formula see original document page 2</formula>
3. 根据权利要求1所述的正弦振动台失真度控制方法，其特征在于所述振动信号为 速度信号时，所述谐波陷波补偿公式为，<formula>formula see original document page 2</formula>
4. 根据权利要求1所述的正弦振动台失真度控制方法，其特征在于所述振动信号为 加速度信号时，所述谐波陷波补偿公式为，<formula>formula see original document page 2</formula>
5. 根据权利要求1至4任一项所述的正弦振动台失真度控制方法，其特征在于所述步骤2进一步包括同步采集振动激光多普勒信号与振动台输入信号，对激光干涉多普勒信 号进行相位展开，获得振动信号，将测得的振动信号与振动台输入信号进行相位比较，从而 计算出振动台振动基频相位延迟。
6. 根据权利要求5所述的正弦振动台失真度控制方法，其特征在于所述步骤4中谐波陷波算法计算信号反馈量的过程为首先判断振动信号主要谐波分量的谐波次数，其次 通过查找表获得该谐波频率的振动台相位延迟，然后，采用谐波陷波补偿公式计算反馈信 号的相位值，从而确定反馈信号量。
全文摘要
本发明涉及一种失真度控制方法，具体地说是一种针对正弦振动台激励信号进行闭环反馈控制的方法，应用于测量与控制领域。本发明正弦振动台失真度控制方法通过分析振动台振动过程中的基频相位延迟和振动加速度谐波，利用谐波陷波算法计算反馈控制数字信号，将该反馈信号与振动台的输入数字信号相加作为振动台新的输入信号，并通过D/A转换后经功率放大器作用于振动台。本发明无需对振动台建模，控制方法简单可行，能快速抑制振动台的失真度，降低非线性因素的影响，提高波形精度，具有控制速度快、可移植性强的特点。
文档编号H03F1/32GK101741318SQ20091026617
公开日2010年6月16日 申请日期2009年12月31日 优先权日2009年12月31日
发明者李新良, 胡春艳, 陈璐 申请人:中国航空工业集团公司北京长城计量测试技术研究所