waveleteeg.m

Decomposition EEG to its 5 Sub-bands using Wavelet

%% Decomposition
close all
wavelet='coif5';
b=wavelet_decompose(y,5,wavelet);
[C,L] = wavedec(y,5,wavelet);
% figure;
% subplot(5,1,1);plot(b(:,1));title('EEG 0-60 Hz');
% subplot(5,1,2);plot(b(:,2));title('Betta');
% subplot(5,1,3);plot(b(:,3));title('Alpha');
% subplot(5,1,4);plot(b(:,4));title('Theta');
% subplot(5,1,5);plot(b(:,5));title('Delta');

b1=wavelet_decompose(b(:,2),2,wavelet);
figure;
subplot(2,1,1);plot(b1(:,1));title('Betta low');
Betta_low_power=sum(abs(b1(:,1)).^2/length(b1(:,1)));
subplot(2,1,2);plot(b1(:,2));title('Betta up');
Betta_up_power=sum(abs(b1(:,2)).^2/length(b1(:,2)));

b2=wavelet_decompose(b(:,3),2,wavelet);
figure;
subplot(2,1,1);plot(b2(:,1));title('Alpha low');
Alpha_low_power=sum(abs(b2(:,1)).^2/length(b2(:,1)));
subplot(2,1,2);plot(b2(:,2));title('Alpha up');
Alpha_up_power=sum(abs(b2(:,2)).^2/length(b2(:,2)));

b3=wavelet_decompose(b(:,4),2,wavelet);
figure;
subplot(2,1,1);plot(b3(:,1));title('Theta low');
Theta_low_power=sum(abs(b3(:,1)).^2/length(b3(:,1)));
subplot(2,1,2);plot(b3(:,2));title('Theta up');
Theta_up_power=sum(abs(b3(:,2)).^2/length(b3(:,2)));

b4=wavelet_decompose(b(:,5),2,wavelet);
figure;
subplot(2,1,1);plot(b4(:,1));title('Delta low');
Delta_low_power=sum(abs(b4(:,1)).^2/length(b4(:,1)));
subplot(2,1,2);plot(b4(:,2));title('Delta up');
Delta_up_power=sum(abs(b4(:,2)).^2/length(b4(:,2)));

%% Total Power in each band
X=[Betta_low_power,Betta_up_power,Alpha_low_power,Alpha_up_power,Theta_low_power,...
    Theta_up_power,Delta_low_power,Delta_up_power];

[Ea,Ed] = wenergy(C,L);
q19=[Ea,Ed];
save q19