This bug is now fixed. Thanks for your report and (correct) diagnosis, Peter
On Thu, Mar 20, 2008 at 12:26:37PM +0100, Michiel Koot wrote: > Dear Peter Lanspeary, > > First of all, thank you for writing and sharing the pwelch routine for > Octave. > > However, I do think that there is an error in how the transfer function is > computed, because now it appears to be the complex conjugate, with the > phase being inversed. > > The way it is computed now is roughly as follows: > > Txy = Pxy./Pxx = Fx.*conj(Fy)./(Fx.*conj(Fx)) = conj(Fy)./conj(Fx) = > conj(Fy./Fx) > > whereas it should be: > > Txy = Fy./Fx > > Could you please verify and correct this? > > With kind regards, > > Michiel Koot ------------------------------------------------------------------------- This SF.net email is sponsored by: Microsoft Defy all challenges. Microsoft(R) Visual Studio 2008. http://clk.atdmt.com/MRT/go/vse0120000070mrt/direct/01/ _______________________________________________ Octave-dev mailing list Octave-dev@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/octave-dev
