Updated header to show it with help.

applications/apertif/matlab/FIR_LPF_ApertifCF.m

@@ -18,8 +18,12 @@
 % along with this program.  If not, see <http://www.gnu.org/licenses/>.
 %
 %-----------------------------------------------------------------------------
-% Author: R. de Wild, 2015
+% Author: R. de Wild, 2015 (Original)
 %         E. Kooistra, 2016
+%
+% Purpose : Calculate FIR filter coefficients for Apertif channel
+%           filterbank using fircls1.
+% Description :
 
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 % workspace initiation

applications/apertif/matlab/bsn_source.m

@@ -19,11 +19,12 @@
 %
 %-----------------------------------------------------------------------------
 % Author: E. Kooistra, 2016
-
+%
 % Purpose : BSN source per block of data
 % Description :
 
 function [state] = bsn_source(ctrl)
+
 % Increment BSN for next call
 state = ctrl;
 state.bsn = ctrl.bsn + 1;

applications/apertif/matlab/dt.m

@@ -19,7 +19,7 @@
 %
 %-----------------------------------------------------------------------------
 % Author: E. Kooistra, 2016
-
+%
 % Purpose : Delay tracking per block of data
 % Description :
 %
@@ -41,6 +41,7 @@
 % * Index = Time + 1, because for Matlab index starts at 1 instead of at 0.
 
 function [state, out_data] = dt(ctrl, in_data)
+
 % Data processing
 N = ctrl.block_size;
 

applications/apertif/matlab/fringe.m

@@ -18,7 +18,8 @@
 % along with this program.  If not, see <http://www.gnu.org/licenses/>.
 %
 %-----------------------------------------------------------------------------
-
+% Author: E. Kooistra, 2016
+%
 % Purpose:
 %   Plot the fringe due to the varying difference in geometrical delay.
 % Description:

applications/apertif/matlab/fringe_stopping.m

@@ -19,8 +19,8 @@
 %
 %-----------------------------------------------------------------------------
 % Author: W. Lubberhuizen, 2015 (top.m for LOFAR Station, uses classes)
-%         Eric Kooistra, 2016
-
+%         Eric Kooistra, 2016 (for Apertif)
+%
 % Purpose : Model data path for delay tracking with polyphase filterbank
 % Description :
 % * General

applications/apertif/matlab/pfft.m

@@ -19,11 +19,12 @@
 %
 %-----------------------------------------------------------------------------
 % Author: E. Kooistra, 2016
-
+%
 % Purpose : FFT per block of data
 % Description :
 
 function [out_data] = pfft(ctrl, in_data)
+
 % Data processing
 out_data = fft(in_data, ctrl.fft_size);
 

applications/apertif/matlab/pfir.m

@@ -19,11 +19,12 @@
 %
 %--------------------------------------------------------------------------
 % Author: E. Kooistra, 2016
-
+%
 % Purpose : Polyphase FIR filter per block of data
 % Description :
 
 function [state, out_data] = pfir(ctrl, in_data)
+
 % Data processing
 for t = 1:ctrl.downsample_factor
     [out_data(t), ctrl.Zdelays(t,:)] = filter(ctrl.coeff(t,:), 1, in_data(t), ctrl.Zdelays(t,:));   % range ctrl.nof_taps coefficients and range ctrl.nof_taps-1 Zdelays

applications/apertif/matlab/pfir_coeff.m

@@ -20,8 +20,7 @@
 %-----------------------------------------------------------------------------
 % Author: W. Lubberhuizen, 2005
 %         E. Kooistra, 2016
-
-function coeff = pfir_coeff(N, L, BWchan, r_pass, r_stop, nof_bits, config)
+%
 % pfir_coeff - compute polyphase filterbank coefficients
 %   Computes the coefficients for a polyphase filterbank with N fft
 %   channels, and subfilter length L and channel bandwidth BWchan.
@@ -38,6 +37,8 @@ function coeff = pfir_coeff(N, L, BWchan, r_pass, r_stop, nof_bits, config)
 %   The returned value is a matrix of dimension (NxL)
 %   The run_pfir_coeff.m shows how this pfir_coeff.m can be used.
 
+function coeff = pfir_coeff(N, L, BWchan, r_pass, r_stop, nof_bits, config)
+
 % requested filter length
 M2=N*L;
 if strcmp(config.design, 'fir1')

applications/apertif/matlab/phase_tracking.m

@@ -18,7 +18,8 @@
 % along with this program.  If not, see <http://www.gnu.org/licenses/>.
 %
 %-----------------------------------------------------------------------------
-
+% Author: E. Kooistra, 2016
+%
 % Purpose:
 %   Plot the phase tracking (PT) phase and phase erroras function of the
 %   varying geometrical delay.

applications/apertif/matlab/quantize.m

@@ -19,8 +19,7 @@
 %
 %-----------------------------------------------------------------------------
 % Author: E. Kooistra, 2016
-
-function out_data = quantize(in_data, in_max, nof_bits, direction, overflow)
+%
 % quantize - quantize signed input data
 %     The output has nof_bits.
 %     Internally in_max maps to 2^(nof_bits-1)-1.
@@ -33,6 +32,8 @@ function out_data = quantize(in_data, in_max, nof_bits, direction, overflow)
 %     . 'wrap' : Wrap within nof_bits signed integer range
 %     The output data range is scaled back to the original in_max.
 
+function out_data = quantize(in_data, in_max, nof_bits, direction, overflow)
+
 if ~exist('direction', 'var'); direction = 'half_away'; end;
 if ~exist('overflow', 'var'); overflow = 'clip'; end;
     

applications/apertif/matlab/run_pfir_coeff.m

@@ -19,7 +19,7 @@
 %
 %-----------------------------------------------------------------------------
 % Author: E. Kooistra, 2016
-
+%
 % Purpose : Run pfir_coeff
 % Description :
 %   Plot FIR filter coefficients and transfer function.

applications/apertif/matlab/spec.m

@@ -19,7 +19,7 @@
 %
 %-----------------------------------------------------------------------------
 % Author: E. Kooistra, 2016
-
+%
 % Purpose : Image spectrum per block of data
 % Description :
 

applications/apertif/matlab/try_quantize.m

@@ -19,7 +19,7 @@
 %
 %-----------------------------------------------------------------------------
 % Author: E. Kooistra, 2016
-
+%
 % Purpose : Try quantize.m
 % Description :
 

applications/apertif/matlab/try_uencode.m

@@ -19,7 +19,7 @@
 %
 %-----------------------------------------------------------------------------
 % Author: E. Kooistra, 2016
-
+%
 % Purpose : Try uencode to determine optimum quantisation
 % Description :
 % Remark : It seems more appropriate to use round() instead of uencode(),

applications/apertif/matlab/wg.m

@@ -19,13 +19,14 @@
 %
 %-----------------------------------------------------------------------------
 % Author: E. Kooistra, 2016
-
+%
 % Purpose : Waveform generator per block of data
 % Description :
 %   The WG data is normalized to -1:1. Data outside this range is clipped,
 %   so use ampl >> 1 to create square wave.
 
 function [state, data] = wg(ctrl)
+
 % Data processing
 t = 0:ctrl.block_size-1;