Clear Filters
Clear Filters

Set cells in matrix equal to sqrt(10^-7), when the cells of the fft of the matrix obey 20.log10(a​bs(mat_tau​))<=-20

1 view (last 30 days)
I have a complex matrix mat_tau. I want for the cells of mat_tau which obey 20.log10(abs(mat_tau))<=-20 to equal to -70, and the same cells in ifft (mat_tau) to equal sqrt(10^-7). How can I do that?
demos serghiou
demos serghiou on 23 Oct 2022
Thanks for the reply. I just want the noise floor components that are less or equal to -20 to be removed from the ifft and can be set as zero instead of sqrt(10^-7). Any suggestions on how to code it?
dpb on 23 Oct 2022
Doesn't matter what the value is, destroying the symmetry and values of the FFT will muck up the IFFT drastically.
Unless the input signal is composed of simple sinusoids at the precise frequencies of the FFT bins, then zero-ing out bins in the frequency domain is equivalent to convolving the signal in the time domain with a series of overlapping rectangular windows at those frequencies. This will end up causing "ringing" in the inverse FFT.
Another way to think of it is that removing a given frequency bin is the same as adding the component at that frequency but with opposite phase. This also can be seen to be a bad idea in the time domain.

Sign in to comment.

Answers (1)

Nikhil on 27 Oct 2022
You can use logical indexing.
mat_tau(20.*log10(abs(mat_tau))<=-20) = -70
Similarly for any matrix.


Find more on Fourier Analysis and Filtering in Help Center and File Exchange

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!