257/CH9/EX9.26/example_9_26.sce


1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16

s=%s
sys=syslin('c',((k)*(s^2-2*s+5))/((s+2)*(s-0.5)))
clf
evans(sys)

//stability

[Ki,s]=kpure(sys)
disp(Ki)

//damping ratio 0.5

f=(-0.3+%i*0.55)
disp("k for damping ratio 0.5 is")
disp(abs(((f+2)*(f-0.5)))/(f^2-2*f+5))