# Identification for Active Vibration Control of Flexible Structure Based on Prony Algorithm

Authors: Xianjun Sheng; Yuanli Kong; Fengyun Zhang; Rui Yang
Journal: Shock and Vibration
 (2016)

---

## Abstract

Flexible structures have been widely used in many fields due to the advantages of light quality, small damping, and strong flexibility. However, flexible structures exhibit the vibration in the process of manipulation, which reduces the pointing precision of the system and causes fatigue of the machine. So, this paper focuses on the identification method for active vibration control of flexible structure. The modal parameters and transfer function of the system are identified from the step response signal based on Prony algorithm, while the vibration is attenuated by using the input shaping technique designed according to the parameters identified from the Prony algorithm. Eventually, the proposed approach is applied to the most common flexible structure, a piezoelectric cantilever beam actuated by Macro Fiber Composite (MFC). The experimental results demonstrate that the Prony algorithm is very effective and accurate on the dynamic modeling of flexible structure and input shaper could significantly reduce the vibration and improve the response speed of system.


---

## Body

## 1. Introduction

Flexible structures experience undesirable vibration in response to environmental and operational forces. The very existence of vibrations not only can limit the accuracy of sensitive instruments or cause significant error in applications, where high positioning precision is essential, but also plays an important role in causing component fatigue and shortening the life of the system [1–3]. Therefore, vibration control of flexible structures has attracted a considerable amount of investigations in many scientific and engineering applications.The vibration of flexible structure is a quite complicated issue concerning many factors and the key technique is to obtain an accurate dynamic model of the system. Therefore, a variety of studies have been done on dynamic modeling for vibration control of flexible structures. An analytic modeling method based on Hamilton’s principles is presented in [4] to procure the dynamic equation of flexible beam and the independent modal method by using piezo-patches as actuators is utilized to realize the vibration suppression. However, when using the analytic modeling, the detailed physical information regarding the whole system is required and the descriptive partial differential equations must be solved, which makes it very difficult to obtain model parameters. It seems that finite element is a better approach compared with analytic model because there is no requirement to solve a partial differential equation [5–7]. Unfortunately, detailed information regarding the materials of structure and boundary conditions is still required. Many efforts were directed towards system identification which is a straightforward method to obtain a structural model from experimental data. An identification method on procuring transfer function from time domain data can be found in [8, 9] and literature [10] put forward a frequency identification approach by swept frequency analysis. And, recently, artificial intelligence techniques are used in system identification for flexible structure, such as genetic algorithm [11] and neural networks [12, 13]. The previous identification methods exhibit good performances on dynamic modeling for flexible structure, but they failed to identify both modal parameters and transfer function simultaneously. Therefore, in this paper, a new identification method based on Prony algorithm is proposed to establish the system dynamic model and combines the input shaping technique to achieve an active vibration control for flexible structure. Prony algorithm uses a linear combination of complex exponential functions to describe equally spaced sampling data and calculates the frequency, amplitude, damping factor, and initial phase of the oscillation system from the response signal. Prony algorithm has succeeded in estimating the modal parameters of power system oscillations and hydro turbine generating sets, referencing articles [14–17], but Prony algorithm has never been used in vibration analysis for flexible structures. So, this paper identifies the modal parameters and transfer function for flexible structure based on Prony algorithm. Then the identification results combine the input shaping techniques to achieve an active vibration control for a piezoelectric cantilever beam.The organization of this paper is as follows. Section2 presents a piezoelectric cantilever beam which is considered as flexible structure in detail and analyzes the system’s characteristic. Section 3 introduces Prony algorithm on identifying the modal parameters and transfer function for flexible structure. Section 4 gives details on control method and the design of input shaper based on identification results. Section 5 presents experimental results on the proposed identification and vibration control. Section 6 summarizes the paper.


## 2. Description and Analysis of the Flexible Structure Behavior

A piezoelectric cantilever beam actuated by Macro Fiber Composites (MFC) is used as flexible structure in this paper. The beam is cantilevered at one end and a MFC patch used as an actuator is bonded on the surface of aluminum layer (see Figure1). The MFC composed of resin and a piezoelectric ceramic with 0° fiber orientation which can offer higher performance, durability, and flexibility compared with piezoelectric ceramics. The MFC used in this paper is developed by NASA and its model is M-8507-P1. The workspace of MFC is 85 mm × 7 mm and free strain epsilon is 1380 ppm. The main characteristics of the piezoelectric cantilever beam are presented in Table 1.Table 1
The main parameters of piezoelectric cantilever beam.
Material
Young’s modulus (GPa)
Density (Kg/m3)
Poisson’s ratio
Length (mm)
Width (mm)
Thickness (mm)
Aluminum layer
69
2800
0.34
150
13
1
Bonding layer
2.5
2700
0.36
101
13
0.1
MFC
30
7450
0.31
101
13
0.3Figure 1
The structure of a piezoelectric cantilever actuated by MFC.The whole experimental framework is shown in Figure2. The computer with software LabVIEW is used as the controller and the data acquisition (DAQ) card whose model is NI USB- 6229 is used for digital-to-analogue (D/A) and analogue-to-digital (A/D) conversions. A voltage amplifier with a fixed gain 200 is used to provide driving voltage for the MFC. In order to detect the displacement of the piezocantilever and validate the experimental performances, a laser sensor from KEYENCE with an accuracy of 1 μm and a high passband is utilized.Figure 2
The schematic of experiment platform.When a step input voltage is applied, the oscillation emerges prominently at the free end of piezoelectric cantilever and attenuates gradually with the time. The existence of oscillation severely degrades the positioning accuracy and increases the settling time. In this paper, the dynamic modeling of piezoelectric cantilever is assumed as linear time-invariant, without consideration of hysteresis and creep. In this paper, the behavior of piezoelectric cantilever is modeled based on the Prony algorithm and the vibration is suppressed by input shaping technique. The whole schematic of identification for active vibration control is depicted in Figure3. The modeling and vibration control are performed step by step. These steps are as follows. Firstly, the transfer function and modal parameters of the system are identified through Prony algorithm by applying the excitation signal to the controlled object. Then the input shaper is designed based on the previously obtained parameters. At last, the reference input signal is convoluted with the impulse sequences gained from the input shaper to realize the vibration suppression.Figure 3
The schematic of identification for active vibration control.


## 3. Modal Parameters and Transfer Function Identification Based on Prony Algorithm

In this section, an identification method based on Prony algorithm is proposed for a flexible structure withn modes. After using the proposed scheme to the flexible system, the modal parameters and transfer function of the system can be obtained effectively.


### 3.1. The Prony Algorithm

Prony analysis supposes that the model is a linear combination of a series of exponential functions which have arbitrary amplitude, phase, frequency, and damping factor and then uses this model to fit equally spaced sampling data to determine corresponding parameters. The general expression of the discrete time signal is as follows:(1)
x
^
n
=
∑
i
=
1
p
b
i
z
i
n
,
n
=
0,1
,
…
,
N
-
1
,where b
i
=
A
i
exp
⁡
(
j
θ
i
), z
i
=
exp
⁡
a
i
+
j
2
π
f
i
Δ
t, p is the number of exponential functions, A
i is the ith amplitude of the exponential function, θ
i is the primary phase, f
i is the frequency, α
i is the damping coefficient, and Δ
t is the sampling interval. The process and steps of solving b
i and z
i and the values of the amplitude, phase, frequency, and damping factor by Prony algorithm were detailed in literature [18, 19]. In order to extract more characteristics of the system, this paper employs the multiple Prony models that differ from the single Prony model in the way of determining the sample function matrix, amplitude, and phase of each model. The main steps of the Prony algorithm are given as follows.Step 1.
Calculate the sample functionr
(
i
,
j
) and construct the matrix R
e by the following:(2)
R
e
=
r
1,0
r
1,1
⋯
r
1
,
p
r
2,0
r
2,1
⋯
r
2
,
p
⋮
⋮
⋮
⋮
r
p
e
,
0
r
p
e
,
1
⋯
r
p
e
,
p
e
p
e
≥
p
,where r
i
,
j
=
∑
n
=
p
e
N
-
1
x
n
-
j
x
∗
n
-
i, i
,
j
=
0,1
,
…
,
p
e.Step 2.
Use the Singular Value Decomposition and Total Least Squares (SVD-TLS) algorithm to determine the rank ofR
e.Step 3.
Solve the linear equation (3) to obtain a
1
,
…
,
a
p by the least squares estimation: (3)
r
0,0
r
0,1
⋯
r
0
,
p
r
1,0
r
1,1
⋯
r
1
,
p
⋮
⋮
⋯
⋮
r
p
,
0
r
p
,
1
⋯
r
p
,
p
1
a
1
⋮
a
p
=
ε
p
0
⋮
0
.Step 4.
Calculate the roots of characteristic polynomial1
+
a
1
z
-
1
+
⋯
+
a
p
z
-
p
=
0 and compute x
^
(
n
) according to recursive formula x
^
(
n
)
=
-
∑
i
=
1
p
a
i
x
^
(
n
-
i
).Step 5.
Calculate parameterb
i according to the following:(4)
1
1
⋯
1
z
1
z
2
⋯
z
p
⋮
⋮
⋮
⋮
z
1
N
-
1
z
2
N
-
1
⋯
z
p
N
-
1
b
1
b
2
⋮
b
p
=
x
^
0
x
^
1
⋮
x
^
N
-
1
.Step 6.
Calculate amplitudeA
i, phase θ
i, frequency f
i, and damping coefficient α
i according to the following:(5)
A
i
=
b
i
,
θ
i
=
arctan
⁡
Im
⁡
b
i
/
Re
⁡
b
i
2
π
Δ
t
,
α
i
=
ln
⁡
z
i
Δ
t
,
f
i
=
arctan
⁡
Im
⁡
z
i
/
Re
⁡
z
i
2
π
Δ
t
.


### 3.2. Transfer Function Identifying Based on Prony Algorithm

In this section, considering a Single Input Single Output (SISO) linear system,I
(
s
) and Y
(
s
), respectively, stand for input signal and output signal in Laplace transform and Y
(
s
)
=
I
(
s
)
G
(
s
); then we suppose that the system transfer function G
(
s
) can be expressed as(6)
G
s
=
∑
i
=
1
n
R
i
s
-
λ
i
,where n is the order of the model and R
i is the residue of transfer function in pole λ
i. Assuming that the input signal consists of a series of delayed signals, then I
(
s
) can be expressed as(7)
I
s
=
c
0
+
c
1
e
-
s
D
1
+
c
2
e
-
s
D
2
+
⋯
+
c
k
e
-
s
D
k
s
-
λ
n
+
1
,where c
i (i
=
0,1
,
…
,
k) is delay factor and D
j (j
=
1
,
…
,
k) is characteristic value.When an input signalI
(
s
) is applied to the system, the output signal Y
(
s
) can be expressed as (8)
Y
s
=
c
0
+
c
1
e
-
s
D
1
+
⋯
+
c
k
e
-
s
D
k
∑
i
=
1
n
R
i
s
-
λ
i
s
-
λ
n
+
1
.By partial fraction method, (8) can be transformed into (9)
Y
s
=
c
0
+
c
1
e
-
s
D
1
+
⋯
+
c
k
e
-
s
D
k
·
Q
n
+
1
s
-
λ
n
+
1
+
∑
i
=
1
n
Q
i
s
-
λ
i
,where Q
i
=
R
i
/
(
λ
i
-
λ
n
+
1
), i
=
1,2
,
…
,
n, and Q
n
+
1
=
-
∑
i
=
1
n
Q
i. Then the Laplace inverse transform of Y
(
s
) can be shown as(10)
Y
t
=
c
0
Q
n
+
1
e
λ
n
+
1
t
+
∑
i
=
1
n
Q
i
e
λ
i
t
u
t
+
∑
i
=
1
k
c
i
Q
n
+
1
e
λ
n
+
1
t
-
D
i
+
∑
i
=
1
n
Q
i
e
λ
i
t
-
D
i
u
t
-
D
i
.In order to apply the Prony method, let τ
=
t
-
D
k and v
(
τ
)
=
y
(
t
+
D
k
); then (10) can be expressed in the form(11)
v
τ
=
∑
j
=
1
n
+
1
B
j
e
λ
i
τ
,where (12)
B
j
=
Q
j
∑
i
=
0
k
c
i
e
λ
j
D
k
-
D
i
j
=
1,2
,
…
,
n
+
1
.Whenn, λ
i, and R
i in (6) are available, the system transfer function can be obtained. Residues are shown, respectively, as(13)
R
j
=
B
j
λ
j
-
λ
n
+
1
∑
i
=
0
k
c
i
e
λ
j
D
k
-
D
i
,
j
=
1,2
,
3
,
…
,
n
.


## 3.1. The Prony Algorithm

Prony analysis supposes that the model is a linear combination of a series of exponential functions which have arbitrary amplitude, phase, frequency, and damping factor and then uses this model to fit equally spaced sampling data to determine corresponding parameters. The general expression of the discrete time signal is as follows:(1)
x
^
n
=
∑
i
=
1
p
b
i
z
i
n
,
n
=
0,1
,
…
,
N
-
1
,where b
i
=
A
i
exp
⁡
(
j
θ
i
), z
i
=
exp
⁡
a
i
+
j
2
π
f
i
Δ
t, p is the number of exponential functions, A
i is the ith amplitude of the exponential function, θ
i is the primary phase, f
i is the frequency, α
i is the damping coefficient, and Δ
t is the sampling interval. The process and steps of solving b
i and z
i and the values of the amplitude, phase, frequency, and damping factor by Prony algorithm were detailed in literature [18, 19]. In order to extract more characteristics of the system, this paper employs the multiple Prony models that differ from the single Prony model in the way of determining the sample function matrix, amplitude, and phase of each model. The main steps of the Prony algorithm are given as follows.Step 1.
Calculate the sample functionr
(
i
,
j
) and construct the matrix R
e by the following:(2)
R
e
=
r
1,0
r
1,1
⋯
r
1
,
p
r
2,0
r
2,1
⋯
r
2
,
p
⋮
⋮
⋮
⋮
r
p
e
,
0
r
p
e
,
1
⋯
r
p
e
,
p
e
p
e
≥
p
,where r
i
,
j
=
∑
n
=
p
e
N
-
1
x
n
-
j
x
∗
n
-
i, i
,
j
=
0,1
,
…
,
p
e.Step 2.
Use the Singular Value Decomposition and Total Least Squares (SVD-TLS) algorithm to determine the rank ofR
e.Step 3.
Solve the linear equation (3) to obtain a
1
,
…
,
a
p by the least squares estimation: (3)
r
0,0
r
0,1
⋯
r
0
,
p
r
1,0
r
1,1
⋯
r
1
,
p
⋮
⋮
⋯
⋮
r
p
,
0
r
p
,
1
⋯
r
p
,
p
1
a
1
⋮
a
p
=
ε
p
0
⋮
0
.Step 4.
Calculate the roots of characteristic polynomial1
+
a
1
z
-
1
+
⋯
+
a
p
z
-
p
=
0 and compute x
^
(
n
) according to recursive formula x
^
(
n
)
=
-
∑
i
=
1
p
a
i
x
^
(
n
-
i
).Step 5.
Calculate parameterb
i according to the following:(4)
1
1
⋯
1
z
1
z
2
⋯
z
p
⋮
⋮
⋮
⋮
z
1
N
-
1
z
2
N
-
1
⋯
z
p
N
-
1
b
1
b
2
⋮
b
p
=
x
^
0
x
^
1
⋮
x
^
N
-
1
.Step 6.
Calculate amplitudeA
i, phase θ
i, frequency f
i, and damping coefficient α
i according to the following:(5)
A
i
=
b
i
,
θ
i
=
arctan
⁡
Im
⁡
b
i
/
Re
⁡
b
i
2
π
Δ
t
,
α
i
=
ln
⁡
z
i
Δ
t
,
f
i
=
arctan
⁡
Im
⁡
z
i
/
Re
⁡
z
i
2
π
Δ
t
.


## 3.2. Transfer Function Identifying Based on Prony Algorithm

In this section, considering a Single Input Single Output (SISO) linear system,I
(
s
) and Y
(
s
), respectively, stand for input signal and output signal in Laplace transform and Y
(
s
)
=
I
(
s
)
G
(
s
); then we suppose that the system transfer function G
(
s
) can be expressed as(6)
G
s
=
∑
i
=
1
n
R
i
s
-
λ
i
,where n is the order of the model and R
i is the residue of transfer function in pole λ
i. Assuming that the input signal consists of a series of delayed signals, then I
(
s
) can be expressed as(7)
I
s
=
c
0
+
c
1
e
-
s
D
1
+
c
2
e
-
s
D
2
+
⋯
+
c
k
e
-
s
D
k
s
-
λ
n
+
1
,where c
i (i
=
0,1
,
…
,
k) is delay factor and D
j (j
=
1
,
…
,
k) is characteristic value.When an input signalI
(
s
) is applied to the system, the output signal Y
(
s
) can be expressed as (8)
Y
s
=
c
0
+
c
1
e
-
s
D
1
+
⋯
+
c
k
e
-
s
D
k
∑
i
=
1
n
R
i
s
-
λ
i
s
-
λ
n
+
1
.By partial fraction method, (8) can be transformed into (9)
Y
s
=
c
0
+
c
1
e
-
s
D
1
+
⋯
+
c
k
e
-
s
D
k
·
Q
n
+
1
s
-
λ
n
+
1
+
∑
i
=
1
n
Q
i
s
-
λ
i
,where Q
i
=
R
i
/
(
λ
i
-
λ
n
+
1
), i
=
1,2
,
…
,
n, and Q
n
+
1
=
-
∑
i
=
1
n
Q
i. Then the Laplace inverse transform of Y
(
s
) can be shown as(10)
Y
t
=
c
0
Q
n
+
1
e
λ
n
+
1
t
+
∑
i
=
1
n
Q
i
e
λ
i
t
u
t
+
∑
i
=
1
k
c
i
Q
n
+
1
e
λ
n
+
1
t
-
D
i
+
∑
i
=
1
n
Q
i
e
λ
i
t
-
D
i
u
t
-
D
i
.In order to apply the Prony method, let τ
=
t
-
D
k and v
(
τ
)
=
y
(
t
+
D
k
); then (10) can be expressed in the form(11)
v
τ
=
∑
j
=
1
n
+
1
B
j
e
λ
i
τ
,where (12)
B
j
=
Q
j
∑
i
=
0
k
c
i
e
λ
j
D
k
-
D
i
j
=
1,2
,
…
,
n
+
1
.Whenn, λ
i, and R
i in (6) are available, the system transfer function can be obtained. Residues are shown, respectively, as(13)
R
j
=
B
j
λ
j
-
λ
n
+
1
∑
i
=
0
k
c
i
e
λ
j
D
k
-
D
i
,
j
=
1,2
,
3
,
…
,
n
.


## 4. Input Shaper Design for Vibration Control

A number of techniques have been proposed to suppress the vibration, such as feedback control [20], positive position feedback (PPF) [21], and input shaper [22]. Both closed-loop feedback technique and PPF use the system’s state to reduce the vibration. However, the input shaping is implemented by convolving a sequence of impulses with a desired system command to produce a shaped input to drive the actuator to reduce the system oscillations [23]. Furthermore, only the modal parameters are required for the design of input shaper and these parameters can be obtained from Prony algorithms presented in Section 3. So, the input shaping technique is utilized to attenuate the oscillations in this paper.Input shaper is a time-delay filter that convolves the command input signal with a series of impulse sequences based on the natural frequencies and damping coefficients of the system. If the timing and amplitudes of impulses are chosen to match modal parameters of the controlled system, vibration can be eliminated [24]. The basic structure of input shaper can be expressed as (14)
C
s
=
∑
i
n
A
i
e
-
s
t
i
,where A
i is the amplitude of the impulse, t
i is the time of the impulse input, and n is the number of impulse sequences.The derivation of the input shaper is based on the linear system theory. An arbitrary linear vibrating system can be specified as a cascaded multiple second-order oscillation system, so the analysis of the original system can be realized by studying a second-order system. A second-order oscillation system with natural frequencyω
n and damping coefficient ξ
(
0
<
ξ
<
1
) can be expressed as(15)
G
s
=
ω
n
2
s
2
+
2
ξ
ω
n
s
+
ω
n
2
.The unit impulse response of a second-order system can be written as(16)
y
t
=
ω
n
1
-
ξ
2
e
-
ξ
ω
n
sin
⁡
ω
d
t
-
t
i
,where ω
d
=
ω
n
1
-
ξ
2. When input shaper is applied to a controlled system, the impulse response of this system can be written as(17)
y
t
=
∑
i
=
1
n
ω
n
A
i
1
-
ξ
2
e
-
ξ
ω
n
sin
⁡
ω
d
t
-
t
i
.Define the ratio of the amplitude of (17) and (16) as vibration of the system: (18)
V
ω
n
,
ξ
=
e
-
ξ
ω
n
t
n
C
2
ω
n
,
ξ
+
S
2
ω
n
,
ξ
,where(19)
C
ω
n
,
ξ
=
∑
i
=
1
n
A
i
e
-
ξ
ω
n
t
n
cos
⁡
ω
n
1
-
ξ
2
t
i
,
S
ω
n
,
ξ
=
∑
i
=
1
n
A
i
e
-
ξ
ω
n
t
n
sin
⁡
ω
n
1
-
ξ
2
t
i
.In order to make the vibration equal to zero,C
(
ω
n
,
ξ
)
=
0 and S
(
ω
n
,
ξ
)
=
0 are prerequisite. When adding different constraints, diverse input shapers can be designed. The zero vibration input shaper (ZV shaper) is the simplest shaper, consisting of two impulses. In order to shorten the time of shaping, we can set t
1
=
0 and let A
i
>
0 and ∑
i
=
1
n
A
i
=
1 to achieve the set point of system. Hence, the constraints of ZV input shaper are expressed as A
1
+
A
2
=
1, A
1
>
0, and A
2
>
0; the solutions can be written as(20)
A
1
=
1
1
+
K
,
t
1
=
0
,
A
2
=
K
1
+
K
,
t
2
=
T
,where K
=
e
-
ξ
π
/
1
-
ξ
2 and T
=
π
/
ω
n
1
-
ξ
2. When the identified modal parameters cannot fit the actual system parameters well, ZV shaper cannot eliminate the vibration effectively. In order to improve the robustness of the system, a constraint (
d
/
d
ω
)
V
(
ω
,
ξ
)
=
0 is applied to configure multiple poles for the system and then a zero vibration and derivative (ZVD) shaper with three impulses can be solved. The parameters of ZVD shaper can be expressed as(21)
A
1
=
1
1
+
2
K
+
K
2
,
t
1
=
0
,
A
2
=
2
K
1
+
2
K
+
K
2
,
t
2
=
T
,
A
3
=
K
2
1
+
2
K
+
K
2
,
t
3
=
2
T
.


## 5. Experiment Results

In order to verify the feasibility of the identification algorithm and test the effect of the vibration control strategy, a piezoelectric cantilever beam described in Section2 is used to establish the experimental platform. The excitation signal is programmed by software LabVIEW. This generated signal is output through the D/A conversion module of data acquisition card whose analog output range is −10 V~+10 V and then the signal is amplified by a voltage amplifier to actuate the MFC. Finally, the displacement of piezoelectric cantilever beam is monitored by a laser sensor put at the tip vertex point of the beam at the free end and recorded at 1 kHz via an A/D converter. The experiment setup is shown in Figure 4.Figure 4
The experiment setup.


### 5.1. The Identification Results Based on Prony Algorithm

A 1 V step signal as excitation signal is applied to piezoelectric cantilever beam system and output displacement of the beam is acquired as identification signal data with the sampling interval of 5 ms in this section. The 1 V step signal is expressed asc
0
=
1, λ
n
+
1
=
0, and c
i
=
0 (i
=
1
,
…
,
k) according to (7). When the step signal lasts for a long time, the creep will engender, so the whole sampling time is chosen as 1 s to avoid the influence of creep effect on the identification. In general, it is inevitable that the noise exists in the measurement of displacement via a sensor, but, unfortunately, the Prony algorithm is sensitive to the data nosing. Hence, the wavelet threshold denoising is used to preprocess the measuring response signal and improve signal-to-noise ratio before using Prony algorithm. The result of filtering is shown in Figure 5. Furthermore, direct current component could lead to poor identification precision when using the Prony algorithm, so it is imperative to extract direct current component from the signal data. The direct current component is 0.195 mm via calculation. Then the response signal with filtering and removing the direct current component is used as identification data and the identification results based on Prony algorithm are shown in Figure 6(a). The identification error is shown in Figure 6(b) and the root mean square error (RMSE) of identification is 0.005 mm. The main parameters obtained from identification results are shown in Table 2.Table 2
Identification results on system parameters.
Order
A
i
f
i/Hz
α
i
λ
i
R
i
1
0.00672
44.209
−45.898
−45.898 + 277.754i
−0.066 − 1.891i
2
0.00672
44.209
−45.898
−45.898 − 277.754i
−0.066 − 1.891i
3
0.03093
0
−4.536
−4.536
0.140
4
0.08045
19.112
−1.4967
−1.497 + 120.087i
0.506 − 9.649i
5
0.08045
19.112
−1.4967
−1.497 − 120.087i
0.506 + 9.649iFigure 5
The result of signal filtering preprocessing.Figure 6
The results of Prony identification.
(a)
The fitting result by Prony
(b)
The error of identificationFrom Figure6, the fitting signal by Prony algorithm can describe the actual experimental data exactly, which could illustrate that Prony algorithm has active effects on dynamic modeling for flexible structure. As shown in Table 2, the order of the system is 5, including a first-order inertial element and two second-order oscillation subsystems. Oscillation mode contains coneigenvalues 1 and 2 (44.209 Hz) and coneigenvalues 4 and 5 (19.112 Hz). It is obvious that the oscillation mode 1 and mode 2 attenuate rapidly because the damping factor of them is relatively large and oscillation amplitude is relatively small; therefore the dominant oscillation frequency is 19.112 Hz. The system transfer function can be expressed as (22) and frequency response curve is shown in Figure 7.(22)
G
s
=
1.022
s
4
+
3474
s
3
+
3.233
×
1
0
5
s
2
+
2.004
×
1
0
8
s
+
1.063
×
1
0
9
s
5
+
99.33
s
4
+
9.439
×
1
0
4
s
3
+
1.987
×
1
0
6
s
2
+
1.15
×
1
0
9
s
+
5.185
×
1
0
9
.Figure 7
The bode diagram of identified system.


### 5.2. Experimental Validation of Performance of Input Shaper

The transfer function and modal parameters obtained from the Prony algorithm are used directly to design the input shaper in this section. The dominant modal parameters are shown as natural frequencyω
n
=
120.10 rad/s and damping ratio ξ
=
0.01246. According to (20), the parameters of ZV input shaper for dominant modal are as follows:(23)
A
1
=
0.5098
,
t
1
=
0
,
A
2
=
0.4902
,
t
2
=
0.02616
.And, according to (21), the parameters of ZVD input shaper for dominant modal are as follows:(24)
A
1
=
0.2599
,
t
1
=
0
,
A
2
=
0.4998
,
t
2
=
0.02616
,
A
3
=
0.2403
,
t
3
=
0.05232
.Figure8 shows the experimental results of the vibration control effect under a 1 V step signal by using ZV shaper and ZVD shaper compared with the one without shaper. It can be concluded that the input shaper has good effect on suppressing vibration and the ZVD shaper is superior to ZV shaper. Under the effect of ZVD input shaper, the settling time is shortened distinctly from 2 s to 0.5 s and the overshoot is null. Furthermore, it can also adequately illustrate that the identification results based on Prony algorithm are effective and accurate.Figure 8
Step response of the system with vibration suppression.


## 5.1. The Identification Results Based on Prony Algorithm

A 1 V step signal as excitation signal is applied to piezoelectric cantilever beam system and output displacement of the beam is acquired as identification signal data with the sampling interval of 5 ms in this section. The 1 V step signal is expressed asc
0
=
1, λ
n
+
1
=
0, and c
i
=
0 (i
=
1
,
…
,
k) according to (7). When the step signal lasts for a long time, the creep will engender, so the whole sampling time is chosen as 1 s to avoid the influence of creep effect on the identification. In general, it is inevitable that the noise exists in the measurement of displacement via a sensor, but, unfortunately, the Prony algorithm is sensitive to the data nosing. Hence, the wavelet threshold denoising is used to preprocess the measuring response signal and improve signal-to-noise ratio before using Prony algorithm. The result of filtering is shown in Figure 5. Furthermore, direct current component could lead to poor identification precision when using the Prony algorithm, so it is imperative to extract direct current component from the signal data. The direct current component is 0.195 mm via calculation. Then the response signal with filtering and removing the direct current component is used as identification data and the identification results based on Prony algorithm are shown in Figure 6(a). The identification error is shown in Figure 6(b) and the root mean square error (RMSE) of identification is 0.005 mm. The main parameters obtained from identification results are shown in Table 2.Table 2
Identification results on system parameters.
Order
A
i
f
i/Hz
α
i
λ
i
R
i
1
0.00672
44.209
−45.898
−45.898 + 277.754i
−0.066 − 1.891i
2
0.00672
44.209
−45.898
−45.898 − 277.754i
−0.066 − 1.891i
3
0.03093
0
−4.536
−4.536
0.140
4
0.08045
19.112
−1.4967
−1.497 + 120.087i
0.506 − 9.649i
5
0.08045
19.112
−1.4967
−1.497 − 120.087i
0.506 + 9.649iFigure 5
The result of signal filtering preprocessing.Figure 6
The results of Prony identification.
(a)
The fitting result by Prony
(b)
The error of identificationFrom Figure6, the fitting signal by Prony algorithm can describe the actual experimental data exactly, which could illustrate that Prony algorithm has active effects on dynamic modeling for flexible structure. As shown in Table 2, the order of the system is 5, including a first-order inertial element and two second-order oscillation subsystems. Oscillation mode contains coneigenvalues 1 and 2 (44.209 Hz) and coneigenvalues 4 and 5 (19.112 Hz). It is obvious that the oscillation mode 1 and mode 2 attenuate rapidly because the damping factor of them is relatively large and oscillation amplitude is relatively small; therefore the dominant oscillation frequency is 19.112 Hz. The system transfer function can be expressed as (22) and frequency response curve is shown in Figure 7.(22)
G
s
=
1.022
s
4
+
3474
s
3
+
3.233
×
1
0
5
s
2
+
2.004
×
1
0
8
s
+
1.063
×
1
0
9
s
5
+
99.33
s
4
+
9.439
×
1
0
4
s
3
+
1.987
×
1
0
6
s
2
+
1.15
×
1
0
9
s
+
5.185
×
1
0
9
.Figure 7
The bode diagram of identified system.


## 5.2. Experimental Validation of Performance of Input Shaper

The transfer function and modal parameters obtained from the Prony algorithm are used directly to design the input shaper in this section. The dominant modal parameters are shown as natural frequencyω
n
=
120.10 rad/s and damping ratio ξ
=
0.01246. According to (20), the parameters of ZV input shaper for dominant modal are as follows:(23)
A
1
=
0.5098
,
t
1
=
0
,
A
2
=
0.4902
,
t
2
=
0.02616
.And, according to (21), the parameters of ZVD input shaper for dominant modal are as follows:(24)
A
1
=
0.2599
,
t
1
=
0
,
A
2
=
0.4998
,
t
2
=
0.02616
,
A
3
=
0.2403
,
t
3
=
0.05232
.Figure8 shows the experimental results of the vibration control effect under a 1 V step signal by using ZV shaper and ZVD shaper compared with the one without shaper. It can be concluded that the input shaper has good effect on suppressing vibration and the ZVD shaper is superior to ZV shaper. Under the effect of ZVD input shaper, the settling time is shortened distinctly from 2 s to 0.5 s and the overshoot is null. Furthermore, it can also adequately illustrate that the identification results based on Prony algorithm are effective and accurate.Figure 8
Step response of the system with vibration suppression.


## 6. Conclusion

This paper proposes an identification method combining with input shaping technique to suppress the vibrations of the flexible system for precise positioning. The dynamic modeling and vibration control are performed in sequence. The dynamic model of the system is considered as a linear time-invariant system, and the Prony algorithm is used to identify the modal parameters and transfer function of the system from the step response signal. In order to improve the accuracy of identification, the data nosing is filtered and the direct current component is removed. The identification results show that the Prony algorithm is quite accurate and effective. The RMSE of identification error is only 0.005 mm and the identified transfer function with five modes can describe the system’s characteristic exactly. The first-order vibration mode is analyzed as 19.112 Hz from the identification results and then the input shaper is designed to succeed an active vibration control. The experimental results show that the vibration is reduced distinctly and the settling time is also shortened significantly after applying the input shaper.The works presented in this paper have confirmed that the Prony algorithm combining with input shaping technique can exhibit admirable performance on identification for active vibration control. Such approach is of great interest because it is easy to implement and the vibration is reduced in open loop. What is more, the identification method presented in this paper is well adopted not only for flexible structure but also for other similar manipulation mechanisms and the identified transfer function can also be used to simulate the system and design other vibration controllers. Finally, the future work will focus on the robustness of the controller considering the parameter variations and environmental disturbance.



---

*Source: 9891054-2016-09-20.xml*