Added tuckeywin function

Author: ppizarror

check_order.m

@@ -0,0 +1,36 @@
+function [n_out, w, trivalwin] = check_order(n_in)
+%CHECK_ORDER Checks the order passed to the window functions.
+% [N,W,TRIVALWIN] = CHECK_ORDER(N_ESTIMATE) will round N_ESTIMATE to the
+% nearest integer if it is not already an integer. In special cases (N is
+% [], 0, or 1), TRIVALWIN will be set to flag that W has been modified.
+
+w = [];
+trivalwin = 0;
+
+if ~(isnumeric(n_in) && isfinite(n_in))
+    error(generatemsgid('InvalidOrder'), 'The order N must be finite.');
+end
+
+% Special case of negative orders:
+if n_in < 0
+    error(generatemsgid('InvalidOrder'), 'Order cannot be less than zero.');
+end
+
+% Check if order is already an integer or empty
+% If not, round to nearest integer.
+if isempty(n_in) || n_in == floor(n_in)
+    n_out = n_in;
+else
+    n_out = round(n_in);
+    warning(generatemsgid('InvalidOrder'), 'Rounding order to nearest integer.');
+end
+
+% Special cases:
+if isempty(n_out) || n_out == 0
+    w = zeros(0, 1); % Empty matrix: 0-by-1
+    trivalwin = 1;
+elseif n_out == 1
+    w = 1;
+    trivalwin = 1;
+end
+

plot_norm_matrix.m

@@ -35,7 +35,8 @@ function plot_norm_matrix(m, mt, mf, t, acc, regname)
 end
 
 %% Plot matrix
-figure();
+f = figure();
+set(f, 'Visible', 'off');
 
 %% Plot matrix
 subplot(2, 1, 1);
@@ -58,5 +59,9 @@ function plot_norm_matrix(m, mt, mf, t, acc, regname)
 xlim([0, t(end)]);
 grid on;
 
+%% Show plot
+set(f, 'Visible', 'on');
+movegui(f, 'center');
+
 end
 

tukeywin.m

@@ -0,0 +1,25 @@
+function w = tukeywin(n, r)
+%TUKEYWIN Tukey window.
+error(nargchk(1, 2, nargin, 'struct')); %#ok<NCHKN>
+
+% Default value for R parameter.
+if nargin < 2 || isempty(r)
+    r = 0.500;
+end
+
+[n, w, trivialwin] = check_order(n);
+if trivialwin, return, end
+
+if r <= 0
+    w = ones(n, 1);
+elseif r >= 1
+    w = hann(n);
+else
+    t = linspace(0, 1, n)';
+    % Defines period of the taper as 1/2 period of a sine wave.
+    per = r / 2;
+    tl = floor(per*(n - 1)) + 1;
+    th = n - tl + 1;
+    % Window is defined in three sections: taper, constant, taper
+    w = [((1 + cos(pi/per*(t(1:tl) - per))) / 2); ones(th-tl-1, 1); ((1 + cos(pi/per*(t(th:end) - 1 + per))) / 2)];
+end