blob: 979ba13bd6d2c0c61c55acab3817c78aa534b2a0 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
|
//given characteristic equation we get GH
s=%s
sys=syslin('c',k/(s*(s^2+8*s+20)))
clf
evans(sys)
//stability
[Ki,s]=kpure(sys)
disp(Ki)
//damping ratio is 0.95
f=(-1.8+%i*0.55)
disp("k for damping ratio 0.95 is")
disp(abs((f*(f^2+8*f+20))))
f=(-03.6+%i*1.1)
disp("k for damping ratio 0.95 is")
disp(abs((f*(f^2+8*f+20))))
|
