What does it mean to deconvolve the impulse response - Signal Processing Stack Exchange most recent 30 from dsp.stackexchange.com 2019-07-22T21:57:24Z https://dsp.stackexchange.com/feeds/question/22111 http://www.creativecommons.org/licenses/by-sa/3.0/rdf https://dsp.stackexchange.com/q/22111 1 What does it mean to deconvolve the impulse response Celdor https://dsp.stackexchange.com/users/4721 2015-03-16T11:28:35Z 2015-03-16T12:26:24Z <p>I am learning DPS and I came across the problem of deconvolution and removing the impulse response from a signal? This still does little sense to me.</p> <p>My understanding of the impulse response is to find out the response, usually in frequency, of a device, a filter, a room etc. If I don't mess things here, to me, the impulse response in signal processing is equivalent to numeric 1. Finding the impulse response is convolving a signal with something that is equivalent to 1 that is the impulse response. So if I need to find what the value of x is, I simply multiply x by 1:</p> <p>$$y = x\times1$$</p> <p>Therefore if $y = 0.8$, I know $x = 0.8$ as well.</p> <p>In DSP, we have the equation</p> <p>$$Y = XH$$</p> <p>So to know the impulse response of $X$, I multiply $X$ by $H$ and obtain $Y$ which is equal to $X$. I've gone that far.</p> <p>I don't understand what deconvolution means. I though by applying the impulse response we have already an answer about the frequency of anything we need. What does deconvolution mean here and why would we look for it?</p> <p>Thanks.</p> https://dsp.stackexchange.com/questions/22111/-/22112#22112 4 Answer by Matt L. for What does it mean to deconvolve the impulse response Matt L. https://dsp.stackexchange.com/users/4298 2015-03-16T11:47:33Z 2015-03-16T12:01:00Z <p>A linear time-invariant (LTI) system is completely characterized by its impulse response $h(t)$. The output signal $y(t)$ given an arbitrary input signal $x(t)$ is given by the convolution of this input signal with the system's impulse response:</p> <p>$$y(t)=(x\star h)(t)\tag{1}$$</p> <p>where $\star$ denotes convolution. In the frequency domain, Eq. (1) becomes a simple multiplication:</p> <p>$$Y(f)=X(f)H(f)$$</p> <p>where $Y(f)$, $X(f)$, and $H(f)$ are the Fourier transforms of the signals $y(t)$, $x(t)$, and $h(t)$, respectively.</p> <p><em>Deconvolution</em> now refers to the process of "undoing" the effect of an LTI system. I.e. you're given the signal $y(t)$ and you want to recover the original signal $x(t)$. E.g., you might want to remove the reverberation added to a signal by a given room. In this case you need to undo the effect of the convolution of the room's impulse response with the original signal.</p> <p>Naive deconvolution can in principle be done by applying an inverse filter to the signal $y(t)$. The inverse filter ideally has a frequency response $G(f)=1/H(f)$. Note that such a filter does not necessarily exist, e.g. if $H(f)$ has zeros at certain frequencies. Another disadvantage of this naive deconvolution is the effect it has on measurement noise. In frequency regions where $H(f)$ has a small magnitude, the magnitude of $G(f)$ becomes large, and the noise is highly amplified. This is why applying an inverse filter is usually not a good solution. If the system's impulse response is known and if the noise can be modeled as (wide-sense) stationary, <a href="http://en.wikipedia.org/wiki/Wiener_deconvolution" rel="nofollow">Wiener deconvolution</a> can be used. In that case, the original signal is estimated using a mean square error criterion, which takes into account the impulse response as well as the noise. If the impulse response is unknown (or can't be estimated), there are methods available referred to as <a href="http://en.wikipedia.org/wiki/Blind_deconvolution" rel="nofollow">blind deconvolution</a>.</p>