Extraction of Correlated Sparse Sources from Signal Mixtures

Pseudocode 1

Step 1

𝑠₁̃,{𝑠₂̃},…,𝑠_𝑁̃ = MHCForward({𝐳}, 𝑣^th) Section 2.1

Step 2

for 𝑚=1∶𝑀

    𝑠̃_𝑁⁽⁰⁾[𝑚]=𝑠_𝑁̃[𝑚] (A.6)

end

𝑠̃_𝑁−1⁽¹⁾,Δ𝐻,NoSub=MHCReverse𝑠̃_𝑁−1,𝑠_𝑁̃,𝜆^th,Corr_min,𝑃 (A.8)

for 𝑖=2∶𝑁−1

    for 𝑗=1∶𝑁

        flagsub(𝑗)=0

    end

    𝑘=0

    for 𝑚=1∶𝑀

        𝑠̃_𝑁−𝑖⁽⁰⁾[𝑚]=𝑠̃_𝑁−𝑖[𝑚]

   end

   while 𝑘≤𝑖

      𝑘=𝑘+1

      for 𝑗=0∶𝑖−1

          Δ𝐻(𝑗)=1𝑒10

      end

      for𝑗=0∶𝑖−1

          if flagsub(𝑗) = 0

             𝑡^(𝑗)_𝑁−𝑖,Δ𝐻𝑗,NoSub𝑗=

             MHCReverse𝑠̃_𝑁−𝑖^(𝑘−1),𝑠̃_𝑁−𝑗^(𝑗),𝜆^th,Corr_min,𝑃

          end

      end

Step 3

      if min(NoSub(𝑗)) > 0

        𝑗_min = value of 𝑗∶Δ𝐻(𝑗) is a minimum

        for 𝑚=1∶𝑀

            𝑠̃_𝑁−𝑖^(𝑘)[𝑚]=𝑡^(𝑗min)_𝑁−𝑖[𝑚]

        end

        flagsub𝑗_min=1

      else

        𝑘=𝑖

        for 𝑚=1∶𝑀

            𝑠̃_𝑁−𝑖^(𝑖)[𝑚]=𝑠̃_𝑁−𝑖[𝑚]

        end

      end

   end

end