On Sun, 17 Dec 2006 07:50:21 -0500, Greg Berchin <firstname.lastname@example.org>
>And here's a nice 48 kHz set:
Unfortunately, I just discovered that it's unstable. It happens,
sometimes, with approximation functions. I'll keep working on it and
post a stable solution if I find one.
Reply by Philip Martel●December 17, 20062006-12-17
"Jerry Avins" <email@example.com> wrote in message
> Al Clark wrote:
>> Jerry Avins <firstname.lastname@example.org> wrote in
>>> Greg Berchin wrote:
>>>> On 15 Dec 2006 22:23:53 -0800, "robert bristow-johnson"
>>>> <email@example.com> wrote:
>>>>> it oughta be okay to just apply the BLT* to the poles/zeros to get a
>>>>> digital filter.
>>>> Depends upon the sampling rate. That 12200 Hz pole pair makes life
>>>> difficult for the Bilinear Transform method if the sampling is 44100.
>>> What about upsampling before the filter and downsampling after it?
>>> Anyhow, is there a phase spec too? I see no reason for one in a
>>> weighting scheme.
>> This has always been one of my solution approaches. It costs a bit more
>> in computation. You still need to over sample by quite a bit if you want
>> the curve to be very accurate since warping is still significant at
> Prewarping will put the corner frequencies in the right place. The problem
> lies in the curve between them. It probably won't depart too severely even
> with a 3:1 oversampling rate, and 2:1 might be good enough. Has anyone
> tried? Experience beats speculation every time.
I would speculate that you are right, but...
> Engineering is the art of making what you want from things you can get.
Reply by Gerold Schrutz●December 18, 20062006-12-18
"robert bristow-johnson" <firstname.lastname@example.org> schrieb im Newsbeitrag
> Gerold Schrutz wrote:
>> Hello Bob,
>> I'll send it to you via email.
> why not post it? c'mon! spread the joy. it's the christmas season.
> r b-j
I've got an C++ Object with coefs, filter structure, ..., so its hard to
If you want to have it, just give me a note.
Reply by Robert Adams●December 18, 20062006-12-18
Presumably the order of this filter is low enough that one of the
factoring programs would work on it (to get cascaded biquads); if so,
then I would guess that some of the large coefficients causing rbj's
incontinence issues would be dramatically lower.
Greg Berchin wrote:
As promised, here's a STABLE 48kHz set. See other thread if interested
in explanation as to why the previous set was unstable.
This set is designed for best performance between 20 Hz and 20 kHz.
Reply by robert bristow-johnson●December 21, 20062006-12-21
Greg Berchin wrote:
> As promised, here's a STABLE 48kHz set. See other thread if interested
> in explanation as to why the previous set was unstable.
> This set is designed for best performance between 20 Hz and 20 kHz.
so what did you have to do? reflect the unstable poles into the
interior of the unit circle?
Reply by Greg Berchin●December 21, 20062006-12-21
On 21 Dec 2006 14:10:13 -0800, "robert bristow-johnson"
>so what did you have to do? reflect the unstable poles into the
>interior of the unit circle?
Naw, I just implemented the algorithm CORRECTLY.
Reply by Randy Yates●December 21, 20062006-12-21
Greg Berchin <email@example.com> writes:
> On 21 Dec 2006 14:10:13 -0800, "robert bristow-johnson"
> <firstname.lastname@example.org> wrote:
>>so what did you have to do? reflect the unstable poles into the
>>interior of the unit circle?
> Naw, I just implemented the algorithm CORRECTLY.
What I hear you saying is that, for a specific design, you had
a problem in your implementation that resulted in an unstable
set of coefficients.
What I'm wondering is if it can be proven that the algorithm in
general will produce a stable IIR transfer function, i.e., how can you
guarantee the poles it generates are always inside the unit circle?
% Randy Yates % "And all that I can do
%% Fuquay-Varina, NC % is say I'm sorry,
%%% 919-577-9882 % that's the way it goes..."
%%%% <email@example.com> % Getting To The Point', *Balance of Power*, ELO