Hello,
I am currently working with the attenComp.m function and noticed on line 510 that the lower triangle coefficients are set to zeros twice.
% zero out lower and upper triangles
ones_mat = ones(N, N);
tv_filter(tril(ones_mat, -ceil(N/2) + 1) == 1) = 0;
tv_filter(tril(ones_mat, -ceil(N/2) + 1) == 1) = 0;
Shouldn't it be something like:
% zero out lower and upper triangles
ones_mat = ones(N, N);
tv_filter(tril(ones_mat, -ceil(N/2) + 1) == 1) = 0;
tv_filter(triu(ones_mat, ceil(N/2) + 1) == 1) = 0;