请写一下用脑电mat文件求DE特征的代码

以下是用MATLAB实现用脑电mat文件求DE特征的代码：

% 读取mat文件
load('eeg_data.mat');

% 设置参数
fs = 128; % 采样频率
window_size = 2*fs; % 窗口大小
window_shift = fs/2; % 窗口移动步长

% 计算DE特征
eeg_data = eeg_data'; % 转置，使得每一行代表一个通道
num_channels = size(eeg_data, 1); % 通道数
num_windows = floor((size(eeg_data, 2) - window_size) / window_shift) + 1; % 窗口数
de_feature = zeros(num_windows, num_channels); % 存储DE特征

for i = 1:num_windows
    % 提取当前窗口的数据
    window_data = eeg_data(:, (i-1)*window_shift+1 : (i-1)*window_shift+window_size);
    
    % 对每个通道进行DE特征计算
    for j = 1:num_channels
        channel_data = window_data(j,:);
        de_feature(i,j) = sum(channel_data.^2) / window_size;
    end
end

% 可视化DE特征
figure;
plot(de_feature(:,1)); % 取第一个通道的DE特征进行可视化
title('DE Feature of Channel 1');
xlabel('Window Index');
ylabel('DE Value');

注意：这里的DE特征计算方法是简单地将每个窗口内的每个通道的平方和除以窗口大小，实际应用中可能需要更加复杂的方法来提取DE特征。此外，这里的代码只用了单个通道的数据，实际应用中可能需要对多个通道进行联合处理。