PNES-Project/freqCheck.m at master · roeysc/PNES-Project · 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
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
clear; clc;

signal = load('C:\Users\roeysc\Desktop\PNES\queries\sinus.mat');

t = signal.data.time{1};

% %% Single Trial
% % y = sin(4.5*2*pi.*t);
% y = sin(40*2*pi.*t) + sin(41*2*pi.*t) + 5*sin(10*2*pi.*t);
%
% for i=1:19
%     signal.data.trial{1}(i,:) = i.*y;
% end
%
% data = signal.data;
%
% save(['C:\Users\roeysc\Desktop\PNES\queries\sinus3.mat'], 'data');
%
% % freq1 = freqAnalysisWarp('C:\Users\roeysc\Desktop\PNES\queries\sinus3.mat', [0:40], 1);
% %
% % frequencies = freq1.freq;
% % pow = freq1.powspctrm(1,19,:);
% % pow = pow(1,:);
% % bar(frequencies, pow


%%

%  Multiple trials with small variance
for j=1:100

y = sin((10+2*randn-1)*2*pi.*t);
% y = sin((40+4*randn)*2*pi.*t) + sin((30+5*randn)*2*pi.*t) + 5*sin((10+6*randn)*2*pi.*t);

            for i=1:19
                signal.data.trial{j}(i,:) = 10.*y;
            end
            signal.data.time{j} = signal.data.time{1};

end

data = signal.data;
save(['C:\Users\roeysc\Desktop\PNES\queries\sinus1.mat'], 'data');

data = [];
signal = load('C:\Users\roeysc\Desktop\PNES\queries\sinus.mat');

for j=1:100

y = sin((15+2*randn-1)*2*pi.*t);
% y = sin((40+4*randn)*2*pi.*t) + sin((30+5*randn)*2*pi.*t) + 5*sin((10+6*randn)*2*pi.*t);

            for i=1:19
                signal.data.trial{j}(i,:) = 10.*y;
            end
            signal.data.time{j} = signal.data.time{1};

end

data = signal.data;
save(['C:\Users\roeysc\Desktop\PNES\queries\sinus2.mat'], 'data');

%% Plot the freq1, freq2 plots
% y = mean( freq1.powspctrm(:,1,:),1 );
% y = y(:);
% x = freq2.freq;
% bar(x,y);
% y = mean( freq2.powspctrm(:,1,:),1 );
% y = y(:);
% x = freq2.freq;
% hold on
% bar(x,y,'r');
% title('PowerSpectrum');
% legend('y1', 'y2');
%