Data_Analysis/plotPhotodiode.m at master · ryantmaloney/Data_Analysis · 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
qfunction plotPhotodiode(unfilteredPhotodiode)


figure(1)
smtlb = sgolayfilt(double(unfilteredPhotodiode),1,91);
sr=30000; %Sample rate
t=1:length(smtlb);
t=t./sr;

startTime=0;
endTime=160;

% Start Threshold
% startTime=45.8295;
% endTime=45.83;

% % End Time
% startTime=79.823;
% endTime=79.831;

disp(1/(endTime-startTime))

subplot(4,1,1)
plot(t, (unfilteredPhotodiode))
axis([startTime, endTime, 1.5*10^4, 1.7*10^4])
subplot(4,1,2)

plot(t, smtlb);
axis([startTime, endTime, 1.56*10^4, 1.63*10^4])

subplot(4,1,3)
extendedChirp=zeros(size(t));
chirp=createChirpStimulus();


stimulusLength=length(chirp)-1;
% offset=floor(40.8298*sr);

offset=calculateChirpStartTime(unfilteredPhotodiode)*sr;

size(chirp)
size(extendedChirp(offset:(offset+stimulusLength)))
extendedChirp(offset:offset+stimulusLength)=extendedChirp(offset:(offset+stimulusLength))+chirp';
plot(t, extendedChirp);
axis([startTime, endTime, 0, 1]);

subplot(4,1,4)
disp(length(t))
disp(length(chirp))
plot(t(1:length(chirp)), chirp)


end