Input:
Measurement matrix
๐ด
, measurements
y
, sparsity level
๐
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
)
๐ผ
(
๐ก
)
=
{
๐
โถ
๐ฅ
๐
(
๐ก
)
โ
0
}
(
9
)
๐ฅ
๐ผ
(
๐ก
)
(
๐ก
)
=
๐ด
โ
๐ผ
(
๐ก
)
๐ฆ
(
1
0
)
๐ฅ
(
๐ก
)
[
1
,
๐
]
โงต
๐ผ
(
๐ก
)
=
0
(
1
1
)
๐
(
๐ก
)
=
๐ฆ
โ
๐ด
๐ฅ
(
๐ก
)
(
1
2
)
๐ก
=
๐ก
+
1
(
1
3
)
end while
(
1
4
)
return
๐ฅ
Algorithm 1:
EMTP-
๐
algorithm.