버터 워스 하이 패스 필터 적용

Question

데이터에 버터 워스 하이 패스 필터를 적용하는 데 문제가 있습니다. Bx 및 By 행렬에 대한 필터를 인쇄하고 싶습니다. 당신은

방법 Bx에 math.fabs()을 적용하고 By는 긍정적 인 값을 얻기 위해, 나는 모두 양과 음의 값을 가질 볼 수 있듯이?

된 Fc = 2 Hz의

나는 100 pT를 아래 값을 차단하고 싶은 내 하이 패스 필터의 경우

나는 그 요구 사항이 있습니다.

plt.ylabel('Pico Tesle [pT]') 
plt.xlabel('Time [ms]') 
plt.grid() 
plt.plot(time[51:-14], Bx[51:-14], label='Canal 1', color='r', linewidth=0.1, linestyle="-") 
plt.plot(time2[1:-14], By[1:-14], label='Canal 3', color='b', linewidth=0.1, linestyle="-") 
plt.legend(loc='upper right', frameon=False,) 

그리고 차트이다 :

내 현재 코드의 일부이다

UPDATE : I는 Bx로 및 매트릭스에 의해 절대 값을 생성하기 위해이 기능을 사용했다. (- picoTesla PT) 나는 그 덕분에받은 무엇

plt.subplot(413) 
np.absolute(fft1) 
plt.plot(time[51:-14], np.absolute(fft1), color='r', linewidth=0.1, linestyle='-') 
plt.grid() 

plt.subplot(414) 
np.absolute(fft2) 
plt.plot(time2[1:-14], np.absolute(fft2), color='b', linewidth=0.1, linestyle='-') 
plt.grid() 

내 모든 측정 (여기가 대기 방전입니다)에만 양의 값에 있다는 것입니다. 첫 번째 줄거리는 채널 1을 보여주고 두 번째 줄거리는 채널 3을 보여주고 세 번째 줄거리는 두 가지 채널 (채널)을 결합한 줄을 보여줍니다.

아니요, 100pT 미만의 모든 측정을 차단하려면 하이 패스 필터를 사용해야합니다. 어떤 아이디어?

Answer 1

이 문제를 해결하기 위해이 방법을 사용했습니다.

Wn = float(1)/HalfSampling 
b, a = signal.butter(3, Wn, 'high', analog=False) 
BxHPF = signal.filtfilt(b, a, Bxfft) 
ByHPF = signal.filtfilt(b, a, Byfft) 
plt.plot(BxTime, BxHPF, label='Canal 1', color='r', linewidth=0.5, linestyle="-") 
plt.plot(ByTime, ByHPF, label='Canal 3', color='b', linewidth=0.5, linestyle="-")