blob: 41b8ad9ab4d3fa1d2787c163005e7f5b489420ee (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
|
//CHAPTER 4_ Motion and Dimensional Measurement
//Caption : Accelerometer
// Example 11// Page 240
disp("fn=20000")
disp("tou=0.6")
disp("f=10000 ")
fn=20000 //('enter the natural frequency of the accelerometer =:')
tou=0.6 //('enter the daping ratio of the accelerometer=:')
f=10000 //('enter the frequency at which transfer function is to be calculated=:')
r=f/fn;
H_mag=1/sqrt((1-r^2)^2+(2*tou*r)^2);
H_phase=atan((2*tou*r)/(1-r^2))*360/(2*%pi);
printf(' The magnitude is %fd and phase is %fd deg\n',H_mag,H_phase)
error=(H_mag-1)*100/1;
printf(' Error at %fd Hz is %d percent\n',f,error)
|
