summaryrefslogtreecommitdiff
path: root/257/CH7/EX7.28/example_7_28.sce
blob: 6d40eb32d3b95bf6842e921eca65be3a086f902e (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
q=poly([0 2 1],'s','coeff');
G=8/q 
H=1

omegaN=sqrt(8)
zeta=2/(2*omegaN)
disp(omegaN,"omegaN = ")
disp(zeta,"zeta = ")

syms s a
T=8/(s^2+s*(2+8*a)+8)
disp((0.7*2*sqrt(8)-2)/8, " a for zeta=0.7 is ")
disp(%e^(-%pi*zeta/sqrt(1-zeta^2))," Mp for zeta=0.3535 is  ")