肌电信号特征提取matlab代码

这里给出一个简单的肌电信号特征提取的MATLAB代码示例，包括时域和频域特征的提取：

% 载入数据，假设数据存储在名为emg_signal的变量中
load('emg_signal.mat');

% 采样频率、信号长度和时间向量
fs = 1000; % 采样频率
L = length(emg_signal); % 信号长度
t = (0:L-1)/fs; % 时间向量

%% 时域特征提取
% 计算信号的均值、方差、标准差、最大值和最小值
mean_emg_signal = mean(emg_signal);
var_emg_signal = var(emg_signal);
std_emg_signal = std(emg_signal);
max_emg_signal = max(emg_signal);
min_emg_signal = min(emg_signal);

% 计算信号的斜率
slope_emg_signal = diff(emg_signal)./diff(t);

%% 频域特征提取
% 计算功率谱密度
[Pxx, f] = pwelch(emg_signal, [], [], [], fs);

% 计算信号的能量
energy_emg_signal = sum(emg_signal.^2)/L;

% 计算信号的频率特征
mean_freq_emg_signal = sum(f.*Pxx)/sum(Pxx);
median_freq_emg_signal = median(f(Pxx == max(Pxx)));
peak_freq_emg_signal = f(Pxx == max(Pxx));

% 输出特征值
disp(['Mean EMG signal: ', num2str(mean_emg_signal)]);
disp(['Variance of EMG signal: ', num2str(var_emg_signal)]);
disp(['Standard deviation of EMG signal: ', num2str(std_emg_signal)]);
disp(['Maximum value of EMG signal: ', num2str(max_emg_signal)]);
disp(['Minimum value of EMG signal: ', num2str(min_emg_signal)]);
disp(['Slope of EMG signal: ', num2str(slope_emg_signal)]);
disp(['Energy of EMG signal: ', num2str(energy_emg_signal)]);
disp(['Mean frequency of EMG signal: ', num2str(mean_freq_emg_signal)]);
disp(['Median frequency of EMG signal: ', num2str(median_freq_emg_signal)]);
disp(['Peak frequency of EMG signal: ', num2str(peak_freq_emg_signal)]);

这个代码示例展示了如何计算肌电信号的均值、方差、标准差、最大值、最小值、斜率、功率谱密度、能量和频率特征。当然，这只是一部分特征，具体的特征选择需要根据具体应用场景和目的进行选择。