EEG_Analysis_Scripts/examples_for_asymmetry.m at master · mdtdev/EEG_Analysis_Scripts · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
ls
EDB.Fs = 128;
EDB.data = EEG.data';
EDB
help averagePowerSpectrum
[avgS, eegFreqs, pileOfSpectra] = averagePowerSpectrum(EDB.data, EDB.Fs, 4, 0.250);
plotColumns(eegFreqs, pileOfSpectra(:,1,:))
hold on
[avgS, eegFreqs, pileOfSpectra] = averagePowerSpectrum(EDB.data, EDB.Fs, 8, 0.250);
hold off
close all
plotColumns(eegFreqs, pileOfSpectra(:,1,:))
hold on
plot(eegFreqs, avgS, 'r', 'lineWidth', 4);
close all
plotColumns(eegFreqs, pileOfSpectra(:,1,:))
plot(eegFreqs, avgS(:, 1), 'r', 'lineWidth', 2);
plot(eegFreqs, avgS(:, 14), 'g', 'lineWidth', 2);
close all
plotColumns(eegFreqs, pileOfSpectra(:,3,:))
hold on
plot(eegFreqs, avgS(:, 3), 'r', 'lineWidth', 2);
plot(eegFreqs, avgS(:, 12), 'g', 'lineWidth', 2);
figure;
plot(eegFreqs, avgS(:, 4), 'r', 'lineWidth', 2);
hold on
plot(eegFreqs, avgS(:, 11), 'g', 'lineWidth', 2);
close all

figure;
plot(eegFreqs, avgS(:, 3), 'r', 'lineWidth', 2);
hold on
plot(eegFreqs, avgS(:, 12), 'g', 'lineWidth', 2);
title('F34');
xlabel('Frequency');
ylabel('Power (Square Units)');


figure;
plot(eegFreqs, avgS(:, 4), 'r', 'lineWidth', 2);
hold on
plot(eegFreqs, avgS(:, 11), 'g', 'lineWidth', 2);
title('FC56');
xlabel('Frequency');
ylabel('Power (Square Units)');