用fieldtrip工具箱中ft_freqstatistics进行聚类检验的代码

FieldTrip是一个开源的神经科学数据分析工具箱，主要用于脑电图和功能性磁共振成像的数据分析。`ft_freqstatistics`函数主要用于对时间序列数据进行频域统计分析，例如计算功率谱密度（PSD）。如果你想在FieldTrip中进行频率范围内的聚类检验，通常涉及两个步骤：首先，对感兴趣的时间窗或通道的功率谱进行统计；然后，对这些统计结果进行聚类。

以下是一个简单的示例，假设你已经准备好了一个包含单个条件下的脑电图信号的mat文件：

% 加载数据
cfg = [];
data = ft_read_raw(cfg, 'your_file.mat');

% 设置时间窗口和频率范围
window = 50; % 时间窗长度，单位为采样点
fmin = 5;   % 聚类开始的最低频率 (Hz)
fmax = 50;  % 聚类结束的最高频率 (Hz)

% 提取功率谱并计算平均值
cfg.method = 'welch'; % 使用Welch's method估计功率谱
cfg.window = window;
cfg.nperseg = window;
cfg.taper = 'hanning';
psd_data = ft_freqstatistics(data, cfg);

% 对每个通道的功率谱进行归一化处理
psd_data.normed = psd_data.values ./ sum(psd_data.values, 2)';

% 创建频率范围索引
freqs = data.Fs * (0:fmax-fmin+1)./length(data.Fs);
[~, idx] = ismember(freqs(ones(1, window), :), [fmin fmax]);

% 进行聚类
cluster_idx = kmeans(psd_data.normed(idx,:), num_clusters); % 将特定频率范围内的数据分为num_clusters个簇