发明名称 一种电液加速度伺服系统随机振动控制方法
摘要 本发明公开了一种电液加速度伺服系统随机振动控制方法,包括以下步骤:利用加速度功率谱密度参考信号设计滤波器;利用白噪声发生器生成白噪声信号;利用滤波器对白噪声信号进行滤波,将滤波器的输出信号作为线性卷积的输入信号;计算电液加速度伺服系统的阻抗;进行逆傅里叶变换;利用重叠保留法计算两个信号的线性卷积;将线性卷积的输出信号作为电液加速度伺服系统和阻抗计算的输入信号。本发明的所有步骤可通过软件编程实现,具有更高的实时性。本发明通过G<sub>2</sub>(f)除以G<sub>1</sub>(f)直接在线计算电液加速度伺服系统的阻抗,不需要通过频率响应函数求取电液加速度伺服系统的阻抗。本发明采用重叠保留法实现实时在线补偿,提高了随机振动模拟实验的控制精度。
申请公布号 CN103399589B 申请公布日期 2015.05.13
申请号 CN201310331386.X 申请日期 2013.08.01
申请人 大连海事大学 发明人 关广丰;王海涛;熊伟;马文琦;孙长乐
分类号 G05D19/02(2006.01)I 主分类号 G05D19/02(2006.01)I
代理机构 大连东方专利代理有限责任公司 21212 代理人 李洪福
主权项 一种电液加速度伺服系统随机振动控制方法,其特征在于:包括以下步骤:A、利用加速度功率谱密度参考信号(1)设计滤波器(2):<img file="FDA0000663343900000011.GIF" wi="420" he="231" />式中,f为频率,|H(f)|为所设计滤波器(2)的幅频特性,R(f)为加速度功率谱密度参考信号(1),a为白噪声信号(5)的功率,<img file="FDA0000663343900000013.GIF" wi="124" he="76" />为所设计滤波器(2)的相频特性,P为滤波器(2)阶数;B、利用白噪声发生器生成白噪声信号(5);C、将白噪声信号(5)作为滤波器(2)的输入信号,利用滤波器(2)对白噪声信号(5)进行滤波,将滤波器(2)的输出信号作为线性卷积(3)的输入信号;D、将线性卷积(3)的输出信号和电液加速度伺服系统(4)的输出信号作为阻抗计算(6)的输入信号,计算电液加速度伺服系统(4)的阻抗<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>Z</mi><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><msub><mi>G</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow></mrow><mrow><msub><mi>G</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow></mrow></mfrac></mrow>]]></math><img file="FDA0000663343900000012.GIF" wi="318" he="145" /></maths>式中,Z(f)为电液加速度伺服系统(4)的阻抗,G<sub>2</sub>(f)为线性卷积(3)的输出信号的自功率谱密度的平均值,G<sub>1</sub>(f)为线性卷积(3)的输出信号和电液加速度伺服系统(4)的输出信号的互功率谱密度的平均值;E、利用阻抗计算(6)的输出信号,进行逆傅里叶变换(7),得到序列z(m):z(m)=IFFT[Z(f)]F、基于滤波器(2)的输出信号和逆傅里叶变换(7)的输出信号,利用重叠保留法计算两个信号的线性卷积(3);设序列z(m)的长度为M点,滤波器(2)输出信号为序列r(n),利用重叠保留法计算r(n)*z(m),式中,*表示序列r(n)与序列z(m)的线性卷积(3);F1、在序列r(n)前面添加M‑1个零点,形成序列r’(n):<img file="FDA0000663343900000021.GIF" wi="524" he="183" />F2、将r’(n)分割为N段序列{r'<sub>0</sub>(n),r'<sub>1</sub>(n),r'<sub>2</sub>(n),...,r'<sub>N‑1</sub>(n)},N由加速度功率谱密度再现实验的持续时间确定;序列r’(n)的分割方法用公式表示为:k=0时,<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><msup><mi>r</mi><mo>&prime;</mo></msup><mi>k</mi></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>=</mo><mfenced open='{' close=''><mtable><mtr><mtd><msup><mi>r</mi><mo>&prime;</mo></msup><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow></mtd><mtd><mn>0</mn><mo>&le;</mo><mi>n</mi><mo>&le;</mo><mi>L</mi><mo>+</mo><mi>M</mi><mo>-</mo><mn>2</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mi>n</mi><mo>></mo><mo>+</mo><mi>M</mi><mo>-</mo><mn>2</mn></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0000663343900000022.GIF" wi="775" he="158" /></maths>k&gt;0时,<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><msup><mi>r</mi><mo>&prime;</mo></msup><mi>k</mi></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>=</mo><mfenced open='{' close=''><mtable><mtr><mtd><msup><mi>r</mi><mo>&prime;</mo></msup><mrow><mo>(</mo><mi>n</mi><mo>+</mo><mi>kL</mi><mo>)</mo></mrow></mtd><mtd><mn>0</mn><mo>&le;</mo><mi>n</mi><mo>&le;</mo><mi>L</mi><mo>+</mo><mi>M</mi><mo>-</mo><mn>2</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mi>n</mi><mo>></mo><mo>+</mo><mi>M</mi><mo>-</mo><mn>2</mn></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0000663343900000023.GIF" wi="798" he="160" /></maths>式中,L是大于1的正整数;F3、在序列z(m)后面添加L‑1个零点,得到序列z'(m),即:<img file="FDA0000663343900000024.GIF" wi="576" he="189" />F4、设y'<sub>k</sub>(l)=r<sub>k</sub>(n)*z(m),利用傅里叶变换计算y'<sub>k</sub>(l):y'<sub>k</sub>(l)=IFFT[FFT(r'<sub>k</sub>(n))FFT(z'(m))]F5、计算r(n)*z(m):设y(i)=r(n)*z(m),则有:<img file="FDA0000663343900000025.GIF" wi="409" he="148" />式中,<img file="FDA0000663343900000026.GIF" wi="898" he="158" />G、将线性卷积(3)的输出信号作为电液加速度伺服系统(4)和阻抗计算(6)的输入信号。
地址 116026 辽宁省大连市高新园区凌海路1号