|
From: | Paul Kienzle |
Subject: | Re: complex output of pwelch |
Date: | Sun, 11 Mar 2007 13:29:51 -0400 |
On Mar 11, 2007, at 1:56 AM, Dmitri A. Sergatskov wrote:
I noticed that at some moment I started getting a complex output while doing b = pwelch(a) (which must be real). The imaginary part is order of eps, but it is still an annoyance. Should we put Pxx = real(Pxx) or there a better fix?
I modified pwelch.m to use real(x.*conj(x)) instead of conj(x). It should now be available on octave-forge.
- Paul
[Prev in Thread] | Current Thread | [Next in Thread] |