blob: ca0e2964c271bc7638efd3414efc2e36d909137c (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
|
s=%s
F=s^3+4*s^2+13*s+50
disp(routh_t(F))
r=coeff(F)
routh=routh_t(F)
n=length(r)
c=0;
for i=1:n
if (routh(i,1)<0)
c=c+1;
end
end
if(c>=1)
printf("system is unstable")
else printf("there are no roots on RHS")
end
syms s k
G=5*k/(s*(1+s/3)*(1+s/6)*18)
H=1
Kv=limit(s*G*H,s,0)
disp(Kv, " Kv = ")
s=%s
F=s^3+9*s^2+18*s+180
disp(routh_t(F))
r=coeff(F)
routh=routh_t(F)
n=length(r)
c=0;
for i=1:n
if (routh(i,1)<0)
c=c+1;
end
end
if(c>=1)
printf("system is unstable")
else printf("there are no roots on RHS")
end
|
