HL2028_Project/plotHeartRateAnalysis.m at master · GunnarI/HL2028_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
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128

%% Subject 1 RR plots
figure
annotation('textbox', [0 0.9 1 0.1], ...
    'String', 'HR of subject 1 before and after Biofeedback', ...
    'EdgeColor', 'none', ...
    'FontWeight', 'bold', ...
    'HorizontalAlignment', 'center')

sub1 = subplot(4,1,1);
hr_rloc_S1B = r_location_S1B(3:end);
plot(t_S1B(hr_rloc_S1B), hr_S1B(2:end))
title('HR before biofeedback (time domain)')
xlabel('Time [s]')
ylabel('HR [Hz]')
sub2 = subplot(4,1,2);
plot(hr_FFT_f_S1B.*1000, hr_FFT_p_S1B);
title('HR before biofeedback (frequency domain, FFT)')
xlabel('Freq. [Hz]')
ylabel('Amplitude')
xlim([40 400])
sub3 = subplot(4,1,3);
hr_rloc_S1A = r_location_S1A(2:end);
plot(t_S1B(hr_rloc_S1A), hr_S1A)
title('HR after biofeedback (time domain)')
xlabel('Time [s]')
ylabel('HR [Hz]')
sub4 = subplot(4,1,4);
plot(hr_FFT_f_S1A.*1000, hr_FFT_p_S1A);
title('HR after biofeedback (frequency domain, FFT)')
xlabel('Freq. [Hz]')
ylabel('Amplitude')
xlim([40 400])

%% Subject 1: QT plots
figure
annotation('textbox', [0 0.9 1 0.1], ...
    'String', 'HR of subject 1 before and after Biofeedback', ...
    'EdgeColor', 'none', ...
    'FontWeight', 'bold', ...
    'HorizontalAlignment', 'center')

sub1 = subplot(4,1,1);
plot(t_S1B(peakLocs_S1B), QTInterval_S1B)
title('QT Interval before biofeedback (time domain)')
xlabel('Time [s]')
ylabel('Interval [s]')
sub2 = subplot(4,1,2);
plot(qt_FFT_f_S1B.*1000, qt_FFT_p_S1B);
title('QT Interval before biofeedback (frequency domain, FFT)')
xlabel('Freq. [Hz]')
ylabel('Amplitude')
xlim([4 15])
sub3 = subplot(4,1,3);
plot(t_S1A(peakLocs_S1A), QTInterval_S1A)
title('QT Interval after biofeedback (time domain)')
xlabel('Time [s]')
ylabel('Interval [s]')
sub4 = subplot(4,1,4);
plot(qt_FFT_f_S1A.*1000, qt_FFT_p_S1A);
title('QT Interval after biofeedback (frequency domain, FFT)')
xlabel('Freq. [Hz]')
ylabel('Amplitude')
xlim([4 15])

%% Subject 2 RR plots
figure
annotation('textbox', [0 0.9 1 0.1], ...
    'String', 'HR of subject 2 before and after Biofeedback', ...
    'EdgeColor', 'none', ...
    'FontWeight', 'bold', ...
    'HorizontalAlignment', 'center')

sub1 = subplot(4,1,1);
hr_rloc_S2B = r_location_S2B(3:end);
plot(t_S2B(hr_rloc_S2B), hr_S2B(2:end))
title('HR before biofeedback (time domain)')
xlabel('Time [s]')
ylabel('HR [Hz]')
sub2 = subplot(4,1,2);
plot(hr_FFT_f_S2B.*1000, hr_FFT_p_S2B);
title('HR before biofeedback (frequency domain, FFT)')
xlabel('Freq. [Hz]')
ylabel('Amplitude')
xlim([4 15])
sub3 = subplot(4,1,3);
hr_rloc_S2A = r_location_S2A(2:end);
plot(t_S2A(hr_rloc_S2A), hr_S2A)
title('HR after biofeedback (time domain)')
xlabel('Time [s]')
ylabel('HR [Hz]')
sub4 = subplot(4,1,4);
plot(hr_FFT_f_S2A.*1000, hr_FFT_p_S2A);
title('HR after biofeedback (frequency domain, FFT)')
xlabel('Freq. [Hz]')
ylabel('Amplitude')
xlim([4 15])

%% Subject 2: QT plots
figure
annotation('textbox', [0 0.9 1 0.1], ...
    'String', 'HR of subject 2 before and after Biofeedback', ...
    'EdgeColor', 'none', ...
    'FontWeight', 'bold', ...
    'HorizontalAlignment', 'center')

sub1 = subplot(4,1,1);
plot(t_S2B(peakLocs_S2B), QTInterval_S2B)
title('QT Interval before biofeedback (time domain)')
xlabel('Time [s]')
ylabel('Interval [s]')
sub2 = subplot(4,1,2);
plot(qt_FFT_f_S2B.*1000, qt_FFT_p_S2B);
title('QT Interval before biofeedback (frequency domain, FFT)')
xlabel('Freq. [Hz]')
ylabel('Amplitude')
xlim([4 15])
sub3 = subplot(4,1,3);
plot(t_S2A(peakLocs_S2A), QTInterval_S2A)
title('QT Interval after biofeedback (time domain)')
xlabel('Time [s]')
ylabel('Interval [s]')
sub4 = subplot(4,1,4);
plot(qt_FFT_f_S2A.*1000, qt_FFT_p_S2A);
title('QT Interval after biofeedback (frequency domain, FFT)')
xlabel('Freq. [Hz]')
ylabel('Amplitude')
xlim([4 15])