Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- X_vect1 = [0.06, 0.10, 0.04, -0.14, -0.04, 0.14, 0.01, -0.15, -0.01, 0.29, 0.05, -0.25, -0.06, 0.25, 0.08, -0.19, -0.16, 0.24, 0.02, -0.03, -0.03, 0.07, 0.01, 0.01, 0.09, -0.01, -0.11, 0.00, -0.01, -0.31, -0.19, 0.14, 0.38, -0.10, -0.37, 0.02, 0.28, -0.04, -0.15, 0.07, 0.08, -0.05, -0.05, -0.05, -0.09, 0.28, 0.02, -0.19, 0.11, -0.08, -0.02, 0.05, -0.04, -0.11, 0.17, 0.11, 0.04, -0.01, 0.11, 0.01, -0.05, 0.03, -0.04, -0.09, 0.06, 0.04, -0.12, -0.02, -0.05, -0.03, 0.04, 0.17, 0.08, 0.01, 0.02, 0.16, -0.33, -0.07, 0.20, 0.14, -0.19, -0.06, 0.05, 0.09, -0.09, -0.05, 0.04, 0.08, -0.03, -0.05, 0.05, 0.00, -0.05, -0.02, 0.08, -0.09, -0.09, 0.16, 0.04, 0.06, -0.03, -0.04, 0.06, 0.08, 0.10, -0.07, -0.00, 0.08, 0.05, -0.05, -0.09, 0.13, 0.05, -0.14, -0.08, -0.06, 0.10, -0.13, -0.18, -0.22, -0.10, 0.05, -0.02, 0.02, 0.04, -0.03, -0.30, -0.06, 0.12, 0.08, -0.12, -0.20, -0.10, 0.00, -0.07, -0.16, 0.08, 0.08, 0.04, -0.09, -0.13, -0.19, 0.12, 0.21, -0.19, -0.40, -0.05, 0.41, -0.03, -0.52, 0.00, 0.16, 0.02, -0.33, -0.14, 0.29, -0.07, -0.35, -0.07, 0.17, -0.06, -0.30, -0.17, 0.24, 0.09, -0.33, -0.19, 0.14, 0.03, -0.29, 0.01, 0.06, -0.12, 0.01, -0.06, -0.03, -0.02, -0.14, -0.07, 0.15, -0.11, -0.26, -0.10, 0.20, -0.05, -0.35, -0.07, 0.14, -0.07, -0.23, -0.00, 0.02, -0.04, -0.05, 0.08, -0.02, 0.04, -0.08, -0.08, -0.03, 0.03, -0.09, -0.13, 0.06, -0.06, -0.12, -0.03, 0.01, 0.12, -0.06, 0.12, 0.10, -0.08, 0.23, -0.04, 0.01, -0.14, 0.07, 0.03, -0.17, -0.20, 0.04, 0.21, -0.32, -0.24, 0.16, 0.05, -0.26, -0.07, -0.01, -0.02, -0.07, -0.11, 0.08, -0.20, 0.11, -0.09, -0.01, 0.13, 0.09, -0.17, -0.12, 0.01, -0.01, -0.13, -0.04, 0.04, -0.11, -0.08, -0.03, -0.02, -0.05, -0.05, -0.09, 0.05, 0.00, 0.14, 0.08, 0.22, -0.05, -0.05, -0.11, 0.08, 0.00, -0.08, -0.05, 0.01, -0.09, 0.06, 0.05, -0.25, 0.13, 0.16, -0.27, 0.00, 0.07, -0.04, -0.14, 0.01, 0.04, -0.09, -0.11, 0.17, -0.03, 0.03, 0.09, -0.28, 0.03, -0.02, -0.01, -0.01, -0.02, -0.13, -0.03, 0.04, -0.06, -0.07, -0.00, -0.07, -0.09];
- Y_vect1 = [-0.31, 0.38, -0.94, -0.56, 0.14, 0.17, -0.75, -0.40, -0.08, -0.08, -0.23, -0.39, -0.52, -0.14, 0.29, -0.59, -1.28, 0.05, 0.60, -0.84, -0.77, 0.21, 0.50, -0.75, -0.67, 0.03, 0.07, -0.41, -0.39, -0.35, -0.28, -0.04, -0.11, -1.03, -0.84, 0.61, -0.18, -1.13, -0.26, 0.07, -0.14, -0.72, -0.92, -0.01, 1.01, -1.64, -1.07, -0.17, 0.63, -0.78, -0.67, -0.35, 0.13, -0.06, -0.85, -0.86, 0.00, 0.50, -1.07, -0.60, -0.02, 0.11, -0.18, -0.51, -0.64, 0.22, -0.10, -0.71, -0.30, -0.21, -0.15, -0.08, -0.30, -0.71, -0.51, 0.76, -0.42, -1.26, -0.14, 0.78, -1.07, -0.71, 0.20, -0.08, -0.37, -0.45, -0.35, -0.34, -0.04, 0.08, -1.06, -0.82, 0.83, -0.09, -1.24, -0.46, 0.46, -0.18, -0.86, -0.19, 0.00, 0.16, -0.67, -0.47, -0.49, 0.02, 0.10, -0.49, -0.65, -0.48, 0.88, -0.65, -1.65, -0.40, 0.67, -0.49, -0.98, -0.54, 0.19, 0.07, -0.94, -0.77, -0.27, 0.74, -0.63, -0.87, -0.06, -0.02, -0.15, -0.42, -0.53, -0.53, 0.42, -0.41, -0.76, -0.23, 0.10, -0.12, -0.52, -0.45, -0.35, 0.08, -0.08, -0.44, -0.61, -0.20, 0.53, -0.65, -1.28, -0.17, 0.60, -0.78, -0.62, -0.25, -0.21, 0.13, -0.20, -0.66, -0.62, 0.48, -0.11, -1.13, -0.31, 0.22, -0.32, -0.64, -0.20, -0.18, -0.03, -0.02, -0.40, -0.73, 0.04, 0.20, -0.66, -0.44, -0.21, -0.03, -0.19, -0.32, -0.46, -0.10, -0.05, -0.12, -0.56, -0.36, 0.01, 0.11, -0.76, -0.72, 0.34, -0.14, -0.58, -0.23, -0.20, -0.07, -0.10, -0.38, -0.35, -0.04, -0.38, -0.06, -0.21, -0.29, -0.14, -0.17, -0.29, -0.38, -0.08, -0.04, -0.40, -0.59, 0.02, 0.37, -1.10, -0.64, 0.10, 0.20, -0.52, -0.51, -0.15, 0.18, -0.33, -0.46, -0.31, 0.03, -0.07, -0.53, -0.34, 0.05, -0.21, -0.20, -0.36, -0.20, -0.09, -0.30, -0.44, -0.23, 0.06, -0.16, -0.62, -0.37, -0.09, 0.10, -0.23, -0.83, -0.56, 0.31, -0.03, -1.02, -0.29, -0.02, 0.12, -0.60, -0.56, -0.36, 0.27, -0.20, -0.77, -0.30, -0.10, 0.07, -0.37, -0.72, -0.36, 0.16, -0.01, -0.70, -0.44, -0.17, 0.20, -0.21, -1.03, -0.38, 0.67, -0.73, -0.77, -0.13, 0.04, -0.13, -0.34, -0.29, -0.26, 0.09, -0.30, -0.43, -0.27, -0.03, -0.10, -0.33, -0.48, -0.23, 0.09, -0.09, -0.75, -0.71, 0.11, -0.01, -0.70, -0.28];
- Z_vect1 = [ -0.32, -0.09, -0.54, -0.40, -0.53, -0.26, -0.39, -0.55, -1.67, -0.30, -0.35, -0.63, -0.56, -0.41, -0.23, -0.35, -0.19, -0.86, -0.22, -0.33, -0.82, -0.44, -0.28, -0.23, -0.64, -0.22, -0.60, -0.12, -0.33, -0.48, -1.49, -0.16, -0.25, -0.82, -0.10, -0.24, -0.21, -0.59, -0.30, -0.56, -0.21, -0.49, -0.69, -1.37, -0.25, -0.34, -0.73, -0.36, -0.46, -0.13, -0.46, -0.30, -1.53, -0.18, -0.34, -0.68, -0.50, -0.37, -0.17, -0.64, 0.01, -1.06, -0.19, -0.46, -0.41, -0.92, -0.31, -0.33, -0.69, 0.02, -0.69, -0.32, -0.52, -0.60, -0.64, -0.28, -0.28, -0.58, -0.02, -1.18, -0.21, -0.49, -0.56, -1.40, -0.35, -0.33, -0.62, 0.12, -1.06, -0.21, -0.47, -0.67, -0.77, -0.54, -0.29, -0.54, -0.35, -1.25, -0.39, -0.29, -0.55, 0.14, -1.20, -0.21, -0.35, -0.49, -0.89, -0.16, -0.23, -0.51, 0.18, -1.04, -0.29, -0.44, -0.68, -0.45, -0.51, -0.17, -0.49, -0.09, -0.91, -0.34, -0.30, -0.68, 0.04, -0.64, -0.18, -0.44, -0.83, -0.44, -0.60, -0.13, -0.47, -0.38, -1.34, -0.31, -0.28, -0.49, 0.11, -1.25, -0.05, -0.33, -0.56, -0.55, -0.46, -0.17, -0.40, -0.30, -1.10, -0.32, -0.32, -0.44, -0.08, -0.95, -0.21, -0.35, -0.53, -0.29, -0.46, -0.27, -0.46, -0.30, -0.96, -0.15, -0.41, -0.68, -0.33, -0.54, -0.32, -0.56, -0.24, -1.08, -0.35, -0.47, -0.62, -0.12, -0.56, -0.27, -0.42, -0.83, -0.66, -0.31, -0.32, -0.46, -0.69, -1.35, -0.22, -0.41, -0.52, -0.08, -0.59, -0.33, -0.44, -0.69, -0.11, -0.43, -0.26, -0.37, -0.86, -0.59, -0.31, -0.20, -0.42, -0.71, -1.05, -0.33, -0.33, -0.47, -0.02, -0.49, -0.29, -0.32, -0.69, -0.21, -0.19, -0.36, -0.39, -0.86, -0.95, -0.25, -0.37, -0.55, -0.18, -0.28, -0.29, -0.42, -0.77, -1.02, -0.26, -0.35, -0.40, -0.84, -1.11, -0.34, -0.36, -0.42, -0.59, -0.50, -0.32, -0.40, -0.67, -0.22, -0.34, -0.21, -0.36, -0.77, -0.38, -0.16, -0.20, -0.43, -0.96, -0.45, -0.36, -0.09, -0.48, -0.76, -0.53, -0.20, -0.44, -0.95, -0.43, -0.24, -0.20, -0.40, -0.91, -0.55, -0.40, -0.35, -0.60, -0.47, -0.44, -0.24, -0.41, -0.88, -0.58, -0.33, -0.25, -0.44, -0.38, -0.65, -0.23, -0.37, -0.81, -0.37, -0.44, -0.26, -0.40, -0.90, -0.79, -0.30, -0.19, -0.44, -0.41, -0.85, -0.15, -0.34, -0.71, -0.58, -0.34, -0.26];
- fs = 1/0.050; % sampling frequency
- LL = length(X_vect1);
- N=1024; % Block size
- f = fs/2*linspace(0, 1, N/2+1); % Frequency values
- dt = 1/fs;
- t = 0:50:14950;
- subplot(311)
- X = fft(X_vect1, N)/LL; % Accel. along X axis
- plot(f, 2*abs(X(1:N/2+1))); grid on;
- title('Espectro de frequência de vibrações') , shg
- subplot(312)
- Y = fft(Y_vect1, N)/LL; % Accel. along Y axis
- plot(f, 2*abs(Y(1:N/2+1))); grid on;
- subplot(313)
- Z = fft(Z_vect1, N)/LL; % Accel. along Z axis
- plot(f, 2*abs(Z(1:N/2+1))); grid on;
- figure(2)
- plot(f, 2*abs(X(1:N/2+1)), 'r', f, 2*abs(Y(1:N/2+1)), 'g' , f, 2*abs(Z(1:N/2+1)), 'b'); grid on;
- xlabel('f (Hz)')
- ylabel('Magnitude')
- legend('Acc X ','Acc Y','Acc Z')
- title('Espectro de frequência de vibrações') , shg
- figure(3)
- plot(t, (X_vect1), 'r', t, (Y_vect1), 'g' , t, (Z_vect1), 'b'); grid on;
- xlabel('t (ms)')
- ylabel('Aceleração (g)')
- legend('Acc X ','Acc Y','Acc Z')
Advertisement
Add Comment
Please, Sign In to add comment