Extraction of Correlated Sparse Sources from Signal Mixtures

Pseudocode 2

function[({𝑠₁̃},{𝑠₂̃},…,{𝑠_𝑁̃})] = MHCForward({𝐳},𝑣^th)

Step 1

for 𝑘=1∶𝑁

    for 𝑚=2∶𝑀

        𝐯[𝑚]=𝐳[𝑚]−𝐳[𝑚−1]

        𝐯[𝑚]=𝐯[𝑚]/|𝐯[𝑚]|

    end

Step 2

    for 𝑚=2∶𝑀

        𝑉_max[𝑚]=max𝑣₁[𝑚],𝑣₂[𝑚],…,𝑣_𝑁[𝑚] (19)

    end

    𝑣_max=max{|𝐯[1]|,|𝐯[2]|,…,|𝐯[𝑀]|} (20)

Step3

    for 𝑖=1∶𝑀−1

        𝑗=𝑖+1;

        if 𝑉_max𝑗>𝑣^th𝑣_max

           𝐝⁻𝑗=𝐯[𝑖]−𝐯𝑗

            𝐝⁺𝑗=𝐯[𝑖]+𝐯𝑗

            Δ𝑗=min𝐝⁻𝑗,𝐝⁺𝑗

        end

    end

    𝑗_min = value of 𝑗 where Δ𝑗 is a minimum

    𝐑_𝑘=𝐯𝑗_min

Step 4

    for 𝑚=1∶𝑀

        𝑠_𝑘̃[𝑚]=𝐳[𝑚]⋅𝐑_𝑘 (6)

        𝐬_𝑘̃[𝑚]=𝑠_𝑘̃[𝑚]𝐑_𝑘 (7)

        𝐳[𝑚]←𝐳[𝑚]−𝐬_𝑘̃[𝑚] (10)

    end

end