diff --git a/applications/apertif/matlab/delay_tracking_pfb.m b/applications/apertif/matlab/delay_tracking_pfb.m
index 87924e047b67a09e3289d841981dcfa80be6d025..18da1b8abc5ac041a0a1b635245bb82431571a6e 100644
--- a/applications/apertif/matlab/delay_tracking_pfb.m
+++ b/applications/apertif/matlab/delay_tracking_pfb.m
@@ -311,6 +311,7 @@ sub_phase(x) = 0;   % force phase of too small signals to 0
 fig=fig+1;
 figure('position', [xfig+fig*dfig yfig-fig*dfig xfigw yfigw]);
 figure(fig);
+
 subplot(2,1,1);
 data = sub_ampl';
 data = data(:);
@@ -319,12 +320,13 @@ title(sprintf('Subband data - amplitude  (o,x = subband %d,%d for WG subband = %
 xlabel(sprintf('Subbands 0:%d at time 0:%d [Tsub]', tb.nof_subbands-1, tb.nof_tsub-1));
 ylabel('Voltage');
 grid on;
+
 subplot(2,1,2);
 data = sub_phase';
 data = data(:);
 plot(tsub_all, data, 'k', tsub_all(sub_I), data(sub_I), 'ko', tsub_all(sub_Iplus1), data(sub_Iplus1), 'kx');
 ylim([-180 180])
-title(sprintf('Subband data - amplitude  (o,x = subband %d,%d for WG subband = %6.3f)', tb.subband_i, tb.subband_iplus1, tb.subband_wg));
+title(sprintf('Subband data - phase  (o,x = subband %d,%d for WG subband = %6.3f)', tb.subband_i, tb.subband_iplus1, tb.subband_wg));
 xlabel(sprintf('Subbands 0:%d at time 0:%d [Tsub]', tb.nof_subbands-1, tb.nof_tsub-1));
 ylabel('Phase [degrees]');
 grid on;
diff --git a/applications/apertif/matlab/one_pfb.m b/applications/apertif/matlab/one_pfb.m
index b71b41ed3a1cb49512463a680dd8d382c5202f8c..d57c29de29a2e6b72dc24a8a5525abb9a794f05e 100644
--- a/applications/apertif/matlab/one_pfb.m
+++ b/applications/apertif/matlab/one_pfb.m
@@ -429,6 +429,7 @@ sub_phase(x) = 0;   % force phase of too small signals to 0
 fig=fig+1;
 figure('position', [xfig+fig*dfig yfig-fig*dfig xfigw yfigw]);
 figure(fig);
+
 subplot(2,1,1);
 data = sub_ampl.';
 data = data(:);
@@ -437,12 +438,13 @@ title(sprintf('Subband data - amplitude  (o,x = subband %d,%d for WG subband = %
 xlabel(sprintf('Subbands 0:%d at time 0:%d [Tsub]', tb.nof_subbands-1, tb.nof_tsub-1));
 ylabel('Voltage');
 grid on;
+
 subplot(2,1,2);
 data = sub_phase.';
 data = data(:);
 plot(tsub_all, data, 'k', tsub_all(sub_I), data(sub_I), 'ko', tsub_all(sub_Iplus1), data(sub_Iplus1), 'kx');
 ylim([-180 180])
-title(sprintf('Subband data - amplitude  (o,x = subband %d,%d for WG subband = %6.3f)', tb.subband_i, tb.subband_iplus1, tb.subband_wg));
+title(sprintf('Subband data - phase  (o,x = subband %d,%d for WG subband = %6.3f)', tb.subband_i, tb.subband_iplus1, tb.subband_wg));
 xlabel(sprintf('Subbands 0:%d at time 0:%d [Tsub]', tb.nof_subbands-1, tb.nof_tsub-1));
 ylabel('Phase [degrees]');
 grid on;
diff --git a/applications/apertif/matlab/psk_pfb.m b/applications/apertif/matlab/psk_pfb.m
index a3d89d4c88f760c16e096d5a3f7fe589b5124dcb..64c4d19f5c35b2893daedc69f3299af1e080bb8e 100644
--- a/applications/apertif/matlab/psk_pfb.m
+++ b/applications/apertif/matlab/psk_pfb.m
@@ -320,6 +320,7 @@ sub_phase(x) = 0;   % force phase of too small signals to 0
 fig=fig+1;
 figure('position', [xfig+fig*dfig yfig-fig*dfig xfigw yfigw]);
 figure(fig);
+
 subplot(2,1,1);
 data = sub_ampl';
 data = data(:);
@@ -328,12 +329,13 @@ title(sprintf('Subband data - amplitude  (o,x = subband %d,%d for WG subband = %
 xlabel(sprintf('Subbands 0:%d at time 0:%d [Tsub]', tb.nof_subbands-1, tb.nof_tsub-1));
 ylabel('Voltage');
 grid on;
+
 subplot(2,1,2);
 data = sub_phase';
 data = data(:);
 plot(tsub_all, data, 'k', tsub_all(sub_I), data(sub_I), 'ko', tsub_all(sub_Iplus1), data(sub_Iplus1), 'kx');
 ylim([-180 180])
-title(sprintf('Subband data - amplitude  (o,x = subband %d,%d for WG subband = %6.3f)', tb.subband_i, tb.subband_iplus1, tb.subband_wg));
+title(sprintf('Subband data - phase  (o,x = subband %d,%d for WG subband = %6.3f)', tb.subband_i, tb.subband_iplus1, tb.subband_wg));
 xlabel(sprintf('Subbands 0:%d at time 0:%d [Tsub]', tb.nof_subbands-1, tb.nof_tsub-1));
 ylabel('Phase [degrees]');
 grid on;
diff --git a/applications/apertif/matlab/run_pfb.m b/applications/apertif/matlab/run_pfb.m
index 5be728f9cf00f6737184953b8720db3f192423fb..67357297540eec2d8a36faa0e8d2cb2d5f87bd83 100644
--- a/applications/apertif/matlab/run_pfb.m
+++ b/applications/apertif/matlab/run_pfb.m
@@ -419,6 +419,7 @@ sub_phase(x) = 0;   % force phase of too small signals to 0
 fig=fig+1;
 figure('position', [xfig+fig*dfig yfig-fig*dfig xfigw yfigw]);
 figure(fig);
+
 subplot(2,1,1);
 data = sub_ampl.';
 data = data(:);
@@ -427,12 +428,13 @@ title(sprintf('Subband data - amplitude  (o,x = subband %d,%d for WG subband = %
 xlabel(sprintf('Subbands 0:%d at time 0:%d [Tsub]', tb.nof_subbands-1, tb.nof_tsub-1));
 ylabel('Voltage');
 grid on;
+
 subplot(2,1,2);
 data = sub_phase.';
 data = data(:);
 plot(tsub_all, data, 'k', tsub_all(sub_I), data(sub_I), 'ko', tsub_all(sub_Iplus1), data(sub_Iplus1), 'kx');
 ylim([-180 180])
-title(sprintf('Subband data - amplitude  (o,x = subband %d,%d for WG subband = %6.3f)', tb.subband_i, tb.subband_iplus1, tb.subband_wg));
+title(sprintf('Subband data - phase  (o,x = subband %d,%d for WG subband = %6.3f)', tb.subband_i, tb.subband_iplus1, tb.subband_wg));
 xlabel(sprintf('Subbands 0:%d at time 0:%d [Tsub]', tb.nof_subbands-1, tb.nof_tsub-1));
 ylabel('Phase [degrees]');
 grid on;
diff --git a/applications/apertif/matlab/run_pfft.m b/applications/apertif/matlab/run_pfft.m
index 4f4bfbcfe5251194d5df1d7cabd44bcad733ae97..fbc1bfe6de6d9ac385b5910e2bc799739d21b090 100644
--- a/applications/apertif/matlab/run_pfft.m
+++ b/applications/apertif/matlab/run_pfft.m
@@ -289,7 +289,7 @@ data = sub_phase.';
 data = data(:);
 plot(tsub_all, data, 'k', tsub_all(sub_I), data(sub_I), 'ko', tsub_all(sub_Iplus1), data(sub_Iplus1), 'kx');
 ylim([-180 180])
-title(sprintf('Subband data - amplitude  (o,x = subband %d,%d for WG subband = %6.3f)', tb.subband_i, tb.subband_iplus1, tb.subband_wg));
+title(sprintf('Subband data - phase  (o,x = subband %d,%d for WG subband = %6.3f)', tb.subband_i, tb.subband_iplus1, tb.subband_wg));
 xlabel(sprintf('Subbands 0:%d at time 0:%d [Tsub]', tb.nof_subbands-1, tb.nof_tsub-1));
 ylabel('Phase [degrees]');
 grid on;