blob: 16f41fe5104ccd7cd436e7a7f75068716f68c8ea (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
|
ss=%s
G=64/(s*(s+9.6))
H=1
T=G/(1+G*H)
omegaN=sqrt(64)
disp(omegaN,"omegaN = ")
zeta=9.6/(2*omegaN)
disp(zeta,"zeta = ")
omegaD=omegaN*(sqrt(1-zeta^2))
theta=atan(sqrt(1-zeta^2)/zeta)
disp(theta," theta = ")
c=((1-(%e^(-zeta*omegaN*t))/sqrt(1-zeta^2)*sin(omegaD*t+theta)))
disp(c, "c = ")
Tp=%pi/omegaD
disp(Tp," Tp = ")
disp(4/(zeta*omegaN)," Ts = ")
|
