پروژه آموزشی استخراج امواج EEG با متلب + دیتا

تعدادی داده EEG داریم، 22 کاناله و هر کدام حدود 2 دقیقه هستند. برای این 30 داده ، یک (Cell(1*30 نوشته شده است که هر درایه اش یک Cell(1*6) دارد که در درایه اولش هر کدام از داده ها وجود دارد، در درایه دوم فرکانس گامای امواج، در درایه سوم فرکانس بتا، درایه چهارم فرکانس الفا ، درایه پنجم فرکانس تتا و درایه ششم فرکانس دلتای کانال های اصلی وجود دارد. دستورات فراخوانی و plot هم موجود است. با متلب کدنویسی کرده ایم. پروژه متلب

با توجه به محدوده فرکانسی زیر اطلاعات مورد نظر از داده EEG استخراج و برای یک مورد رسم میشود.

 

نام کانال	محدوده فرکانسی
گاما	بالاتر از 32 هرتز
بتا	از 16 تا 31 هرتز
آلفا	از 8 تا 15 هرتز
تتا	از 4 تا 7 هرتز
دلتا	پایینتر از 4 هرتز

 

توضیح کد:

% loading data

load F001.mat

% change variable name

x1 = F001;

 

داده خوانده میشود. شامل 20 کانال که هر کدام به طول تقریبی 16000 نمونه است.

 

% cell initialization

x = {};

% sampling frequency estimation

fs = floor(size(x1,1)/120);

 

سلول مقداردهی اولیه میشود و با توجه به طول زمانی 120 ثانیه ای، فرکانس نمونه برداری تخمین زده میشود.

 

% loop over data

for i=1:size(x1,2)

    % cell initialization

    tmp = {};

 

حلقه جهت ذخیره سازی کانالهای مختلف در سلول شروع میشود.

    % fill missed data and choosing 1st cell for EEG data

    tmp{1} = fillmissing(x1(:,i),’linear’);

 

داده اصلی به صورت خطی از فرآیند پرکردن نقاط خالی عبور داده میشود. این کار برای جبران نقاط NaN در داده ها است. این داده سپس در سلول اول داخلی قرار میگیرد. دقت کنید که سلول اصلی تحت عنوان x خود شامل سلولهایی داخلی است که هر کدام شامل 6 عضو است.