.


:




:

































 

 

 

 


. ᒺ




1

 

- ᒺ, .

 

 

 

A= , B = , C = , D = [0].

 

 

, , .

1. . MATLAB :

A = [ 0 0 0 1; -9.7 0 2.2 -7.1; 6.1*exp(-3) 0 -0.12 0.99;0 0 -0.25 -0.96 ];

 

B = [0; 0; 0; -0.5];

 

C = [1 0 0 0; 0 1 0 0; 0 0 1 0; 0 0 0 1];

 

D = [0];

 

SYS = ss(A, B, C, D),

 

 

a =

x1 x2 x3 x4

x1 0 0 0 1

x2 -9.7 0 2.2 -7.1

x3 0.3037 0 -0.12 0.99

x4 0 0 -0.25 -0.96

 

 

b =

u1

x1 0

x2 0

x3 0

x4 -0.5

 

 

c =

x1 x2 x3 x4

y1 1 0 0 0

y2 0 1 0 0

y3 0 0 1 0

y4 0 0 0 1

 

 

d =

u1

y1 0

y2 0

y3 0

y4 0

 

Continuous-time model.

 

sys - .

2. . :

W=tf (SYS).

 

 

Transfer function from input to output...

-0.5 s - 0.06

#1: -----------------------------------

s^3 + 1.08 s^2 + 0.3627 s + 0.07593

 

3.55 s^2 + 4.187 s + 0.2479

#2: ---------------------------------------

s^4 + 1.08 s^3 + 0.3627 s^2 + 0.07593 s

 

-0.495 s - 0.1519

#3: -----------------------------------

s^3 + 1.08 s^2 + 0.3627 s + 0.07593

 

-0.5 s^2 - 0.06 s

#4: -----------------------------------

s^3 + 1.08 s^2 + 0.3627 s + 0.07593

 

 

.

3. .

) :

pole (SYS).

 

ans =

-0.1781 + 0.2705i

-0.1781 - 0.2705i

-0.7238

 

Matlab

 

: Matlab ().

 

1. LTI- W:

 

>> W=tf ([-.5 -0.06], [1 1.08 0.3627 0.001525])

 

 

Transfer function:

-0.5 s - 0.06

------------------------------------

s^3 + 1.08 s^2 + 0.3627 s + 0.001525

 

2. pole, zero.

 

>> pole(W)

 

ans =

-0.5379 + 0.2623i

-0.5379 - 0.2623i

-0.0043

 

 

>> zero(W)

 

ans =

-0.1200

 

 

3. step(W).

 

>>step(W)

 

step(W,12),grid,title('W')

 

 

 

matlab

 

4. impulse(w).

 

>> impulse(W)

 

5. , bode(W)

 

>>bode(W)

 

 

6. , nyquist(w)

 

>> nyquist(W)

 

solve

 

n=solve('s^3 + 1.08*s^2 + 0.3627*s + 0.001525=0')

 

n =

[ -.53787082010075255328590205528788-.26232439065009344961818894504406*i]

[ -.53787082010075255328590205528788+.26232439065009344961818894504406*i]

[ -.42583597984948934281958894242347e-2]

 

 

s=solve('s^4 + 1.08*s^3 + 0.3627*s^2 + 0.0001525*s=0')

 

s =

[ 0]

[ -.53978950739976351444531001809595-.26621946005720329112758771583220*i]

[ -.53978950739976351444531001809595+.26621946005720329112758771583220*i]

[ -.42098520047297110937996380809810e-3]

 

m=('s^4 + 1.08*s^3 + 0.3627*s^2 + 0.0001525*s=0')

 

m =

s^4 + 1.08*s^3 + 0.3627*s^2 + 0.0001525*s=0

 

p=solve(m)

 

p =

[ 0]

[ -.53978950739976351444531001809595-.26621946005720329112758771583220*i]

[ -.53978950739976351444531001809595+.26621946005720329112758771583220*i]

[ -.42098520047297110937996380809810e-3]

 

 

, :

-0,0043

-0,5379 + 0,2623/

-0,5379-0,2623/

) ,

 

damp ('sys').

 

 

Eigenvalue Damping Freq. (rad/s)

 

-5.26e-001 + 8.50e-001i 5.26e-001 1.00e+000

-5.26e-001 - 8.50e-001i 5.26e-001 1.00e+000

 

 

, . (Eigenvalue) , {Damping) - , (Freq. (rad/s)) - (/).

 

4. . :

step (SYS).

.step (SYS)

 

 





:


: 2015-05-05; !; : 508 |


:

:

, .
==> ...

1578 - | 1472 -


© 2015-2024 lektsii.org - -

: 0.022 .