From: Pete Fraser on 1 Jul 2010 21:48 "Vladimir Vassilevsky" <nospam(a)nowhere.com> wrote in message news:OP-dnf44bO7hJ7HRnZ2dnUVZ_t-dnZ2d(a)giganews.com... > Check the FDLS method of Greg Berchin. That is one of the most elegant > methods of designing to a prototype. I think I must be doing something dumb. I started off with a really simple half-band FIR, just as a test. The input file was: 11 100000 0 1 0 1 5000 1 -72 1 10000 1 -144 1 15000 1 -216 1 20000 1 -288 1 25000 .5 -360 1 30000 0 -432 1 35000 0 -504 1 40000 0 -576 1 45000 0 -648 1 50000 0 -720 1 N=8 D=0 delta=0 and I get B = -3.55068459956424e-018 -0.0454056666613377 2.42666919829597e-016 0.292936161872951 0.5 0.326331656363899 -2.00804673259988e-016 -0.0874864446451108 -1.76118166519034e-017 so it's doing a lot of things right. The center is 0.5, four are close to 0, but 0.29 is not 0.32, and -0.04 is not -0.08. What am I doing wrong? Thanks Pete
|
Pages: 1 Prev: Catalytic Fixed-Point Toolbox Next: Using mpcmove from embedded matlab R14 |