# A New Approach for the Stability Analysis of Wave Networks

Authors: Ya Xuan Zhang; Gen Qi Xu
Journal: Abstract and Applied Analysis
 (2014)

---

## Abstract

We introduce a new approach to investigate the stability of controlled tree-shaped wave networks and subtrees of complex wave networks. It is motivated by regarding the network as branching out from a single edge. We present the recursive relations of the Laplacian transforms of adjacent edges of the system in its branching order, which form the characteristic equation. In the stability analysis, we estimate the infimums of the recursive expressions in the inverse order based on the spectral analysis. It is a feasible way to check whether the system is exponentially stable under any control strategy or parameter choice. As an application we design the control law and study the stability of a 12-edge tree-shaped wave network.


---

## Body

## 1. Introduction

The stability analysis of networks of flexible elements, such as vibrating strings, beams, membranes, and plates has been intensively and extensively studied in the past decades (see [1–3] for some general works and [4–17] for wave networks). A wave network can be asymptotically stabilized by some appropriate feedback control strategy. However, its exponential stability is difficult to check. The exponential stability of a dynamic system is usually analyzed by multiplier method ([18, 19], etc.), resolvent estimation ([20, 21], etc.), or spectral analysis and Riesz basis approach ([4, 8–12, 14], etc.). For wave networks, an appropriate multiplier has not been found so far and the resolvent is too complicated to estimate; thus the spectral analysis becomes a possible technique. In fact, suppose that the spectral determined growth assumption (SDGA) holds for a network, which means that the energy decay rate of the network is exactly the supremum of the real parts of its spectra. Then a dissipative system is exponentially stable if and only if the imaginary axis is not an asymptote of its spectra. That is, the characteristic equation D
(
λ
) of the network satisfies
(1)
inf
⁡
σ
∈
ℝ
|
D
(
i
σ
)
|
>
0
.
Generally speaking, most controlled wave networks satisfy SDGA under certain conditions. However, as for the spectrum of these networks, it is only known that the spectra lie in a vertical strip in the complex plane. More precise properties, such as the asymptote, of the spectra remain unknown in general. Besides, the characteristic equation has been derived for wave networks of some special types (e.g., [2] presented nice results about the characteristic equation of generic trees), but its expression form does not seem very convenient for checking (1). These are challenging problems in both stability analysis of wave networks and spectral theory of nonadjoint operators.In this paper, we will introduce a relatively easy approach to study the exponential stability of a controlled wave network. It is valid for tree-shaped networks and subtrees of complex networks. The basic idea is motivated by the natural growth of a tree. Regarding a tree-shaped network as branching out from a single edge, we can obtain a simple form of its characteristic equation by recurrence, which leads to an easy verification of (1). Our idea is illustrated below in detail.First we consider a tree of only one edge, the dynamic behavior of which is governed by the following wave equation and boundary conditions:(2)
T
1
y
1
,
x
x
(
x
,
t
)
=
m
1
y
1
,
t
t
(
x
,
t
)
,
x
∈
(
0,1
)
,
t
>
0
y
1
(
0
,
t
)
=
0
T
1
y
1
,
x
(
1
,
t
)
=
R
1
(
t
)
,
where R
1
(
t
) denotes the extra force acting at the boundary end of the edge.Then we consider another tree ofN edges, whose dynamic behavior is governed by
(3)
T
k
y
k
,
x
x
(
x
,
t
)
=
m
k
y
k
,
t
t
(
x
,
t
)
,
x
∈
(
0,1
)
,
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
t
>
0
,
k
=
1,2
,
…
,
N
y
1
(
0
,
t
)
=
0
y
1
(
1
,
t
)
=
y
k
(
0
,
t
)
,
k
=
2
,
…
,
N
T
1
y
1
,
x
(
1
,
t
)
=
∑
k
=
2
N
T
k
y
k
,
x
(
0
,
t
)
T
k
y
k
,
x
(
1
,
t
)
=
R
k
(
t
)
,
k
=
2
,
…
,
N
,
where R
k
(
t
), k
=
2
,
…
,
N are the extra forces acting at the corresponding boundary ends.Comparing (2) and (3), we can regard the tree-shaped wave network (3) as branching out from the single-edge tree (2), based on the “branching force principle”:
(4)
R
1
(
t
)
=
∑
k
=
2
N
‍
T
k
y
k
,
x
(
0
,
t
)
.We can see this “branching” process clearly in Figure1.Tree-shaped wave network (3) branching out from the single-edge tree (2).
(a)
(b)According to this principle, we can describe more complicated tree-shaped wave networks. Moreover, we can deduce from it the characteristic equation of the networks, which has a simple form to check (1). This can be done by the Laplacian transform of (3):
(5)
T
k
f
k
,
x
x
(
x
,
λ
)
=
m
k
λ
2
f
k
(
x
,
λ
)
,
k
k
x
∈
(
0,1
)
,
k
=
1,2
,
…
,
N
f
1
(
0
,
λ
)
=
0
f
1
(
1
,
λ
)
=
f
k
(
0
,
λ
)
,
k
=
2
,
…
,
N
T
1
f
1
,
x
(
1
,
λ
)
=
R
^
1
(
λ
)
R
^
1
(
λ
)
=
∑
k
=
2
N
T
k
f
k
,
x
(
0
,
λ
)
T
k
f
k
,
x
(
1
,
λ
)
=
R
^
k
(
λ
)
,
k
=
2
,
…
,
N
,
where f
k
(
x
,
λ
) and R
^
k
(
λ
) are the transforms of y
k
(
x
,
t
) and R
k
(
t
), respectively.Note that forλ
≠
0, the solutions f
k
(
x
,
λ
), k
=
1,2
,
…
,
N to the differential equations in (5) have the general form
(6)
f
k
(
x
,
λ
)
=
c
k
[
F
k
(
λ
)
cosh
⁡
λ
ρ
k
x
-
G
k
(
λ
)
sinh
⁡
λ
ρ
k
x
]
,
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
ρ
k
=
m
k
T
k
,
k
=
1,2
,
…
,
N
,
where F
k
(
λ
)
,
G
k
(
λ
)
∈
ℂ are arbitrary constants related to λ, and c
k are nonzero auxiliary constants.Substituting (6) into the first boundary condition in (5), we obtain that
(7)
F
1
(
λ
)
=
0
.Substituting (6) into f
k
(
1
,
λ
) and the last boundary condition in (5), we have that
(8)
f
k
(
1
,
λ
)
=
c
k
[
F
k
(
λ
)
cosh
⁡
λ
ρ
k
-
G
k
(
λ
)
sinh
⁡
λ
ρ
k
]
,
R
^
k
(
λ
)
=
λ
T
k
ρ
k
c
k
[
F
k
(
λ
)
sinh
⁡
λ
ρ
k
-
G
k
(
λ
)
cosh
⁡
λ
ρ
k
]
,
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
=
2
,
…
,
N
.
From the above we can express the coefficients F
k
(
λ
) and G
k
(
λ
) by f
k
(
1
,
λ
) and R
^
k
(
λ
) for k
=
2
,
…
,
N as follows:
(9)
F
k
(
λ
)
=
1
c
k
T
k
ρ
k
[
T
k
ρ
k
f
k
(
1
,
λ
)
cosh
⁡
λ
ρ
k
-
1
λ
R
^
k
(
λ
)
sinh
⁡
λ
ρ
k
]
,
G
k
(
λ
)
=
1
c
k
T
k
ρ
k
[
T
k
ρ
k
f
k
(
1
,
λ
)
sinh
⁡
λ
ρ
k
-
1
λ
R
^
k
(
λ
)
cosh
⁡
λ
ρ
k
]
,
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
=
2
,
…
,
N
,
where R
^
k
(
λ
), k
=
2
,
…
,
N are to be determined.Substituting (6) into other connective conditions in (5), we get that
(10)
R
^
1
(
λ
)
=
λ
c
1
T
1
ρ
1
[
F
1
(
λ
)
sinh
⁡
λ
ρ
1
-
G
1
(
λ
)
cosh
⁡
λ
ρ
1
]
=
-
λ
∑
k
=
2
N
c
k
T
k
ρ
k
G
k
(
λ
)
f
1
(
1
,
λ
)
=
c
1
[
F
1
(
λ
)
sinh
⁡
λ
ρ
1
-
G
1
(
λ
)
cosh
⁡
λ
ρ
1
]
=
c
k
F
k
(
λ
)
,
k
=
2
,
…
,
N
.
Solving the above equations, we can express F
1
(
λ
) and G
1
(
λ
) by F
k
(
λ
), G
k
(
λ
), k
=
2
,
…
,
N with appropriately chosen constants c
k, k
=
2
,
…
,
N. It remains to be determined the expressions of F
k
(
λ
)
,
G
k
(
λ
), k
=
2
,
…
,
N, that is, the expressions of R
^
k
(
λ
) (or equivalently R
k
(
t
)), k
=
2
,
…
,
N. This can be done in the following two cases.Case 1.
The network under consideration is described as (3).
LetR
k
(
t
), k
=
2
,
…
,
N be the feedback controls, such as the collocated velocity feedback control R
k
(
t
)
=
-
α
y
k
,
x
(
1
,
t
), k
=
2
,
…
,
N. Then we can directly get the expressions of F
k
(
λ
)
,
G
k
(
λ
), k
=
2
,
…
,
N by (9).Case 2.
The network under consideration is a bigger tree branching out from (3).
We proceed similarly to the boundary of the system by recurrence.In both cases, we can finally obtain the recursive expressions ofF
1
(
λ
), and hence F
1
(
λ
)
=
0 is the characteristic equation of the system. The process is in accordance with the natural growth of a tree. In the stability analysis, we only need to estimate the infimums of each recursive expression in the inverse order to assert whether (1) holds or not. Thus the exponential stability analysis can be carried out in an easy manner. This can be similarly carried out for subtrees of complex wave networks.The content is arranged as follows. In Section2, we state our main results for general tree-shaped wave network and present a complete proof. In Section 3, we investigate a 12-edge tree-shaped wave network as an application to the main results. In Section 4, we give a conclusion.


## 2. Main Results

In this section, we will give a rigorous statement of the main results in general. For this aim, we deal with a general subtree𝒜 of a complex wave network G. We adopt the notations introduced in [1, pp. 104] to describe a tree or a subtree.Let𝒜 be a subtree (or a tree). By the degree of a vertex of 𝒜 we mean the number of the edges that branch out from it. A vertex is called a boundary vertex if its degree is one and an interior vertex otherwise. Let N be the number of the edges of 𝒜.Choose a boundary vertex as the root of𝒜, denoted by ℛ. The other edges and vertices are denoted by e
α
¯ and v
α
¯, respectively, where α
¯
=
(
α
1
,
…
,
α
k
) is a multi-index (possibly empty) defined by recurrence in the following way. For the edge containing ℛ we choose the empty index; that is, it is denoted by e and its vertex different from ℛ is denoted by v. Assume that the interior vertex v
α
¯, contained in the edge e
α
¯, has multiplicity m
α
¯
+
1. Then there are m
α
¯ edges, different from e
α
¯, branching out from v
α
¯. Denote them by e
α
¯
∘
β, where β
=
1,2
,
…
,
m
α
¯, α
¯
∘
β
=
(
α
1
,
…
,
α
k
,
β
), and the other vertex of the edge e
α
¯
∘
β is denoted by v
α
¯
∘
β.Let𝒜
M and 𝒜
S be the set of all the interior and boundary vertices, respectively, ℛ being excepted. Set 𝒥
M
=
{
α
¯
∣
v
α
¯
∈
𝒜
M
}, 𝒥
S
=
{
α
¯
∣
v
α
¯
∈
𝒜
S
}, and 𝒥
=
𝒥
M
∪
𝒥
S.Suppose without loss of generality that the length of the edgee
α
¯ is 1. Then e
α
¯ can be parameterized by its arc length by means of the function π
α
¯ defined as π
α
¯
:
[
0,1
]
→
e
α
¯
,
α
¯
∈
𝒥 such that π
α
¯
(
0
)
=
v
α
¯ and π
α
¯
(
1
) is the other vertex of e
α
¯. In this way, e
α
¯ and its end points can be identified with the interval [
0,1
].Assume that a tree-shaped subnetwork coincides with𝒜 at rest and performs a small perpendicular vibration. Let y
α
¯
=
y
α
¯
(
x
,
t
)
:
[
0,1
]
×
ℝ
+
→
ℝ describe the transversal displacement of e
α
¯ at position π
α
¯
(
x
) and time t. Suppose that all the interior vertices satisfy the geometric continuity condition and Kirchhoff law. Then the dynamical behavior of the subtree 𝒜 is governed by
(11)
T
α
¯
y
x
x
α
¯
(
x
,
t
)
=
m
α
¯
y
t
t
α
¯
(
x
,
t
)
,
x
∈
(
0,1
)
,
t
>
0
,
α
¯
∈
𝒥
,
y
α
¯
(
1
,
t
)
=
y
α
¯
∘
β
(
0
,
t
)
,
α
¯
∈
𝒥
M
,
β
=
1,2
,
…
,
m
α
¯
,
T
α
¯
y
x
α
¯
(
1
,
t
)
=
R
α
¯
(
t
)
,
α
¯
∈
𝒥
M
,
R
α
¯
(
t
)
=
∑
β
=
1
m
α
¯
‍
T
α
¯
∘
β
y
x
α
¯
∘
β
(
0
,
t
)
,
α
¯
∈
𝒥
M
.The Laplacian transform of the subtree (11) is
(12)
T
α
¯
f
x
x
α
¯
(
x
,
λ
)
=
m
α
¯
λ
2
f
α
¯
(
x
,
λ
)
,
x
∈
(
0,1
)
,
α
¯
∈
𝒥
,
f
α
¯
(
1
,
λ
)
=
f
α
¯
∘
β
(
0
,
λ
)
,
α
∈
𝒥
M
,
β
=
1,2
,
…
,
m
α
¯
,
T
α
¯
f
x
α
¯
(
1
,
λ
)
=
R
^
α
¯
(
λ
)
,
α
¯
∈
𝒥
M
,
R
^
α
¯
(
λ
)
=
∑
β
=
1
m
α
¯
‍
T
α
¯
∘
β
f
x
α
¯
∘
β
(
0
,
λ
)
,
α
¯
∈
𝒥
M
,
where f
α
¯
(
x
,
λ
) and R
^
α
¯
(
λ
) are the transforms of y
α
¯
(
x
,
t
) and R
α
¯
(
t
), respectively.According to the differential equations in (12), we can set
(13)
f
α
¯
(
x
,
λ
)
=
c
α
¯
[
F
α
¯
(
λ
)
cosh
⁡
λ
ρ
α
¯
x
-
G
α
¯
(
λ
)
sinh
⁡
λ
ρ
α
¯
x
]
,
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
ρ
α
¯
=
m
α
¯
T
α
¯
,
α
¯
∈
𝒥
,
where the coefficients c
α
¯ are nonzero auxiliary coefficients; F
α
¯
(
λ
)
,
G
α
¯
(
λ
) are arbitrary constants related to λ, which are to be determined as recursive expressions.For any fixedα
¯
∈
𝒥
M, we consider the edge e
α
¯ and the m
α
¯ edges e
α
¯
∘
β
,
β
=
1,2
,
…
,
m
α
¯, which branch out from e
α
¯. Their vibration is governed by
(14)
T
α
¯
y
x
x
α
¯
(
x
,
t
)
=
m
α
¯
y
t
t
α
¯
(
x
,
t
)
,
x
∈
(
0,1
)
,
t
>
0
T
α
¯
∘
β
y
x
x
α
¯
∘
β
(
x
,
t
)
=
m
α
¯
∘
β
y
t
t
α
¯
∘
β
(
x
,
t
)
,
x
∈
(
0,1
)
,
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
t
>
0
,
β
=
1,2
,
…
,
m
α
¯
y
α
¯
(
1
,
t
)
=
y
α
¯
∘
β
(
0
,
t
)
,
β
=
1,2
,
…
,
m
α
¯
T
α
¯
y
x
α
¯
(
1
,
t
)
=
R
α
¯
(
t
)
R
α
¯
(
t
)
=
∑
β
=
1
m
α
¯
‍
T
α
¯
∘
β
y
x
α
¯
∘
β
(
0
,
t
)
.According to (12) and (13), the Laplacian transform of (14) yields
(15)
f
α
¯
(
x
,
λ
)
=
c
α
¯
[
F
α
¯
(
λ
)
cosh
⁡
λ
ρ
α
¯
x
-
G
α
¯
(
λ
)
sinh
⁡
λ
ρ
α
¯
x
]
f
α
¯
∘
β
(
x
,
λ
)
=
c
α
¯
∘
β
[
F
α
¯
∘
β
(
λ
)
cosh
⁡
λ
ρ
α
¯
∘
β
x
-
G
α
¯
∘
β
(
λ
)
sinh
⁡
λ
ρ
α
¯
∘
β
x
]
,
k
k
k
k
k
k
k
k
k
k
k
k
k
i
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
β
=
1,2
,
…
,
m
α
¯
f
α
¯
(
1
,
λ
)
=
f
α
¯
∘
β
(
0
,
λ
)
,
β
=
1,2
,
…
,
m
α
¯
T
α
¯
f
x
α
¯
(
1
,
λ
)
=
R
^
α
¯
(
λ
)
R
^
α
¯
(
λ
)
=
∑
β
=
1
m
α
¯
‍
T
α
¯
∘
β
f
x
α
¯
∘
β
(
0
,
λ
)
.The following theorem indicates the recursive relations between the coefficientsF
α
¯
(
λ
), G
α
¯
(
λ
), and F
α
¯
∘
β
(
λ
), G
α
¯
∘
β
(
λ
), β
=
1,2
,
…
,
m
α
¯ of the adjacent edges. It provides us with the easy manner for the exponential stability analysis.Theorem 1.
LetG be a complex wave network and 𝒜 be a subtree of G described by (11). Assume that λ
≠
0. Then for any α
¯
∈
𝒥
M, the recursive expressions of F
α
¯
(
λ
) and G
α
¯
(
λ
) in F
α
¯
∘
β
(
λ
), G
α
¯
∘
β
(
λ
), β
=
1,2
,
…
,
m
α
¯ are given by
(16)
F
α
¯
(
λ
)
=
1
T
α
¯
ρ
α
¯
[
T
α
¯
ρ
α
¯
∏
β
=
1
m
α
¯
F
α
¯
∘
β
(
λ
)
cosh
⁡
λ
ρ
α
¯
k
k
k
k
k
k
k
k
+
sinh
⁡
λ
ρ
α
¯
∑
β
=
1
m
α
¯
T
α
¯
∘
β
ρ
α
¯
∘
β
G
α
¯
∘
β
(
λ
)
k
k
k
k
k
k
k
k
×
∏
γ
=
1
,
γ
≠
β
m
α
¯
F
α
¯
∘
γ
(
λ
)
]
G
α
¯
(
λ
)
=
1
T
α
¯
ρ
α
¯
[
T
α
¯
ρ
α
¯
∏
β
=
1
m
α
¯
F
α
¯
∘
β
(
λ
)
sinh
⁡
λ
ρ
α
¯
k
k
k
k
k
k
k
k
+
cosh
⁡
λ
ρ
α
¯
∑
β
=
1
m
α
¯
T
α
¯
∘
β
ρ
α
¯
∘
β
G
α
¯
∘
β
(
λ
)
k
k
k
k
k
k
k
k
×
∏
γ
=
1
,
γ
≠
β
m
α
¯
F
α
¯
∘
γ
(
λ
)
]
.
Moreover, if
(17)
inf
⁡
σ
∈
ℝ
|
F
α
¯
∘
β
(
i
σ
)
|
>
0
,
ℜ
(
G
α
¯
∘
β
(
i
σ
)
F
α
¯
∘
β
(
i
σ
)
¯
)
>
0
hold for every β
=
1,2
,
…
,
m
α
¯, then
(18)
inf
⁡
σ
∈
ℝ
|
F
α
¯
(
i
σ
)
|
>
0
,
ℜ
(
G
α
¯
(
i
σ
)
F
α
¯
(
i
σ
)
¯
)
>
0
.Proof.
According to (15), we obtain that
(19)
f
α
¯
(
1
,
λ
)
=
c
α
¯
[
F
α
¯
(
λ
)
cosh
⁡
λ
ρ
α
¯
-
G
α
¯
(
λ
)
sinh
⁡
λ
ρ
α
¯
]
R
^
α
¯
(
λ
)
=
T
α
¯
f
x
α
¯
(
1
,
λ
)
=
λ
c
α
¯
T
α
¯
ρ
α
¯
[
F
α
¯
(
λ
)
sinh
⁡
λ
ρ
α
¯
-
G
α
¯
(
λ
)
cosh
⁡
λ
ρ
α
¯
]
,
which indicates
(20)
F
α
¯
(
λ
)
=
1
c
α
¯
T
α
¯
ρ
α
¯
[
T
α
¯
ρ
α
¯
f
α
¯
(
1
,
λ
)
cosh
⁡
λ
ρ
α
¯
-
1
λ
R
^
α
¯
(
λ
)
sinh
⁡
λ
ρ
α
¯
]
G
α
¯
(
λ
)
=
1
c
α
¯
T
α
¯
ρ
α
¯
[
T
α
¯
ρ
α
¯
f
α
¯
(
1
,
λ
)
sinh
⁡
λ
ρ
α
¯
-
1
λ
R
^
α
¯
(
λ
)
cosh
⁡
λ
ρ
α
¯
]
.
Note that(21)
f
α
¯
(
1
,
λ
)
=
f
α
¯
∘
β
(
0
,
λ
)
=
c
α
¯
∘
β
F
α
¯
∘
β
(
λ
)
,
β
=
1,2
,
…
,
m
α
¯
,
R
^
α
¯
(
λ
)
=
∑
β
=
1
m
α
¯
T
α
¯
∘
β
f
x
α
¯
∘
β
(
0
,
λ
)
=
-
λ
∑
β
=
1
m
α
¯
T
α
¯
∘
β
ρ
α
¯
∘
β
c
α
¯
∘
β
G
α
¯
∘
β
(
λ
)
.
Substituting the above into (20), we have
(22)
F
α
¯
(
λ
)
=
1
c
α
¯
T
α
¯
ρ
α
¯
[
∑
β
=
1
m
α
¯
T
α
¯
∘
β
ρ
α
¯
∘
β
c
α
¯
∘
β
G
α
¯
∘
β
(
λ
)
T
α
¯
ρ
α
¯
c
α
¯
∘
β
F
α
¯
∘
β
(
λ
)
cosh
⁡
λ
ρ
α
¯
k
k
k
k
k
k
k
k
+
sinh
⁡
λ
ρ
α
¯
∑
β
=
1
m
α
¯
T
α
¯
∘
β
ρ
α
¯
∘
β
c
α
¯
∘
β
G
α
¯
∘
β
(
λ
)
]
G
α
¯
(
λ
)
=
1
c
α
¯
T
α
¯
ρ
α
¯
[
∑
β
=
1
m
α
¯
T
α
¯
∘
β
ρ
α
¯
∘
β
c
α
¯
∘
β
G
α
¯
∘
β
(
λ
)
T
α
¯
ρ
α
¯
c
α
¯
∘
β
F
α
¯
∘
β
(
λ
)
sinh
⁡
λ
ρ
α
¯
k
k
k
k
k
k
k
k
+
cosh
⁡
λ
ρ
α
¯
∑
β
=
1
m
α
¯
T
α
¯
∘
β
ρ
α
¯
∘
β
c
α
¯
∘
β
G
α
¯
∘
β
(
λ
)
]
.
Take(23)
c
α
¯
∘
β
=
c
α
¯
∏
γ
=
1
,
γ
≠
β
m
α
¯
‍
F
α
¯
∘
γ
(
λ
)
,
β
=
1,2
,
…
,
m
α
¯
,
then (22) yields
(24)
F
α
¯
(
λ
)
=
1
T
α
¯
ρ
α
¯
[
T
α
¯
ρ
α
¯
∏
β
=
1
m
α
¯
F
α
¯
∘
β
(
λ
)
cosh
⁡
λ
ρ
α
¯
k
k
k
k
k
k
k
k
+
sinh
⁡
λ
ρ
α
¯
∑
β
=
1
m
α
¯
T
α
¯
∘
β
ρ
α
¯
∘
β
G
α
¯
∘
β
(
λ
)
k
k
k
k
k
k
k
k
×
∏
γ
=
1
,
γ
≠
β
m
α
¯
F
α
¯
∘
γ
(
λ
)
]
G
α
¯
(
λ
)
=
1
T
α
¯
ρ
α
¯
[
T
α
¯
ρ
α
¯
∏
β
=
1
m
α
¯
F
α
¯
∘
β
(
λ
)
sinh
⁡
λ
ρ
α
¯
k
k
k
k
k
k
k
k
+
cosh
⁡
λ
ρ
α
¯
∑
β
=
1
m
α
¯
T
α
¯
∘
β
ρ
α
¯
∘
β
G
α
¯
∘
β
(
λ
)
k
k
k
k
k
k
k
k
×
∏
γ
=
1
,
γ
≠
β
m
α
¯
F
α
¯
∘
γ
(
λ
)
]
.
Thus we have proved (16) in the theorem.
Sinceinf
⁡
σ
∈
ℝ
|
F
α
¯
∘
β
(
i
σ
)
|
>
0, β
=
1,2
,
…
,
m
α
¯, we derive from (16) that
(25)
T
α
¯
ρ
α
¯
F
α
¯
(
i
σ
)
∏
β
=
1
m
α
¯
F
α
¯
∘
β
(
i
σ
)
=
T
α
¯
ρ
α
¯
cos
⁡
σ
ρ
α
¯
+
i
sin
σ
ρ
α
¯
∑
β
=
1
m
α
¯
T
α
¯
∘
β
ρ
α
¯
∘
β
G
α
¯
∘
β
(
i
σ
)
F
α
¯
∘
β
(
i
σ
)
=
T
α
¯
ρ
α
¯
cos
⁡
σ
ρ
α
¯
-
sin
σ
ρ
α
¯
∑
β
=
1
m
α
¯
T
α
¯
∘
β
ρ
α
¯
∘
β
ℑ
(
G
α
¯
∘
β
(
i
σ
)
F
α
¯
∘
β
(
i
σ
)
¯
)
|
F
α
¯
∘
β
(
i
σ
)
|
2
+
i
sin
σ
ρ
α
¯
∑
β
=
1
m
α
¯
T
α
¯
∘
β
ρ
α
¯
∘
β
ℜ
(
G
α
¯
∘
β
(
i
σ
)
F
α
¯
∘
β
(
i
σ
)
¯
)
|
F
α
¯
∘
β
(
i
σ
)
|
2
.
Thus the assumption (17) reads that
(26)
inf
⁡
σ
∈
ℝ
|
F
α
¯
(
i
σ
)
|
>
0
.
To complete the proof, we deduce from (16) that
(27)
(
T
α
¯
ρ
α
¯
)
2
G
α
¯
(
i
σ
)
F
α
¯
(
i
σ
)
¯
=
[
i
T
α
¯
ρ
α
¯
∏
β
=
1
m
α
¯
F
α
¯
∘
β
(
i
σ
)
sin
σ
ρ
α
¯
k
k
k
+
cos
⁡
λ
ρ
α
¯
∑
β
=
1
m
α
¯
T
α
¯
∘
β
ρ
α
¯
∘
β
G
α
¯
∘
β
(
i
σ
)
∏
γ
=
1
,
γ
≠
β
m
α
¯
F
α
¯
∘
γ
(
i
σ
)
]
·
[
T
α
¯
ρ
α
¯
∏
β
=
1
m
α
¯
F
α
¯
∘
β
(
i
σ
)
¯
cos
⁡
σ
ρ
α
¯
k
k
k
k
-
i
sin
σ
ρ
α
¯
∑
β
=
1
m
α
¯
T
α
¯
∘
β
ρ
α
¯
∘
β
G
α
¯
∘
β
(
i
σ
)
¯
∏
γ
=
1
,
γ
≠
β
m
α
¯
F
α
¯
∘
γ
(
i
σ
)
¯
]
=
i
2
(
T
α
¯
ρ
α
¯
)
2
∏
β
=
1
m
α
¯
|
F
α
¯
∘
β
(
i
σ
)
|
2
sin
2
σ
ρ
α
¯
-
i
2
sin
2
σ
ρ
α
¯
∑
β
=
1
m
α
¯
|
T
α
¯
∘
β
ρ
α
¯
∘
β
G
α
¯
∘
β
(
i
σ
)
∏
γ
=
1
,
γ
≠
β
m
α
¯
F
α
¯
∘
γ
(
i
σ
)
|
2
+
T
α
¯
ρ
α
¯
sin
2
σ
ρ
α
¯
∏
β
=
1
m
α
¯
F
α
¯
∘
β
(
i
σ
)
×
∑
β
=
1
m
α
¯
T
α
¯
∘
β
ρ
α
¯
∘
β
G
α
¯
∘
β
(
i
σ
)
¯
∏
γ
=
1
,
γ
≠
β
m
α
¯
F
α
¯
∘
γ
(
i
σ
)
¯
+
T
α
¯
ρ
α
¯
cos
⁡
2
σ
ρ
α
¯
∏
β
=
1
m
α
¯
F
α
¯
∘
β
(
i
σ
)
¯
×
∑
β
=
1
m
α
¯
T
α
¯
∘
β
ρ
α
¯
∘
β
G
α
¯
∘
β
(
i
σ
)
∏
γ
=
1
,
γ
≠
β
m
α
¯
F
α
¯
∘
γ
(
i
σ
)
:
=
Θ
1
+
Θ
2
+
Θ
3
+
Θ
4
,
where ℜ
Θ
1
=
ℜ
Θ
2
=
0, and
(28)
ℜ
(
Θ
3
+
Θ
4
)
=
T
α
¯
ρ
α
¯
∑
β
=
1
m
α
¯
T
α
¯
∘
β
ρ
α
¯
∘
β
×
ℜ
(
G
α
¯
∘
β
(
i
σ
)
¯
F
α
¯
∘
β
(
i
σ
)
sin
2
σ
ρ
α
¯
k
+
G
α
¯
∘
β
(
i
σ
)
F
α
¯
∘
β
(
i
σ
)
¯
cos
⁡
2
σ
ρ
α
¯
)
×
∏
γ
=
1
,
γ
≠
β
m
α
¯
|
F
α
¯
∘
γ
(
i
σ
)
|
2
=
T
α
¯
ρ
α
¯
∑
β
=
1
m
α
¯
T
α
¯
∘
β
ρ
α
¯
∘
β
ℜ
(
G
α
¯
∘
β
(
i
σ
)
F
α
¯
∘
β
(
i
σ
)
¯
)
×
∏
γ
=
1
,
γ
≠
β
m
α
¯
|
F
α
¯
∘
γ
(
i
σ
)
|
2
.
Thus(29)
ℜ
(
G
α
¯
(
i
σ
)
F
α
¯
(
i
σ
)
¯
)
=
1
T
α
¯
ρ
α
¯
∑
β
=
1
m
α
¯
T
α
¯
∘
β
ρ
α
¯
∘
β
×
ℜ
(
G
α
¯
∘
β
(
i
σ
)
F
α
¯
∘
β
(
i
σ
)
¯
)
∏
γ
=
1
,
γ
≠
β
m
α
¯
|
F
α
¯
∘
γ
(
i
σ
)
|
2
.
Therefore the assumption (17) together with the equality above implies that
(30)
ℜ
(
G
α
¯
(
i
σ
)
F
α
¯
(
i
σ
)
¯
)
>
0
.
The proof is then complete.Remark 2.
In (11), there is no interior controller. If a velocity feedback controller is equipped at some interior vertex v
α
¯, α
¯
∈
𝒥
M, then we only have to add a term -
a
α
¯
y
t
α
¯
(
1
,
t
) to R
α
¯
(
t
), where a
α
¯ is the nonnegative feedback gain constant, and -
a
α
¯
y
t
α
¯
(
1
,
t
) can be regarded as a branching damping. In this situation, the results in Theorem 1 can be proved similarly. The only difference is that the recursive expression (16) is modified by
(31)
F
α
¯
(
λ
)
=
1
T
α
¯
ρ
α
¯
[
w
α
¯
(
λ
)
∏
β
=
1
m
α
¯
F
α
¯
∘
β
(
λ
)
k
k
k
k
k
k
+
sinh
⁡
λ
ρ
α
¯
∑
β
=
1
m
α
¯
T
α
¯
∘
β
ρ
α
¯
∘
β
G
α
¯
∘
β
(
λ
)
k
k
k
k
k
k
×
∏
γ
=
1
,
γ
≠
β
m
α
¯
F
α
¯
∘
γ
(
λ
)
]
G
α
¯
(
λ
)
=
1
T
α
¯
ρ
α
¯
[
v
α
¯
(
λ
)
∏
β
=
1
m
α
¯
F
α
¯
∘
β
(
λ
)
k
k
k
k
k
k
+
cosh
⁡
λ
ρ
α
¯
∑
β
=
1
m
α
¯
T
α
¯
∘
β
ρ
α
¯
∘
β
G
α
¯
∘
β
(
λ
)
k
k
k
k
k
k
×
∏
γ
=
1
,
γ
≠
β
m
α
¯
F
α
¯
∘
γ
(
λ
)
]
,
where
(32)
w
α
¯
(
λ
)
=
T
α
¯
ρ
α
¯
cosh
⁡
λ
ρ
α
¯
+
a
α
¯
sinh
⁡
λ
ρ
α
¯
v
α
¯
(
λ
)
=
T
α
¯
ρ
α
¯
sinh
⁡
λ
ρ
α
¯
+
a
α
¯
cosh
⁡
λ
ρ
α
¯
,
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
a
α
¯
≥
0
.Now we consider the special case that𝒜 is a tree-shaped network. We apply Theorem 1 to its stability analysis. Herein, the boundary conditions are set to be
(33)
T
α
¯
y
x
α
¯
(
1
,
t
)
=
R
α
¯
(
t
)
,
α
¯
∈
𝒥
S
.If the rootℛ is clamped, that is, y
(
0
,
t
)
=
0, then we get from its Laplacian transform f
(
0
,
λ
)
=
0 that F
(
λ
)
=
0, which is the characteristic equation of 𝒜. Using (16) in Theorem 1, we obtain the recursive expression of F
(
λ
). Moreover, if we choose appropriate R
α
¯
(
t
)
,
α
¯
∈
𝒥
S such that
(34)
inf
⁡
σ
∈
ℝ
|
F
α
¯
(
i
σ
)
|
>
0
,
ℜ
(
G
α
¯
(
i
σ
)
F
α
¯
(
i
σ
)
¯
)
>
0
,
∀
α
¯
∈
𝒥
S
,
then the second assertion of Theorem 1 says that
(35)
inf
⁡
σ
∈
ℝ
|
F
(
i
σ
)
|
>
0
.
The inequality (35) guarantees the exponential stability of 𝒜, provided that 𝒜 satisfies SDGA. If (34)-(35) are not true, the exponential stability of 𝒜 no longer holds. However, we can still obtain the asymptotic stability of it as long as F
(
i
σ
)
≠
0, σ
∈
ℝ.If the rootℛ is free, we can similarly analyze the stability. In fact, we can prove similarly that if
(36)
inf
⁡
σ
∈
ℝ
|
G
α
¯
∘
β
(
i
σ
)
|
>
0
,
ℜ
(
G
α
¯
∘
β
(
i
σ
)
F
α
¯
∘
β
(
i
σ
)
¯
)
>
0
hold for every β
=
1,2
,
…
,
m
α
¯, then
(37)
inf
⁡
σ
∈
ℝ
|
G
α
¯
(
i
σ
)
|
>
0
,
ℜ
(
G
α
¯
(
i
σ
)
F
α
¯
(
i
σ
)
¯
)
>
0
.In the same way, Theorem1 can also be applied for the stability analysis of subtrees of complex networks.


## 3. Application Example

In this section, we will give an example of 12-edge tree-shaped wave network as an application to Theorem1. We will see that the stability analysis of the controlled system can be simply carried out by our approach.The structure of the network is shown in the rightmost of Figure2. It is viewed as branching out from a single edge; the process is divided into 4 steps (also see Figure 2).Figure 2
The “branching” process of the 12-edge tree-shaped wave network.The network is described by(38)
T
α
¯
y
x
x
α
¯
(
x
,
t
)
=
m
α
¯
y
t
t
α
¯
(
x
,
t
)
,
x
∈
(
0,1
)
,
t
>
0
,
α
¯
∈
𝒥
y
(
0
,
t
)
=
0
y
α
¯
(
1
,
t
)
=
y
α
¯
∘
β
(
0
,
t
)
,
α
¯
∈
𝒥
M
,
β
=
1,2
T
α
¯
y
x
α
¯
(
1
,
t
)
=
R
α
¯
(
t
)
,
α
¯
∈
𝒥
M
R
α
¯
(
t
)
=
∑
β
=
1
2
T
α
¯
∘
β
y
x
α
¯
∘
β
(
0
,
t
)
,
α
¯
∈
𝒥
M
T
α
¯
y
x
α
¯
(
1
,
t
)
=
R
α
¯
(
t
)
,
α
¯
∈
𝒥
S
,
where 𝒥
M
=
{
empty
index
,
(
1
)
,
(
2
)
,
(
3
)
,
(
1,1
)
}, 𝒥
S
=
{
(
1,2
)
,
(
2,1
)
,
(
2,2
)
,
(
3,1
)
,
(
3,2
)
,
(
1,1
,
1
)
,
(
1,1
,
2
)
}, and 𝒥
=
𝒥
M
∪
𝒥
S.Since no edge branches out fromv
α
¯
,
α
¯
∈
𝒥
S, we take the extra forces to be the collocated velocity feedback control:
(39)
R
α
¯
(
t
)
=
-
a
α
¯
y
t
α
¯
(
1
,
t
)
,
a
α
¯
>
0
,
α
¯
∈
𝒥
S
.Thus the 12-edge tree-shaped wave network is described by (38)-(39). A standard regulation of (38)-(39) results in its vector-valued form, which falls into the model discussed in [14]. Thus all the theoretic results in [14] hold, among which we emphasize that(1)
the system is well-posed;
(2)
λ
≠
0; and
(3)
ifa
α
¯
≠
T
α
¯
/
m
α
¯, α
¯
∈
𝒥
S, then the SDGA satisfies.Moreover, we have the following assertion.Proposition 3.
LetD
(
λ
)
,
λ
∈
ℂ be the characteristic equation of (38)-(39). Suppose that a
α
¯
≠
T
α
¯
/
m
α
¯
,
α
¯
∈
𝒥
S. If inf
⁡
σ
∈
ℝ
⁡
|
D
(
i
σ
)
|
≠
0, then the system (38)-(39) is exponentially stable.In the sequel, we will deduce the characteristic equation by our approach shown in Section2, prove inf
⁡
σ
∈
ℝ
⁡
|
D
(
i
σ
)
|
≠
0 by Theorem 1, and obtain the exponential stability of system (38)-(39) according to Proposition 3.To this end, we first derive the recursive expressions of the coefficientsF
α
¯
(
λ
) and G
α
¯
(
λ
)
,
α
¯
∈
𝒥 in its branching order. Taking the Laplacian transform of (38)-(39), one has
(40)
f
α
¯
(
x
,
λ
)
=
c
α
¯
[
F
α
¯
(
λ
)
cosh
⁡
λ
ρ
α
¯
x
-
G
α
¯
(
λ
)
sinh
⁡
λ
ρ
α
¯
x
]
,
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
ρ
α
¯
=
m
α
¯
T
α
¯
,
α
¯
∈
𝒥
,
(41)
f
(
0
,
λ
)
=
0
,
(42)
f
(
1
,
λ
)
=
f
(
k
)
(
0
,
λ
)
T
f
x
(
1
,
λ
)
=
R
^
(
λ
)
=
∑
k
=
1
3
T
(
k
)
f
x
(
k
)
(
0
,
λ
)
f
(
k
)
(
1
,
λ
)
=
f
(
k
,
j
)
(
0
,
λ
)
T
(
k
)
f
x
(
k
)
(
1
,
λ
)
=
R
^
(
k
)
(
λ
)
=
∑
j
=
1
2
T
(
k
,
j
)
f
x
(
k
,
j
)
(
0
,
λ
)
f
(
1,1
)
(
1
,
λ
)
=
f
(
1,1
,
l
)
(
0
,
λ
)
T
(
1,1
)
f
x
(
1,1
)
(
1
,
λ
)
=
R
^
(
1,1
)
(
λ
)
=
∑
l
=
1
2
T
(
1,1
,
l
)
f
x
(
1,1
,
l
)
(
0
,
λ
)
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
=
1,2
,
3
,
j
,
l
=
1,2
,
(43)
T
α
¯
f
x
α
¯
(
1
,
λ
)
=
R
^
α
¯
(
λ
)
=
-
λ
a
α
¯
f
α
¯
(
1
,
λ
)
,
α
¯
∈
𝒥
S
.Firstly, substituting (40) into (41) indicates
(44)
F
(
λ
)
=
0
which is the characteristic equation of the system (38)-(39) to be determined.Then, according to (16) in Theorem 1, we get the recursive expressions of F
α
¯
(
λ
) and G
α
¯
(
λ
)
,
α
¯
∈
𝒥
M:
(45)
F
(
λ
)
=
1
T
ρ
[
∏
j
=
1
,
j
≠
k
3
F
(
j
)
(
λ
)
T
ρ
cosh
⁡
λ
ρ
∏
k
=
1
3
F
(
k
)
(
λ
)
+
sinh
⁡
λ
ρ
∑
k
=
1
3
T
(
k
)
ρ
(
k
)
G
(
k
)
(
λ
)
×
∏
j
=
1
,
j
≠
k
3
F
(
j
)
(
λ
)
]
G
(
λ
)
=
1
T
ρ
[
∏
j
=
1
,
j
≠
k
3
F
(
j
)
(
λ
)
T
ρ
sinh
⁡
λ
ρ
∏
k
=
1
3
F
(
k
)
(
λ
)
+
cosh
⁡
λ
ρ
∑
k
=
1
3
T
(
k
)
ρ
(
k
)
G
(
k
)
(
λ
)
×
∏
j
=
1
,
j
≠
k
3
F
(
j
)
(
λ
)
]
,
(46)
F
(
k
)
(
λ
)
=
1
T
(
k
)
ρ
(
k
)
×
[
T
(
k
)
ρ
(
k
)
cosh
⁡
λ
ρ
(
k
)
∏
j
=
1
2
F
(
k
,
j
)
(
λ
)
k
k
k
+
sinh
⁡
λ
ρ
(
k
)
∑
j
=
1
2
T
(
k
,
j
)
ρ
(
k
,
j
)
G
(
k
,
j
)
(
λ
)
k
k
k
×
∏
l
=
1
,
l
≠
j
2
F
(
k
,
l
)
(
λ
)
]
G
(
k
)
(
λ
)
=
1
T
(
k
)
ρ
(
k
)
×
[
T
(
k
)
ρ
(
k
)
sinh
⁡
λ
ρ
(
k
)
∏
j
=
1
2
F
(
k
,
j
)
(
λ
)
k
k
k
+
cosh
⁡
λ
ρ
(
k
)
∑
j
=
1
2
T
(
k
,
j
)
ρ
(
k
,
j
)
G
(
k
,
j
)
(
λ
)
k
k
k
×
∏
l
=
1
,
l
≠
j
2
F
(
k
,
l
)
(
λ
)
]
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
=
1,2
,
3
,
(47)
F
(
1,1
)
(
λ
)
=
1
T
(
1,1
)
ρ
(
1,1
)
×
[
∏
j
=
1
,
j
≠
k
2
F
(
1,1
,
j
)
(
λ
)
T
(
1,1
)
ρ
(
1,1
)
cosh
⁡
λ
ρ
(
1,1
)
∏
k
=
1
2
F
(
1,1
,
k
)
(
λ
)
k
k
k
+
sinh
⁡
λ
ρ
(
1,1
)
∑
k
=
1
2
T
(
1,1
,
k
)
ρ
(
1,1
,
k
)
G
(
1,1
,
k
)
(
λ
)
k
k
k
×
∏
j
=
1
,
j
≠
k
2
F
(
1,1
,
j
)
(
λ
)
]
G
(
1,1
)
(
λ
)
=
1
T
(
1,1
)
ρ
(
1,1
)
×
[
∏
j
=
1
,
j
≠
k
2
F
(
1,1
,
j
)
(
λ
)
T
(
1,1
)
ρ
(
1,1
)
sinh
⁡
λ
ρ
(
1,1
)
∏
k
=
1
2
F
(
1,1
,
k
)
(
λ
)
k
k
k
+
cosh
⁡
λ
ρ
(
1,1
)
∑
k
=
1
2
T
(
1,1
,
k
)
ρ
(
1,1
,
k
)
G
(
1,1
,
k
)
(
λ
)
k
k
k
×
∏
j
=
1
,
j
≠
k
2
F
(
1,1
,
j
)
(
λ
)
]
.Finally, we substitute (40) into (43) to get
(48)
F
α
¯
(
λ
)
(
T
α
¯
ρ
α
¯
sinh
⁡
λ
ρ
α
¯
+
α
α
¯
cosh
⁡
λ
ρ
α
¯
)
-
G
α
¯
(
λ
)
(
T
α
¯
ρ
α
¯
cosh
⁡
λ
ρ
α
¯
+
α
α
¯
sinh
⁡
λ
ρ
α
¯
)
=
0
,
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
α
¯
∈
𝒥
S
.
Thus we can take
(49)
F
α
¯
(
λ
)
=
T
α
¯
ρ
α
¯
cosh
⁡
λ
ρ
α
¯
+
a
α
¯
sinh
⁡
λ
ρ
α
¯
G
α
¯
(
λ
)
=
T
α
¯
ρ
α
¯
sinh
⁡
λ
ρ
α
¯
+
a
α
¯
cosh
⁡
λ
ρ
α
¯
,
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
α
¯
∈
𝒥
S
.By now, one has obtained all the recursive expressions (44)–(49) of the coefficients F
α
¯
(
λ
) and G
α
¯
(
λ
)
,
α
¯
∈
𝒥, which form the characteristic equation of (38)-(39). Now we will prove the exponential stability of the system (38)-(39). Proposition 3 says it suffices to prove that F
(
λ
) defined in (45) satisfies inf
⁡
σ
∈
ℝ
⁡
|
F
(
i
σ
)
|
>
0. We will verify this by easily estimating the infimums of (44)–(49) in the inverse order.By (49), a direct calculation leads to
(50)
F
α
¯
(
i
σ
)
=
T
α
¯
ρ
α
¯
cos
⁡
σ
ρ
α
¯
+
i
a
α
¯
sin
σ
ρ
α
¯
,
G
α
¯
(
i
σ
)
=
i
T
α
¯
ρ
α
¯
sin
σ
ρ
α
¯
+
a
α
¯
cos
⁡
σ
ρ
α
¯
.
Thus one has
(51)
inf
⁡
σ
∈
ℝ
|
F
α
¯
(
i
σ
)
|
=
inf
⁡
σ
∈
ℝ
|
G
α
¯
(
i
σ
)
|
=
min
⁡
{
a
α
¯
,
T
α
¯
ρ
α
¯
}
>
0
for α
¯
∈
𝒥
S.Moreover, it holds that(52)
ℜ
(
G
α
¯
(
i
σ
)
F
α
¯
(
i
σ
)
¯
)
=
ℜ
(
(
T
α
¯
ρ
α
¯
cos
⁡
σ
ρ
α
¯
-
i
a
α
¯
sin
σ
ρ
α
¯
)
k
k
k
k
×
(
i
T
α
¯
ρ
α
¯
sin
σ
ρ
α
¯
+
a
α
¯
cos
⁡
σ
ρ
α
¯
)
)
=
a
α
¯
T
α
¯
ρ
α
¯
>
0
,
α
¯
∈
𝒥
S
.Thus Theorem1, (34)-(35), and (51)-(52) read that
(53)
inf
⁡
σ
∈
ℝ
|
F
(
λ
)
|
>
0
.Therefore Proposition3 indicates the theorem below.Theorem 4.
Let the 12-edge tree-shaped wave network be given by (38)-(39). If a
α
¯
≠
T
α
¯
/
m
α
¯, α
¯
∈
𝒥
S, then it is exponentially stable.Remark 5.
Theorem4 implies that the tree-shaped wave network is exponentially stabilized by boundary collocated velocity feedback control. We can also study the stability of the system under other control laws. As stated in Remark 2, by adding the term -
a
α
¯
y
t
α
¯
(
1
,
t
) to R
α
¯
(
t
), α
¯
∈
𝒥
M, the characteristic equation of the modified controlled system can be obtained similarly and its stability can be analyzed in the same way. For instance, if T
α
¯
=
ρ
α
¯
=
1, α
¯
∈
𝒥, fewer controllers are needed for the system to be exponentially stable. This is shown in Figure 3, where the dots stand for the controllers. We omit the details.Figure 3
Fewer controllers are needed for the exponential stability.


## 4. Conclusion

In this paper, we introduce a new approach for the stability analysis of wave networks. Viewing the network as branching out from a single edge, we divide this process into several steps. Then we get the recursive expressions of the Laplacian transform of the adjacent edges of the system in its branching order. These recursive expressions together form the characteristic equation of the system. In the stability analysis, we estimate the infimums of these recursive expressions in the inverse order. Then the exponential stability of the system can be finished in an easy manner. This approach is valid for tree-shaped networks and subtrees of complex networks. For the stability analysis of networks with circuits, it needs some improvements, which will be studied in the future.



---

*Source: 724512-2014-05-11.xml*