matlab脑电信号频域分析butterworth滤波器代码

以下是 Matlab 脑电信号频域分析 butterworth 滤波器代码：

% 导入数据
data = load('EEG_data.mat');
eeg = data.EEG;
fs = data.fs;

% 设计滤波器
low_freq = 1; % 低通截止频率
high_freq = 50; % 高通截止频率
order = 4; % 滤波器阶数
[b, a] = butter(order, [low_freq high_freq]/(fs/2));

% 滤波
eeg_filtered = filtfilt(b, a, eeg);

% 绘制原始信号和滤波后的信号
t = (:length(eeg)-1)/fs;
figure;
subplot(2,1,1);
plot(t, eeg);
title('原始信号');
xlabel('时间 (s)');
ylabel('电压 (uV)');
subplot(2,1,2);
plot(t, eeg_filtered);
title('滤波后的信号');
xlabel('时间 (s)');
ylabel('电压 (uV)');