summaryrefslogtreecommitdiff
path: root/257/CH7/EX7.34/example_7_34.sce
blob: a04cdfd7fada2daabf5a91c064124694266ea022 (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
syms s k
q=poly([0 0 1],'s','coeff');
G=100/q //gain FACTOR=100
H=1+k*s
y=G*H 

T=G/(1+G*H)

omegaN=sqrt(100)           //comparing denominator with standard form
disp(omegaN,"omegaN = ")
zeta=100*k/(2*omegaN)
disp(zeta,"zeta=")

Mp=4.32;
zeta=0.7071;
disp(zeta/5," k = ")