blob: 6f69e6e584edb4ef1f3366acf7b3028ee7154023 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
|
s=%s
Mr=2 //given
omegaR=3 //given
zeta=roots(16*s^4 - 16*s^2 + 1) //Mr=1/(2*zeta*sqrt(1-zeta^2))
zeta(3,1)=0.933
zeta(2,1)=0.0669
disp(zeta)
omegaN=omegaR/(sqrt(1-2*(0.2588)^2))
disp(omegaN,"omegaN = ")
TF= (omegaN^2)/poly([omegaN^2 2*0.2588*omegaN 1],'s',"coeff")
disp(TF," transfer function = ")
omegaD=omegaN*sqrt(1-(0.2588))
Tr=(%pi-(atan(sqrt(1-(0.2588)^2)/0.2588)))/(omegaD)
disp(Tr,"Tr = ")
Tp=%pi/omegaD;
disp(Tp,"Tp= ")
Ts=4/(0.2588*omegaN)
disp(Ts,"Ts = ")
Tosc=2*%pi/omegaD
disp(Tosc,"Tosc = ")
N=Ts/Tosc;
disp(N,"number of oscillations = ")
Mp=%e^(-0.2588*%pi/(sqrt(1-(0.2588)^2)))
disp(Mp,"Mp = ")
|
