使用MATLAB进行PSD特征提取EEG可视化的示例代码:
时间: 2024-03-06 13:51:36 浏览: 283
基于MATLAB和FFT(快速傅里叶变换)的脑电EEG信号的功率谱密度(PSD)提取方法试验
以下是使用MATLAB进行PSD特征提取EEG可视化的示例代码:
```matlab
% 加载EEG数据(使用EEGLAB工具箱)
EEG = pop_loadset('mydata.set');
% PSD特征提取
[PSD,f] = pwelch(EEG.data',[],[],EEG.srate);
% 可视化PSD
figure;
plot(f,10*log10(PSD));
xlabel('Frequency (Hz)');
ylabel('Power/Frequency (dB/Hz)');
title('Power Spectral Density');
% 提取alpha频带PSD特征
alphaIdx = (f >= 8) & (f <= 13);
alphaPSD = mean(PSD(alphaIdx,:));
% 提取beta频带PSD特征
betaIdx = (f >= 13) & (f <= 30);
betaPSD = mean(PSD(betaIdx,:));
% 提取theta频带PSD特征
thetaIdx = (f >= 4) & (f <= 8);
thetaPSD = mean(PSD(thetaIdx,:));
```
此代码加载EEG数据,使用`pwelch`函数提取PSD特征,并可视化PSD。然后,提取alpha、beta和theta频带的PSD特征。
阅读全文