帮我写一个可以对表面肌电信号分段并进行特征提取的matlab代码

好的，以下是一个可以对表面肌电信号分段并进行特征提取的 Matlab 代码：

%% 导入数据
data = load('emg_data.mat');
emg = data.emg;

%% 滤波
fs = 1000; % 采样率
fc = 20; % 截止频率
[b, a] = butter(4, fc/(fs/2), 'high'); % 高通滤波器
emg_filtered = filtfilt(b, a, emg); % 滤波

%% 信号分段
window_size = 250; % 窗口大小，单位：ms
overlap = 0.5; % 窗口重叠比例
step_size = window_size * (1 - overlap); % 步长
num_windows = floor((size(emg_filtered, 1) - window_size) / step_size) + 1; % 分段个数
windows = zeros(num_windows, window_size, size(emg_filtered, 2)); % 存储分段数据
for i = 1:num_windows
    start_idx = (i - 1) * step_size + 1;
    end_idx = start_idx + window_size - 1;
    windows(i, :, :) = emg_filtered(start_idx:end_idx, :);
end

%% 特征提取
features = zeros(num_windows, 4); % 存储特征
for i = 1:num_windows
    window = squeeze(windows(i, :, :));
    features(i, :) = rms(window); % 计算 RMS 特征
end

上述代码中，`emg_data.mat` 包含一个名为 `emg` 的矩阵，大小为 1250x4。代码中采用了高通滤波器将低频噪声滤除，然后将信号分段。每个窗口的大小为 250ms，重叠比例为 50%。最后，对每个窗口计算了 RMS 特征，将特征保存到了一个名为 `features` 的矩阵中，大小为 `num_windows`x4。