blob: b6b0e522703167a0c377c0dbdbf617d3ea6ce966 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
|
Mp=30 //given
Ts=5
zeta=0.358
omegaN=4/(zeta*Ts)
disp(omegaN,"omegaN = ")
omegaD=omegaN*(sqrt(1-zeta^2))
Tp=%pi/omegaD
disp(Tp," Tp = ")
TF=omegaN^2/(s^2+2*zeta*omegaN*s+omegaN^2)
disp(TF, "transfer function = ")
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 = ")
|
