2
$\begingroup$

I am given $|H(\omega)|$, I wonder if minimum phase stable causal filter is unique and how to calculate it.

$\endgroup$
  • $\begingroup$ if you're wanting to design a finite-order IIR (determine the order $N$ and the $2N+1$ coefficients, you have to use something like the Prony method or Greg Berchin's FDLS. $\endgroup$ – robert bristow-johnson Mar 19 '15 at 21:04
  • $\begingroup$ @robertbristow-johnson: Prony's method designs an IIR filter given a prescribed impulse response. How would you use it here when the magnitude of the frequency response is given? $\endgroup$ – Matt L. Mar 20 '15 at 8:08
  • $\begingroup$ @MattL., i guess i would first create a minimum-phase complex frequency response, and compute an impulse response from that. or use Greg's FDLS, since it's frequency domain. $\endgroup$ – robert bristow-johnson Mar 20 '15 at 19:40
2
$\begingroup$

If $H(\omega)=e^{\alpha(\omega)+j\phi(\omega)}$ is a minimum phase frequency response, then the attenuation $\alpha(\omega)$ and the phase $\phi(\omega)$ are related by the following Hilbert transform relationship:

$$\phi(\omega)=-\frac{1}{\pi}\int_{-\infty}^{\infty}\frac{\alpha(\Omega)}{\omega-\Omega}d\Omega$$

So $\phi(\omega)$ is uniquely determined by $\alpha(\omega)$. This is the corresponding wikipedia entry.

$\endgroup$
  • $\begingroup$ @user40129: I just looked around but couldn't find anything convincing. In 'The Fourier Integral and Its Applications' by Papoulis there's proof of a similar relationship using contour integration. And in Oppenheim and Schafer's book (Discrete-Time Signal Processing) there's something like a proof of the equivalent relationship for discrete-time systems, but part of it is left as a 'problem'. $\endgroup$ – Matt L. Mar 19 '15 at 18:31
  • $\begingroup$ Proof is at Page 206 of 'The Fourier Integral and Its Applications' by Papoulis. Thanks! $\endgroup$ – user40129 Mar 19 '15 at 19:50

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Not the answer you're looking for? Browse other questions tagged or ask your own question.