Solutions for Exercises in A First Course in Fuzzy and Neural Control by Nguyen & Prasad
Solutions for Exercises in A First Course in Fuzzy and Neural Control by Nguyen & Prasad
mulate a high-level system description. Given this breadth for problem solving,
students are expected to develop the hierarchical structure of any system.
The fact is that there is a diversity of answers, and the advice is to let
students brainstorm their ideas. These exercises are meant to expand student
thinking about sensors and possible feedback architectures that will improve
overall controllability and observability of systems.
ail
dynamics of water discharged from the kettle to increase in mass of the bucket.
MATHEMATICAL
MODELS IN CONTROL
3
4 CHAPTER 2. MATHEMATICAL MODELS IN CONTROL
Exercise 2
For a linear system the input v(t) produces an output y(t). Let y(t) =
Bej( t+θ) . Expressing the output in terms of real and imaginary quanti-
ties yields
y(t) = Re B cos( t + θ) + Im B sin( t + θ)
We can conclude therefore that an input Re v(t) = Re A cos( t + φ) pro-
duces an output Re y(t) = Re B cos( t + θ).
b. The response to an input v(t) is y(t). For a linear system this can be
expressed as y(t) = Av(t), where A is a linear operator. Therefore
dy(t) dv(t)
=A
dt dt
Exercise 3
or that Rank [O0 ] = 1 < 2, proves the system is not fully observable.
Note: Finding the rank of the matrix does not by itself explain why
the system is not observable. To understand the reason why, we need
to examine the transfer function as requested in (c).
(c) The transfer function of a linear time-invariant system is given by
where A, B, C, and D are the matrices obtained from the state and
output equations, and I is the identity matrix. In this system, as
well as in the others, D = 0.
½∙ ¸ ¸¾−1 ∙ ¸
∙
£ ¤ s 0 1 0 1
Y (s) = 2 1 − +0
0 s 2 2 0
½∙ ¸¾−1 ∙ ¸
£ ¤ s−1 0 1
= 2 1
−2 s−2 0
⎧⎡ ⎤⎫
⎪ (s − 2) ⎪
⎪
⎨ 0 ⎪∙ ¸
£ ¤ ⎢ (s − 1)(s − 2) ⎥⎬ 1
= 2 1 ⎢ ⎥
⎪ ⎣ 2 (s − 1) ⎦⎪ 0
⎪
⎩ ⎪
⎭
(s − 1)(s − 2) (s − 1)(s − 2)
⎧⎡ ⎤⎫
⎪ (s − 2) ⎪
⎪ ⎪
£ ¤ ⎨⎢ (s − 1)(s − 2) ⎥⎬
= 2 1 ⎢ ⎥
⎪ ⎣ 2 ⎦⎪
⎪
⎩ ⎪
⎭
(s − 1)(s − 2)
2(s − 2) 2 2(s − 2) + 2
= + =
(s − 1)(s − 2) (s − 1)(s − 2) (s − 1)(s − 2)
2(s − 1) 2
= =
(s − 1)(s − 2) (s − 2)
det [λI − A] = 0
Now
½∙ ¸ ∙ ¸¾
λ 0 1 0
det [λI − A] = det −
0 λ 2 2
½∙ ¸¾
λ−1 0
= det = (λ − 1) (λ − 2)
−2 λ−2
(λ − 1) (λ − 2) = 0
tCO2=rank(CO);
% Set up the observability matrix
OB=[C’ A’*C’ (A’^2)*C’];
% Test for full observability
% You can either use the determinant or directly compute the rank
tOB1=det(OB);
tOB2=rank(OB);
% Obtain the numerator and denominator coefficients of the transfer function
[num,den]=ss2tf(A,B,C,D,1);
% To test the stability compute the eigenvalues.
% Eigenvalues are the roots of the characteristic equation
e=eig(A);