where denotes the thresholding operator for the singular-values, and denotes the threshold of the operator. Defining an optimal threshold , however, is often challenging because a constant threshold is inadequate to handle an inhomogeneous distribution of noise energy. Huang et al. (2016) developed a more elegant way to shrinkage the singular values by deriving a damping operator:

where is called the damping matrix, is an identity matrix, and is referred to as the damping threshold matrix. is the th singular value in the un-truncated singular matrix . in equation 9 is called the damping factor, which is used to control the strength of the damping operator. In a special case, when , , indicating that the damped rank-reduction method reverts to the traditional rank-reduction method. The mathematical details to derive the damping operator can be found in Huang et al. (2016) and Huang et al. (2017). The damping operator is used to shrinkage the singular values in the TSVD formula to reduce the residual noise:

Compared with the thresholding method in equation 7, where a rigid threshold is used to shrinkage the singular values, the damping operator applies variable thresholds to different singular values. According to equation 8, the threshold decreases as the residual noise becomes less dominant, i.e., the singular value becomes larger. The algorithm workflow for the damped rank-reduction method (DRR) is outlined in Algorithm 2.

2020-12-06