function signal_win = blackmanWin(signal, signal_duration, sampling_freq, pad)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%function signal_win = blackmanWin(signal)
%
% Inputs:
%	- signal: signal of interest
%
% Output:
%	- signal_win: signal of interest on which a blackman window was applied
%
% Author: Guillaume Gibert, guillaume.gibert@ecam.fr
% Date: 15/03/2024
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% create the  temporal array
t=-signal_duration/2:1/sampling_freq:(signal_duration/2)-1/sampling_freq;
% window duration is half of signal duration
windowDuration = signal_duration/2;


blackmanWin = zeros(1, length(t));

for l_sample=1:length(signal)
	blackmanWin(l_sample) = (0.42 - 0.5 * cos(2*pi*(l_sample)/length(signal)) + 0/08*cos(4*pi*(l_sample)/length(signal)));
end

% plot Blackman window
figure;
title("Blackman Window"); hold on;
plot(t, blackmanWin);

% apply the Blackman window
for l_sample=1:length(signal)
	signal_win(l_sample) = signal(l_sample) * blackmanWin(l_sample);
end

figure;
title("Original and Windowed signals"); hold on;
plot(signal); hold on;
plot(signal_win);

frequencySpectrum(signal_win, sampling_freq/2,pad);