计算脑电信号相位锁定值 Phase Locking Value

最新推荐文章于 2024-08-29 20:03:21 发布

油泼西红柿

最新推荐文章于 2024-08-29 20:03:21 发布

阅读量7.8k

点赞数 8

分类专栏：脑机接口

本文链接：https://blog.csdn.net/qq_37813206/article/details/117923269

版权

本文介绍了计算脑电信号相位锁定值（Phase Locking Value, PLV）的方法，涉及希尔伯特变换、相位差计算及指数平均。PLV值在0-1之间，0表示相位无同步。提供了Matlab函数的链接及使用示例，展示如何处理16通道x840数据点x80次试验的脑电信号，滤波范围为8-12Hz，并绘制了特定通道的PLV结果。" 132312011,7337247,Kafka在大数据处理中的应用与实践,"['大数据处理', '实时计算', '消息队列', 'Kafka', '复杂事件处理']

摘要由CSDN通过智能技术生成

PLV：相位同步信息在时域的计算值

分别计算两个通道数据的希尔伯特变换，然后计算数据的瞬时相位差值，然后根据所有点的相位差计算自然常数为底的指数，将指数的绝对值求和平均，得到PLV的值

值的区间是0-1，0说明相位完全没同步

Matlab函数：文末有附
https://www.mathworks.com/matlabcentral/fileexchange/31600-phase-locking-value

用法：
脑电数据格式是16个通道x840个数据点x80个trial
采样率1200Hz
滤波区间8-12Hz
绘制3通道和4通道的PLV值

  eegData=samples(:,:,find(labels==0));
  srate = 1200; %Hz
  filtSpec.order = 50;
  filtSpec.range = [8 12]; %Hz
  [plv] = pn_eegPLV(eegData, srate, filtSpec);
  figure;
  plot(squeeze(plv(:,3,4)));
  xlabel('Time (s)'); 
  ylabel('Plase Locking Value');

结果：
在这里插入图片描述
函数：

function [plv] = pn_eegPLV(eegData, srate, filtSpec, dataSelectArr)
% Computes the Phase Locking Value (PLV) for an EEG dataset.
%
% Input parameters:
%   eegData is a 3D matrix numChannels x numTimePoints x numTrials
%   srate is the sampling rate of the EEG data
%   filtSpec is the filter specification to filter the EEG signal in the
%     desired frequency band of interest. It is a structure with two
%     fields, order and range. 
%       Range specifies the limits of the frequency
%     band, for example, put filtSpec.range = [35 45] for gamma band.
%       Specify the order of the FIR filter in filtSpec.order. A useful
%     rule of thumb can be to include about 4 to 5 cycles of the desired
%     signal. For example, filtSpec.order = 50 for eeg data sampled at
%     500 Hz corresponds to 100 ms and contains ~4 cycles of gamma band
%     (40 Hz).
%   dataSelectArr (OPTIONAL) is a logical 2D matrix of size