Main_2_compute_mp_atoms.m

% Copyright (c) 2021 Mohammad Fathi Al-Sa'd
% 
% Permission is hereby granted, free of charge, to any person obtaining a
% copy of this software and associated documentation files (the "Software"),
% to deal in the Software without restriction, including without limitation
% the rights to use, copy, modify, merge, publish, distribute, sublicense,
% and/or sell copies of the Software, and to permit persons to whom the
% Software is furnished to do so, subject to the following conditions:
% 
% The above copyright notice and this permission notice shall be included
% in all copies or substantial portions of the Software.
% 
% THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS
% OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
% FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL
% THE AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
% LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING
% FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER
% DEALINGS IN THE SOFTWARE.
%
% Email: mohammad.al-sad@tuni.fi, alsad.mohamed@gmail.com
%
% The following reference should be cited whenever this script is used:
% M. Al-Sa'd, B. Boashash, M. Gabbouj, "Design of an Optimal Piece-Wise
% Spline Wigner-Ville Distribution for TFD Performance Evaluation and
% Comparison", IEEE Transactions on Signal Processing, 2021.
%
% Last Modification: 21-May-2021
%
% Description:
% This main script decomposes the database signals using the matching
% pursuit approach with chirplet atoms and saves the results in
% MP_atoms.mat. It generates a chirplet dictionary with approximately 6
% million atoms and decomposes each signal in the database to reach a 1%
% minimum relative L2 error.

%% Initialization
clear; close all; clc;
addpath(genpath('Functions'));

%% Loading database
load(['Data' filesep 'Database' filesep 'pw_wvd_database.mat'],'z','fs','N','K');

%% Generate chirplet dictionary
dict = chirplet_dictionary(N,fs);

%% Main
ep  = 1;
Atoms = cell(1,K);
for k = 1:K
    disp(100*k/K)
    Atoms{k} = mp_atoms(z{k},dict,ep);
end
save(['Data' filesep 'MP' filesep 'MP_atoms.mat'],'Atoms');