Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- import matplotlib.pyplot as mpl
- import numpy as np
- from scipy import signal
- f_sig = list(map(float, input().split(" ")))
- f_sig = list(map(float, f_sig))
- f_sig = np.fft.fft(f_sig)
- b, a = signal.butter(3, 0.05)
- zi = signal.lfilter_zi(b, a)
- zi = zi * int(f_sig[0])
- av_sig, _ = signal.lfilter(b, a,f_sig,-1,zi)
- zi=zi*int(av_sig[0])
- av_sig2, _ = signal.lfilter(b,a,av_sig,-1,zi)
- av_sig2 = list(map(float, av_sig2))
- for i in range(len(av_sig2)):
- if (av_sig2[i] < 0):
- av_sig2[i] = 0 - abs(av_sig2[i])
- else:
- av_sig2[i] = abs(av_sig2[i])
- mpl.plot(av_sig2)
- mpl.show()
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement