Re: Remez Octave Signal vs Scipi.signal

[Dmitri A. Sergatskov]

On Sun, Sep 3, 2017 at 10:49 PM, Thomas D. Dean <address@hidden> wrote:

figure()
hold on
subplot(3,1,1)
pwelch(x2,4,0,2048,Fs)
title('pwelch x2')
subplot(3,1,2)
pwelch(x4_lpf,4,0,2048,Fs_y)
title('pwelch x4_lpf')
subplot(3,1,3)
pwelch(x4_filt,4,0,2048,Fs_y)
title('pwelch x4_filt')
hold off

​I am not a digital filter guy, but

​

pwelch(x4_lpf,4,0,2048,Fs_y)

​does not make any sense to me.

​Try

​pwelch(x4_lpf,[],[],2048,Fs_y, "loglog", "linear")

(you pick up 4 data-point window, pad it with 0 to 2048, do fft on that,

make a zero overlap -- what is the point of pwelch then?)

Also there is a decimate() function in signal pkg.

​Dmitri.

--

​