Subversion Repositories svnkaklik

Rev

Go to most recent revision | Details | Last modification | View Log

Rev Author Line No. Line
1040 kaklik 1
\documentclass[12pt,notitlepage,fleqn]{article}
2
 
3
\usepackage[czech]{babel}
4
\usepackage[pdftex]{graphicx}
5
\usepackage{fancyhdr,multicol} %nastavení češtiny, fancy, grafiky, sloupce
6
\usepackage[utf8]{inputenc} %vstupni soubory v kodovani UTF-8
7
\usepackage[a4paper,text={17cm,25cm},centering]{geometry} %nastavení okrajů
8
\usepackage{rotating}
9
 
10
\begin{document}
11
 
12
\section*{Řešení 2. zadané úlohy		-	 Jakub Kákona}
13
 
14
\begin{enumerate}
15
\item Výpočet pomocí Jordanova kanonického tvaru
16
 
17
Víme, že vlastní čísla matice A jsou kořeny charakteristického polynomu
18
 
19
$det (\lambda I - A) = (\lambda - 1) (\lambda - 2)(\lambda - 2)$
20
 
21
Vlastní čísla potom jsou $\lambda _1 = 1, \lambda _2 = 2, \lambda _3 = 2$ .
22
 
23
Dvojka je násobným kořenem, proto může tvořit více bloků. 
24
 
25
$C \lambda _23 = dim A - h(\lambda I - A) = 3 -h 
26
\left[ \begin{array}{ccc}
27
1 & -4 & -10 \\
28
 
29
 
30
\end{array}
31
\right]
32
= 2
33
\Rightarrow
34
J =  
35
\left[ \begin{array}{ccc}
36
1 & 0 & 0 \\
37
 
38
 
39
\end{array}
40
\right]
41
 $
42
 
43
Výpočet vlastních vektorů náležících vlastním číslům.
44
 
45
$(\lambda _1 I - A) v _1 = 0
46
\left[ \begin{array}{ccc}
47
 
48
 
49
 
50
\end{array}
51
\right] v_1
52
= 0
53
\Rightarrow
54
v_1 =  
55
\left[ \begin{array}{c}
56
1 \\
57
 
58
 
59
\end{array}
60
\right]
61
 $
62
 
63
 
64
$(\lambda _2 I - A) v _2 = 0
65
\left[ \begin{array}{ccc}
66
1 & -4 & -10 \\
67
 
68
 
69
\end{array}
70
\right] v_2
71
= 0
72
\Rightarrow
73
v_2 =  
74
\left[ \begin{array}{c}
75
10 \\
76
 
77
1 \\
78
\end{array}
79
\right]
80
 $
81
 
82
 
83
$(\lambda _3 I - A) v _3 = 0
84
\left[ \begin{array}{ccc}
85
1 & -4 & -10 \\
86
 
87
 
88
\end{array}
89
\right] v_3
90
= 0
91
\Rightarrow
92
v_3 =  
93
\left[ \begin{array}{c}
94
4 \\
95
1 \\
96
 
97
\end{array}
98
\right]
99
 $
100
 
101
Transformační matice vytvořená z vlastních vektorů 
102
 
103
 
104
$P = \left[v_1, v_2, v_3 \right]
105
\left[ \begin{array}{ccc}
106
1 & 10 & 4 \\
107
 
108
 
109
\end{array}
110
\right]
111
 $
112
 
113
Její inverze je 
114
 
115
$P^{-1} = \left[v_1, v_2, v_3 \right]
116
= \left[ \begin{array}{ccc}
117
1 & -4 & -10 \\
118
 
119
 
120
\end{array}
121
\right]
122
$
123
 
124
$e^{At} = P e^{At} P^{-1}
125
= \left[ 
126
\begin{array}{ccc}
127
1 & 10 & 4 \\
128
 
129
 
130
\end{array}
131
\right]
132
\left[ 
133
\begin{array}{ccc}
134
e^{t} & 0 & 0 \\
135
 
136
 
137
\end{array}
138
\right]
139
\left[ 
140
\begin{array}{ccc}
141
1 & -4 & -10 \\
142
 
143
 
144
\end{array}
145
\right]
146
= \\
147
=
148
\left[ 
149
\begin{array}{ccc}
150
e^{t} & -4e^{t}+4e^{2t} & -10e^{t}+10e^{2t} \\
151
 
152
 
153
\end{array}
154
\right]
155
=
156
\left[ 
157
\begin{array}{ccc}
158
1 & -4 & -10 \\
159
 
160
 
161
\end{array}
162
\right] 
163
e^{t} +
164
\left[ 
165
\begin{array}{ccc}
166
 
167
 
168
 
169
\end{array}
170
\right] 
171
e^{2t}
172
$
173
 
174
Metoda s použitím Laplaceovy transformace:
175
 
176
$e^{At} = L^{-1}{(sI - A)^{-1}}$
177
 
178
 
179
$(sI - A)^{-1} =
180
\left[ \begin{array}{ccc}
181
s-1 & -4 & -10 \\
182
 
183
 
184
\end{array}
185
\right]^{-1}
186
= \frac{1}{(s-1)(s-2)^2}
187
\left[ \begin{array}{ccc}
188
(s-2)^2 & -4(s-2) & -10(s-2) \\
189
 
190
 
191
\end{array}
192
\right]^{-1}=\\
193
=
194
\left[ \begin{array}{ccc}
195
1 & -4 & 10 \\
196
 
197
 
198
\end{array}
199
\right]^{-1}
200
\frac{1}{s-1}
201
+
202
\left[ \begin{array}{ccc}
203
 
204
 
205
 
206
\end{array}
207
\right]^{-1}
208
\frac{1}{s-2}
209
$
210
 
211
Z toho pak 
212
 
213
$
214
e^{At} = L^{-1} \left\{
215
\left[ \begin{array}{ccc}
216
1 & -4 & 10 \\
217
 
218
 
219
\end{array}
220
\right]^{-1}
221
\frac{1}{s-1}
222
+
223
\left[ \begin{array}{ccc}
224
 
225
 
226
 
227
\end{array}
228
\right]^{-1}
229
\frac{1}{s-2}
230
\right\}
231
=\\
232
=
233
\left[ 
234
\begin{array}{ccc}
235
1 & -4 & -10 \\
236
 
237
 
238
\end{array}
239
\right] 
240
e^{t} +
241
\left[ 
242
\begin{array}{ccc}
243
 
244
 
245
 
246
\end{array}
247
\right] 
248
e^{2t}
249
$
250
 
251
\item Dosazením definičních matic systému získáme rovnice tvaru: $\dot x =Ax, y =Cx$
252
Provedením Laplaceovy transformace pak $sX -x_0 =AX, Y =CX$
253
 
254
Z toho vyjádříme $X = (sI - A)^{-1} x_0, Y =C(sI - A)^{-1}x_0$.
255
 
256
$Y =
257
\left[ \begin{array}{ccc}
258
1 & -1 & 1 \\
259
\end{array}
260
\right]
261
\left[ \begin{array}{ccc}
262
\frac{1}{s+1} & \frac{1}{(s+1)^2} & 0 \\
263
 
264
 
265
\end{array}
266
\right] x_0 = 
267
\left[ \begin{array}{ccc}
268
\frac{1}{s+1} & \frac{-s}{(s+1)^2} & \frac{1}{s-2} \\
269
\end{array}
270
\right] x_0 $
271
 
272
$ L\{(te^{-t})\} =  \frac{-s}{(s+1)^2}$
273
 
274
$X =
275
\left[ \begin{array}{ccc}
276
\frac{1}{s+1} & \frac{1}{(s+1)^2} & 0 \\
277
 
278
 
279
\end{array}
280
\right] x_0 
281
$
282
 
283
 Stav $x(0)$ pro který bude platit požadovaná podmínka je $
284
\left[ \begin{array}{c}
285
1 \\
286
1 \\
287
 
288
\end{array}
289
\right]
290
$
291
 
292
\item Požadavkem je konstantní stav $x_0$ proto dosadíme: $ x_0 = Ax_0 + Bu(k)$ z toho $(I - A)x_0 = Bu(k)$
293
 
294
$X =
295
\left(
296
\left[ \begin{array}{cc}
297
1 & 0 \\
298
 
299
\end{array}
300
\right] -
301
\left[ \begin{array}{cc}
302
1 & 0 \\
303
 
304
\end{array}
305
\right] \right)
306
\left[ \begin{array}{c}
307
-2 \\
308
 1 \\
309
\end{array}
310
\right]=
311
\left[ \begin{array}{c}
312
1 \\
313
 2 \\
314
\end{array}
315
\right] u(k)
316
$
317
 
318
$
319
\left[ \begin{array}{c}
320
1 \\
321
 2 \\
322
\end{array}
323
\right]
324
=
325
\left[ \begin{array}{c}
326
1 \\
327
2 \\
328
\end{array}
329
\right] u(k) \rightarrow u(k) = 1(k) 
330
$
331
 
332
Potřebná sekvence vstupů pro splnění podmínek odpovídá jednotkovému skoku.  
333
 
334
\end{enumerate} 
335
 
336
 
337
\end{document}