

I also asked remez to work harder by increasing maxiter. hope yall beings safe out here and doing well. I relaxed the transition bandwidth to 2 Hz and increased the order to 300. I gave it a shot with the bandwidth you need (20 Hz). You may need to try a few combinations before you get a filter you like. The degrees of freedom you have available are: filter bandwidth, filter order, and transistion bandwidth.

# Order hard set to 200 as in MATLAB codeĪnd then to take a look at what’s going on with a fft filterweights = (firls(200,ffrequencies,idealresponse)) įft_filtkern = fft_filtkern./maximum(fft_filtkern) # normalized to 1.0 for visual comparison ease (1+transition_width)*(center_freq+filter_frequency_spread), (1-transition_width)*(center_freq-filter_frequency_spread), center_freq = 20 # in Hzįilter_frequency_spread = 6 # Hz +/- the center frequency
#Filter designer matlab code#
I checked out the documentation, and did the following to reproduce something similar to the MATLAB code I was using. What exactly is this value? And is it the same as the filter order for butterworth and the “order- n FIR filter” (sorry if this latter part may be beyond the scope of this discussion) Design of FIR Filter Using Window method in Matlab - Stack Overflow Step 2: Designing a Filter :: SPTool: A. Is there a way to make a smoother transition (assuming this is the reason) like this in Julia?Īlso the filter order here comes out to something like 192, this doesn’t work with Butterworth(192). And then the shape of the filter is defined by idealresponse. So ffrequencies defines the frequencies of interest, normalized to the nyquist frequency. In this episode Brady Volpe of Nimble This and The Volpe Firm, and John Downey (of Cisco) answer viewers questions about upstream SNR / Diplex filters. One of the things that author does is define a specific shape of the filter (I think to make a smoother transition?) with the following Matlab code: nyquist = EEG.srate/2 įilter_order = round(3*(EEG.srate/lower_filter_bound)) įfrequencies = /nyquist įilterweights = firls(filter_order,ffrequencies,idealresponse) Thank you for your replies! They really helped me getting started.
