Extraction of Correlated Sparse Sources from Signal Mixtures

Pseudocode 3

function𝑠̃_A⁽¹⁾,Δ𝐻,NoSub= MHCReverse𝑠_Ã,𝑠_B̃,𝜆^th,Corr_min,𝑃

Step 1

for 𝑚=1∶𝑀

    𝐳[𝑚]=𝑠_Ã[𝑚],𝑠_B̃[𝑚]

end

for 𝑚=2∶𝑀

    𝐯[𝑚]=𝐳[𝑚]−𝐳[𝑚−1] Compute velocity vector, (11)

end

Step 2

for 𝑚=2∶𝑀

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

end

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

Step 3

for 𝑖=𝑃∶𝑀

    for 𝑗=1∶𝑃

        𝐳^local𝑗=𝐳[𝑖−𝑃+𝑗]

    end

    𝐶[𝑖] = corrcoef(𝑧^local₁,𝑧^local₂)

end

Step 4

for 𝑖=𝑃∶𝑀−1

    𝑗=𝑖+1;

    if 𝑣₁[𝑖]>𝜆^th𝑣_max & 𝑣₂[𝑖]>𝜆^th𝑣_max & 𝑣₁𝑗>𝜆^th𝑣_max&𝑣₂𝑗>𝜆^th𝑣_max

    &  |𝐶[𝑖]|>Corr_min

        Flag = 1

         𝐝⁻𝑗=𝐯[𝑖]−𝐯𝑗

         𝐝⁺𝑗=𝐯[𝑖]+𝐯𝑗

         Δ𝑗=min𝐝⁻𝑗,𝐝⁺𝑗

    else

        Flag = 0

    end

end

if Flag = 1

    Δ𝐻=minΔ𝑗

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

    𝐑_B=𝐯𝑗_min

end

Step 5

if Flag = 1

  NoSub = 1

   for 𝑚=1∶𝑀

      𝑠̃_A⁽¹⁾[𝑚]=𝑠_Ã[𝑚]−(𝐳[𝑚]⋅𝐑_B)𝑅_B⁽¹⁾  (35)

   end

else

  NoSub = 0

  for 𝑚=1∶𝑀

      𝑠̃_A⁽¹⁾[𝑚]=𝑠_Ã[𝑚]

  end

end