blob: f36eda6363804e57323057e4742c1762f820d2f7 (
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
|
//CHAPTER 4_ Motion and Dimensional Measurement
//Caption : Seismic pickup
// Example 10// Page 238
disp("r1=0.2;")
disp("r2=0.6 ")
disp("tou=0.05")
r1=0.2;
r2=0.6 // given
tou=0.05; //given
wn=1600 //(' enter the natural frequency=:')
disp("H1=1/sqrt((1-r1^2)^2+(2*tou*r1)^2)")
H1=1/sqrt((1-r1^2)^2+(2*tou*r1)^2);
H1_phase=-atan((2*tou*r1)/(1-r1^2))*360/(2*%pi);
disp("H1_phase=-atan((2*tou*r1)/(1-r1^2))*360/(2*%pi)")
H2=1/sqrt((1-r2^2)^2+(2*tou*r2)^2);
H2_phase=-atan((2*tou*r2)/(1-r2^2))*360/(2*%pi);
//In order to obtain the amplitude of relative displacement, transfer function must be multiplied by amplitude of the input signal and the static sensitivty of the pickup (1/wn^2) for each frequency
//amp1=H1/wn^2;
//amp2=H2/wn^2;
tou2=0.6; // given
H11=1/sqrt((1-r1^2)^2+(2*tou2*r1)^2);
H11_phase=-atan((2*tou2*r1)/(1-r1^2))*360/(2*%pi);
H22=1/sqrt((1-r2^2)^2+(2*tou2*r2)^2);
H22_phase=-atan((2*tou2*r2)/(1-r2^2))*360/(2*%pi);
//amp11=H11/wn^2;
//amp22=H22/wn^2;
printf('the magnitude of the transfer function will be %fd and %fd while the phases will shift by %fd and %fd for tou=0.05\n',H1,H2,H1_phase,H2_phase)
printf('the magnitude of the transfer function will be %fd and %fd while the phases will shift by %fd and %fd for tou=0.6\n',H11,H22,H11_phase,H22_phase)
|
