用matlab计算mat文件脑电爆发抑制比的完整代码及各参数意义

以下是MATLAB计算脑电爆发抑制比的完整代码，以及各参数的解释：

% Load data from .mat file
load('data.mat');
eeg_data = data.eeg; % EEG data
fs = data.fs; % sampling frequency
win_len = 4; % window length in seconds
win_ovlap = 0.5; % window overlap in seconds

% Calculate power spectral density (PSD) using Welch's method
nfft = 2^(nextpow2(win_len*fs)); % number of FFT points
window = hamming(win_len*fs); % window function
noverlap = win_ovlap*fs; % number of overlapping samples
[S, f] = pwelch(eeg_data, window, noverlap, nfft, fs);

% Calculate power in the alpha and beta frequency bands
alpha_freq = [8 13]; % alpha frequency range
beta_freq = [14 30]; % beta frequency range
alpha_idx = (f >= alpha_freq(1) & f <= alpha_freq(2));
beta_idx = (f >= beta_freq(1) & f <= beta_freq(2));
alpha_power = trapz(S(alpha_idx,:), 1);
beta_power = trapz(S(beta_idx,:), 1);

% Calculate burst suppression ratio (BSR)
bsr = (beta_power - alpha_power) ./ (beta_power + alpha_power);

% Plot PSD and BSR
figure;
subplot(2,1,1);
semilogx(f, 10*log10(S));
xlabel('Frequency (Hz)');
ylabel('Power/Frequency (dB/Hz)');
title('Power Spectral Density');
legend('EEG');
subplot(2,1,2);
plot(bsr);
ylim([-1 1]);
xlabel('Time (s)');
ylabel('BSR');
title('Burst Suppression Ratio');

各参数的解释如下：

- `data.mat`：需要计算脑电爆发抑制比的EEG数据，存储在MATLAB数据文件中。
- `eeg_data`：EEG数据，是一个一维向量。
- `fs`：采样频率。
- `win_len`：窗口长度，单位为秒。
- `win_ovlap`：窗口重叠长度，单位为秒。
- `nfft`：FFT点数，是2的幂次方。
- `window`：窗口函数，这里使用的是汉明窗。
- `noverlap`：重叠样本数。
- `S`：PSD，是一个矩阵，每一列代表一个时间窗口的PSD。
- `f`：频率向量，与PSD矩阵S的行数一致。
- `alpha_freq`：alpha波频率范围。
- `beta_freq`：beta波频率范围。
- `alpha_idx`：alpha波频率范围对应的索引。
- `beta_idx`：beta波频率范围对应的索引。
- `alpha_power`：alpha波功率，是一个向量，每个元素代表一个时间窗口的alpha波功率。
- `beta_power`：beta波功率，是一个向量，每个元素代表一个时间窗口的beta波功率。
- `bsr`：脑电爆发抑制比，是一个向量，每个元素代表一个时间窗口的BSR值。