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)