Guest User

Untitled

a guest
Mar 19th, 2018
84
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 1.54 KB | None | 0 0
  1. [b,a] = butter(5, [0.75*2/30, 5.0*2/30], 'bandpass');
  2. y = filtfilt(b, a, input_signal)
  3.  
  4. y = signal.filtfilt(b, a, input_signal, padtype = None)
  5.  
  6. passband = [0.75*2/30, 5.0*2/30]
  7. b, a = scipy.signal.butter(5, passband, 'bandpass')
  8.  
  9. b = array([ 5.49209388e-03, 0.00000000e+00, -2.74604694e-02,
  10. -1.97776791e-17, 5.49209388e-02, 2.47220989e-17,
  11. -5.49209388e-02, -1.97776791e-17, 2.74604694e-02,
  12. 0.00000000e+00, -5.49209388e-03])
  13. a = array([ 1. , -6.52098852, 19.50384534, -35.47189804,
  14. 43.65758795, -38.07760914, 23.83047021, -10.55670367,
  15. 3.16710078, -0.58122912, 0.04945954])
  16.  
  17. b = 0.0055, 0, -0.0275, 0, 0.0549, 0, -0.0549, 0, 0.0275, 0, -0.0055
  18. a = 1.0000, -6.5210, 19.5038, -35.4719, 43.6576, -38.0776, 23.8305,
  19. -10.5567, 3.1671, -0.5812, 0.0495
  20.  
  21. from scipy.io import loadmat
  22. from scipy.signal import butter, filtfilt
  23. from matplotlib.pyplot import plot
  24.  
  25. signaldata = loadmat('signaldata.mat')
  26.  
  27. input_signal = signaldata['input_signal'][0]
  28.  
  29. passband = [0.75*2/30, 5.0*2/30]
  30. b, a = butter(5, passband, 'bandpass')
  31.  
  32. y = filtfilt(b, a, input_signal)
  33. plot(y)
  34.  
  35. input_signal = [1,2,3,4,5,1,2,3,4,5,1,2,3,4,1,2,3,4,5,1,2,3,4,5,1,2,3,4,5,1,2,3,4,5]
  36. [b,a] = butter(5, [0.75*2/30, 5.0*2/30], 'bandpass');
  37. y = filtfilt(b, a, input_signal)
  38.  
  39. from scipy.signal import butter, filtfilt
  40. input_signal = [1,2,3,4,5,1,2,3,4,5,1,2,3,4,1,2,3,4,5,1,2,3,4,5,1,2,3,4,5,1,2,3,4,5]
  41. [b,a] = butter(5, [0.75*2/30, 5.0*2/30], 'bandpass');
  42. y = filtfilt(b, a, input_signal, padtype = 'odd', padlen=3*(max(len(b),len(a))-1))
Add Comment
Please, Sign In to add comment