blob: e2f1bae77e874e492bf7a0c1a3d7086e2886e6d1 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
|
s=poly(0,'s');
omegaR=7
zeta=abs(roots(poly([0.043766 0 -1 0 1],'s',"coeff")))
disp(zeta)
disp("but for zeta>0.7 , Mr does not exist, so neglect higher value")
zeta=0.2141
disp(zeta)
omegaN=omegaR/(sqrt(1-2*(zeta)^2))
disp(omegaN,"omegaN = ")
disp(omegaN*sqrt(1-2*zeta^2 + sqrt(2-4*zeta^2+4*zeta^4))," B.W. = ")
|
