单向HP滤波 python

单向HP滤波（One-sided High Pass Filtering）是一种信号处理技术，用于去除信号中的低频成分，保留高频成分。在Python中，可以使用不同的库和方法来实现单向HP滤波。

一种常用的方法是使用SciPy库中的`scipy.signal`模块。该模块提供了丰富的信号处理函数，包括单向HP滤波函数`scipy.signal.filtfilt`。

下面是一个使用`scipy.signal.filtfilt`函数进行单向HP滤波的示例代码：

import numpy as np
from scipy import signal

# 生成示例信号
t = np.linspace(0, 1, 1000)
x = np.sin(2 * np.pi * 5 * t) + np.sin(2 * np.pi * 50 * t)

# 设计滤波器
cutoff_freq = 10  # 截止频率
nyquist_freq = 0.5 * 1000  # 奈奎斯特频率
b, a = signal.butter(4, cutoff_freq / nyquist_freq, btype='high', analog=False)

# 进行滤波
filtered_x = signal.filtfilt(b, a, x)

# 绘制原始信号和滤波后的信号
import matplotlib.pyplot as plt

plt.figure()
plt.plot(t, x, label='Original Signal')
plt.plot(t, filtered_x, label='Filtered Signal')
plt.xlabel('Time')
plt.ylabel('Amplitude')
plt.legend()
plt.show()

在上述代码中，首先生成了一个示例信号`x`，然后使用`signal.butter`函数设计了一个4阶Butterworth滤波器，指定了截止频率为10Hz。最后，使用`signal.filtfilt`函数对信号进行滤波，得到滤波后的信号`filtered_x`。