Scribd
Upload
56 views2 pages

Observer Design: System Matrices

This document describes the design of an observer for a system with system matrices A, B, C. It provides the open-loop and desired closed-loop characteristic polynomials, po and pc. It then transforms the system into observer form with matrices Aof and Cof using a transformation matrix T. It calculates the observability matrices O and Oof and uses them to find the observer gain matrix L, which designs an observer with the desired closed-loop dynamics specified by pc. Verification is done by comparing L to using the place command.

Uploaded by

yiğit_Öksüz
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
56 views2 pages

Observer Design: System Matrices

This document describes the design of an observer for a system with system matrices A, B, C. It provides the open-loop and desired closed-loop characteristic polynomials, po and pc. It then transforms the system into observer form with matrices Aof and Cof using a transformation matrix T. It calculates the observability matrices O and Oof and uses them to find the observer gain matrix L, which designs an observer with the desired closed-loop dynamics specified by pc. Verification is done by comparing L to using the place command.

Uploaded by

yiğit_Öksüz
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 2

Observer Design

System Matrices
A = [zeros(2,2),eye(2); -10, 10, -2, 2; 60, -660, 12, -12]
B = [0; 0; 1/300; -0.02]
C = [1,1,0,0]
A =
0
0
-10
60
B =

0
0
10
-660

1
0
-2
12

0
1
2
-12

0
0
0.0033
-0.0200
C =
1
1

Observer Design (basic principles)


po=poly(A) % Open-loop characteristic polynomial
a=po(2:length(po) )
po =
1.0e+003 *
0.0010
a =
1.0e+003 *
0.0140

0.0140

0.6700

1.2000

0.6700

1.2000

6.0000

6.0000

Lob=[-15,-16,-18,-20]
pc=poly(Lob) % Desired closed-loop characteristic polynomial
ad=pc(2:length(pc))
Lob =
-15
pc =

-16

-18

-20

69

1778

20280

69

1778

20280

86400

86400

ad =

Aof=[-a',[eye(3);0,0,0]]
Cof=[1,0,0,0]
Aof =
1.0e+003 *
-0.0140
-0.6700
-1.2000
-6.0000

0.0010
0
0
0

% State matrix for observer form


% Input matrix for observer form

0
0.0010
0
0

0
0
0.0010
0

Cof =
1

O=obsv(A,C)
% Observability matrix for the given realization.
Oof=obsv(Aof,Cof) % Observability matrix for the observer form.
T =O\Oof
% The transformation matrix.
O =
1
0
50
-700

1
0
-650
6700

0
1
10
-90

Oof =
1.0e+004 *
0.0001
0
-0.0014
0.0001
-0.0474
-0.0014
1.4816
-0.0474

0
0
0.0001
-0.0014

0
0
0
0.0001

T =
-0.0100
1.0100
2.1484
-16.1484

0.0016
-0.0016
-0.0054
0.0054

-0.0000
0.0000
0.0016
-0.0016

Lo=(ad-a)'
L=T*Lo

-0.0054
0.0054
-0.0100
1.0100

0
1
-10
-510

% Observer gain matrix for observer form


% Observer gain matrix for the given realization

Lo =
1.0e+004 *
0.0055
0.1108
1.9080
8.0400
L =
20.0154
34.9846
130.8615
207.1385
Using the MATLAB command place gives the same answer
place(A',C',Lob)'
place: ndigits= 15
ans =
20.0154
34.9846
130.8615
207.1385

You might also like