Input:
Measurement matrix
๐ด
, measurements
๐ฆ
, thresholding parameter
๐ฝ
Output:
The reconstructed signal
๐ฅ
(
1
)
Initialization:
(
2
)
๐ก
=
1
//iteration number
(
3
)
๐ฅ
(
0
)
=
๐
//initial signal
(
4
)
๐
(
0
)
=
๐ฆ
//initial residual
(
5
)
๐ผ
(
0
)
=
โ
//initial support set
(
6
)
while
halting criterion false
do
(
7
)
๐ฅ
(
๐ก
)
=
๐ฅ
(
๐ก
โ
1
)
+
๐ด
โ
๐
(
๐ก
โ
1
)
(
8
)
๐ผ
(
๐ก
)
=
{
๐
โถ
|
๐ฅ
๐
(
๐ก
)
|
>
๐ฝ
๐ก
m
a
x
|
๐ฅ
|
}
(
9
)
๐ฅ
๐ผ
(
๐ก
)
(
๐ก
)
=
๐ด
โ
๐ผ
(
๐ก
)
๐ฆ
(
1
0
)
๐ฅ
(
๐ก
)
[
1
,
๐
]
โงต
๐ผ
(
๐ก
)
=
0
(
1
1
)
๐
(
๐ก
)
=
๐ฆ
โ
๐ด
๐ฅ
(
๐ก
)
(
1
2
)
๐ก
=
๐ก
+
1
(
1
3
)
end while
(
1
4
)
return
๐ฅ
Algorithm 2:
EMTP-
๐ฝ
algorithm.