% http://zafarrafii.com

% https://github.com/zafarrafii

% https://www.linkedin.com/in/zafarrafii/

% Define the properties

properties (Access = private, Constant = true)

% Get the number of samples and channels in the audio signal

[number_samples,number_channels] = size(audio_signal);

% Initialize the STFT

audio_stft = zeros(window_length,number_times,number_channels);

% Loop over the channels

end

audio_spectrogram = abs(audio_stft(1:window_length/2+1,:,:));

beat_spectrum = repet.beatspectrum(mean(audio_spectrogram,3).^2);

period_range = round(repet.period_range*sampling_frequency/step_length);

% Repeating period in time frames given the period range

repeating_period = repet.periods(beat_spectrum,period_range);

cutoff_frequency = ceil(repet.cutoff_frequency*(window_length-1)/sampling_frequency);

% Initialize the background signal

% Loop over the channels

for channel_index = 1:number_channels

repeating_mask = repet.mask(audio_spectrogram(:,:,channel_index),repeating_period);

repeating_mask(2:cutoff_frequency+1,:) = 1;

repeating_mask = cat(1,repeating_mask,repeating_mask(end-1:-1:2,:));

background_signal1 ...

= repet.istft(repeating_mask.*audio_stft(:,:,channel_index),window_function,step_length);