Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- [b,a] = butter(5, [0.75*2/30, 5.0*2/30], 'bandpass');
- y = filtfilt(b, a, input_signal)
- y = signal.filtfilt(b, a, input_signal, padtype = None)
- passband = [0.75*2/30, 5.0*2/30]
- b, a = scipy.signal.butter(5, passband, 'bandpass')
- b = array([ 5.49209388e-03, 0.00000000e+00, -2.74604694e-02,
- -1.97776791e-17, 5.49209388e-02, 2.47220989e-17,
- -5.49209388e-02, -1.97776791e-17, 2.74604694e-02,
- 0.00000000e+00, -5.49209388e-03])
- a = array([ 1. , -6.52098852, 19.50384534, -35.47189804,
- 43.65758795, -38.07760914, 23.83047021, -10.55670367,
- 3.16710078, -0.58122912, 0.04945954])
- b = 0.0055, 0, -0.0275, 0, 0.0549, 0, -0.0549, 0, 0.0275, 0, -0.0055
- a = 1.0000, -6.5210, 19.5038, -35.4719, 43.6576, -38.0776, 23.8305,
- -10.5567, 3.1671, -0.5812, 0.0495
- from scipy.io import loadmat
- from scipy.signal import butter, filtfilt
- from matplotlib.pyplot import plot
- signaldata = loadmat('signaldata.mat')
- input_signal = signaldata['input_signal'][0]
- passband = [0.75*2/30, 5.0*2/30]
- b, a = butter(5, passband, 'bandpass')
- y = filtfilt(b, a, input_signal)
- plot(y)
- 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]
- [b,a] = butter(5, [0.75*2/30, 5.0*2/30], 'bandpass');
- y = filtfilt(b, a, input_signal)
- from scipy.signal import butter, filtfilt
- 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]
- [b,a] = butter(5, [0.75*2/30, 5.0*2/30], 'bandpass');
- y = filtfilt(b, a, input_signal, padtype = 'odd', padlen=3*(max(len(b),len(a))-1))
Add Comment
Please, Sign In to add comment