blob: c84912b9d855d90b20f72d7211ce0e04456760d4 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
|
q=poly([0 10 1],'s','coeff');
G=k/q
H=1
T=G/(1+G*H)
zeta=0.5
k=100/(4*zeta^2)
disp(k," k = ")
omegaN=sqrt(k)
disp(omegaN,"omegaN = ")
Ts=4/(zeta*omegaN)
disp(Ts," Ts = ")
omegaD=omegaN*(sqrt(1-zeta^2))
Tp=%pi/omegaD
disp(Tp," Tp = ")
disp(%e^(-%pi*zeta/sqrt(1-zeta^2))," Mp = ")
|
