Re: Signal Package Filter Design

[Dmitri A. Sergatskov]

On Sun, Sep 3, 2017 at 4:52 PM, Doug Stewart <address@hidden> wrote:

On Sun, Sep 3, 2017 at 3:50 PM, Mike Miller <address@hidden> wrote:

On Sun, Sep 03, 2017 at 03:12:42 -0700, Thomas D. Dean wrote:
> I am attempting to duplicate a python application in octave.
>
> octave:112> pkg list signal
> Package Name  | Version | Installation directory
> --------------+---------+-----------------------
>       signal *|   1.3.2 | /home/tomdean/octave/signal-1.3.2
>
>
> remez(n_taps, [0, f_bw, f_bw+(Fs/2-f_bw)/4,Fs/2],[1,0],Hz=Fs)
>
> > pwd
> /home/tomdean/octave/signal-1.3.2
> > ls x86_64-pc-linux-gnu-api-v51/
> __fwht__.oct     cl2bp.oct  medfilt1.oct  sosfilt.oct
> __ultrwin__.oct  doc-cache  remez.oct     upfirdn.oct
>
> Note, no PKG_ADD file
>
> The remez.oct file is there, but, there does not seem to be any way to call
> it.
>
> How do I do this?

Calling remez works for me.

--
mike

_

It does not work for me.

Tried with latest stable and dev.

--------------+---------+-----------------------

     control *|   3.0.0 | /home/doug/octave/control-3.0.0

      signal *|   1.3.2 | /home/doug/octave/signal-1.3.2

    symbolic *|   2.5.0 | /home/doug/octave/symbolic-2.5.0

>> help remez

error: get_help_text: /home/doug/octave/signal-1.3.2/x86_64-pc-linux-gnu-api-v51/remez.oct: failed to load: /home/doug/octave/signal-1.3.2/x86_64-pc-linux-gnu-api-v51/remez.oct: undefined symbol: _ZTIN6octave12chunk_bufferE

error: called from

    help at line 95 column 18

>>

Doug

​Works for me (with dev):

​octave:5>  help remez
'remez' is a function from the file /home/dima/octave/signal-1.3.2/x86_64-pc-linux-gnu-api-v51/remez.oct

 -- Loadable Function: B = remez (N, F, A)
 -- Loadable Function: B = remez (N, F, A, W)
 -- Loadable Function: B = remez (N, F, A, W, FTYPE)
 -- Loadable Function: B = remez (N, F, A, W, FTYPE, GRIDDENSITY)
     Parks-McClellan optimal FIR filter design.
     N
          gives the number of taps in the returned filter
     F
          gives frequency at the band edges [b1 e1 b2 e2 b3 e3 ...]
     A
          gives amplitude at the band edges [a(b1) a(e1) a(b2) a(e2)
          ...]
     W
          gives weighting applied to each band
     FTYPE
          is "bandpass", "hilbert" or "differentiator"
     GRIDDENSITY
          determines how accurately the filter will be constructed.  The
          minimum value is 16, but higher numbers are slower to compute.

     Frequency is in the range (0, 1), with 1 being the Nyquist
     frequency.

octave:6> remez(n_taps, [0, f_bw, f_bw+(Fs/2-f_bw)/4,Fs/2],[1,0],Hz=Fs)
error: 'n_taps' undefined near line 1 column 7